Dynamic pictures of membrane proteins in two-dimensional crystal, lipid bilayer and detergent as revealed by site-directed solid-state 13C NMR

We have compared site-directed 13C solid-state NMR spectra of [3-13C]Ala- and/or [1-13C]Val-labeled membrane proteins, including bacteriorhodopsin (bR), pharaonis phoborhodopin (ppR), its cognate transducer (pHtrII) and Escherichia coli diacylglycerol kinase (DGK), in two-dimensional (2D) crystal, lipid bilayers, and detergent. Restricted fluctuation motions of these membrane proteins due to oligomerization of bR by specific protein-protein interactions in the 2D crystalline lattice or protein complex between ppR and pHtrII provide the most favorable environment to yield well-resolved, fully visible 13C NMR signals for [3-13C]Ala-labeled proteins. In contrast, several signals from such membrane proteins were broadened or lost owing to interference of inherent fluctuation frequencies (10(4)-10(5)Hz) with frequency of either proton decoupling or magic angle spinning, if their 13C NMR spectra were recorded as a monomer in lipid bilayers at ambient temperature. The presence of such protein dynamics is essential for the respective proteins to achieve their own biological functions. Finally, spectral broadening found for bR and DGK in detergents were discussed.     

Structure and plasticity of the human immunodeficiency virus gp41 fusion domain in lipid micelles and bilayers

A thorough understanding of the structure of fusion domains of enveloped viruses in changing lipid environments helps us to formulate mechanistic models on how they might function in mediating viral entry by membrane fusion. We have expressed the N-terminal fusion domain of HIV-1 gp41 as a construct that is water-soluble in the absence of membranes, but that also binds with high affinity to lipid micelles and bilayers in their presence. We have solved the structure and studied the dynamics of this domain bound to dodecylphosphocholine micelles by homo- and heteronuclear NMR spectroscopy. The fusion peptide forms a stable hydrophobic helix from Ile(4) to Ala(14), but is increasingly more disordered and dynamic in a segment of intermediate polarity that stretches from Ala(15) to Ser(23). When bound to lipid bilayers at low concentration, the HIV fusion domain is also largely alpha-helical, as determined by CD and FTIR spectroscopy. However, at higher protein/lipid ratios, the domain is partially converted to form beta-structures in lipid bilayers. Controlled lipid mixing occurs at concentrations that support the alpha-helical, but not the beta-strand conformation.     

Structure and orientation of sarcolipin in lipid environments.

Sarcolipin (SLN) is a 31 amino acid integral membrane protein that regulates Ca-ATPase activity in skeletal muscle. Here, we report the three-dimensional structure and topology of synthetic SLN in lipid environments, as determined by solution and solid-state NMR spectroscopy. 2D solution NMR experiments were performed on SLN solubilized in sodium dodecyl sulfate (SDS) micelles. We found that SLN adopts a highly defined alpha-helical conformation from F9 through R27, with a backbone RMSD of 0.65 A and a side chain RMSD of 1.66 A. The N-terminus (M1 through L8) and the C-terminus (S28 through Y31) are mostly unstructured. The orientation of the SLN was determined using one-dimensional (15)N NMR solid-state spectroscopy. The protein was incorporated into phospholipid bilayers prepared from a mixture of 1,2-dioleoyl-sn-glycero-3-phosphocholine and 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine. The (15)N chemical shift solid-state spectra from selectively labeled SLN samples indicate that SLN orients perpendicularly to the plane of the membrane bilayers. These results support the proposed mechanism of Ca-ATPase regulation of SLN via protein-protein intramembranous interactions between the highly conserved transmembrane domains of the Ca-ATPase and the conserved transmembrane domain of SLN.     

Correlation of gene and protein structures in the FXYD family proteins

The FXYD family proteins are auxiliary subunits of the Na,K-ATPase, expressed primarily in tissues that specialize in fluid or solute transport, or that are electrically excitable. These proteins range in size from about 60 to 160 amino acid residues, and share a core homology of 35 amino acid residues in and around a single transmembrane segment. Despite their relatively small sizes, they are all encoded by genes with six to nine small exons. We show that the helical secondary structures of three FXYD family members, FXYD1, FXYD3, and FXYD4, determined in micelles by NMR spectroscopy, reflect the structures of their corresponding genes. The coincidence of helical regions, and connecting segments, with the positions of intron-exon junctions in the genes, support the hypothesis that the FXYD proteins may have been assembled from discrete structural modules through exon shuffling.     

