Hydrophobic mismatch between proteins and lipids in membranes

Expression, purification, and reconstitution of the voltage-sensing domain from Ci-VSP

The voltage-sensing domain (VSD) is the common scaffold responsible for the functional behavior of voltage-gated ion channels, voltage sensitive enzymes, and proton channels. Because of the position of the voltage dependence of the available VSD structures, at present, they all represent the activated state of the sensor. Yet in the absence of a consensus resting state structure, the mechanistic details of voltage sensing remain controversial. The voltage dependence of the VSD from Ci-VSP (Ci-VSD) is dramatically right shifted, so that at 0 mV it presumably populates the putative resting state. Appropriate biochemical methods are an essential prerequisite for generating sufficient amounts of Ci-VSD protein for high-resolution structural studies. Here, we present a simple and robust protocol for the expression of eukaryotic Ci-VSD in Escherichia coli at milligram levels. The protein is pure, homogeneous, monodisperse, and well-folded after solubilization in Anzergent 3-14 at the analyzed concentration (~0.3 mg/mL). Ci-VSD can be reconstituted into liposomes of various compositions, and initial site-directed spin labeling and electron paramagnetic resonance (EPR) spectroscopic measurements indicate its first transmembrane segment folds into an α-helix, in agreement with the homologous region of other VSDs. On the basis of our results and enhanced relaxation EPR spectroscopy measurement, Ci-VSD reconstitutes essentially randomly in proteoliposomes, precluding straightforward application of transmembrane voltages in combination with spectroscopic methods. Nevertheless, these results represent an initial step that makes the resting state of a VSD accessible to a variety of biophysical and structural approaches, including X-ray crystallography, spectroscopic methods, and electrophysiology in lipid bilayers.

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.

Aurein 2.3 functionality is supported by oblique orientated α-helical formation

In this study, an amphibian antimicrobial peptide, aurein 2.3, was predicted to use oblique orientated α-helix formation in its mechanism of membrane destabilisation. Molecular dynamic (MD) simulations and circular dichroism (CD) experimental data suggested that aurein 2.3 exists in solution as unstructured monomers and folds to form predominantly α-helical structures in the presence of a dimyristoylphosphatidylcholine membrane. MD showed that the peptide was highly surface active, which supported monolayer data where the peptide induced surface pressure changes>34 mNm(-1). In the presence of a lipid membrane MD simulations suggested that under hydrophobic mismatch the peptide is seen to insert via oblique orientation with a phenylalanine residue (PHE3) playing a key role in the membrane interaction. There is evidence of snorkelling leucine residues leading to further membrane disruption and supporting the high level of lysis observed using calcein release assays (76%). Simulations performed at higher peptide/lipid ratio show peptide cooperativity is key to increased efficiency leading to pore-formation.

An Extension and Further Validation of an All-Atomistic Force Field for Biological Membranes

Biological membranes are versatile in composition and host intriguing molecular processes. In order to be able to study these systems, an accurate model Hamiltonian or force field (FF) is a necessity. Here, we report the results of our extension of earlier developed all-atomistic FF parameters for fully saturated phospholipids that complements an earlier parameter set for saturated phosphatidylcholine lipids (J. Phys. Chem. B, 2012, 116, 3164-3179). The FF, coined Slipids (Stockholm lipids), now also includes parameters for unsaturated phosphatidylcholine and phosphatidylethanolamine lipids, e.g., POPC, DOPC, SOPC, POPE, and DOPE. As the extended set of parameters is derived with the same philosophy as previously applied, the resulting FF has been developed in a fully consistent manner. The capabilities of Slipids are demonstrated by performing long simulations without applying any surface tension and using the correct isothermal-isobaric (NPT) ensemble for a range of temperatures and carefully comparing a number of properties with experimental findings. Results show that several structural properties are very well reproduced, such as scattering form factors, NMR order parameters, thicknesses, and area per lipid. Thermal dependencies of different thicknesses and area per lipid are reproduced as well. Lipid diffusion is systematically slightly underestimated, whereas the normalized lipid diffusion follows the experimental trends. This is believed to be due to the lack of collective movement in the relatively small bilayer patches used. Furthermore, the compatibility with amino acid FFs from the AMBER family is tested in explicit transmembrane complexes of the WALP23 peptide with DLPC and DOPC bilayers, and this shows that Slipids can be used to study more complex and biologically relevant systems.

