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

Sterol affinity for phospholipid bilayers is influenced by hydrophobic matching between lipids and transmembrane peptides

Lipid self-organization is believed to be essential for shaping the lateral structure of membranes, but it is becoming increasingly clear that also membrane proteins can be involved in the maintenance of membrane architecture. Cholesterol is thought to be important for the lateral organization of eukaryotic cell membranes and has also been implicated to take part in the sorting of cellular transmembrane proteins. Hence, a good starting point for studying the influence of lipid-protein interactions on membrane trafficking is to find out how transmembrane proteins influence the lateral sorting of cholesterol in phospholipid bilayers. By measuring equilibrium partitioning of the fluorescent cholesterol analog cholestatrienol between large unilamellar vesicles and methyl-β-cyclodextrin the effect of hydrophobic matching on the affinity of sterols for phospholipid bilayers was determined. Sterol partitioning was measured in 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC), 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) bilayers with and without WALP19, WALP23 or WALP27 peptides. The results showed that the affinity of the sterol for the bilayers was affected by hydrophobic matching. An increasing positive hydrophobic mismatch led to stronger sterol binding to the bilayers (except in extreme situations), and a large negative hydrophobic mismatch decreased the affinity of the sterol for the bilayer. In addition, peptide insertion into the phospholipid bilayers was observed to depend on hydrophobic matching. In conclusion, the results showed that hydrophobic matching can affect lipid-protein interactions in a way that may facilitate the formation of lateral domains in cell membranes. This could be of importance in membrane trafficking.

Altering hydrophobic sequence lengths shows that hydrophobic mismatch controls affinity for ordered lipid domains (rafts) in the multitransmembrane strand protein perfringolysin O

The hypothesis that mismatch between transmembrane (TM) length and bilayer width controls TM protein affinity for ordered lipid domains (rafts) was tested using perfringolysin O (PFO), a pore-forming cholesterol-dependent cytolysin. PFO forms a multimeric barrel with many TM segments. The properties of PFO mutants with lengthened or shortened TM segments were compared with that of PFO with wild type TM sequences. Both mutant and wild type length PFO exhibited cholesterol-dependent membrane insertion. Maximal PFO-induced pore formation occurred in vesicles with wider bilayers for lengthened TM segments and in thinner bilayers for shortened TM segments. In diC(18:0) phosphatidylcholine (PC)/diC(14:1) PC/cholesterol vesicles, which form ordered domains with a relatively thick bilayer and disordered domains with a relatively thin bilayer, affinity for ordered domains was greatest with lengthened TM segments and least with shortened TM segments as judged by FRET. Similar results were observed by microscopy in giant vesicles containing sphingomyelin in place of diC(18:0) PC. In contrast, in diC(16:0) PC/diC(14:0) PC/diC(20:1) PC/cholesterol vesicles, which should form ordered domains with a relatively thin bilayer and disordered domains with a relatively thick bilayer, relative affinity for ordered domains was greatest with shortened TM segments and least with lengthened TM segments. The inability of multi-TM segment proteins (unlike single TM segment proteins) to adapt to mismatch by tilting may explain the sensitivity of raft affinity to mismatch. The difference in width sensitivity for single and multi-TM helix proteins may link raft affinity to multimeric state and thus control the assembly of multimeric TM complexes in rafts.

The determinants of hydrophobic mismatch response for transmembrane helices

Hydrophobic mismatch arises from a difference in the hydrophobic thickness of a lipid membrane and a transmembrane protein segment, and is thought to play an important role in the folding, stability and function of membrane proteins. We have investigated the possible adaptations that lipid bilayers and transmembrane α-helices undergo in response to mismatch, using fully-atomistic molecular dynamics simulations totaling 1.4 μs. We have created 25 different tryptophan-alanine-leucine transmembrane α-helical peptide systems, each composed of a hydrophobic alanine-leucine stretch, flanked by 1-4 tryptophan side chains, as well as the β-helical peptide dimer, gramicidin A. Membrane responses to mismatch include changes in local bilayer thickness and lipid order, varying systematically with peptide length. Adding more flanking tryptophan side chains led to an increase in bilayer thinning for negatively mismatched peptides, though it was also associated with a spreading of the bilayer interface. Peptide tilting, bending and stretching were systematic, with tilting dominating the responses, with values of up to ~45° for the most positively mismatched peptides. Peptide responses were modulated by the number of tryptophan side chains due to their anchoring roles and distributions around the helices. Potential of mean force calculations for local membrane thickness changes, helix tilting, bending and stretching revealed that membrane deformation is the least energetically costly of all mismatch responses, except for positively mismatched peptides where helix tilting also contributes substantially. This comparison of energetic driving forces of mismatch responses allows for deeper insight into protein stability and conformational changes in lipid membranes.

The role of tryptophan side chains in membrane protein anchoring and hydrophobic mismatch

Tryptophan (Trp) is abundant in membrane proteins, preferentially residing near the lipid-water interface where it is thought to play a significant anchoring role. Using a total of 3 μs of molecular dynamics simulations for a library of hydrophobic WALP-like peptides, a long poly-Leu α-helix, and the methyl-indole analog, we explore the thermodynamics of the Trp movement in membranes that governs the stability and orientation of transmembrane protein segments. We examine the dominant hydrogen-bonding interactions between the Trp and lipid carbonyl and phosphate moieties, cation-π interactions to lipid choline moieties, and elucidate the contributions to the thermodynamics that serve to localize the Trp, by ~4 kcal/mol, near the membrane glycerol backbone region. We show a striking similarity between the free energy to move an isolated Trp side chain to that found from a wide range of WALP peptides, suggesting that the location of this side chain is nearly independent of the host transmembrane segment. Our calculations provide quantitative measures that explain Trp's role as a modulator of responses to hydrophobic mismatch, providing a deeper understanding of how lipid composition may control a range of membrane active peptides and proteins.

