Surface binding of alamethicin stabilizes its helical structure: molecular dynamics simulations.

Alamethicin is an amphipathic alpha-helical peptide that forms ion channels. An early event in channel formation is believed to be the binding of alamethicin to the surface of a lipid bilayer. Molecular dynamics simulations are used to compare the structural and dynamic properties of alamethicin in water and alamethicin bound to the surface of a phosphatidylcholine bilayer. The bilayer surface simulation corresponded to a loosely bound alamethicin molecule that interacted with lipid headgroups but did not penetrate the hydrophobic core of the bilayer. Both simulations started with the peptide molecule in an alpha-helical conformation and lasted 2 ns. In water, the helix started to unfold after approximately 300 ps and by the end of the simulation only the N-terminal region of the peptide remained alpha-helical and the molecule had collapsed into a more compact form. At the surface of the bilayer, loss of helicity was restricted to the C-terminal third of the molecule and the rod-shaped structure of the peptide was retained. In the surface simulation about 10% of the peptide/water H-bonds were replaced by peptide/lipid H-bonds. These simulations suggest that some degree of stabilization of an amphipathic alpha-helix occurs at a bilayer surface even without interactions between hydrophobic side chains and the acyl chain core of the bilayer.     

Analysis and evaluation of channel models: simulations of alamethicin

Alamethicin is an antimicrobial peptide that forms stable channels with well-defined conductance levels. We have used extended molecular dynamics simulations of alamethicin bundles consisting of 4, 5, 6, 7, and 8 helices in a palmitoyl-oleolyl-phosphatidylcholine bilayer to evaluate and analyze channel models and to link the models to the experimentally measured conductance levels. Our results suggest that four helices do not form a stable water-filled channel and might not even form a stable intermediate. The lowest measurable conductance level is likely to correspond to the pentamer. At higher aggregation numbers the bundles become less symmetrical. Water properties inside the different-sized bundles are similar. The hexamer is the most stable model with a stability comparable with simulations based on crystal structures. The simulation was extended from 4 to 20 ns or several times the mean passage time of an ion. Essential dynamics analyses were used to test the hypothesis that correlated motions of the helical bundles account for high-frequency noise observed in open channel measurements. In a 20-ns simulation of a hexameric alamethicin bundle, the main motions are those of individual helices, not of the bundle as a whole. A detailed comparison of simulations using different methods to treat long-range electrostatic interactions (a twin range cutoff, Particle Mesh Ewald, and a twin range cutoff combined with a reaction field correction) shows that water orientation inside the alamethicin channels is sensitive to the algorithms used. In all cases, water ordering due to the protein structure is strong, although the exact profile changes somewhat. Adding an extra 4-nm layer of water only changes the water ordering slightly in the case of particle mesh Ewald, suggesting that periodicity artifacts for this system are not serious.     

Conformation of peptides in lipid membranes studied by x-ray grazing incidence scattering

Although the antimicrobial, fungal peptide alamethicin has been extensively studied, the conformation of the peptide and the interaction with lipid bilayers as well as the mechanism of channel gating are still not completely clear. As opposed to studies of the crystalline state, the polypeptide structures in the environment of fluid bilayers are difficult to probe. We have investigated the conformation of alamethicin in highly aligned stacks of model lipid membranes by synchrotron-based x-ray scattering. The (wide-angle) scattering distribution has been measured by reciprocal space mappings. A pronounced scattering signal is observed in samples of high molar peptide/lipid ratio which is distinctly different from the scattering distribution of an ideal helix in the transmembrane state. Beyond simple models of ideal helices, the data is analyzed in terms of models based on atomic coordinates from the Brookhaven Protein Data Bank, as well as from published molecular dynamics simulations. The results can be explained by assuming a wide distribution of helix tilt angles with respect to the membrane normal and a partial insertion of the N-terminus into the membrane.     

Peptide Aggregation and Pore Formation in a Lipid Bilayer: A Combined Coarse-Grained and All Atom Molecular Dynamics Study

