Alamethicin channels in a membrane: molecular dynamics simulations

The Dynamic Range of Acidity: Tracking Rules for the Unidirectional Penetration of Cellular Compartments

Labeled ammonium cations with pK_a∼7.4 accumulate in acidic organelles because they can be neutralized transiently to cross the membrane at cytosolic pH 7.2 but not at their internal pH<5.5. Retention in early endosomes with less acidic internal pH was achieved recently using weaker acids of up to pK_a 9.8. We report here that primary ammonium cations with higher pK_a 10.6, label early endosomes more efficiently. This maximized early endosome tracking coincides with increasing labeling of Golgi networks with similarly weak internal acidity. Guanidinium cations with pK_a 13.5 cannot cross the plasma membrane in monomeric form and label the plasma membrane with selectivity for vesicles embarking into endocytosis. Self‐assembled into micelles, guanidinium cations enter cells like arginine‐rich cell‐penetrating peptides and, driven by their membrane potential, penetrate mitochondria unidirectionally despite their high inner pH. The resulting tracking rules with an approximated dynamic range of pK_a change ∼3.5 are expected to be generally valid, thus enabling the design of chemistry tools for biology research in the broadest sense. From a practical point of view, most relevant are two complementary fluorescent flipper probes that can be used to image the mechanics at the very beginning of endocytosis.

Spiers Memorial Lecture: Analysis and de novo design of membrane-interactive peptides

Membrane-peptide interactions play critical roles in many cellular and organismic functions, including protection from infection, remodeling of membranes, signaling, and ion transport. Peptides interact with membranes in a variety of ways: some associate with membrane surfaces in either intrinsically disordered conformations or well-defined secondary structures. Peptides with sufficient hydrophobicity can also insert vertically as transmembrane monomers, and many associate further into membrane-spanning helical bundles. Indeed, some peptides progress through each of these stages in the process of forming oligomeric bundles. In each case, the structure of the peptide and the membrane represent a delicate balance between peptide-membrane and peptide-peptide interactions. We will review this literature from the perspective of several biologically important systems, including antimicrobial peptides and their mimics, α-synuclein, receptor tyrosine kinases, and ion channels. We also discuss the use of de novo design to construct models to test our understanding of the underlying principles and to provide useful leads for pharmaceutical intervention of diseases.

NMR studies on the conformation, stability and dynamics of alamethicin in methanol

Alamethicin is an antibiotic peptide comprising 20 amino acid residues and functions as an ion channel in biological membranes. Natural alamethicins have a variety of amino acid sequences. Two of them, used as a mixed sample in this study, are: UPUAUAQUVUGLUPVUUQQO and UPUAUUQUVUGLUPVUUQQO, where U and O represent α-aminoisobutyric acid and phenylalaninol, respectively. As indicated, only the amino acid at position six differs, and the two alamethicins are referred to as alamethicin-A6 and -U6, respectively. The conformation and thermal stability of alamethicin-A6 and -U6 in methanol were examined using proton nuclear magnetic resonance (NMR) spectroscopy. Both alamethicins form an α-helix between the 2nd and 11th residues. The N-terminal, 19th and C-terminal residues take a non-helical conformation. The structure between the 12th and 18th residues has not been well determined due to the absence of cross peaks in the two-dimensional NMR data. The α-helices are maintained up to 54 °C at least. In contrast to these similarities, it has been found that the length of the α-helix of alamethicin-U6 is somewhat shorter than that of alamethicin-A6, the intra-molecular hydrogen bonds formed by the amide proton of the seventh residue is much more thermally stable for alamethicin-U6 than for alamethicin-A6, and the C-terminal residue of alamethicin-U6 has higher mobility than that of alamethicin-A6. The mobility of the N- and C-terminal residues is discussed on the basis of a model chain which consists of particles connected by rigid links, and the physiological significance of the mobility is emphasized.

Phosphorylation of conserved phosphoinositide binding pocket regulates sorting nexin membrane targeting

Sorting nexins anchor trafficking machines to membranes by binding phospholipids. The paradigm of the superfamily is sorting nexin 3 (SNX3), which localizes to early endosomes by recognizing phosphatidylinositol 3-phosphate (PI3P) to initiate retromer-mediated segregation of cargoes to the trans-Golgi network (TGN). Here we report the solution structure of full length human SNX3, and show that PI3P recognition is accompanied by bilayer insertion of a proximal loop in its extended Phox homology (PX) domain. Phosphoinositide (PIP) binding is completely blocked by cancer-linked phosphorylation of a conserved serine beside the stereospecific PI3P pocket. This “PIP-stop” releases endosomal SNX3 to the cytosol, and reveals how protein kinases control membrane assemblies. It constitutes a widespread regulatory element found across the PX superfamily and throughout evolution including of fungi and plants. This illuminates the mechanism of a biological switch whereby structured PIP sites are phosphorylated to liberate protein machines from organelle surfaces.