Monte Carlo simulations of peptide-membrane interactions with the MCPep web server

The MCPep server (http://bental.tau.ac.il/MCPep/) is designed for non-experts wishing to perform Monte Carlo (MC) simulations of helical peptides in association with lipid membranes. MCPep is a web implementation of a previously developed MC simulation model. The model has been tested on a variety of peptides and protein fragments. The simulations successfully reproduced available empirical data and provided new molecular insights, such as the preferred locations of peptides in the membrane and the contribution of individual amino acids to membrane association. MCPep simulates the peptide in the aqueous phase and membrane environments, both described implicitly. In the former, the peptide is subjected solely to internal conformational changes, and in the latter, each MC cycle includes additional external rigid body rotational and translational motions to allow the peptide to change its location in the membrane. The server can explore the interaction of helical peptides of any amino-acid composition with membranes of various lipid compositions. Given the peptide’s sequence or structure and the natural width and surface charge of the membrane, MCPep reports the main determinants of peptide–membrane interactions, e.g. average location and orientation in the membrane, free energy of membrane association and the peptide’s helical content. Snapshots of example simulations are also provided.

The Transmembrane Helix Tilt May Be Determined by the Balance between Precession Entropy and Lipid Perturbation

Hydrophobic helical peptides interact with lipid bilayers in various modes, determined by the match between the length of the helix’s hydrophobic core and the thickness of the hydrocarbon region of the bilayer. For example, long helices may tilt with respect to the membrane normal to bury their hydrophobic cores in the membrane, and the lipid bilayer may stretch to match the helix length. Recent molecular dynamics simulations and potential of mean force calculations have shown that some TM helices whose lengths are equal to, or even shorter than, the bilayer thickness may also tilt. The tilt is driven by a gain in the helix precession entropy, which compensates for the free energy penalty resulting from membrane deformation. Using this free energy balance, we derived theoretically an equation of state, describing the dependence of the tilt on the helix length and membrane thickness. To this end, we conducted coarse-grained Monte Carlo simulations of the interaction of helices of various lengths with lipid bilayers of various thicknesses, reproducing and expanding the previous molecular dynamics simulations. Insight from the simulations facilitated the derivation of the theoretical model. The tilt angles calculated using the theoretical model agree well with our simulations and with previous calculations and measurements.

Influence of hydrophobic mismatch on structures and dynamics of gramicidin a and lipid bilayers

Gramicidin A (gA) is a 15-amino-acid antibiotic peptide with an alternating L-D sequence, which forms (dimeric) bilayer-spanning, monovalent cation channels in biological membranes and synthetic bilayers. We performed molecular dynamics simulations of gA dimers and monomers in all-atom, explicit dilauroylphosphatidylcholine (DLPC), dimyristoylphosphatidylcholine (DMPC), dioleoylphosphatidylcholine (DOPC), and 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) bilayers. The variation in acyl chain length among these different phospholipids provides a way to alter gA-bilayer interactions by varying the bilayer hydrophobic thickness, and to determine the influence of hydrophobic mismatch on the structure and dynamics of both gA channels (and monomeric subunits) and the host bilayers. The simulations show that the channel structure varied little with changes in hydrophobic mismatch, and that the lipid bilayer adapts to the bilayer-spanning channel to minimize the exposure of hydrophobic residues. The bilayer thickness, however, did not vary monotonically as a function of radial distance from the channel. In all simulations, there was an initial decrease in thickness within 4-5 Å from the channel, which was followed by an increase in DOPC and POPC or a further decrease in DLPC and DMPC bilayers. The bilayer thickness varied little in the monomer simulations-except one of three independent simulations for DMPC and all three DLPC simulations, where the bilayer thinned to allow a single subunit to form a bilayer-spanning water-permeable pore. The radial dependence of local lipid area and bilayer compressibility is also nonmonotonic in the first shell around gA dimers due to gA-phospholipid interactions and the hydrophobic mismatch. Order parameters, acyl chain dynamics, and diffusion constants also differ between the lipids in the first shell and the bulk. The lipid behaviors in the first shell around gA dimers are more complex than predicted from a simple mismatch model, which has implications for understanding the energetics of membrane protein-lipid interactions.

The residue composition of the aromatic anchor of the second transmembrane helix determines the signaling properties of the aspartate/maltose chemoreceptor Tar of Escherichia coli

Repositioning of the tandem aromatic residues (Trp-209 and Tyr-210) at the cytoplasmic end of the second transmembrane helix (TM2) modulates the signal output of the aspartate/maltose chemoreceptor of Escherichia coli (Tar(Ec)). Here, we directly assessed the effect of the residue composition of the aromatic anchor by studying the function of a library of Tar(Ec) variants that possess all possible combinations of Ala, Phe, Tyr, and Trp at positions 209 and 210. We identified three important properties of the aromatic anchor. First, a Trp residue at position 209 was required to maintain clockwise (CW) signal output in the absence of adaptive methylation, but adaptive methylation restored the ability of all of the mutant receptors to generate CW rotation. Second, when the aromatic anchor was replaced with tandem Ala residues, signaling was less compromised than when an Ala residue occupied position 209 and an aromatic residue occupied position 210. Finally, when Trp was present at position 209, the identity of the residue at position 210 had little effect on baseline signal output or aspartate chemotaxis, although maltose taxis was significantly affected by some substitutions at position 210. All of the mutant receptors we constructed supported some level of aspartate and maltose taxis in semisolid agar swim plates, but those without Trp at position 209 were overmethylated in their baseline signaling state. These results show the importance of the cytoplasmic aromatic anchor of TM2 in maintaining the baseline Tar(Ec) signal output and responsiveness to attractant signaling.