We present a simulation study where different resolutions, namely coarse-grained (CG) and all-atom (AA) molecular dynamics simulations, are used sequentially to combine the long timescale reachable by CG simulations with the high resolution of AA simulations, to describe the complete processes of peptide aggregation and pore formation by alamethicin peptides in a hydrated lipid bilayer. In the 1-μs CG simulations the peptides spontaneously aggregate in the lipid bilayer and exhibit occasional transitions between the membrane-spanning and the surface-bound configurations. One of the CG systems at t = 1 μs is reverted to an AA representation and subjected to AA simulation for 50 ns, during which water molecules penetrate the lipid bilayer through interactions with the peptide aggregates, and the membrane starts leaking water. During the AA simulation significant deviations from the α-helical structure of the peptides are observed, however, the size and arrangement of the clusters are not affected within the studied time frame. Solid-state NMR experiments designed to match closely the setup used in the molecular dynamics simulations provide strong support for our finding that alamethicin peptides adopt a diverse set of configurations in a lipid bilayer, which is in sharp contrast to the prevailing view of alamethicin oligomers formed by perfectly aligned helical alamethicin peptides in a lipid bilayer.     

Simulation studies of the interaction of antimicrobial peptides and lipid bilayers

Experimental studies of a number of antimicrobial peptides are sufficiently detailed to allow computer simulations to make a significant contribution to understanding their mechanisms of action at an atomic level. In this review we focus on simulation studies of alamethicin, melittin, dermaseptin and related antimicrobial, membrane-active peptides. All of these peptides form amphipathic alpha-helices. Simulations allow us to explore the interactions of such peptides with lipid bilayers, and to understand the effects of such interactions on the conformational dynamics of the peptides. Mean field methods employ an empirical energy function, such as a simple hydrophobicity potential, to provide an approximation to the membrane. Mean field approaches allow us to predict the optimal orientation of a peptide helix relative to a bilayer. Molecular dynamics simulations that include an atomistic model of the bilayer and surrounding solvent provide a more detailed insight into peptide-bilayer interactions. In the case of alamethicin, all-atom simulations have allowed us to explore several steps along the route from binding to the membrane surface to formation of transbilayer ion channels. For those antimicrobial peptides such as dermaseptin which prefer to remain at the surface of a bilayer, molecular dynamics simulations allow us to explore the favourable interactions between the peptide helix sidechains and the phospholipid headgroups.     

Interactions of alpha-helices with lipid bilayers: a review of simulation studies

Membrane proteins, of which the majority seem to contain one or more alpha-helix, constitute approx. 30% of most genomes. A complete understanding of the nature of helix/bilayer interactions is necessary for an understanding of the structural principles underlying membrane proteins. This review describes computer simulation studies of helix/bilayer interactions. Key experimental studies of the interactions of alpha-helices and lipid bilayers are briefly reviewed. Surface associated helices are found in some membrane-bound enzymes (e.g. prostaglandin synthase), and as stages in the mechanisms of antimicrobial peptides and of pore-forming bacterial toxins. Transmembrane alpha-helices are found in most integral membrane proteins, and also in channels formed by amphipathic peptides or by bacterial toxins. Mean field simulations, in which the lipid bilayer is approximated as a hydrophobic continuum, have been used in studies of membrane-active peptides (e.g. alamethicin, melittin, magainin and dermaseptin) and of simple membrane proteins (e.g. phage Pf1 coat protein). All atom molecular dynamics simulations of fully solvated bilayers with transmembrane helices have been applied to: the constituent helices of bacteriorhodopsin; peptide-16 (a simple model TM helix); and a number of pore-lining helices from ion channels. Surface associated helices (e.g. melittin and dermaseptin) have been simulated, as have alpha-helical bundles such as bacteriorhodopsin and alamethicin. From comparison of the results from the two classes of simulation, it emerges that a major theoretical challenge is to exploit the results of all atom simulations in order to improve the mean field approach.     