Multiple Pathway Relaxation Enhancement in the System Composed of Three Paramagnetic Species: Nitroxide Radical-Ln(3+)-O2

Longitudinal relaxation of nitroxide spin-labels has been measured for a membrane-incorporated α-helical polypeptide in the presence and absence of residual amounts of membrane-dissolved O2 and paramagnetic Dy(3+) ions. Such a model system, containing three different types of paramagnetic species, provides an important example of nonadditivity of two different relaxation channels for the nitroxide spins.

Wax-tear and meibum protein, wax-β-carotene interactions in vitro using infrared spectroscopy

Protein-meibum and terpenoids-meibum lipid interactions could be important in the etiology of meibomian gland dysfunction (MGD) and dry eye symptoms. In the current model studies, attenuated total reflectance (ATR) infrared (IR) spectroscopy was used to determine if the terpenoid β-carotene and the major proteins in tears and meibum affect the hydrocarbon chain conformation and carbonyl environment of wax, an abundant component of meibum. The main finding of these studies is that mucin binding to wax disordered slightly the conformation of the hydrocarbon chains of wax and caused the wax carbonyls to become hydrogen bonded or experience a more hydrophilic environment. Lysozyme and lactoglobulin, two proteins shown to bind to monolayers of meibum, did not have such an effect. Keratin and β-carotene did not affect the fluidity (viscosity) or environment of the carbonyl moieties of wax. Based on these results, tetraterpenoids are not likely to influence the structure of meibum in the meibomian glands. In addition, these findings suggest that it is unlikely that keratin blocks meibomian glands by causing the meibum to become more viscous. Among the tear fluid proteins studied, mucin is the most likely to influence the conformation and carbonyl environment of meibum at the tear film surface.

Fluid membranes can drive linear aggregation of adsorbed spherical nanoparticles

Using computer simulations, we show that lipid membranes can mediate linear aggregation of spherical nanoparticles binding to it for a wide range of biologically relevant bending rigidities. This result is in net contrast with the isotropic aggregation of nanoparticles on fluid interfaces or the expected clustering of isotropic insertions in biological membranes. We present a phase diagram indicating where linear aggregation is expected and compute explicitly the free-energy barriers associated with linear and isotropic aggregation. Finally, we provide simple scaling arguments to explain this phenomenology.

Coupling exo- and endocytosis: an essential role for PIP₂ at the synapse

Chemical synapses are specialist points of contact between two neurons, where information transfer takes place. Communication occurs through the release of neurotransmitter substances from small synaptic vesicles in the presynaptic terminal, which fuse with the presynaptic plasma membrane in response to neuronal stimulation. However, as neurons in the central nervous system typically only possess ~200 vesicles, high levels of release would quickly lead to a depletion in the number of vesicles, as well as leading to an increase in the area of the presynaptic plasma membrane (and possible misalignment with postsynaptic structures). Hence, synaptic vesicle fusion is tightly coupled to a local recycling of synaptic vesicles. For a long time, however, the exact molecular mechanisms coupling fusion and subsequent recycling remained unclear. Recent work now indicates a unique role for the plasma membrane lipid phosphatidylinositol 4,5-bisphosphate (PIP(2)), acting together with the vesicular protein synaptotagmin, in coupling these two processes. In this work, we review the evidence for such a mechanism and discuss both the possible advantages and disadvantages for vesicle recycling (and hence signal transduction) in the nervous system. This article is part of a Special Issue entitled Lipids and Vesicular Transport.

Role of unsaturated lipid and ergosterol in ethanol tolerance of model yeast biomembranes