Probing the lipid-protein interface using model transmembrane peptides with a covalently linked acyl chain

The aim of this study was to gain insight into how interactions between proteins and lipids in membranes are sensed at the protein-lipid interface. As a probe to analyze this interface, we used deuterium-labeled acyl chains that were covalently linked to a model transmembrane peptide. First, a perdeuterated palmitoyl chain was coupled to the Trp-flanked peptide WALP23 (Ac-CGWW(LA)(8)LWWA-NH(2)), and the deuterium NMR spectrum was analyzed in di-C18:1-phosphatidylcholine (PC) bilayers. We found that the chain order of this peptide-linked chain is rather similar to that of a noncovalently coupled perdeuterated palmitoyl chain, except that it exhibits a slightly lower order. Similar results were obtained when site-specific deuterium labels were used and when the palmitoyl chain was attached to the more-hydrophobic model peptide WLP23 (Ac-CGWWL(17)WWA-NH(2)) or to the Lys-flanked peptide KALP23 (Ac-CGKK(LA)(8)LKKA-NH(2)). The experiments showed that the order of both the peptide-linked chains and the noncovalently coupled palmitoyl chains in the phospholipid bilayer increases in the order KALP23 < WALP23 < WLP23. Furthermore, changes in the bulk lipid bilayer thickness caused by varying the lipid composition from di-C14:1-PC to di-C18:1-PC or by including cholesterol were sensed rather similarly by the covalently coupled chain and the noncovalently coupled palmitoyl chains. The results indicate that the properties of lipids adjacent to transmembrane peptides mostly reflect the properties of the surrounding lipid bilayer, and hence that (at least for the single-span model peptides used in this study) annular lipids do not play a highly specific role in protein-lipid interactions.

Green proteorhodopsin reconstituted into nanoscale phospholipid bilayers (nanodiscs) as photoactive monomers

Over 4000 putative proteorhodopsins (PRs) have been identified throughout the oceans and seas of the Earth. The first of these eubacterial rhodopsins was discovered in 2000 and has expanded the family of microbial proton pumps to all three domains of life. With photophysical properties similar to those of bacteriorhodopsin, an archaeal proton pump, PRs are also generating interest for their potential use in various photonic applications. We perform here the first reconstitution of the minimal photoactive PR structure into nanoscale phospholipid bilayers (nanodiscs) to better understand how protein-protein and protein-lipid interactions influence the photophysical properties of PR. Spectral (steady-state and time-resolved UV-visible spectroscopy) and physical (size-exclusion chromatography and electron microscopy) characterization of these complexes confirms the preparation of a photoactive PR monomer within nanodiscs. Specifically, when embedded within a nanodisc, monomeric PR exhibits a titratable pK(a) (6.5-7.1) and photocycle lifetime (∼100-200 ms) that are comparable to the detergent-solubilized protein. These ndPRs also produce a photoactive blue-shifted absorbance, centered at 377 or 416 nm, that indicates that protein-protein interactions from a PR oligomer are required for a fast photocycle. Moreover, we demonstrate how these model membrane systems allow modulation of the PR photocycle by variation of the discoidal diameter (i.e., 10 or 12 nm), bilayer thickness (i.e., 23 or 26.5 Å), and degree of saturation of the lipid acyl chain. Nanodiscs also offer a highly stable environment of relevance to potential device applications.

Oligomerization state of photosynthetic core complexes is correlated with the dimerization affinity of a transmembrane helix

In the Rhodobacter (Rba.) species of photosynthetic purple bacteria, a single transmembrane α-helix, PufX, is found within the core complex, an essential photosynthetic macromolecular assembly that performs the absorption and the initial processing of light energy. Despite its structural simplicity, many unresolved questions surround PufX, the most important of which is its location within the photosynthetic core complex. One proposed placement of PufX is at the center of a core complex dimer, where two PufX helices associate in the membrane and form a homodimer. Inability for PufX of certain Rba. species to form a homodimer is thought to lead to monomeric core complexes. In the present study, we employ a combination of computational and experimental techniques to test the hypothesized homodimerization of PufX. We carry out a systematic investigation to measure the dimerization affinity of PufX from four Rba. species, Rba. blasticus , Rba. capsulatus , Rba. sphaeroides , and Rba. veldkampii , using a molecular dynamics-based free-energy method, as well as experimental TOXCAT assays. We found that the four PufX helices have substantially different dimerization affinities. Both computational and experimental techniques demonstrate that species with dimeric core complexes have PufX that can potentially form a homodimer, whereas the one species with monomeric core complexes has a PufX with little to no dimerization propensity. Our analysis of the helix-helix interface revealed a number of positions that may be important for PufX dimerization and the formation of a hydrogen-bond network between these GxxxG-containing helices. Our results suggest that the different oligomerization states of core complexes in various Rba. species can be attributed, among other factors, to the different propensity of its PufX helix to homodimerize.

Electric field driven changes of a gramicidin containing lipid bilayer supported on a Au(111) surface