Sorting nexin 3 (SNX3) is a phosphatidylinositol 3-phosphate binding protein that localizes to early endosomes. Here the authors use NMR to resolve SNX3′s membrane interactions, revealing that membrane binding is regulated through phosphorylation of a conserved serine by its lipid recognition site.

Conformational and thermodynamic properties of non-canonical α,α-dialkyl glycines in the peptaibol Alamethicin: molecular dynamics studies

In this work, we investigate the structure, dynamic and thermodynamic properties of noncanonical disubstituted amino acids (α,α-dialkyl glycines), also known as non-natural amino acids, in the peptaibol Alamethicin. The amino acids under study are Aib (α-amino isobutyric acid or α-methyl alanine), Deg (α,α-diethyl glycine), Dpg (α,α-dipropyl glycine), Dibg (α,α-di-isobutyl glycine), Dhg (α,α-dihexyl glycine), DΦg (α,α-diphenyl glycine), Dbzg (α,α-dibenzyl glycine), Ac6c (α,α-cyclohexyl glycine), and Dmg (α,α-dihydroxymethyl glycine). It is hypothesized that these amino acids are able to induce well-defined secondary structure in peptidomimetics. To test this hypothesis, new peptidomimetics of Alamethicin were constructed by replacing the native Aib positions of Alamethicin by one or more new α,α-dialkyl glycines. Dhg and Ac6c demonstrated the capacity to induce well-defined α-helical structures. Dhg and Ac6c also promote the thermodynamic stabilization of these peptides in a POPC model membrane and are better alternatives to the Aib in Alamethicin. These noncanonical amino acids also improved secondary structure properties, revealing preorganization in water and maintenance of α helical structure in POPC. We show that it is possible to optimize the helicity and thermodynamic properties of native Alamethicin, and we suggest that these amino acids could be incorporated in other peptides with similar structural effect.

Membrane thickness and the mechanism of action of the short peptaibol trichogin GA IV

Trichogin GA IV (GAIV) is an antimicrobial peptide of the peptaibol family, like the extensively studied alamethicin (Alm). GAIV acts by perturbing membrane permeability. Previous data have shown that pore formation is related to GAIV aggregation and insertion in the hydrophobic core of the membrane. This behavior is similar to that of Alm and in agreement with a barrel-stave mechanism, in which transmembrane oriented peptides aggregate to form a channel. However, while the 19-amino acid long Alm has a length comparable to the membrane thickness, GAIV comprises only 10 amino acids, and its helix is about half the normal bilayer thickness. Here, we report the results of neutron reflectivity measurements, showing that GAIV inserts in the hydrophobic region of the membrane, causing a significant thinning of the bilayer. Molecular dynamics simulations of GAIV/membrane systems were also performed. For these studies we developed a novel approach for constructing the initial configuration, by embedding the short peptide in the hydrophobic core of the bilayer. These calculations indicated that in the transmembrane orientation GAIV interacts strongly with the polar phospholipid headgroups, drawing them towards its N- and C-termini, inducing membrane thinning and becoming able to span the bilayer. Finally, vesicle leakage experiments demonstrated that GAIV activity is significantly higher with thinner membranes, becoming similar to that of Alm when the bilayer thickness is comparable to its size. Overall, these data indicate that a barrel-stave mechanism of pore formation might be possible for GAIV and for similarly short peptaibols despite their relatively small size.

Mechanisms of peptide-induced pore formation in lipid bilayers investigated by oriented 31P solid-state NMR spectroscopy

There is a considerable interest in understanding the function of antimicrobial peptides (AMPs), but the details of their mode of action is not fully understood. This motivates extensive efforts in determining structural and mechanistic parameters for AMP’s interaction with lipid membranes. In this study we show that oriented-sample ³¹P solid-state NMR spectroscopy can be used to probe the membrane perturbations and -disruption by AMPs. For two AMPs, alamethicin and novicidin, we observe that the majority of the lipids remain in a planar bilayer conformation but that a number of lipids are involved in the peptide anchoring. These lipids display reduced dynamics. Our study supports previous studies showing that alamethicin adopts a transmembrane arrangement without significant disturbance of the surrounding lipids, while novicidin forms toroidal pores at high concentrations leading to more extensive membrane disturbance.

Backbone dynamics of alamethicin bound to lipid membranes: spin-echo electron paramagnetic resonance of TOAC-spin labels