We present a combined atomic force microscopy and fluorescence microscopy study of the behavior of a ternary supported lipid bilayer system containing a saturated lipid (DPPC), an unsaturated lipid (DOPC), and ergosterol in the presence of high ethanol (20 vol %). We find that the fluorescent probe Texas Red DHPE preferentially partitions into the ethanol-induced interdigitated phase, which allows the use of fluorescence imaging to investigate the phase behavior of the system. Atomic force microscopy and fluorescence images of samples with the same lipid mixture show good agreement in sample morphology and area fractions of the observed phases. Using area fractions obtained from fluorescence images over a broad range of compositions, we constructed a phase diagram of the DPPC/DOPC/ergosterol system at 20 vol % ethanol. The phase diagram clearly shows that increasing unsaturated lipid and/or ergosterol protects the membrane by preventing the formation of the interdigitated phase. This result supports the hypothesis that yeast cells increase ergosterol and unsaturated lipid content to prevent interdigitation and maintain an optimal membrane thickness as ethanol concentration increases during anaerobic fermentations. Changes in plasma membrane composition provide an important survival factor for yeast cells to deter ethanol toxicity.

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.

Simulation of fusion-mediated nanoemulsion interactions with model lipid bilayers

Perfluorocarbon-based nanoemulsion particles have become promising platforms for the delivery of therapeutic and diagnostic agents to specific target cells in a non-invasive manner. A "contact-facilitated" delivery mechanism has been proposed wherein the emulsifying phospholipid monolayer on the nanoemulsion surface contacts and forms a lipid complex with the outer monolayer of target cell plasma membrane, allowing cargo to diffuse to the surface of target cell. While this mechanism is supported by experimental evidence, its molecular details are unknown. The present study develops a coarse-grained model of nanoemulsion particles that are compatible with the MARTINI force field. Simulations using this coarse-grained model have demonstrated multiple fusion events between the particles and a model vesicular lipid bilayer. The fusion proceeds in the following sequence: dehydration at the interface, close apposition of the particles, protrusion of hydrophobic molecules to the particle surface, transient lipid complex formation, absorption of nanoemulsion into the liposome. The initial monolayer disruption acts as a rate-limiting step and is strongly influenced by particle size as well as by the presence of phospholipids supporting negative spontaneous curvature. The core-forming perfluorocarbons play critical roles in initiating the fusion process by facilitating protrusion of hydrophobic moieties into the interface between the two particles. This study directly supports the hypothesized nanoemulsion delivery mechanism and provides the underlying molecular details that enable engineering of nanoemulsions for a variety of medical applications.

Aggregation of model membrane proteins, modulated by hydrophobic mismatch, membrane curvature, and protein class

Aggregation of transmembrane proteins is important for many biological processes, such as protein sorting and cell signaling, and also for in vitro processes such as two-dimensional crystallization. We have used large-scale simulations to study the lateral organization and dynamics of lipid bilayers containing multiple inserted proteins. Using coarse-grained molecular dynamics simulations, we have studied model membranes comprising ∼7000 lipids and 16 identical copies of model cylindrical proteins of either α-helical or β-barrel types. Through variation of the lipid tail length and hence the degree of hydrophobic mismatch, our simulations display levels of protein aggregation ranging from negligible to extensive. The nature and extent of aggregation are shown to be influenced by membrane curvature and the shape or orientation of the protein. Interestingly, a model β-barrel protein aggregates to form one-dimensional strings within the bilayer plane, whereas a model α-helical bundle forms two-dimensional clusters. Overall, it is clear that the nature and extent of membrane protein aggregation is dependent on several aspects of the proteins and lipids, including hydrophobic mismatch, protein class and shape, and membrane curvature.

The role of the lipid matrix for structure and function of the GPCR rhodopsin

Photoactivation of rhodopsin in lipid bilayers results within milliseconds in a metarhodopsin I (MI)-metarhodopsin II (MII) equilibrium that is very sensitive to the lipid composition. It has been well established that lipid bilayers that are under negative curvature elastic stress from incorporation of lipids like phosphatidylethanolamines (PE) favor formation of MII, the rhodopsin photointermediate that is capable of activating G protein. Furthermore, formation of the MII state is favored by negatively charged lipids like phosphatidylserine and by lipids with longer hydrocarbon chains that yield bilayers with larger membrane hydrophobic thickness. Cholesterol and rhodopsin-rhodopsin interactions from crowding of rhodopsin molecules in lipid bilayers shift the MI-MII equilibrium towards MI. A variety of mechanisms seems to be responsible for the large, lipid-induced shifts between MI and MII: adjustment of the thickness of lipid bilayers to rhodopsin and adjustment of rhodopsin helicity to the thickness of bilayers, curvature elastic deformations in the lipid matrix surrounding the protein, direct interactions of PE headgroups and polyunsaturated hydrocarbon chains with rhodopsin, and direct or lipid-mediated interactions between rhodopsin molecules. This article is part of a Special Issue entitled: Membrane protein structure and function.