Langmuir-Blodgett and Langmuir-Schaeffer methods were employed to deposit a mixed bilayer consisting of 90% of 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 10% of gramicidin (GD), a short 15 residue ion channel forming peptide, onto a Au(111) electrode surface. This architecture allowed us to investigate the effect of the electrostatic potential applied to the electrode on the orientation and conformation of DMPC molecules in the bilayer containing the ion channel. The charge density data were determined from chronocoulometry experiments. The electric field and the potential across the membrane were determined through the use of charge density curves. The magnitudes of potentials across the gold-supported biomimetic membrane were comparable to the transmembrane potential acting on a natural membrane. The information regarding the orientation and conformation of DMPC and GD molecules in the bilayer was obtained from photon polarization modulation infrared reflection absorption spectroscopy (PMIRRAS) measurements. The results show that the bilayer is adsorbed, in direct contact with the metal surface, when the potential across the interface is more positive than -0.4 V and is lifted from the gold surface when the potential across the interface is more negative than -0.4 V. This change in the state of the bilayer has a significant impact on the orientation and conformation of the phospholipid and gramicidin molecules. The potential induced changes in the membrane containing peptide were compared to the changes in the structure of the pure DMPC bilayer determined in earlier studies.

Response of GWALP transmembrane peptides to changes in the tryptophan anchor positions

While the interfacial partitioning of charged or aromatic anchor residues may determine the preferred orientations of transmembrane peptide helices, the dependence of helix orientation on anchor residue position is not well understood. When anchor residue locations are changed systematically, some adaptations of the peptide-lipid interactions may be required to compensate for the altered interfacial interactions. Recently, we have developed a novel transmembrane peptide, termed GW(5,19)ALP23 (acetyl-GGALW(5)LALALALALALALW(19)LAGA-ethanolamide), which proves to be a well-behaved sequence for an orderly investigation of protein-lipid interactions. Its roughly symmetric nature allows for shifting the anchoring Trp residues by one Leu-Ala pair inward (GW(7,17)ALP23) or outward (GW(3,21)ALP23), thus providing fine adjustments of the formal distance between the tryptophan residues. With no other obvious anchoring features present, we postulate that the inter-Trp distance may be crucial for aspects of the peptide-lipid interaction. Importantly, the amino acid composition is identical for each of the resulting related GWALP23 sequences, and the radial separation between the pairs of Trp residues on each side of the transmembrane α-helix remains similar. Here we address the adaptation of the aforementioned peptides to the varying Trp locations by means of solid-state (2)H nuclear magnetic resonance experiments in varying lipid bilayer membrane environments. All of the GW(x,y)ALP23 sequence isomers adopt transmembrane orientations in DOPC, DMPC, and DLPC environments, even when the Trp residues are quite closely spaced, in GW(7,17)ALP23. Furthermore, the dynamics for each peptide isomer are less extensive than for peptides possessing additional interfacial Trp residues. The helical secondary structure is maintained more strongly within the Trp-flanked core region than outside of the Trp boundaries. Deuterium-labeled tryptophan indole rings in the GW(x,y)ALP23 peptides provide additional insights into the behavior of the Trp side chains. A Trp side chain near the C-terminus adopts a different orientation and undergoes somewhat faster dynamics than a corresponding Trp side chain located an equivalent distance from the N-terminus. In contrast, as the inter-Trp distance changes, the variations among the average orientations of the Trp indole rings at either terminus are systematic yet fairly small. We conclude that subtle adjustments to the peptide tilt, and to the N- and C-terminal Trp side chain torsion angles, permit the GW(x,y)ALP23 peptides to maintain preferred transmembrane orientations while adapting to lipid bilayers with differing hydrophobic thicknesses.

Exploring peptide-membrane interactions with coarse-grained MD simulations

The interaction of α-helical peptides with lipid bilayers is central to our understanding of the physicochemical principles of biological membrane organization and stability. Mutations that alter the position or orientation of an α-helix within a membrane, or that change the probability that the α-helix will insert into the membrane, can alter a range of membrane protein functions. We describe a comparative coarse-grained molecular dynamics simulation methodology, based on self-assembly of a lipid bilayer in the presence of an α-helical peptide, which allows us to model membrane transmembrane helix insertion. We validate this methodology against available experimental data for synthetic model peptides (WALP23 and LS3). Simulation-based estimates of apparent free energies of insertion into a bilayer of cystic fibrosis transmembrane regulator-derived helices correlate well with published data for translocon-mediated insertion. Comparison of values of the apparent free energy of insertion from self-assembly simulations with those from coarse-grained molecular dynamics potentials of mean force for model peptides, and with translocon-mediated insertion of cystic fibrosis transmembrane regulator-derived peptides suggests a nonequilibrium model of helix insertion into bilayers.

An AFM study of solid-phase bilayers of unsaturated PC lipids and the lateral distribution of the transmembrane model peptide WALP23 in these bilayers

An altered lipid packing can have a large influence on the properties of the membrane and the lateral distribution of proteins and/or peptides that are associated with the bilayer. Here, it is shown by contact-mode atomic force microscopy that the surface topography of solid-phase bilayers of PC lipids with an unsaturated cis bond in their acyl chains shows surfaces with a large number of line-type packing defects, in contrast to the much smoother surfaces observed for saturated PC lipids. Di-n:1-PC (n = 20, 22, 24) and (16:0,18:1)-PC (POPC) were used. Next, the influence of an altered lipid environment on the lateral distribution of the single α-helical model peptide WALP23 was studied by incorporating the peptide in the bilayers of di-n:1-PC (n = 20, 22, 24) and (16:0,18:1)-PC unsaturated lipids. The presence of WALP23 leads to an increase in the number of packing defects but does not lead to the formation of the striated domains that were previously observed in bilayers of saturated PC lipids and WALP. This is ascribed to the less efficient lateral lipid packing of the unsaturated lipids, while the increase in packing defects is probably an indirect effect of the peptide. Finally, the fact that an altered lipid packing affects the distribution of WALP23 is also confirmed in an additional experiment where the solvent TFE (2,2,2-trifluorethanol) is added to bilayers of di-16:0-PC/WALP23. At 3.5 vol% TFE, the previous striated ordering of the peptide is abolished and replaced by loose lines.