Structural dynamics and topology of phospholamban in oriented lipid bilayers using multidimensional solid-state NMR

Phospholamban (PLN), a single-pass membrane protein, regulates heart muscle contraction and relaxation by reversible inhibition of the sarco(endo)plasmic reticulum Ca-ATPase (SERCA). Studies in detergent micelles and oriented lipid bilayers have shown that in its monomeric form PLN adopts a dynamic L shape (bent or T state) that is in conformational equilibrium with a more dynamic R state. In this paper, we use solid-state NMR on both uniformly and selectively labeled PLN to refine our initial studies, describing the topology and dynamics of PLN in oriented lipid bilayers. Two-dimensional PISEMA (polarization inversion spin exchange at the magic angle) experiments carried out in DOPC/DOPE mixed lipid bilayers reveal a tilt angle of the transmembrane domain with respect to the static magnetic field, of 21 +/- 2 degrees and, at the same time, map the rotation angle of the transmembrane domain with respect to the bilayer. PISEMA spectra obtained with selectively labeled samples show that the cytoplasmic domain of PLN is helical and makes an angle of 93 +/- 6 degrees with respect to the bilayer normal. In addition, using samples tilted by 90 degrees , we find that the transmembrane domain of PLN undergoes fast long-axial rotational diffusion about the bilayer normal with the cytoplasmic domain undergoing this motion and other complex dynamics, scaling the values of chemical shift anisotropy. While this dynamic was anticipated by previous solution NMR relaxation studies in micelles, these measurements in the anisotropic lipid environment reveal new dynamic and conformational features encoded in the free protein that might be crucial for SERCA recognition and subsequent inhibition.     

Structures of the FXYD regulatory proteins in lipid micelles and membranes

The FXYD membrane proteins constitute a family of conserved auxiliary subunits of the Na,K-ATPase, and have been the focus of recent attention due to their ability to finely regulate the activity of the enzyme complex in various physiological settings. In this review we describe the structures of the proteins, as well as their dynamics and their associations with the lipid bilayer membrane, which we have recently determined by NMR spectroscopy. Although the proteins are relatively small, their genes contain as many as six to nine small exons, and the coincidence of structured protein segments with their genetic elements suggests assembly from discrete structural modules through exon shuffling. The three-dimensional structures and backbone dynamics provide the foundation for understanding their intra-membrane association with the Na,K-ATPase α subunit, and the structure of FXYD1 suggests a mechanism whereby the phosphorylation of conserved Ser residues, by protein kinases A and C, could induce a conformational change in the cytoplasmic domain of the protein, to modulate its interaction with the α subunit.     

Thermal, dynamic and structural properties of drug AT1 antagonist olmesartan in lipid bilayers

It is proposed that AT1 antagonists (ARBs) exert their biological action by inserting into the lipid membrane and then diffuse to the active site of AT1 receptor. Thus, lipid bilayers are expected to be actively involved and play a critical role in drug action. For this reason, the thermal, dynamic and structural effects of olmesartan alone and together with cholesterol were studied using differential scanning calorimetry (DSC), 13C magic-angle spinning (MAS) nuclear magnetic resonance (NMR), cross-polarization (CP) MAS NMR, and Raman spectroscopy as well as small- and wide angle X-ray scattering (SAXS and WAXS) on dipalmitoyl-phosphatidylcholine (DPPC) multilamellar vesicles. 13C CP/MAS spectra provided direct evidence for the incorporation of olmesartan and cholesterol in lipid bilayers. Raman and X-ray data revealed how both molecules modify the bilayer's properties. Olmesartan locates itself at the head-group region and upper segment of the lipid bilayers as 13C CP/MAS spectra show that its presence causes significant chemical shift changes mainly in the A ring of the steroidal part of cholesterol. The influence of olmesartan on DPPC/cholesterol bilayers is less pronounced. Although, olmesartan and cholesterol are residing at the same region of the lipid bilayers, due to their different sizes, display distinct impacts on the bilayer's properties. Cholesterol broadens significantly the main transition, abolishes the pre-transition, and decreases the membrane fluidity above the main transition. Olmesartan is the only so far studied ARB that increases the gauche:trans ratio in the liquid crystalline phase. These significant differences of olmesartan may in part explain its distinct pharmacological profile.     

Effects of glycosylation on the structure and dynamics of eel calcitonin in micelles and lipid bilayers determined by nuclear magnetic resonance spectroscopy.