Proline-induced hinges in transmembrane helices: possible roles in ion channel gating

A number of ion channels contain transmembrane (TM) alpha-helices that contain proline-induced molecular hinges. These TM helices include the channel-forming peptide alamethicin (Alm), the S6 helix from voltage-gated potassium (Kv) channels, and the D5 helix from voltage-gated chloride (CLC) channels. For both Alm and KvS6, experimental data implicate hinge-bending motions of the helix in an aspect of channel gating. We have compared the hinge-bending motions of these TM helices in bilayer-like environments by multi-nanosecond MD simulations in an attempt to describe motions of these helices that may underlie possible modes of channel gating. Alm is an alpha-helical channel-forming peptide, which contains a central kink associated with a Gly-x-x-Pro motif in its sequence. Simulations of Alm in a TM orientation for 10 ns in an octane slab indicate that the Gly-x-x-Pro motif acts as a molecular hinge. The S6 helix from Shaker Kv channels contains a Pro-Val-Pro motif. Modeling studies and recent experimental data suggest that the KvS6 helix may be kinked in the vicinity of this motif. Simulations (10 ns) of an isolated KvS6 helix in an octane slab and in a POPC bilayer reveal hinge-bending motions. A pattern-matching approach was used to search for possible hinge-bending motifs in the TM helices of other ion channel proteins. This uncovered a conserved Gly-x-Pro motif in TM helix D5 of CLC channels. MD simulations of a model of hCLC1-D5 spanning an octane slab suggest that this channel also contains a TM helix that undergoes hinge-bending motion. In conclusion, our simulations suggest a model in which hinge-bending motions of TM helices may play a functional role in the gating mechanisms of several different families of ion channels.     

Asymmetrical ion-channel model inferred from two-dimensional crystallization of a peptide antibiotic

The structural organization of ion channels formed in lipid membranes by amphiphilic alpha-helical peptides is deduced by applying direct structural methods to different lipid/alamethicin systems. Alamethicin represents a hydrophobic alpha-helical peptide antibiotic forming voltage-gated ion channels in lipid membranes. Here the first direct evidence for the existence of large-scale two-dimensional crystalline domains of alamethicin helices, oriented parallel to the air/water interface, is presented using synchrotron x-ray diffraction, fluorescence microscopy, and surface pressure/area isotherms. Proofs are obtained that the antibiotic peptide injected into the aqueous phase under phospholipid monolayers penetrates these monolayers, phase separates, and forms domains within the lipid environment, keeping the same, parallel orientation of the alpha-helices with respect to the phospholipid/water interface. A new asymmetrical, "lipid-covered ring" model of the voltage-gated ion channel of alamethicin is inferred from the structural results presented, and the mechanism of ion-channel formation is discussed.     

Folding is not required for bilayer insertion: replica exchange simulations of an alpha-helical peptide with an explicit lipid bilayer

We implement the replica exchange molecular dynamics algorithm to study the interactions of a model peptide (WALP-16) with an explicitly represented DPPC membrane bilayer. We observe the spontaneous, unbiased insertion of WALP-16 into the DPPC bilayer and its folding into an alpha-helix with a transbilayer orientation. The free energy surface suggests that the insertion of the peptide into the DPPC bilayer precedes secondary structure formation. Although the peptide has some propensity to form a partially helical structure in the interfacial region of the DPPC/water system, this state is not a productive intermediate but rather an off-pathway trap for WALP-16 insertion. Equilibrium simulations show that the observed insertion/folding pathway mirrors the potential of mean force (PMF). Calculation of the enthalpic and entropic contributions to this PMF show that the surface bound conformation of WALP-16 is significantly lower in energy than other conformations, and that the insertion of WALP-16 into the bilayer without regular secondary structure is enthalpically unfavorable by 5-10 kcal/mol/residue. The observed insertion/folding pathway disagrees with the dominant conceptual model, which is that a surface-bound helix is an obligatory intermediate for the insertion of alpha-helical peptides into lipid bilayers. In our simulations, the observed insertion/folding pathway is favored because of a large (>100 kcal/mol) increase in system entropy that occurs when the unstructured WALP-16 peptide enters the lipid bilayer interior. The insertion/folding pathway that is lowest in free energy depends sensitively on the near cancellation of large enthalpic and entropic terms. This suggests the possibility that intrinsic membrane peptides may have a diversity of insertion/folding behaviors depending on the exact system of peptide and lipid under consideration.     