Peptide model helices in lipid membranes: insertion, positioning, and lipid response on aggregation studied by X-ray scattering

Studying membrane active peptides or protein fragments within the lipid bilayer environment is particularly challenging in the case of synthetically modified, labeled, artificial, or recently discovered native structures. For such samples the localization and orientation of the molecular species or probe within the lipid bilayer environment is the focus of research prior to an evaluation of their dynamic or mechanistic behavior. X-ray scattering is a powerful method to study peptide/lipid interactions in the fluid, fully hydrated state of a lipid bilayer. For one, the lipid response can be revealed by observing membrane thickening and thinning as well as packing in the membrane plane; at the same time, the distinct positions of peptide moieties within lipid membranes can be elucidated at resolutions of up to several angstroms by applying heavy-atom labeling techniques. In this study, we describe a generally applicable X-ray scattering approach that provides robust and quantitative information about peptide insertion and localization as well as peptide/lipid interaction within highly oriented, hydrated multilamellar membrane stacks. To this end, we have studied an artificial, designed β-helical peptide motif in its homodimeric and hairpin variants adopting different states of oligomerization. These peptide lipid complexes were analyzed by grazing incidence diffraction (GID) to monitor changes in the lateral lipid packing and ordering. In addition, we have applied anomalous reflectivity using synchrotron radiation as well as in-house X-ray reflectivity in combination with iodine-labeling in order to determine the electron density distribution ρ(z) along the membrane normal (z axis), and thereby reveal the hydrophobic mismatch situation as well as the position of certain amino acid side chains within the lipid bilayer. In the case of multiple labeling, the latter technique is not only applicable to demonstrate the peptide's reconstitution but also to generate evidence about the relative peptide orientation with respect to the lipid bilayer.

Membrane anchoring and interaction between transmembrane domains are crucial for K+ channel function

The small viral channel Kcv is a Kir-like K(+) channel of only 94 amino acids. With this simple structure, the tetramer of Kcv represents the pore module of all complex K(+) channels. To examine the structural contribution of the transmembrane domains (TMDs) to channel function, we performed Ala scanning mutagenesis of the two domains and tested the functionality of the mutants in a yeast complementation assay. The data reveal, in combination with computational models, that the upper halves of both TMDs, which face toward the external medium, are rather rigid, whereas the inner parts are more flexible. The rigidity of the outer TMD is conferred by a number of essential aromatic amino acids that face the membrane and probably anchor this domain in the bilayer. The inner TMD is intimately connected with the rigid part of the outer TMD via π···π interactions between a pair of aromatic amino acids. This structural principle is conserved within the viral K(+) channels and also present in Kir2.2, implying a general importance of this architecture for K(+) channel function.

Assembly of the m2 tetramer is strongly modulated by lipid chain length

The influenza virus matrix protein 2 (M2) assembles into a tetramer in the host membrane during viral uncoating and maturation. It has been used as a model system to understand the relative contributions of protein-lipid and protein-protein interactions to membrane protein structure and association. Here we investigate the effect of lipid chain length on the association of the M2 transmembrane domain into tetramers using Förster resonance energy transfer. We observe that the interactions between the M2 helices are much stronger in 1,2-dilauroyl-sn-glycero-3-phosphocholine than in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine bilayers. Thus, lipid chain length and bilayer thickness not only modulate peptide interactions, but could also be a major determinant of the association of transmembrane helices into functional membrane protein oligomers.

Lipid packing drives the segregation of transmembrane helices into disordered lipid domains in model membranes

Cell membranes are comprised of multicomponent lipid and protein mixtures that exhibit a complex partitioning behavior. Regions of structural and compositional heterogeneity play a major role in the sorting and self-assembly of proteins, and their clustering into higher-order oligomers. Here, we use computer simulations and optical microscopy to study the sorting of transmembrane helices into the liquid-disordered domains of phase-separated model membranes, irrespective of peptide-lipid hydrophobic mismatch. Free energy calculations show that the enthalpic contribution due to the packing of the lipids drives the lateral sorting of the helices. Hydrophobic mismatch regulates the clustering into either small dynamic or large static aggregates. These results reveal important molecular driving forces for the lateral organization and self-assembly of transmembrane helices in heterogeneous model membranes, with implications for the formation of functional protein complexes in real cells.

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-mus 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 approximately 10 degrees 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.

Transmembrane peptides influence the affinity of sterols for phospholipid bilayers