The three-dimensional structures of eel calcitonin (CT) and two glycosylated CT derivatives, [Asn(GlcNAc)3]-CT (CT-GlcNAc) and [Asn(Man6-GlcNAc2)3]-CT (CT-M6), in micelles were determined by solution NMR spectroscopy. The topologies of these peptides associated with oriented lipid bilayers were determined with solid-state NMR. All of the peptides were found to have an identical conformation in micelles characterized by an amphipathic alpha-helix consisting of residues Ser5 through Leu19 followed by an unstructured region at the C-terminus. The overall conformation of the peptide moiety was not affected by the glycosylation. Nevertheless, comparison of the relative exchange rates of the Leu12 amide proton might suggest the possibility that fluctuations of the alpha-helix are reduced by glycosylation. The presence of NOEs between the carbohydrate and the peptide moieties of CT-GlcNAc and CT-M6 and the amide proton chemical shift data suggested that the carbohydrate interacted with the peptide, and this might account for the conformational stabilization of the alpha-helix. Both the unmodified CT and the glycosylated CT were found to have orientations with their helix axes parallel to the plane of the lipid bilayers by solid-state NMR spectroscopy.     

The interaction of amino-deuteromethylated melittin with phospholipid membranes studied by deuterium NMR

Melittin, deuteromethylated on each of the four amino groups (Gly-1 N alpha and Lys-7, 21, and 23 N epsilon), was prepared by reductive methylation using deuteroformaldehyde and NaBD3CN. Deuterium NMR spectra were obtained for the modified peptide (D-melittin) bound to phospholipid bilayers and erythrocyte ghosts. D-Melittin at 4 mol% (peptide:lipid) induced reversible transitions between extended bilayers and micelles at the phase-transition temperature in dimyristoylphosphatidylcholine (DMPC) bilayers. These changes in lipid morphology did not occur at 1 mol% D-melittin: DMPC and the peptide was highly motionally restricted in gel in gel-phase lipid.     

Atomic view of calcium-induced clustering of phosphatidylserine in mixed lipid bilayers

Membranes play key regulatory roles in biological processes, with bilayer composition exerting marked effects on binding affinities and catalytic activities of a number of membrane-associated proteins. In particular, proteins involved in diverse processes such as vesicle fusion, intracellular signaling cascades, and blood coagulation interact specifically with anionic lipids such as phosphatidylserine (PS) in the presence of Ca(2+) ions. While Ca(2+) is suspected to induce PS clustering in mixed phospholipid bilayers, the detailed structural effects of this ion on anionic lipids are not established. In this study, combining magic angle spinning (MAS) solid-state NMR (SSNMR) measurements of isotopically labeled serine headgroups in mixed lipid bilayers with molecular dynamics (MD) simulations of PS lipid bilayers in the presence of different counterions, we provide site-resolved insights into the effects of Ca(2+) on the structure and dynamics of lipid bilayers. Ca(2+)-induced conformational changes of PS in mixed bilayers are observed in both liposomes and Nanodiscs, a nanoscale membrane mimetic of bilayer patches. Site-resolved multidimensional correlation SSNMR spectra of bilayers containing (13)C,(15)N-labeled PS demonstrate that Ca(2+) ions promote two major PS headgroup conformations, which are well resolved in two-dimensional (13)C-(13)C, (15)N-(13)C, and (31)P-(13)C spectra. The results of MD simulations performed on PS lipid bilayers in the presence or absence of Ca(2+) provide an atomic view of the conformational effects underlying the observed spectra.     

Membrane composition determines pardaxin's mechanism of lipid bilayer disruption

