D.P Tieleman

An alamethicin channel in a lipid bilayer: molecular dynamics simulations.

We present the results of 2-ns molecular dynamics (MD) simulations of a hexameric bundle of Alm helices in a 1-palmitoyl-2-oleoylphosphatidylcholine bilayer. These simulations explore the dynamic properties of a model of a helix bundle channel in a complete phospholipid bilayer in an aqueous environment. We explore the stability and conformational dynamics of the bundle in a phospholipid bilayer. We also investigate the effect on bundle stability of the ionization state of the ring of Glu18 side chains. If all of the Glu18 side chains are ionised, the bundle is unstable; if none of the Glu18 side chains are ionized, the bundle is stable. pKA calculations suggest that either zero or one ionized Glu18 is present at neutral pH, correlating with the stable form of the helix bundle. The structural and dynamic properties of water in this model channel were examined. As in earlier in vacuo simulations (Breed et al., 1996 .Biophys. J. 70:1643-1661), the dipole moments of water molecules within the pore were aligned antiparallel to the helix dipoles. This contributes to the stability of the helix bundle.

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.

Voltage-Dependent Insertion of Alamethicin at Phospholipid/Water and Octane/Water Interfaces

Understanding the binding and insertion of peptides in lipid bilayers is a prerequisite for understanding phenomena such as antimicrobial activity and membrane-protein folding. We describe molecular dynamics simulations of the antimicrobial peptide alamethicin in lipid/water and octane/water environments, taking into account an external electric field to mimic the membrane potential. At cis-positive potentials, alamethicin does not insert into a phospholipid bilayer in 10 ns of simulation, due to the slow dynamics of the peptide and lipids. However, in octane N-terminal insertion occurs at field strengths from 0.33 V/nm and higher, in simulations of up to 100 ns duration. Insertion of alamethicin occurs in two steps, corresponding to desolvation of the Gln7 side chain, and the backbone of Aib10 and Gly11. The proline induced helix kink angle does not change significantly during insertion. Polyalanine and alamethicin form stable helices both when inserted in octane and at the water/octane interface, where they partition in the same location. In water, both polyalanine and alamethicin partially unfold in multiple simulations. We present a detailed analysis of the insertion of alamethicin into the octane slab and the influence of the external field on the peptide structure. Our findings give new insight into the mechanism of channel formation by alamethicin and the structure and dynamics of membrane-associated helices.

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.

Defining the Transmembrane Helix of M2 Protein from Influenza A by Molecular Dynamics Simulations in a Lipid Bilayer

Integral membrane proteins containing at least one transmembrane (TM) alpha-helix are believed to account for between 20% and 30% of most genomes. There are several algorithms that accurately predict the number and position of TM helices within a membrane protein sequence. However, these methods tend to disagree over the beginning and end residues of TM helices, posing problems for subsequent modeling and simulation studies. Molecular dynamics (MD) simulations in an explicit lipid and water environment are used to help define the TM helix of the M2 protein from influenza A virus. Based on a comparison of the results of five different secondary structure prediction algorithms, three different helix lengths (an 18mer, a 26mer, and a 34mer) were simulated. Each simulation system contained 127 POPC molecules plus approximately 3500-4700 waters, giving a total of approximately 18,000-21,000 atoms. Two simulations, each of 2 ns duration, were run for the 18mer and 26mer, and five separate simulations were run for the 34mer, using different starting models generated by restrained in vacuo MD simulations. The total simulation time amounted to 11 ns. Analysis of the time-dependent secondary structure of the TM segments was used to define the regions that adopted a stable alpha-helical conformation throughout the simulation. This analysis indicates a core TM region of approximately 20 residues (from residue 22 to residue 43) that remained in an alpha-helical conformation. Analysis of atomic density profiles suggested that the 18mer helix revealed a local perturbation of the lipid bilayer. Polar side chains on either side of this region form relatively long-lived H-bonds to lipid headgroups and water molecules.

Alamethicin channels in a membrane: molecular dynamics simulations

Alamethicin (Alm) is a 20 residue peptide which forms a kinked alpha-helix in membrane and membrane-mimetic environments. Ion channels formed by intramembraneous aggregates of Alm are thought to be formed by bundles of approximately parallel Alm helices surrounding a central bilayer pore. Different channel conductance levels correspond to different numbers of helices per bundle, ranging from N = 5 to N > 8. Calculation of the predicted pKA values of the ring of Glu18 sidechains at the C-terminal mouth of the pore suggests that at neutral pH most or all of these sidechains will remain protonated. Nanosecond molecular dynamics (MD) simulations of N = 5, 6, 7 and 8 bundles of Alm helices in a POPC bilayer have been run, corresponding to a total simulation time of 4 ns. These simulations explore the stability and conformational dynamics of these helix bundle channels when embedded in a full phospholipid bilayer in an aqueous environment. The structural and dynamic properties of water in these model channels are examined. As in earlier in vacuo simulations (J. Breed, R. Sankararamakrishnan, I. D. Kerr and M. S. P. Sansom, Biophys. J., 1996, 70, 1643) the dipole moments of water molecules within the pores are aligned antiparallel to the helix dipoles. This helps to contribute to the stability of the helix bundles.