Cholesterol is distributed unevenly between different cellular membrane compartments, and the cholesterol content increases from the inner bilayers toward the plasma membrane. It has been suggested that this cholesterol gradient is important in the sorting of transmembrane proteins. Cholesterol has also been to shown play an important role in lateral organization of eukaryotic cell membranes. In this study the aim was to determine how transmembrane proteins influence the lateral distribution of cholesterol in phospholipid bilayers. Insight into this can be obtained by studying how cholesterol interacts with bilayer membranes of different composition in the presence of designed peptides that mimic the transmembrane helices of proteins. For this purpose we developed an assay in which the partitioning of the fluorescent cholesterol analog CTL between LUVs and mbetaCD can be measured. Comparison of how cholesterol and CTL partitioning between mbetaCD and phospholipid bilayers with different composition suggests that CTL sensed changes in bilayer composition similarly as cholesterol. Therefore, the results obtained with CTL can be used to understand cholesterol distribution in lipid bilayers. The effect of WALP23 on CTL partitioning between DMPC bilayers and mbetaCD was measured. From the results it was clear that WALP23 increased both the order in the bilayers (as seen from CTL and DPH anisotropy) and the affinity of the sterol for the bilayer in a concentration dependent way. Although WALP23 also increased the order in DLPC and POPC bilayers the effects on CTL partitioning was much smaller with these lipids. This indicates that proteins have the largest effect on sterol interactions with phospholipids that have longer and saturated acyl chains. KALP23 did not significantly affect the acyl chain order in the phospholipid bilayers, and inclusion of KALP23 into DMPC bilayers slightly decreased CTL partitioning into the bilayer. This shows that transmembrane proteins can both decrease and increase the affinity of sterols for the lipid bilayers surrounding proteins. This is likely to affect the sterol distribution within the bilayer and thereby the lateral organization in biomembranes.

Molecular recognition at the membrane-water interface: controlling integral peptide helices by off-membrane nucleobase pairing

The aggregation and organization of membrane proteins and transmembrane peptides is related to the interacting molecular species itself and strongly depends on the lipid environment. Because of the complexity and dynamics of these interactions, they are often hardly traceable and nearly impossible to predict. For this reason, peptide model systems are a valuable tool in studying membrane associated processes since they are synthetically accessible and can be readily modified. To control and study the aggregation of peptide transmembrane domains (TMDs) the interacting interfaces of the TMDs themselves can be altered. A second less extensively studied approach targets the TMD assembly by using interaction and recognition of domains at the membrane outside as frequently found in the membrane protein interplay and protein assembly. In the present study, double helical transmembrane domains were designed and synthesized on the basis of a recently reported d,l-alternating peptide pore motif derived from gramicidin A. The highly hydrophobic and aromatic transmembrane peptide was covalently functionalized with a short peptide nucleic acid (PNA) used as specific outer-membrane recognition unit. The PNA sequences were chosen with high polarity to ensure localization within the aqueous phase. To estimate the impact of the membrane adjacent recognition on the TMD assembly by Förster resonance energy transfer (FRET), fluorescence probes were covalently attached to the side chains of the membrane spanning peptide helices. Dimerization of the TMD-peptide/PNA conjugates within unilamellar lipid vesicles was observed. The dimer/monomer ratio of TMDs can be controlled by temperature variation.

Interpretation of 2H-NMR experiments on the orientation of the transmembrane helix WALP23 by computer simulations

Orientation, dynamics, and packing of transmembrane helical peptides are important determinants of membrane protein structure, dynamics, and function. Because it is difficult to investigate these aspects by studying real membrane proteins, model transmembrane helical peptides are widely used. NMR experiments provide information on both orientation and dynamics of peptides, but they require that motional models be interpreted. Different motional models yield different interpretations of quadrupolar splittings (QS) in terms of helix orientation and dynamics. Here, we use coarse-grained (CG) molecular dynamics (MD) simulations to investigate the behavior of a well-known model transmembrane peptide, WALP23, under different hydrophobic matching/mismatching conditions. We compare experimental (2)H-NMR QS (directly measured in experiments), as well as helix tilt angle and azimuthal rotation (not directly measured), with CG MD simulation results. For QS, the agreement is significantly better than previously obtained with atomistic simulations, indicating that equilibrium sampling is more important than atomistic details for reproducing experimental QS. Calculations of helix orientation confirm that the interpretation of QS depends on the motional model used. Our simulations suggest that WALP23 can form dimers, which are more stable in an antiparallel arrangement. The origin of the preference for the antiparallel orientation lies not only in electrostatic interactions but also in better surface complementarity. In most cases, a mixture of monomers and antiparallel dimers provides better agreement with NMR data compared to the monomer and the parallel dimer. CG MD simulations allow predictions of helix orientation and dynamics and interpretation of QS data without requiring any assumption about the motional model.

Variation of the lateral mobility of transmembrane peptides with hydrophobic mismatch

A hydrophobic mismatch between protein length and membrane thickness can lead to a modification of protein conformation, function, and oligomerization. To study the role of hydrophobic mismatch, we have measured the change in mobility of transmembrane peptides possessing a hydrophobic helix of various length d(pi) in lipid membranes of giant vesicles. We also used a model system where the hydrophobic thickness of the bilayers, h, can be tuned at will. We precisely measured the diffusion coefficient of the embedded peptides and gained access to the apparent size of diffusing objects. For bilayers thinner than d(pi), the diffusion coefficient decreases, and the derived characteristic sizes are larger than the peptide radii. Previous studies suggest that peptides accommodate by tilting. This scenario was confirmed by ATR-FTIR spectroscopy. As the membrane thickness increases, the value of the diffusion coefficient increases to reach a maximum at h approximately = d(pi). We show that this variation in diffusion coefficient is consistent with a decrease in peptide tilt. To do so, we have derived a relation between the diffusion coefficient and the tilt angle, and we used this relation to derive the peptide tilt from our diffusion measurements. As the membrane thickness increases, the peptides raise (i.e., their tilt is reduced) and reach an upright position and a maximal mobility for h approximately = d(pi). Using accessibility measurements, we show that when the membrane becomes too thick, the peptide polar heads sink into the interfacial region. Surprisingly, this "pinching" behavior does not hinder the lateral diffusion of the transmembrane peptides. Ultimately, a break in the peptide transmembrane anchorage is observed and is revealed by a "jump" in the D values.

Reconstruction of atomistic details from coarse-grained structures