Determining Peptide Partitioning Properties via Computer Simulation

The transfer of polypeptide segments into lipid bilayers to form transmembrane helices represents the crucial first step in cellular membrane protein folding and assembly. This process is driven by complex and poorly understood atomic interactions of peptides with the lipid bilayer environment. The lack of suitable experimental techniques that can resolve these processes both at atomic resolution and nanosecond timescales has spurred the development of computational techniques. In this review, we summarize the significant progress achieved in the last few years in elucidating the partitioning of peptides into lipid bilayer membranes using atomic detail molecular dynamics simulations. Indeed, partitioning simulations can now provide a wealth of structural and dynamic information. Furthermore, we show that peptide-induced bilayer distortions, insertion pathways, transfer free energies, and kinetic insertion barriers are now accurate enough to complement experiments. Further advances in simulation methods and force field parameter accuracy promise to turn molecular dynamics simulations into a powerful tool for investigating a wide range of membrane active peptide phenomena.     

MD simulations of Mistic: conformational stability in detergent micelles and water

Mistic is an unusual membrane protein from Bacillus subtilis. It appears to fold and insert autonomously into a lipid bilayer and has been suggested as a tool that aids the targeting of eukaryotic membrane proteins to bacterial membranes. The NMR structure of Mistic in detergent (LDAO) micelles has revealed it to be a four alpha-helix bundle. From a structural perspective, Mistic does not resemble other membrane proteins. Its external surface is not very hydrophobic, and standard methods do not predict any of its helices to be in the transmembrane orientation. Molecular dynamics simulations (simulation times approximately 30 ns) in water and in detergent micelles have been used to explore the conformational stability of Mistic as a function of its environment. In water, the protein is stable, exhibiting no significant change in fold on a 30 ns time scale. In contrast, in three simulations in detergent micelles, the partial unfolding of Mistic occurred, whereby the H4 helix drifted away from the H1-H3 core. This was due to the penetration of detergent molecules between H4 and the remainder of the protein. This is unlike the behavior of several other membrane proteins, both alpha-helix bundles and beta-barrels, in comparable detergent micelle simulations. The unfolding of H4 from the H1-H3 core of Mistic could be partially reversed by a simulation in which the detergent molecules were removed, and the unfolded protein was simulated in water. These results suggest that Mistic may not be a stable integrated membrane protein but rather that it may undergo a conformational change upon interaction with a membrane or membrane-like environment.     

The energetics of transmembrane helix insertion into a lipid bilayer

Free energy profiles for insertion of a hydrophobic transmembrane protein α-helix (M2 from CFTR) into a lipid bilayer have been calculated using coarse-grained molecular dynamics simulations and umbrella sampling to yield potentials of mean force along a reaction path corresponding to translation of a helix across a lipid bilayer. The calculated free energy of insertion is smaller when a bilayer with a thinner hydrophobic region is used. The free energies of insertion from the potentials of mean force are compared with those derived from a number of hydrophobicity scales and with those derived from translocon-mediated insertion. This comparison supports recent models of translocon-mediated insertion and in particular suggests that: 1), helices in an about-to-be-inserted state may be located in a hydrophobic region somewhat thinner than the core of a lipid bilayer; and/or 2), helices in a not-to-be-inserted state may experience an environment more akin (e.g., in polarity/hydrophobicity) to the bilayer/water interface than to bulk water.     

