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.     

A theoretical model for the association of amphiphilic transmembrane peptides in lipid bilayers

A theoretical model is proposed for the association of trans-bilayer peptides in lipid bilayers. The model is based on a lattice model for the pure lipid bilayer, which accounts accurately for the most important conformational states of the lipids and their mutual interactions and statistics. Within the lattice formulation the bilayer is formed by two independent monolayers, each represented by a triangular lattice, on which sites the lipid chains are arrayed. The peptides are represented by regular objects, with no internal flexibility, and with a projected area on the bilayer plane corresponding to a hexagon with seven lattice sites. In addition, it is assumed that each peptide surface at the interface with the lipid chains is partially hydrophilic, and therefore interacts with the surrounding lipid matrix via selective anisotropic forces. The peptides would therefore assemble in order to shield their hydrophilic residues from the hydrophobic surroundings. The model describes the self-association of peptides in lipid bilayers via lateral and rotational diffusion, anisotropic lipid-peptide interactions, and peptide-peptide interactions involving the peptide hydrophilic regions. The intent of this model study is to analyse the conditions under which the association of trans-bilayer and partially hydrophilic peptides (or their dispersion in the lipid matrix) is lipid-mediated, and to what extent it is induced by direct interactions between the hydrophilic regions of the peptides. The model properties are calculated by a Monte Carlo computer simulation technique within the canonical ensemble. The results from the model study indicate that direct interactions between the hydrophilic regions of the peptides are necessary to induce peptide association in the lipid bilayer in the fluid phase. Furthermore, peptides within each aggregate are oriented in such a way as to shield their hydrophilic regions from the hydrophobic environment. The average number of peptides present in the aggregates formed depends on the degree of mismatch between the peptide hydrophobic length and the lipid bilayer hydrophobic thickness: The lower the degree of mismatch is the higher this number is. Received: 30 December 1996 / Accepted: 9 May 1997     

A structural dissection of amino acid substitutions in helical transmembrane proteins

The evolution of protein folds is under strong constraints from their surrounding environment. Although folding in water‐soluble proteins is driven primarily by hydrophobic forces, the nature of the forces that determine the folding and stability of transmembrane proteins are still not fully understood. Furthermore, the chemically heterogeneous lipid bilayer has a non‐uniform effect on protein structure. In this article, we attempt to get an insight into the nature of this effect by examining the impact of various types of local structure environment on amino acid substitution, based on alignments of high‐resolution structures of polytopic helical transmembrane proteins combined with sequences of close homologs. Compared to globular proteins, burying amino acid sidechains, especially hydrophilic ones, led to a lower increase in conservation in both the lipid‐water interface region and the hydrocarbon core region. This observation is due to surface residues in HTM proteins especially in the HC region being relatively highly conserved, suggesting higher evolutionary constraints from their specific interactions with the surrounding lipid molecules. Polar and small residues, particularly Pro and Gly, show a noticeable increase in conservation as they are positioned more towards the centre of the membrane, which is consistent with their recognized key roles in structural stability. In addition, the examination of hydrogen bonds in the membrane environment identified some exposed hydrophilic residues being better conserved when not hydrogen‐bonded to other residues, supporting the importance of lipid‐protein sidechain interactions. The conclusions presented in this study highlight the distinct features of substitution matrices that take into account the membrane environment, and their potential role in improving sequence‐structure alignments of transmembrane proteins. Proteins 2010; © 2010 Wiley‐Liss, Inc.     

Revisiting Hydrophobic Mismatch with Free Energy Simulation Studies of Transmembrane Helix Tilt and Rotation