We present an algorithm to reconstruct atomistic structures from their corresponding coarse-grained (CG) representations and its implementation into the freely available molecular dynamics (MD) program package GROMACS. The central part of the algorithm is a simulated annealing MD simulation in which the CG and atomistic structures are coupled via restraints. A number of examples demonstrate the application of the reconstruction procedure to obtain low-energy atomistic structural ensembles from their CG counterparts. We reconstructed individual molecules in vacuo (NCQ tripeptide, dipalmitoylphosphatidylcholine, and cholesterol), bulk water, and a WALP transmembrane peptide embedded in a solvated lipid bilayer. The first examples serve to optimize the parameters for the reconstruction procedure, whereas the latter examples illustrate the applicability to condensed-phase biomolecular systems.

Changes in transmembrane helix alignment by arginine residues revealed by solid-state NMR experiments and coarse-grained MD simulations

Independent experimental and computational approaches show agreement concerning arginine/membrane interactions when a single arginine is introduced at selected positions within the membrane-spanning region of acetyl-GGALW(5)LALALAL(12)AL(14)ALALW(19)LAGA-ethanolamide, designated GWALP23. Peptide sequence isomers having Arg in position 12 or position 14 display markedly different behaviors, as deduced by both solid-state NMR experiments and coarse-grained molecular dynamics (CG-MD) simulations. With respect to the membrane normal of DOPC or DPPC lipid bilayer membranes, GWALP23-R14 shows one major state whose apparent average tilt is approximately 10 degrees greater than that of GWALP23. The presence of R14 furthermore induces bilayer thinning and peptide displacement to "lift" the charged guanidinium toward the bilayer surface. By contrast, GWALP23-R12 exhibits multiple states that are in slow exchange on the NMR time scale, with CG-MD simulations indicating two distinct positions with different screw rotation angles in the membrane, along with an increased tendency to exit the lipid bilayer.

Structural features of fusogenic model transmembrane domains that differentially regulate inner and outer leaflet mixing in membrane fusion

The transmembrane domains of fusion proteins are known to be important for their fusogenic activity. In an effort to systematically investigate the structure/function relationships of transmembrane domains we had previously designed LV-peptides that mimic natural fusion protein TMDs in their ability to drive fusion after incorporation into liposomal membranes. Here, we investigate the impact of different structural features of LV-peptide TMDs on inner and outer leaflet mixing. We find that fusion driven by the helical peptides involves a hemifusion intermediate as previously seen for natural fusion proteins. Helix backbone dynamics enhances fusion by selectively promoting outer leaflet mixing. Furthermore, the hydrophobic length of the peptides as well as covalent attachment of long acyl chains affects outer and inner leaflet mixing to different extents. Different structural features of transmembrane domains thus appear to differentially influence the rearrangements of lipids in fusion initiation and the hemifusion-to-fusion transition. The relevance of these findings in respect to the function of natural fusion proteins is discussed.

Characterization of membrane-protein interactions for the leucine transporter from Aquifex aeolicus by molecular dynamics calculations

Multinanosecond molecular dynamics (MD) simulations have been employed to characterize the interaction of an integral membrane protein (IMP), the leucine transmitter from Aquifex aeolicus (Yamashita et al., Nature 2005, 437, 215-223), with hydrated lipid bilayer membranes in their physiologically relevant liquid crystalline phases. Analysis of the MD trajectories for dimyristoyl phosphatidylcholine (DMPC), 1-palmitoyl-2-oleoyl phosphatidylcholine (POPC), and 1-palmitoyl-2-oleoyl phosphatidylethanolamine (POPE) focused on the contacts between aromatic and basic side chains of the IMP with the lipid head groups and water. Structural fluctuations of the IMP were investigated as well as the contact dynamics of neighboring lipids. In characterizing the IMP-membrane systems, the behaviors of the protein's cytoplasmic and periplasmic parts are considered separately. All three lipid membranes show a rather similar overall level of association with the IMP. However, for DMPC there is a better matching of the membrane core to the hydrophobic transmembrane portion of the IMP. The closed cytoplasmic end of the IMP exhibits a higher degree of association with lipids than the more open periplasmic end, an observation which correlates with the more compact structure and a slower dynamics of surrounding lipids.

A glycophorin A-like framework for the dimerization of photosynthetic core complexes

The core complex in photosynthetic bacteria plays a central role in photosynthesis. This molecular assembly is composed of two protein complexes, viz., the light-harvesting complex I (LH1), which absorbs sunlight by means of the protein-bound bacteriochlorophylls, and the reaction center (RC), which uses the light-excitation energy absorbed by the LH complexes to produce a transmembrane (TM) charge gradient, subsequently employed for energy conversion. In Rhodobacter (Rba.) sphaeroides, the core complex contains, in addition, two copies of the single TM alpha-helix protein, PufX, and forms a (RC-LH1-PufX)(2) dimer. To this date, no high-resolution structure has been reported for the entire core complex. In particular, the location of PufX within the (RC-LH1-PufX)(2) dimer is still the subject of much debate. Here, one of the proposed locations for PufX, requiring its dimerization, is examined. The PufX-dimer model on the basis of the glycophorin A (GpA) dimer was constructed, and its robustness was probed through a series of molecular dynamics (MD) simulations. The free-energy change due to the replacement of Gly35 by valine was also determined to assess whether this mutation is responsible for distinct PufX oligomerization states in different Rba. species. The present study shows that PufX helices form a stable GpA-like dimer with a helix-helix crossing angle that could constitute the molecular basis of the reported highly bent and V-shaped structure of the Rba. sphaeroides core complex dimer.

Dynamic transitions of membrane-active peptides