Pardaxin is a membrane-lysing peptide originally isolated from the fish Pardachirus marmoratus. The effect of the carboxy-amide of pardaxin (P1a) on bilayers of varying composition was studied using (15)N and (31)P solid-state NMR of mechanically aligned samples and differential scanning calorimetry (DSC). (15)N NMR spectroscopy of [(15)N-Leu(19)]P1a found that the orientation of the peptide's C-terminal helix depends on membrane composition. It is located on the surface of lipid bilayers composed of 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) and is inserted in lipid bilayers composed of 1,2-dimyristoyl-phosphatidylcholine (DMPC). The former suggests a carpet mechanism for bilayer disruption whereas the latter is consistent with a barrel-stave mechanism. The (31)P chemical shift NMR spectra showed that the peptide significantly disrupts lipid bilayers composed solely of zwitterionic lipids, particularly bilayers composed of POPC, in agreement with a carpet mechanism. P1a caused the formation of an isotropic phase in 1-palmitoyl-2-oleoyl-phosphatidylethanolamine (POPE) lipid bilayers. This, combined with DSC data that found P1a reduced the fluid lamellar-to-inverted hexagonal phase transition temperature at very low concentrations (1:50,000), is interpreted as the formation of a cubic phase and not micellization of the membrane. Experiments exploring the effect of P1a on lipid bilayers composed of 4:1 POPC:cholesterol, 4:1 POPE:cholesterol, 3:1 POPC:1-palmitoyl-2-oleoyl-phosphatidylglycerol (POPG), and 3:1 POPE:POPG were also conducted, and the presence of anionic lipids or cholesterol was found to reduce the peptide's ability to disrupt bilayers. Considered together, these data demonstrate that the mechanism of P1a is dependent on membrane composition.     

Identification of trapped and boundary lipid binding sites in M13 coat protein/lipid complexes by deuterium NMR spectroscopy

The major coat protein of M13 bacteriophage has been incorporated into bilayers of 1,2-dimyristoyl-sn-glycero-3-phosphocholine, deuterated in the trimethyl segments of the choline headgroup (DMPC-d9). Two-component deuterium and phosphorus-31 NMR spectra have been observed from bilayer complexes containing the coat protein, indicating slow exchange (on the deuterium quadrupole anisotropy and phosphorus-31 chemical shift averaging time scales) of lipid molecules of less than 10(3) Hz between two motionally distinct environments in the complexes. The fraction of the isotropic spectral component increases with increasing M13 protein concentration, and this component is attributed to lipid headgroups, which are disordered relative to their order in protein-free bilayers. The activation energy of the fast local motions of the trimethyl groups of the choline residue in the headgroup decreases from 23 kJ mol-1 in the pure lipid bilayers to 20 kJ mol-1 for the protein-associated lipid headgroups. The chemical exchange rate of lipid molecules between the two motionally distinct environments has been estimated to be 20-50 Hz by steady-state line-shape simulations of the deuterium spectra of DMPC-d9/M13 coat protein complexes using exchange-coupled modified Bloch equations. The off-rate was, as expected from one-to-one exchange, independent of the L/P ratio; tau off -1 = 0.23 kHz. It is suggested that the protein-associated lipid may be trapped between closely packed parallel aggregates of M13 coat protein and that the high local concentration of protein in a one-dimensional arrangement in lipid bilayers may be required for the fast reassembly of phage particles before release from an infected cell.     

Mechanism of lipid bilayer disruption by the human antimicrobial peptide,LL-37

LL-37 is an amphipathic, alpha-helical, antimicrobial peptide. (15)N chemical shift and (15)N dipolar-shift spectroscopy of site-specifically labeled LL-37 in oriented lipid bilayers indicate that the amphipathic helix is oriented parallel to the surface of the bilayer. This surface orientation is maintained in both anionic and zwitterionic bilayers and at different temperatures and peptide concentrations, ruling out a barrel-stave mechanism for bilayer disruption by LL-37. In contrast, electrostatic factors, the type of lipid, and the presence of cholesterol do affect the extent to which LL-37 perturbs the lipids in the bilayer as observed with (31)P NMR. The (31)P spectra also show that micelles or other small, rapidly tumbling membrane fragments are not formed in the presence of LL-37, excluding a detergent-like mechanism. LL-37 does increase the lamellar to inverted hexagonal phase transition temperature of both PE model lipid systems and Escherichia coli lipids, demonstrating that it induces positive curvature strain in these environments. These results support a toroidal pore mechanism of lipid bilayer disruption by LL-37.     

Phase structure of liposome in lipid mixtures

Gas microbubbles present in ultrasound imaging contrast agents are stabilized by lipid aggregates that typically contain a mixture of lipids. In this study, the phase structure of the lipid mixtures that contained two or three lipids was investigated using three different methods: dynamic light scattering, ¹H NMR, and microfluidity measurements with fluorescence probes. Three lipids that are commonly present in imaging agents (DPPC, DPPE-PEG, and DPPA) were used. Two types of systems, two-lipid model systems and simulated imaging systems were investigated. The results show that liposomes were the dominant aggregates in all the samples studied. The polar PEG side chains from the PEGylated lipid lead to the formation of micelles and micellar aggregates in small sizes. In the ternary lipid systems, almost all the lipids were present in bilayers with micelles absent and free lipids at very low concentration. These results suggest that liposomes, not micelles, contribute to the stabilization of microbubbles in an ultrasound imaging contrast agent.     