Protein-lipid interaction and bilayer regulation of membrane protein functions are largely controlled by the hydrophobic match between the transmembrane (TM) domain of membrane proteins and the surrounding lipid bilayer. To systematically characterize responses of a TM helix and lipid adaptations to a hydrophobic mismatch, we have performed a total of 5.8-μs umbrella sampling simulations and calculated the potentials of mean force (PMFs) as a function of TM helix tilt angle under various mismatch conditions. Single-pass TM peptides called WALPn (n = 16, 19, 23, and 27) were used in two lipid bilayers with different hydrophobic thicknesses to consider hydrophobic mismatch caused by either the TM length or the bilayer thickness. In addition, different flanking residues, such as alanine, lysine, and arginine, instead of tryptophan in WALP23 were used to examine their influence. The PMFs, their decomposition, and trajectory analysis demonstrate that 1), tilting of a single-pass TM helix is the major response to a hydrophobic mismatch; 2), TM helix tilting up to ∼10° is inherent due to the intrinsic entropic contribution arising from helix precession around the membrane normal even under a negative mismatch; 3), the favorable helix-lipid interaction provides additional driving forces for TM helix tilting under a positive mismatch; 4), the minimum-PMF tilt angle is generally located where there is the hydrophobic match and little lipid perturbation; 5), TM helix rotation is dependent on the specific helix-lipid interaction; and 6), anchoring residues at the hydrophilic/hydrophobic interface can be an important determinant of TM helix orientation.     

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.     

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.     

Effect of lipid structural modifications on their intermolecular hydrogen bonding interactions and membrane functions

The large number of different membrane lipids with various structural modifications and properties and the characteristic lipid composition of different types of membranes suggest that different lipids have specific functions in the membrane. Many of the varying properties of lipids with different polar head groups and in different ionization states can be attributed to the presence of interactive or repulsive forces between the head groups in the bilayer. The interactive forces are hydrogen bonds between hydrogen bond donating groups such as --P--OH,--OH, and--NH3+ and hydrogen bond accepting groups such as --P--O- and --COO-. These interactions increase the lipid phase transition temperature and can account for the tendency of certain lipids to go into the hexagonal phase and the dependence of this tendency on the pH and ionization state of the lipid. The presence or absence of these interactions can also affect the penetration of hydrophobic substances into the bilayer, including hydrophobic residues of membrane proteins. Evidence for this suggestion has been gathered from studies of the myelin basic protein, a water-soluble protein with a number of hydrophobic residues. In this way the lipid composition can affect the conformation and activity of membrane proteins. Since hydrogen-bonding interactions depend on the ionization state of the lipid, they can be altered by changes in the environment which affect the pK of the ionizable groups. The formation of the hexagonal phase or inverted micelles, the conformation and activity of membrane proteins, and other functions mediated by lipids could thus be regulated in this way.     

Electrostatic interactions during acidic phospholipid reactivation of DnaA protein, the Escherichia coli initiator of chromosomal replication.

The initiation of Escherichia coli chromosomal replication by DnaA protein is strongly influenced by the tight binding of the nucleotides ATP and ADP. Anionic phospholipids in a fluid bilayer promote the conversion of inactive ADP-DnaA protein to replicatively active ATP-DnaA protein in vitro, and thus likely play a key role in regulating DnaA activity. Previous studies have revealed that, during this reactivation, a specific region of DnaA protein inserts into the hydrophobic portion of the lipid bilayer in an acidic phospholipid-dependent manner. To elucidate the requirement for acidic phospholipids in the reactivation process, the contribution of electrostatic forces in the interaction of DnaA and lipid was examined. DnaA-lipid binding required anionic phospholipids, and DnaA-lipid binding as well as lipid-mediated release of DnaA-bound nucleotide were inhibited by increased ionic strength, suggesting the involvement of electrostatic interactions in these processes. As the vesicular content of acidic phospholipids was increased, both nucleotide release and DnaA-lipid binding increased in a linear, parallel manner. Given that DnaA-membrane binding, the insertion of DnaA into the membrane, and the consequent nucleotide release all require anionic phospholipids, the acidic headgroup may be necessary to recruit DnaA protein to the membrane for insertion and subsequent reactivation for replication.     

Nature of Interactions between PEO-PPO-PEO Triblock Copolymers and Lipid Membranes: (II) Role of Hydration Dynamics Revealed by Dynamic Nuclear Polarization