Alamethicin F50/5 is a hydrophobic peptide that is devoid of charged residues and that induces voltage-dependent ion channels in lipid membranes. The peptide backbone is likely to be involved in the ion conduction pathway. Electron spin-echo spectroscopy of alamethicin F50/5 analogs in which a selected Aib residue (at position n = 1, 8, or 16) is replaced by the TOAC amino-acid spin label was used to study torsional dynamics of the peptide backbone in association with phosphatidylcholine bilayer membranes. Rapid librational motions of limited angular amplitude were observed at each of the three TOAC sites by recording echo-detected spectra as a function of echo delay time, 2tau. Simulation of the time-resolved spectra, combined with conventional EPR measurements of the librational amplitude, shows that torsional fluctuations of the peptide backbone take place on the subnanosecond to nanosecond timescale, with little temperature dependence. Associated fluctuations in polar fields from the peptide could facilitate ion permeation.

pH modulation of transport properties of alamethicin oligomers inserted in zwitterionic-based artificial lipid membranes

Electric features of biological membranes are major determinants of the function and physiological manifestation of membrane-penetrating peptides, and such features are prone to be modulated by the properties of the surrounding aqueous medium. In this work, we demonstrate that pH plays crucial roles in modulating electric characteristics of zwitterionic-based artificial lipid membranes. The effect of pH on electrical properties of such membranes was probed by evaluating the transport properties of embedded alamethicin oligomers over a wide range of pH values (i.e., 0.65, 2.08, 2.94, 7 and 10.1). Our data strongly support the paradigm of a pH-dependent variation of the surface and membrane dipole potential which, in conjunction with possible lateral pressure effects within the lipid membrane, lead to a non-monotonic modulation of the electrical conductance of alamethicin oligomers. As expected, pH modulation of transport properties through the alamethicin oligomer is more visible for narrower pores (that is, the 1st conductive state) with slightly better cation selectivity as compared to larger oligomers.

Understanding the energetics of helical peptide orientation in membranes

Understanding the energetic factors determining the positioning and orientation of single-helical peptides in membranes is of fundamental interest in structural biology. Here, a simple 5-slab continuum dielectric model for the membrane is examined that distinguishes between the solvent, headgroup, and core regions. An analytical solution for the electrostatic solvation of a single dipole and an all-atom model of N-methylacetamide are used to demonstrate the effect of the dielectric boundaries in the system on peptide dipole orientation. The dipole orientation energy is shown to dominate the electrostatic solvation energy of a polyalanine helix in the membrane. With an additional surface-area-dependent term to account for the cavity formation in the aqueous region, the continuum electrostatics description is used to examine several helical peptides, the atoms of which are explicitly represented with a molecular mechanics force field. The experimentally determined tilt angles of a number of peptides of alternating alanine and leucine residues, and of glycophorin and melittin, are accurately reproduced by the model. The factors determining the tilt angles and their fluctuations are analyzed. The tilt angles of the simpler peptides are found to increase approximately linearly with peptide length; this effect is also rationalized. The analysis and model presented here provide a step toward the prediction of helical membrane protein structure.

Voltage-dependent energetics of alamethicin monomers in the membrane

The implicit membrane model IMM1 is extended to include the effect of transmembrane potential and used to investigate the optimal membrane binding configurations and energies for alamethicin helices. In the absence of voltage, the lowest energy configuration is on the membrane surface with a tilt allowing the N terminus to be fully buried. Slightly higher in energy is an also tilted configuration with the N terminus deeper in the membrane and almost crossing the membrane. In 26A membranes and in the presence of 0.1V voltage, the TM orientation becomes lower in energy. This is consistent with the assumption that voltage induces a transition from the interfacial to the inserted (TM) orientation. This effect of voltage is smaller in thicker membranes. The results are compared to previous experimental and theoretical studies and the findings are discussed in relation to the mechanism of channel formation by alamethicin.

The implementation of slab geometry for membrane-channel molecular dynamics simulations

Slab geometric boundary conditions are applied in the molecular dynamics simulation of a simple membrane-channel system. The results of the simulation were compared to those of an analogous system using normal three-dimensional periodic boundary conditions. Analysis of the dynamics and electrostatics of the system show that slab geometric periodicity eliminates the artificial bulk water orientational polarization that is present while using normal three-dimensional periodicity. Furthermore, even though the water occupancy and volume of our simple channel is the same when using either method, the electrostatic properties are considerably different when using slab geometry. In particular, the orientational polarization of water is seen to be different in the interior of the channel. This gives rise to a markedly different electric field within the channel. We discuss the implications of slab geometry for the future simulation of this type of system and for the study of channel transport properties.

A tethered bilayer sensor containing alamethicin channels and its detection of amiloride based inhibitors

Alamethicin, a small transmembrane peptide, inserts into a tethered bilayer membrane (tBLM) to form ion channels, which we have investigated using electrical impedance spectroscopy. The number of channels formed is dependent on the incubation time, concentration of the alamethicin and the application of DC voltage. The properties of the ion channels when formed in tethered bilayers are similar to those for such channels assembled into black lipid membranes (BLMs). Furthermore, amiloride and certain analogs can inhibit the channel pores, formed in the tBLMs. The potency and concentration of the inhibitors can be determined by measuring the change of impedance. Our work illustrates the possibility of using a synthetic tBLM for the study of small peptide voltage dependent ion channels. A potential application of such a device is as a screening tool in drug discovery processes.

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.