Interaction of the peptide antibiotic alamethicin with bilayer- and non-bilayer-forming lipids: influence of increasing alamethicin concentration on the lipids supramolecular structures

Incorporation of the helical antimicrobial peptide alamethicin from aqueous phase into hydrated phases of dioleoylphosphatidylethanolamine (DOPE) and dioleoylphosphatidylcholine (DOPC) was investigated within a range of peptide concentrations and temperatures by time-resolved synchrotron X-ray diffraction. It was found that alamethicin influences the organizations of the non-bilayer-forming (DOPE) and the bilayer-forming (DOPC) lipids in different ways. In DOPC, only the bilayer thickness was affected, while in DOPE new phases were induced. At low peptide concentrations (<1.10(-4) M), an inverted hexagonal (H(II)) phase was observed as with DOPE dispersions in pure buffer solution. A coexistence of two cubic structures was found at the critical peptide concentration for induction of new lipid/peptide phases. The first one Q224 (space group Pn3m) was identified within the entire temperature region studied (from 1 to 45 degrees C) and was found in coexistence with H(II)-phase domains. The second lipid/peptide cubic structure was present only at temperatures below 16 degrees C and its X-ray reflections were better fitted by a Q212 (P4(3)32) space group, rather than by the expected Q229 (Im3m) space group. At alamethicin concentrations of 1 mM and higher, a nonlamellar phase transition from a Q224 cubic phase into an H(II) phase was observed. Within the investigated range of peptide concentrations, lamellar structures of two different bilayer periods were established with the bilayer-forming lipid DOPC. They correspond to lipid domains of associated and nonassociated helical peptide. The obtained X-ray results suggest that the amphiphilic alamethicin molecules adsorb from the aqueous phase at the lipid head group/water interface of the DOPE and DOPC membranes. At sufficiently high (>1.10(-4) M) solution concentrations, the peptide is probably accommodated in the head group region of the lipids thus inducing structural features of mixed lipid/peptide phases.     

Molecular dynamics investigation of membrane-bound bundles of the channel-forming transmembrane domain of viral protein U from the human immunodeficiency virus HIV-1

Molecular dynamics (MD) simulations have been carried out on bundles of the channel-forming transmembrane (TM) domain of the viral protein U (VPU(1-27) and VPU(6-27)) from the human immunodeficiency virus (HIV-1). Simulations of hexameric and pentameric bundles of VPU(6-27) in an octane/water membrane mimetic system suggested that the pentamer is the preferred oligomer. Accordingly, an unconstrained pentameric helix bundle of VPU(1-27) was then placed in a hydrated palmitoyl-oleyl-3-n-glycero-phosphatidylethanolamine (POPE) lipid bilayer and its structural properties calculated from a 3-ns MD run. Some water molecules, initially inside the channel lumen, were expelled halfway through the simulation and the bundle adopted a conical structure reminiscent of previous MD results obtained for VPU(6-27) in an octane/water system. The pore constriction generated may correspond to a closed state of the channel and underlies the relocation of the W residue toward the pore lumen. The relative positions of the helices with respect to the bilayer and their interactions with the lipids are discussed. The observed structure is stabilized via specific interactions between the VPU helices and the carbonyl oxygen atoms of the lipid molecules, particularly at the Q and S residues.     

High count rates with total internal reflection fluorescence correlation spectroscopy

We have modeled the structure of KirBac1.1 in an open state using as a starting point the structure of KirBac1.1 in its closed conformation (Protein Data Bank 1P7B). To test the validity of the open-state model, molecular dynamics simulations in octane, a lipid bilayer mimetic, were carried out. Simulations of the closed conformer were used for comparison purposes. The total simulation time was approximately 138 ns. The initial open model was refined by using projection maps obtained from electron microscopy experiments on two-dimensional crystals of the inwardly rectifying K+ channel KirBac3.1 from Magentospirillum magnetotacticum captured in its open state (C. Vénien-Bryan, unpublished data). Significant movements of the outer helices take place in going from the closed to the open model in agreement with structural and biochemical data in potassium channels, which suggests that gating is accomplished by a conformational change that takes place in the transmembrane domain upon an external stimulus. The motion of the inner helices is mainly achieved by bending at conserved glycine residues that have been previously reported to act as molecular hinges. Overall, these simulations suggest that the open conformer is stable, providing a plausible all-atom model that will enable the study of potential gating mechanisms in more detail.     