Lipid-controlled peptide topology and interactions in bilayers: structural insights into the synergistic enhancement of the antimicrobial activities of PGLa and magainin 2

To gain further insight into the antimicrobial activities of cationic linear peptides, we investigated the topology of each of two peptides, PGLa and magainin 2, in oriented phospholipid bilayers in the presence and absence of the other peptide and as a function of the membrane lipid composition. Whereas proton-decoupled (15)N solid-state NMR spectroscopy indicates that magainin 2 exhibits stable in-plane alignments under all conditions investigated, PGLa adopts a number of different membrane topologies with considerable variations in tilt angle. Hydrophobic thickness is an important parameter that modulates the alignment of PGLa. In equimolar mixtures of PGLa and magainin 2, the former adopts transmembrane orientations in dimyristoyl-, but not 1-palmitoyl-2-oleoyl-, phospholipid bilayers, whereas magainin 2 remains associated with the surface in all cases. These results have important consequences for the mechanistic models explaining synergistic activities of the peptide mixtures and will be discussed. The ensemble of data suggests that the thinning of the dimyristoyl membranes caused by magainin 2 tips the topological equilibrium of PGLa toward a membrane-inserted configuration. Therefore, lipid-mediated interactions play a fundamental role in determining the topology of membrane peptides and proteins and thereby, possibly, in regulating their activities as well.

Membrane protein integration into the endoplasmic reticulum

Most integral membrane proteins are targeted, inserted and assembled in the endoplasmic reticulum membrane. The sequential and potentially overlapping events necessary for membrane protein integration take place at sites termed translocons, which comprise a specific set of membrane proteins acting in concert with ribosomes and, probably, molecular chaperones to ensure the success of the whole process. In this minireview, we summarize our current understanding of helical membrane protein integration at the endoplasmic reticulum, and highlight specific characteristics that affect the biogenesis of multispanning membrane proteins.

Membrane thickness varies around the circumference of the transmembrane protein BtuB

BtuB is a large outer-membrane β-barrel protein that belongs to a class of active transport proteins that are TonB-dependent. These TonB-dependent transporters are based upon a 22-stranded antiparallel β-barrel, which is notably asymmetric in its length. Here, site-directed spin labeling and simulated annealing were used to locate the membrane lipid interface surrounding BtuB when reconstituted into phosphatidylcholine bilayers. Positions on the outer facing surface of the β-barrel and the periplasmic turns were spin-labeled and distances from the label to the membrane interface estimated by progressive power saturation of the electron paramagnetic resonance spectra. These distances were then used as atom-to-plane distance restraints in a simulated annealing routine, to dock the protein to two independent planes and produce a model representing the average position of the lipid phosphorus atoms at each interface. The model is in good agreement with the experimental data; however, BtuB is mismatched to the bilayer thickness and the resulting planes representing the bilayer interface are not parallel. In the model, the membrane thickness varies by 11 Å around the circumference of the protein, indicating that BtuB distorts the bilayer interface so that it is thinnest on the short side of the protein β-barrel.

Specificity of intramembrane protein-lipid interactions

Our concept of biological membranes has markedly changed, from the fluid mosaic model to the current model that lipids and proteins have the ability to separate into microdomains, differing in their protein and lipid compositions. Since the breakthrough in crystallizing membrane proteins, the most powerful method to define lipid-binding sites on proteins has been X-ray and electron crystallography. More recently, chemical biology approaches have been developed to analyze protein-lipid interactions. Such methods have the advantage of providing highly specific cellular probes. With the advent of novel tools to study functions of individual lipid species in membranes together with structural analysis and simulations at the atomistic resolution, a growing number of specific protein-lipid complexes are defined and their functions explored. In the present article, we discuss the various modes of intramembrane protein-lipid interactions in cellular membranes, including examples for both annular and nonannular bound lipids. Furthermore, we will discuss possible functional roles of such specific protein-lipid interactions as well as roles of lipids as chaperones in protein folding and transport.