Amphiphilic poly(ethylene oxide)-poly(propylene oxide)-poly(ethylene oxide) copolymers, also known as poloxamers, have broad biomembrane activities. To illustrate the nature of these activities, (1)H Overhauser dynamic nuclear polarization NMR spectroscopy was employed to sensitively detect polymer-lipid membrane interactions through the modulation of local hydration dynamics in lipid membranes. Our study shows P188, the most hydrophilic poloxamer that is a known membrane sealant, weakly adsorbs on the membrane surface, yet effectively retards membrane hydration dynamics. Contrarily, P181, the most hydrophobic poloxamer that is a known membrane permeabilizer, initially embeds at lipid headgroups and enhances intrabilayer water diffusivity. Unprecedented resolution for differentiating weak surface adsorption versus translocation of polymers to membranes is obtained by probing local water diffusivity in lipid bilayer systems. Our results illustrate that the relative hydrophilic/hydrophobic ratio of the polymer dictates its functions. These findings gleaned from local hydration dynamics are well supported by a thermodynamics study presented in the accompanying paper (Wang, J.-Y.; Marks, J. M.; Lee, K. Y. C. Biomacromolecules, 2012, DOI: 10.1021/bm300847x).     

Specific interactions of tryptophan with phosphatidylcholine and digalactosyldiacylglycerol in pure and mixed bilayers in the dry and hydrated state

Amphiphilic solutes play an important role in the desiccation tolerance of plant cells, because they can reversibly partition into cellular membranes during dehydration. Their effects on membrane stability depend on their chemical structure, but also on the lipid composition of the host membrane. We have shown recently that tryptophan destabilizes liposomes during freezing. The degree of destabilization depends on the presence of glycolipids in the membranes, but not on the phase preference (bilayer or non-bilayer) of the lipids in mixtures with the bilayer lipid phosphatidylcholine. Here, we have investigated the influence of tryptophan on the phase behavior and intermolecular interactions in dry and hydrated bilayers made from the phospholipid egg phosphatidylcholine and the plant chloroplast glycolipid digalactosyldiacylglycerol, or from a mixture (1:1) of these lipids, using Fourier-transform infrared spectroscopy. To distinguish effects of the hydrophobic ring structure of tryptophan from those of the amino acid moiety, we also performed experiments with the hydrophilic amino acid glycine. Our data show that there are specific interactions between tryptophan and either phospholipid or glycolipid in the dry state, as well as H-bonding interactions between the lipids and both solutes. In the rehydrated state, the H-bonding interactions between amino acids and lipids are mostly replaced by interactions between water and lipids, while the hydrophobic interactions between lipids and tryptophan mostly persist.     

Peptide Nanopores and Lipid Bilayers: Interactions by Coarse-Grained Molecular-Dynamics Simulations

A set of 49 protein nanopore-lipid bilayer systems was explored by means of coarse-grained molecular-dynamics simulations to study the interactions between nanopores and the lipid bilayers in which they are embedded. The seven nanopore species investigated represent the two main structural classes of membrane proteins (α-helical and β-barrel), and the seven different bilayer systems range in thickness from ∼28 to ∼43 Å. The study focuses on the local effects of hydrophobic mismatch between the nanopore and the lipid bilayer. The effects of nanopore insertion on lipid bilayer thickness, the dependence between hydrophobic thickness and the observed nanopore tilt angle, and the local distribution of lipid types around a nanopore in mixed-lipid bilayers are all analyzed. Different behavior for nanopores of similar hydrophobic length but different geometry is observed. The local lipid bilayer perturbation caused by the inserted nanopores suggests possible mechanisms for both lipid bilayer-induced protein sorting and protein-induced lipid sorting. A correlation between smaller lipid bilayer thickness (larger hydrophobic mismatch) and larger nanopore tilt angle is observed and, in the case of larger hydrophobic mismatches, the simulated tilt angle distribution seems to broaden. Furthermore, both nanopore size and key residue types (e.g., tryptophan) seem to influence the level of protein tilt, emphasizing the reciprocal nature of nanopore-lipid bilayer interactions.     

Hydrophobic forces drive spontaneous membrane insertion of the bacteriophage Pf3 coat protein without topological control

Bacterial integral inner membrane proteins are either translocated across the lipid bilayer using an energy-driven enzyme, such as the Sec translocase, or they might interact directly with the membrane due to hydrophobic forces. We report that the single-spanning Pf3 coat protein is spontaneously inserted into the membrane of Escherichia coli and requires the electrical component of the membrane potential (DeltaPsi) to translocate its N-terminal region. This results in a final N(out)C(in) orientation of the protein in the cytoplasmic membrane, due the potential-driven translocation of the aspartyl residue at position 18 in the hydrophilic N-terminal tail. Uncharged protein tails are only translocated when the hydrophobic transmembrane region of the protein has been extended. An extended transmembrane anchor allows membrane insertion in the absence of an electrochemical membrane potential, but also causes the loss of a strict determination of the topology.     