Membrane-active peptides or protein segments play an important role in many biological processes at the cellular interface to the environment. They are involved, e.g., in cellular fusion or host defense, where they can cause not only merging but also the destabilization of cell membranes. Many factors determine how these typically amphipathic peptides interact with the lipid bilayer. For example, the peptide orientation in the membrane determines which parts of the peptide are exposed to the hydrophobic bilayer interior or to the polar lipid/water interface. As another example, oligomerization is required for many activities such as pore formation. Peptides have been often classified according to a single characteristic mode of interaction with the bilayer, but over the years a more versatile picture has emerged. It appears that any single peptide can adopt several different alignments and/or oligomeric states in response to changes in the environment. For instance, many antimicrobial peptides adopt a surface-parallel alignment at low concentration, but they tilt obliquely into or even fully insert transmembrane into the bilayer above a critical peptide-to-lipid ratio, often in the form of oligomeric pores. Similar changes in peptide orientation or oligomeric state have been observed as a function of, e.g., temperature, lipid composition, pH, or induced by a synergistic partner peptide. Such transitions between peptide states can be regarded as the result of a re-adjustment in the balance between peptide-peptide and peptide-lipid interactions, as the environment conditions are changed. Though often studied in model membrane systems, such rich variety of peptide states is even more likely to occur in native biomembranes with their diverse compositions and physicochemical properties. The ability to undergo transitions between different states thus plays a fundamental role for the biological activities of membrane-active peptides.

An intramembrane aromatic network determines pentameric assembly of Cys-loop receptors

Cys-loop receptors are pentameric ligand-gated ion channels (pLGICs) that mediate fast synaptic transmission. Here functional pentameric assembly of truncated fragments comprising the ligand-binding N-terminal ectodomains and the first three transmembrane helices, M1-M3, of both the inhibitory glycine receptor (GlyR) alpha1 and the 5HT(3)A receptor subunits was found to be rescued by coexpressing the complementary fourth transmembrane helix, M4. Alanine scanning identified multiple aromatic residues in M1, M3 and M4 as key determinants of GlyR assembly. Homology modeling and molecular dynamics simulations revealed that these residues define an interhelical aromatic network, which we propose determines the geometry of M1-M4 tetrahelical packing such that nascent pLGIC subunits must adopt a closed fivefold symmetry. Because pLGIC ectodomains form random nonstoichiometric oligomers, proper pentameric assembly apparently depends on intersubunit interactions between extracellular domains and intrasubunit interactions between transmembrane segments.

Lipid rafts as functional heterogeneity in cell membranes

Biological membranes are not structurally passive solvents of amphipathic proteins and lipids. Rather, it appears their constituents have evolved intrinsic characteristics that make homogeneous distribution of components unlikely. As a case in point, the concept of lipid rafts has received considerable attention from biologists and biophysicists since the formalization of the hypothesis more than 10 years ago. Today, it is clear that sphingolipid and cholesterol can self-associate into micron-scaled phases in model membranes and that these lipids are involved in the formation of highly dynamic nanoscale heterogeneity in the plasma membrane of living cells. However, it remains unclear whether these entities are manifestations of the same principle. A powerful means by which the molecular organization of rafts can be assessed is through analysis of their functionalized condition. Raft heterogeneity can be activated to coalesce and laterally reorganize/stabilize bioactivity in cell membranes. Evaluation of this property suggests that functional raft heterogeneity arises through principles of lipid-driven phase segregation coupled to additional chemical specificities, probably involving proteins.

Relevance of fatty acid covalently bound to Escherichia coli alpha-hemolysin and membrane microdomains in the oligomerization process

alpha-Hemolysin (HlyA) is an exotoxin secreted by some pathogenic strains of Escherichia coli that causes lysis of several mammalian cells, including erythrocytes of different species. HlyA is synthesized as a protoxin, pro-HlyA, which is activated by acylation at two internal lysines Lys-563 and Lys-689. It has been proposed that pore formation is the mechanism of cytolytic activity for this toxin, as shown in experiments with whole cells, planar lipid membranes, and liposomes, but these experiments have yielded conflicting results about the structure of the pore. In this study, HlyA cysteine replacement mutant proteins of amino acids have been labeled with Alexa-488 and Alexa-546. Fluorescence resonance energy transfer measurements, employing labeled toxin bound to sheep ghost erythrocytes, have demonstrated that HlyA oligomerizes on erythrocyte membranes. As the cytotoxic activity is absolutely dependent on acylation, we have studied the role of acylation in the oligomerization, demonstrating that fatty acids are essential in this process. On the other hand, fluorescence resonance energy transfer and the hemolytic activity decrease when the erythrocyte ghosts are cholesterol-depleted, hence indicating the role of membrane microdomains in the clustering of HlyA. Simultaneously, HlyA was found in detergent-resistant membranes. Pro-HlyA has also been found in detergent-resistant membranes, thus demonstrating that the importance of acyl chains in toxin oligomerization is the promotion of protein-protein interaction. These results change the concept of the main role assigned to acyl chain in the targeting of proteins to membrane microdomains.

Protein-lipid interactions of bacteriophage M13 gene 9 minor coat protein (Review)

Gene 9 protein is one of the minor coat proteins of bacteriophage M13. The protein plays a role in the assembly process by associating with the host membrane by protein-lipid interactions. The availability of chemically synthesized protein has enabled the biophysical characterization of the membrane-bound state of the protein by using model membrane systems. This paper summarizes, discusses and further interprets this work in the light of the current state of the literature, leading to new possible models of the coat protein in a membrane. The biological implications of these findings related to the membrane-bound phage assembly are indicated.