Toward understanding protocell mechanosensation

Mechanosensitive (MS) channels can prevent bacterial bursting during hypo-osmotic shocks by responding to increases in lateral tension at the membrane level through an integrated and coordinated opening mechanism. Mechanical regulation in protocells could have been one of the first mechanisms to evolve in order to preserve their integrity against changing environmental conditions. How has the rich functional diversity found in present cells been created throughout evolution, and what did the primordial MS channels look like? This review has been written with the aim of identifying which factors may have been important for the appearance of the first osmotic valve in a prebiotic context, and what this valve may have been like. It highlights the mechanical properties of lipid bilayers, the association of peptides as aggregates in membranes, and the conservation of sequence motifs as central aspects to understand the evolution of proteins that gate below the tension required for spontaneous pore formation and membrane rupture. The arguments developed here apply to both MscL and MscS homologs, but could be valid to mechano-susceptible proteins in general.

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.

Modulation of plant mitochondrial VDAC by phytosterols

We have investigated the effect of cholesterol and two abundant phytosterols (sitosterol and stigmasterol) on the voltage-dependent anion-selective channel (VDAC) purified from mitochondria of bean seeds (Phaseolus coccineus). These sterols differ by the degree of freedom of their lateral chain. We show that VDAC displays sensitivity to the lipid-sterol ratio and to the type of sterol found in the membrane. The main findings of this study are: 1), cholesterol and phytosterols modulate the selectivity but only stigmasterol alters the voltage-dependence of the plant VDAC in the range of sterol fraction found in the plant mitochondrial membrane; 2), VDAC unitary conductance is not affected by the addition of sterols; 3), the effect of sterols on the VDAC is reversible upon sterol depletion with 10 μM methyl-β-cyclodextrins; and 4), phytosterols are essential for the channel gating at salt concentration prevailing in vivo. A quantitative analysis of the voltage-dependence indicates that stigmasterol inhibits the transition of the VDAC in the lowest subconductance states.

Structure of rhomboid protease in a lipid environment

Structures of the prokaryotic homologue of rhomboid proteases reveal a core of six transmembrane helices, with the active-site residues residing in a hydrophilic cavity. The native environment of rhomboid protease is a lipid bilayer, yet all the structures determined thus far are in a nonnative detergent environment. There remains a possibility of structural artefacts arising from the use of detergents. In an attempt to address the effect of detergents on the structure of rhomboid protease, crystals of GlpG, an Escherichia coli rhomboid protease in a lipid environment, were obtained using two alternative approaches. The structure of GlpG refined to 1. 7-Å resolution was obtained from crystals grown in the presence of lipid bicelles. This structure reveals well-ordered and partly ordered lipid molecules forming an annulus around the protein. Lipid molecules adapt to the surface features of protein and arrange such that they match the hydrophobic thickness of GlpG. Virtually identical two-dimensional crystals were also obtained after detergent removal by dialysis. A comparison of an equivalent structure determined in a completely delipidated detergent environment provides insights on how detergent substitutes for lipid. A detergent molecule is also observed close to the active site, helping to postulate a model for substrate binding and hydrolysis in rhomboids.

Channel and nonchannel forms of spin-labeled gramicidin in membranes and their equilibria

Channel and nonchannel forms of gramicidin A (GA) were studied by ESR in various lipid environments using new mono- and double-spin-labeled compounds. For GA channels, we demonstrate here how pulse dipolar ESR can be used to determine the orientation of the membrane-traversing molecule relative to the membrane normal and to study subtle effects of lipid environment on the interspin distance in the spin-labeled gramicidin channel. To study nonchannel forms of gramicidin, pulse dipolar ESR was used first to determine interspin distances corresponding to monomers and double-helical dimers of spin-labeled GA molecules in the organic solvents trifluoroethanol and octanol. The same distances were then observed in membranes. Since detection of nonchannel forms in the membrane is complicated by aggregation, we suppressed any dipolar spectra from intermolecular interspin distances arising from the aggregates by using double-labeled GA in a mixture with excess unlabeled GA. In hydrophobic mismatching lipids (L(β) phase of DPPC), gramicidin channels dissociate into free monomers. The backbone structure of the monomeric form is similar to a monomeric unit of the channel dimer. In addition to channels and monomers, the double-helical conformation of gramicidin is present in some membrane environments. In the gel phase of saturated phosphatidylcholines, the fraction of double helices increases in the following order: DLPC < DMPC < DSPC < DPPC. The equilibrium DHD/monomer ratio in DPPC was determined. In membranes, the double-helical form is present only in aggregates. In addition, we studied the effect of N-terminal substitution in the GA molecule upon channel formation. This work demonstrates how pulsed dipolar ESR may be utilized to study complex equilibria of peptides in membranes.