Reversible transition between α-helix and β-sheet conformation of a transmembrane domain

Despite the important functions of protein transmembrane domains, their structure and dynamics are often scarcely known. The SNARE proteins VAMP/synaptobrevin and syntaxin 1 are implicated in membrane fusion. Using different spectroscopic approaches we observed a marked sensitivity of their transmembrane domain structure in regard to the lipid/peptide ratio. In the dilute condition, peptides corresponding to the complete transmembrane domain fold into an α-helix inserted at ∼ 35° to the normal of the membranes, an observation in line with molecular simulations. Upon an increase in the peptide/lipid ratio, the peptides readily exhibited transition to β-sheet structure. Moreover, the insertion angle of these β-sheets increased to 54° and was accompanied by a derangement of lipid acyl chains. For both proteins the transition from α-helix to β-sheet was reversible under certain conditions by increasing the peptide/lipid ratio. This phenomenon was observed in different model systems including multibilayers and small unilamellar vesicles. In addition, differences in peptide structure and transitions were observed when using distinct lipids (DMPC, DPPC or DOPC) thus indicating parameters influencing transmembrane domain structure and conversion from helices to sheets. The putative functional consequences of this unprecedented dynamic behavior of a transmembrane domain are discussed.     

Structure of magainin and alamethicin in model membranes studied by x-ray reflectivity

We have investigated the structure of lipid bilayers containing varied molar ratios of different lipids and the antimicrobial peptides magainin and alamethicin. For this structural study, we have used x-ray reflectivity on highly aligned solid-supported multilamellar lipid membranes. The reflectivity curves have been analyzed by semi-kinematical reflectivity theory modeling the bilayer density profile rho(z). Model simulations of the reflectivity curves cover a large range of vertical momentum transfer q(z), and yield excellent agreement between data and theory. The structural changes observed as a function of the molar peptide/lipid concentration P/L are discussed in a comparative way.     

Location and dynamics of alamethicin in unilamellar vesicles and thylakoids as model systems. A spin label study

Location and dynamics of the voltage-dependent pore-forming icosapeptide alamethicin have been studied using spin labels which were linked directly and via spacers to the C-terminus of the amphiphilic alpha-helix. Ion-transport activities of these derivatives were found to be very similar to those of natural alamethicin in green plant thylakoids chosen as a model system. The shape of the electron spin resonance spectra indicates segmental motion of the nitroxide rather than rotation of the whole peptide. A population of spins showing narrow lines in the presence of thylakoids or lipid vesicles is attributed to alamethicin in the aqueous solution. A second population shows rotational correlation times greater than 10(-9) s and is bound to the membranes, the C-termini residing in an environment with a polarity close to that of water. This population is inaccessible to the hydrophilic, charged line broadening agent chromium oxalate. Since spectral shapes and amplitudes of spectra are unchanged by additions of unlabelled peptide, it is concluded that the ESR detectable spins are bound to peptides essentially in the monomeric state. Alamethicin induced pore formation under flash illumination is demonstrated by measurement of kinetics of proton deposition in the thylakoid interior. When pores are opened by illuminating thylakoids and thus applying a membrane potential, mainly the bound population is affected by a process reversibly suppressing the signal, whereas only limited disappearance of label from the external medium is detected. Apparently, the potential causes a change in the conformation of the peptide which leads to a further immobilisation of the label, possibly due to a deeper insertion of the alpha-helices into the lipid membrane. However, evidence has been presented experimentally that there is no detectable change of potential prior to the opening of the pore.