Effect of head group and curvature on binding of the antimicrobial peptide tritrpticin to lipid membranes

In this work we examine the interaction between the 13-residue cationic antimicrobial peptide (AMP) tritrpticin (VRRFPWWWPFLRR, TRP3) and model membranes of variable lipid composition. The effect on peptide conformational properties was investigated by means of CD (circular dichroism) and fluorescence spectroscopies. Based on the hypothesis that the antibiotic acts through a mechanism involving toroidal pore formation, and taking into account that models of toroidal pores imply the formation of positive curvature, we used large unilamellar vesicles (LUV) to mimic the initial step of peptide-lipid interaction, when the peptide binds to the bilayer membrane, and micelles to mimic the topology of the pore itself, since these aggregates display positive curvature. In order to more faithfully assess the role of curvature, micelles were prepared with lysophospholipids containing (qualitatively and quantitatively) head groups identical to those of bilayer phospholipids. CD and fluorescence spectra showed that, while TRP3 binds to bilayers only when they carry negatively charged phospholipids, binding to micelles occurs irrespective of surface charge, indicating that electrostatic interactions play a less predominant role in the latter case. Moreover, the conformations acquired by the peptide were independent of lipid composition in both bilayers and micelles. However, the conformations were different in bilayers and in micelles, suggesting that curvature has an influence on the secondary structure acquired by the peptide. Fluorescence data pointed to an interfacial location of TRP3 in both types of aggregates. Nevertheless, experiments with a water soluble fluorescence quencher suggested that the tryptophan residues are more accessible to the quencher in micelles than in bilayers. Thus, we propose that bilayers and micelles can be used as models for the two steps of toroidal pore formation.     

The structural topology of wild-type phospholamban in oriented lipid bilayers using 15N solid-state NMR spectroscopy

For the first time, 15N solid-state NMR experiments were conducted on wild-type phospholamban (WT-PLB) embedded inside mechanically oriented phospholipid bilayers to investigate the topology of its cytoplasmic and transmembrane domains. 15N solid-state NMR spectra of site-specific 15N-labeled WT-PLB indicate that the transmembrane domain has a tilt angle of 13 degrees+/-6 degrees with respect to the POPC (1-palmitoyl-2-oleoyl-sn-glycero-phosphocholine) bilayer normal and that the cytoplasmic domain of WT-PLB lies on the surface of the phospholipid bilayers. Comparable results were obtained from site-specific 15N-labeled WT-PLB embedded inside DOPC/DOPE (1,2-dioleoyl-sn-glycero-3-phosphocholine/1,2-dioleoyl-sn-glycero-3-phosphoethanolamine) mechanically oriented phospholipids' bilayers. The new NMR data support a pinwheel geometry of WT-PLB, but disagree with a bellflower structure in micelles, and indicate that the orientation of the cytoplasmic domain of the WT-PLB is similar to that reported for the monomeric AFA-PLB mutant.     

Structure and Dynamics of the Membrane-Bound Form of Pf1 Coat Protein: Implications of Structural Rearrangement for Virus Assembly

The three-dimensional structure of the membrane-bound form of the major coat protein of Pf1 bacteriophage was determined in phospholipid bilayers using orientation restraints derived from both solid-state and solution NMR experiments. In contrast to previous structures determined solely in detergent micelles, the structure in bilayers contains information about the spatial arrangement of the protein within the membrane, and thus provides insights to the bacteriophage assembly process from membrane-inserted to bacteriophage-associated protein. Comparisons between the membrane-bound form of the coat protein and the previously determined structural form found in filamentous bacteriophage particles demonstrate that it undergoes a significant structural rearrangement during the membrane-mediated virus assembly process. The rotation of the transmembrane helix (Q16-A46) around its long axis changes dramatically (by 160°) to obtain the proper alignment for packing in the virus particles. Furthermore, the N-terminal amphipathic helix (V2-G17) tilts away from the membrane surface and becomes parallel with the transmembrane helix to form one nearly continuous long helix. The spectra obtained in glass-aligned planar lipid bilayers, magnetically aligned lipid bilayers (bicelles), and isotropic lipid bicelles reflect the effects of backbone motions and enable the backbone dynamics of the N-terminal helix to be characterized. Only resonances from the mobile N-terminal helix and the C-terminus (A46) are observed in the solution NMR spectra of the protein in isotropic q > 1 bicelles, whereas only resonances from the immobile transmembrane helix are observed in the solid-state ¹H/¹⁵N-separated local field spectra in magnetically aligned bicelles. The N-terminal helix and the hinge that connects it to the transmembrane helix are significantly more dynamic than the rest of the protein, thus facilitating structural rearrangement during bacteriophage assembly.     