Atomistic simulation of lipid and DiI dynamics in membrane bilayers under tension

Membrane tension modulates cellular processes by initiating changes in the dynamics of its molecular constituents. To quantify the precise relationship between tension, structural properties of the membrane, and the dynamics of lipids and a lipophilic reporter dye, we performed atomistic molecular dynamics (MD) simulations of DiI-labeled dipalmitoylphosphatidylcholine (DPPC) lipid bilayers under physiological lateral tensions ranging from -2.6 mN m(-1) to 15.9 mN m(-1). Simulations showed that the bilayer thickness decreased linearly with tension consistent with volume-incompressibility, and this thinning was facilitated by a significant increase in acyl chain interdigitation at the bilayer midplane and spreading of the acyl chains. Tension caused a significant drop in the bilayer's peak electrostatic potential, which correlated with the strong reordering of water and lipid dipoles. For the low tension regime, the DPPC lateral diffusion coefficient increased with increasing tension in accordance with free-area theory. For larger tensions, free area theory broke down due to tension-induced changes in molecular shape and friction. Simulated DiI rotational and lateral diffusion coefficients were lower than those of DPPC but increased with tension in a manner similar to DPPC. Direct correlation of membrane order and viscosity near the DiI chromophore, which was just under the DPPC headgroup, indicated that measured DiI fluorescence lifetime, which is reported to decrease with decreasing lipid order, is likely to be a good reporter of tension-induced decreases in lipid headgroup viscosity. Together, these results offer new molecular-level insights into membrane tension-related mechanotransduction and into the utility of DiI in characterizing tension-induced changes in lipid packing.

A comparative study on the ability of two implicit solvent lipid models to predict transmembrane helix tilt angles

Free-energy profiles describing the relative orientation of membrane proteins along predefined coordinates can be efficiently calculated by means of umbrella simulations. Such simulations generate reliable orientational distributions but are difficult to converge because of the very long equilibration times of the solvent and the lipid bilayer in explicit representation. Two implicit lipid membrane models are here applied in combination with the umbrella sampling strategy to the simulation of the transmembrane (TM) helical segment from virus protein U (Vpu). The models are used to study both orientation and energetics of this α-helical peptide as a function of hydrophobic mismatch. We observe that increasing the degree of positive hydrophobic mismatch increased the tilt angle of Vpu. These findings agree well with experimental data and as such validate the solvation models used in this study.

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.

Morphological changes of supported lipid bilayers induced by lysozyme: planar domain formation vs. multilayer stacking

Total internal reflection fluorescence microscopy (TIRFM) has been utilized to explore the effect of cationic protein lysozyme (Lz) on the morphology of solid-supported lipid bilayers (SLBs) comprised of zwitterionic lipid phosphatidylcholine (PC) and its mixture with anionic lipid cardiolipin (CL). Kinetic TIRFM imaging of different systems revealed subtle interplay between lipid lateral segregation accompanied by exchange of neutral and acidic lipids in the protein-lipid interaction zone, and the formation of lipid multilayer stacks. The switch between these states was shown to be controlled by CL content. In weakly charged SLBs containing 5 mol% CL, assembling of CL molecules into planar domains upon Lz adsorption has been observed while at higher content of anionic lipid (25 mol%) in-plane domains tend to transform into multilayer stacks, thereby ensuring the most thermodynamically-favorable membrane conformation.

Conformational plasticity of the influenza A M2 transmembrane helix in lipid bilayers under varying pH, drug binding, and membrane thickness