Molecular dynamics simulations of wild-type and mutant forms of the Mycobacterium tuberculosis MscL channel

The crystal structure of the Mycobacterium tuberculosis homolog of the bacterial mechanosensitive channel of large conductance (Tb-MscL) provides a unique opportunity to consider mechanosensitive signal transduction at the atomic level. Molecular dynamics simulations of the Tb-MscL channel embedded in an explicit lipid bilayer and of its C-terminal helical bundle alone in aqueous solvent were performed. C-terminal calculations imply that although the helix bundle structure is relatively unstable at physiological pH, it may have been stabilized under low pH conditions such as those used in the crystallization of the channel. Specific mutations to the C-terminal region, which cause a similar conservation of the crystal structure conformation, have also been identified. Full channel simulations were performed for the wild-type channel and two experimentally characterized gain-of-function mutants, V21A and Q51E. The wild-type Tb-MscL trajectory gives insight into regions of relative structural stability and instability in the channel structure. Channel mutations led to observable changes in the trajectories, such as an alteration of intersubunit interactions in the Q51E mutant. In addition, interesting patterns of protein-lipid interactions, such as hydrogen bonding, arose in the simulations. These and other observations from the simulations are relevant to previous and ongoing experimental studies focusing on characterization of the channel.

Characterization of (Na+, K+)-ATPase in gill microsomes of the freshwater shrimp Macrobrachium olfersii

To better understand the adaptive strategies that led to freshwater invasion by hyper-regulating Crustacea, we prepared a microsomal (Na+, K+)-ATPase by differential centrifugation of a gill homogenate from the freshwater shrimp Macrobrachium olfersii. Sucrose gradient centrifugation revealed a light fraction containing most of the (Na+, K+)-ATPase activity, contaminated with other ATPases, and a heavy fraction containing negligible (Na+, K+)-ATPase activity. Western blotting showed that M. olfersii gill contains a single alpha-subunit isoform of about 110 kDa. The (Na+, K+)-ATPase hydrolyzed ATP with Michaelis Menten kinetics with K5, = 165+/-5 microM and Vmax = 686.1+/-24.7 U mg(-1). Stimulation by potassium (K0.5 = 2.4+/-0.1 mM) and magnesium ions (K0.5 = 0.76+/-0.03 mM) also obeyed Michaelis-Menten kinetics, while that by sodium ions (K0.5 = 6.0+/-0.2 mM) exhibited site site interactions (n = 1.6). Ouabain (K0.5 = 61.6+/-2.8 microM) and vanadate (K0.5 = 3.2+/-0.1 microM) inhibited up to 70% of the total ATPase activity, while thapsigargin and ethacrynic acid did not affect activity. The remaining 30% activity was inhibited by oligomycin, sodium azide and bafilomycin A. These data suggest that the (Na+, K+)-ATPase corresponds to about 70% of the total ATPase activity; the remaining 30%, i.e. the ouabain-insensitive ATPase activity, apparently correspond to F0F1- and V-ATPases, but not Ca-stimulated and Na- or K-stimulated ATPases. The data confirm the recent invasion of the freshwater biotope by M. olfersii and suggest that (Na+, K+)-ATPase activity may be regulated by the Na+ concentration of the external medium.

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.

The mechanism of channel formation by alamethicin as viewed by molecular dynamics simulations

Alamethicin is a 20-residue channel-forming peptide that forms a stable amphipathic alpha-helix in membrane and membrane-mimetic environments. This helix contains a kink induced by a central Gly-X-X-Pro sequence motif. Alamethicin channels are activated by a cis positive transbilayer voltage. Channel activation is suggested to correspond to voltage-induced insertion of alamethicin helices in the bilayer. Alamethicin forms multi-conductance channels in lipid bilayers. These channels are formed by parallel bundles of transmembrane helices surrounding a central pore. A change in the number of helices per bundle switches the single channel conductance level. Molecular dynamics simulations of alamethicin in a number of different environments have been used to explore its channel-forming properties. These simulations include: (i) alamethicin in solution in water and in methanol; (ii) a single alamethicin helix at the surface of a phosphatidylcholine bilayer; (iii) single alamethicin helices spanning a phosphatidylcholine bilayer; and (iv) channels formed by bundles of 5, 6, 7 or 8 alamethicin helices spanning a phosphatidylcholine bilayer. The total simulation time is c. 30 ns. Thus, these simulations provide a set of dynamic snapshots of a possible mechanism of channel formation by this peptide.

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.