Correct folding of the beta-barrel of the human membrane protein VDAC requires a lipid bilayer

Spontaneous membrane insertion and folding of beta-barrel membrane proteins from an unfolded state into lipid bilayers has been shown previously only for few outer membrane proteins of Gram-negative bacteria. Here we investigated membrane insertion and folding of a human membrane protein, the isoform 1 of the voltage-dependent anion-selective channel (hVDAC1) of mitochondrial outer membranes. Two classes of transmembrane proteins with either alpha-helical or beta-barrel membrane domains are known from the solved high-resolution structures. VDAC forms a transmembrane beta-barrel with an additional N-terminal alpha-helix. We demonstrate that similar to bacterial OmpA, urea-unfolded hVDAC1 spontaneously inserts and folds into lipid bilayers upon denaturant dilution in the absence of folding assistants or energy sources like ATP. Recordings of the voltage-dependence of the single channel conductance confirmed folding of hVDAC1 to its active form. hVDAC1 developed first beta-sheet secondary structure in aqueous solution, while the alpha-helical structure was formed in the presence of lipid or detergent. In stark contrast to bacterial beta-barrel membrane proteins, hVDAC1 formed different structures in detergent micelles and phospholipid bilayers, with higher content of beta-sheet and lower content of alpha-helix when inserted and folded into lipid bilayers. Experiments with mixtures of lipid and detergent indicated that the content of beta-sheet secondary structure in hVDAC1 decreased at increased detergent content. Unlike bacterial beta-barrel membrane proteins, hVDAC1 was not stable even in mild detergents such as LDAO or dodecylmaltoside. Spontaneous folding of outer membrane proteins into lipid bilayers indicates that in cells, the main purpose of membrane-inserted or associated assembly factors may be to select and target beta-barrel membrane proteins towards the outer membrane instead of actively assembling them under consumption of energy as described for the translocons of cytoplasmic membranes.     

Probing ground and excited states of phospholamban in model and native lipid membranes by magic angle spinning NMR spectroscopy

In this paper, we analyzed the ground and excited states of phospholamban (PLN), a membrane protein that regulates sarcoplasmic reticulum calcium ATPase (SERCA), in different membrane mimetic environments. Previously, we proposed that the conformational equilibria of PLN are central to SERCA regulation. Here, we show that these equilibria detected in micelles and bicelles are also present in native sarcoplasmic reticulum lipid membranes as probed by MAS solid-state NMR. Importantly, we found that the kinetics of conformational exchange and the extent of ground and excited states in detergent micelles and lipid bilayers are different, revealing a possible role of the membrane composition on the allosteric regulation of SERCA. Since the extent of excited states is directly correlated to SERCA inhibition, these findings open up the exciting possibility that calcium transport in the heart can be controlled by the lipid bilayer composition. This article is part of a Special Issue entitled: Membrane protein structure and function.     

Molecular dynamics simulations of indolicidin association with model lipid bilayers

Identifying the mechanisms responsible for the interaction of peptides with cell membranes is critical to the design of new antimicrobial peptides and membrane transporters. We report here the results of a computational simulation of the interaction of the 13-residue peptide indolicidin with single-phase lipid bilayers of dioleoylphosphatidylcholine, distearoylphosphatidylcholine, dioleoylphosphatidylglycerol, and distearoylphosphatidylglycerol. Ensemble analysis of the membrane-bound peptide revealed that, in contrast to the extended, linear backbone structure reported for indolicidin in sodium dodecyl sulphate detergent micelles, the peptide adopts a boat-shaped conformation in both phosphatidylglycerol and phosphatidylcholine lipid bilayers, similar to that reported for dodecylphosphocholine micelles. In agreement with fluorescence and NMR experiments, simulations confirmed that the peptide localizes in the membrane interface, with the distance between phosphate headgroups of each leaflet being reduced in the presence of indolicidin. These data, along with a concomitant decrease in lipid order parameters for the upper-tail region, suggest that indolicidin binding results in membrane thinning, consistent with recent in situ atomic force microscopy studies.