Membrane proteins change their conformations to respond to environmental cues, thus conformational plasticity is important for function. The influenza A M2 protein forms an acid-activated proton channel important for the virus lifecycle. Here we have used solid-state NMR spectroscopy to examine the conformational plasticity of membrane-bound transmembrane domain of M2 (M2TM). (13)C and (15)N chemical shifts indicate coupled conformational changes of several pore-facing residues due to changes in bilayer thickness, drug binding, and pH. The structural changes are attributed to the formation of a well-defined helical kink at G34 in the drug-bound state and in thick lipid bilayers, nonideal backbone conformation of the secondary-gate residue V27 in the presence of drug, and nonideal conformation of the proton-sensing residue H37 at high pH. The chemical shifts constrained the (ϕ, ψ) torsion angles for three "basis" states, the equilibrium among which explains the multiple resonances per site in the NMR spectra under different combinations of bilayer thickness, drug binding, and pH conditions. Thus, conformational plasticity is important for the proton conduction and inhibition of M2TM. The study illustrates the utility of NMR chemical shifts for probing the structural plasticity and folding of membrane proteins.

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.

Protein folding in membranes

Separation of cells and organelles by bilayer membranes is a fundamental principle of life. Cellular membranes contain a baffling variety of proteins, which fulfil vital functions as receptors and signal transducers, channels and transporters, motors and anchors. The vast majority of membrane-bound proteins contain bundles of alpha-helical transmembrane domains. Understanding how these proteins adopt their native, biologically active structures in the complex milieu of a membrane is therefore a major challenge in today's life sciences. Here, we review recent progress in the folding, unfolding and refolding of alpha-helical membrane proteins and compare the molecular interactions that stabilise proteins in lipid bilayers. We also provide a critical discussion of a detergent denaturation assay that is increasingly used to determine membrane-protein stability but is not devoid of conceptual difficulties.

LCP-Tm: an assay to measure and understand stability of membrane proteins in a membrane environment

Structural and functional studies of membrane proteins are limited by their poor stability outside the native membrane environment. The development of novel methods to efficiently stabilize membrane proteins immediately after purification is important for biophysical studies, and is likely to be critical for studying the more challenging human targets. Lipidic cubic phase (LCP) provides a suitable stabilizing matrix for studying membrane proteins by spectroscopic and other biophysical techniques, including obtaining highly ordered membrane protein crystals for structural studies. We have developed a robust and accurate assay, LCP-Tm, for measuring the thermal stability of membrane proteins embedded in an LCP matrix. In its two implementations, protein denaturation is followed either by a change in the intrinsic protein fluorescence on ligand release, or by an increase in the fluorescence of a thiol-binding reporter dye that measures exposure of cysteines buried in the native structure. Application of the LCP-Tm assay to an engineered human beta2-adrenergic receptor and bacteriorhodopsin revealed a number of factors that increased protein stability in LCP. This assay has the potential to guide protein engineering efforts and identify stabilizing conditions that may improve the chances of obtaining high-resolution structures of intrinsically unstable membrane proteins.

Structure and alignment of the membrane-associated peptaibols ampullosporin A and alamethicin by oriented 15N and 31P solid-state NMR spectroscopy

Ampullosporin A and alamethicin are two members of the peptaibol family of antimicrobial peptides. These compounds are produced by fungi and are characterized by a high content of hydrophobic amino acids, and in particular the alpha-tetrasubstituted amino acid residue ?-aminoisobutyric acid. Here ampullosporin A and alamethicin were uniformly labeled with (15)N, purified and reconstituted into oriented phophatidylcholine lipid bilayers and investigated by proton-decoupled (15)N and (31)P solid-state NMR spectroscopy. Whereas alamethicin (20 amino acid residues) adopts transmembrane alignments in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) or 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) membranes the much shorter ampullosporin A (15 residues) exhibits comparable configurations only in thin membranes. In contrast the latter compound is oriented parallel to the membrane surface in 1,2-dimyristoleoyl-sn-glycero-3-phosphocholine and POPC bilayers indicating that hydrophobic mismatch has a decisive effect on the membrane topology of these peptides. Two-dimensional (15)N chemical shift -(1)H-(15)N dipolar coupling solid-state NMR correlation spectroscopy suggests that in their transmembrane configuration both peptides adopt mixed alpha-/3(10)-helical structures which can be explained by the restraints imposed by the membranes and the bulky alpha-aminoisobutyric acid residues. The (15)N solid-state NMR spectra also provide detailed information on the helical tilt angles. The results are discussed with regard to the antimicrobial activities of the peptides.

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.