Probing the interaction forces between hydrophobic peptides and supported lipid bilayers using AFM

Despite the vast body of literature that has accumulated on tilted peptides in the past decade, direct information on the forces that drive their interaction with lipid membranes is lacking. Here, we attempted to use atomic force microscopy (AFM) to explore the interaction forces between the Simian immunodeficiency virus peptide and phase-separated supported bilayers composed of various lipids, i.e. dipalmitoylphosphatidylcholine, dioleoylphosphatidylcholine, dioleoylphosphatidic acid and dipalmitoylphosphatidylethanolamine. Histidine-tagged peptides were attached onto AFM tips terminated with nitrilotriacetate and tri(ethylene glycol) groups, an approach expected to ensure optimal exposure of the C-terminal hydrophobic domain. Force-distance curves recorded between peptide-tips and the different bilayer domains always showed a long-range repulsion upon approach and a lack of adhesion upon retraction, in marked contrast with the hydrophobic nature of the peptide. To explain this unexpected behaviour, we suggest a mechanism in which lipids are pulled out from the bilayer due to strong interactions with the peptide-tip, in agreement with the very low force needed to extract lipids from supported bilayers.     

High-resolution AFM of membrane proteins directly incorporated at high density in planar lipid bilayer

The heterologous expression and purification of membrane proteins represent major limitations for their functional and structural analysis. Here we describe a new method of incorporation of transmembrane proteins in planar lipid bilayer starting from 1 pmol of solubilized proteins. The principle relies on the direct incorporation of solubilized proteins into a preformed planar lipid bilayer destabilized by dodecyl-beta-maltoside or dodecyl-beta-thiomaltoside, two detergents widely used in membrane biochemistry. Successful incorporations are reported at 20 degrees C and at 4 degrees C with three bacterial photosynthetic multi-subunit membrane proteins. Height measurements by atomic force microscopy (AFM) of the extramembraneous domains protruding from the bilayer demonstrate that proteins are unidirectionally incorporated within the lipid bilayer through their more hydrophobic domains. Proteins are incorporated at high density into the bilayer and on incubation diffuse and segregate into protein close-packing areas. The high protein density allows high-resolution AFM topographs to be recorded and protein subunits organization delineated. This approach provides an alternative experimental platform to the classical methods of two-dimensional crystallization of membrane proteins for the structural analysis by AFM. Furthermore, the versatility and simplicity of the method are important intrinsic properties for the conception of biosensors and nanobiomaterials involving membrane proteins.     

Interaction of the Mechanosensitive Channel,MscS, with the Membrane Bilayer through Lipid Intercalation into Grooves and Pockets

All membrane proteins have dynamic and intimate relationships with the lipids of the bilayer that may determine their activity. Mechanosensitive channels sense tension through their interaction with the lipids of the membrane. We have proposed a mechanism for the bacterial channel of small conductance, MscS, that envisages variable occupancy of pockets in the channel by lipid chains. Here, we analyze protein–lipid interactions for MscS by quenching of tryptophan fluorescence with brominated lipids. By this strategy, we define the limits of the bilayer for TM1, which is the most lipid exposed helix of this protein. In addition, we show that residues deep in the pockets, created by the oligomeric assembly, interact with lipid chains. On the cytoplasmic side, lipids penetrate as far as the pore-lining helices and lipid molecules can align along TM3b perpendicular to lipids in the bilayer. Cardiolipin, free fatty acids, and branched lipids can access the pockets where the latter have a distinct effect on function. Cholesterol is excluded from the pockets. We demonstrate that introduction of hydrophilic residues into TM3b severely impairs channel function and that even “conservative” hydrophobic substitutions can modulate the stability of the open pore. The data provide important insights into the interactions between phospholipids and MscS and are discussed in the light of recent developments in the study of Piezo1 and TrpV4.     

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.