Lipids and cholesterol as regulators of traffic in the endomembrane system

The endomembrane system of eukaryotic cells uses membrane-enclosed carriers to move diverse macromolecules among different membrane-bound compartments, a requirement for cells to secrete and take up molecules from their environment. Two recycling pathways-biosynthetic and endocytic, each with specific lipid components-make up this system, with the Golgi apparatus mediating transport between the two. Here, we integrate lipid-based mechanisms into the description of this system. A partitioning model of the Golgi apparatus is discussed as a working hypothesis to explain how membrane lipids and proteins that are segregated based on lateral lipid partitioning support the unique composition of the biosynthetic and endocytic recycling pathways in the face of constant trafficking of molecular constituents. We further discuss how computational modeling can allow for interpretation of experimental findings and provide mechanistic insight into these important cellular pathways.

Sorting of lens aquaporins and connexins into raft and nonraft bilayers: role of protein homo-oligomerization

Two classes of channel-forming proteins in the eye lens, the water channel aquaporin-0 (AQP-0) and the connexins Cx46 and Cx50, are preferentially located in different regions of lens plasma membranes (1,2). Because these membranes contain high concentrations of cholesterol and sphingomyelin, as well as phospholipids such as phosphatidylcholine with unsaturated hydrocarbon chains, microdomains (rafts) form in these membranes. Here we test the hypothesis that sorting into lipid microdomains can play a role in the disposition of AQP-0 and the connexins in the plane of the membrane. For both crude membrane fractions and proteoliposomes composed of lens proteins in phosphatidylcholine/sphingomyelin/cholesterol lipid bilayers, detergent extraction experiments showed that the connexins were located primarily in detergent soluble membrane (DSM) fractions, whereas AQP-0 was found in both detergent resistant membrane and DSM fractions. Analysis of purified AQP-0 reconstituted in raft-containing bilayers showed that the microdomain location of AQP-0 depended on protein/lipid ratio. AQP-0 was located almost exclusively in DSMs at a 1:1200 AQP-0/lipid ratio, whereas approximately 50% of the protein was sequestered into detergent resistant membranes at a 1:100 ratio, where freeze-fracture experiments show that AQP-0 oligomerizes (3). Consistent with these detergent extraction results, confocal microscopy images showed that AQP-0 was sequestered into raft microdomains in the 1:100 protein/lipid membranes. Taken together these results indicate that AQP-0 and connexins can be segregated in the membrane by protein-lipid interactions as modified by AQP-0 homo-oligomerization.

Alamethicin aggregation in lipid membranes

X-ray scattering features induced by aggregates of alamethicin (Alm) were obtained in oriented stacks of model membranes of DOPC(diC18:1PC) and diC22:1PC. The first feature obtained near full hydration was Bragg rod in-plane scattering near 0.11 A(-1) in DOPC and near 0.08 A(-1) in diC22:1PC at a 1:10 Alm:lipid ratio. This feature is interpreted as bundles consisting of n Alm monomers in a barrel-stave configuration surrounding a water pore. Fitting the scattering data to previously published molecular dynamics simulations indicates that the number of peptides per bundle is n = 6 in DOPC and n >or= 9 in diC22:1PC. The larger bundle size in diC22:1PC is explained by hydrophobic mismatch of Alm with the thicker bilayer. A second diffuse scattering peak located at q(r) approximately 0.7 A(-1) is obtained for both DOPC and diC22:1PC at several peptide concentrations. Theoretical calculations indicate that this peak cannot be caused by the Alm bundle structure. Instead, we interpret it as being due to two-dimensional hexagonally packed clusters in equilibrium with Alm bundles. As the relative humidity was reduced, interactions between Alm in neighboring bilayers produced more peaks with three-dimensional crystallographic character that do not index with the conventional hexagonal space groups.