The role of membrane thickness in charged protein–lipid interactions

Charged amino acids are known to be important in controlling the actions of integral and peripheral membrane proteins and cell disrupting peptides. Atomistic molecular dynamics studies have shed much light on the mechanisms of membrane binding and translocation of charged protein groups, yet the impact of the full diversity of membrane physico-chemical properties and topologies has yet to be explored. Here we have performed a systematic study of an arginine (Arg) side chain analog moving across saturated phosphatidylcholine (PC) bilayers of variable hydrocarbon tail length from 10 to 18 carbons. For all bilayers we observe similar ion-induced defects, where Arg draws water molecules and lipid head groups into the bilayers to avoid large dehydration energy costs. The free energy profiles all exhibit sharp climbs with increasing penetration into the hydrocarbon core, with predictable shifts between bilayers of different thickness, leading to barrier reduction from 26 kcal/mol for 18 carbons to 6 kcal/mol for 10 carbons. For lipids of 10 and 12 carbons we observe narrow transmembrane pores and corresponding plateaus in the free energy profiles. Allowing for movements of the protein and side chain snorkeling, we argue that the energetic cost for burying Arg inside a thin bilayer will be small, consistent with recent experiments, also leading to a dramatic reduction in pK_a shifts for Arg. We provide evidence that Arg translocation occurs via an ion-induced defect mechanism, except in thick bilayers (of at least 18 carbons) where solubility-diffusion becomes energetically favored. Our findings shed light on the mechanisms of ion movement through membranes of varying composition, with implications for a range of charged protein–lipid interactions and the actions of cell-perturbing peptides. This article is part of a Special Issue entitled: Membrane protein structure and function.     

On the role of anionic lipids in charged protein interactions with membranes

We investigate the role of anionic lipids in the binding to, and subsequent movement of charged protein groups in lipid membranes, to help understand the role of membrane composition in all membrane-active protein sequences. We demonstrate a small effect of phosphatidylglycerol (PG) lipids on the ability of an arginine (Arg) side chain to bind to, and cross a lipid membrane, despite possessing a neutralizing charge. We observe similar membrane deformations in lipid bilayers composed of phosphatidylcholine (PC) and PC/PG mixtures, with comparable numbers of water and lipid head groups pulled into the bilayer hydrocarbon core, and prohibitively large ~20 kcal/mol barriers for Arg transfer across each bilayer, dropping by just 2-3 kcal/mol due to the binding of PG lipids. We explore the causes of this small effect of introducing PG lipids and offer an explanation in terms of the limited membrane interaction for the choline groups of PC lipids bound to the translocating ion. Our calculations reveal a surprising lack of preference for Arg binding to PG lipids themselves, but a small increase in interfacial binding affinity for lipid bilayers containing PG lipids. These results help to explain the nature of competitive lipid binding to charged protein sequences, with implications for a wide range of membrane binding domains and cell perturbing peptides.     

Phospholipid binding of synthetic talin peptides provides evidence for an intrinsic membrane anchor of talin

Talin, an actin-binding protein, is assumed to anchor at the membrane via an intrinsic amino acid sequence. Three N-terminal talin fragments, 21-39 (S19), 287-304 (H18), and 385-406 (H17) have been proposed as potential membrane anchors. The interaction of the corresponding synthetic peptides with lipid model systems was investigated with CD spectroscopy, isothermal titration calorimetry, and monolayer expansion measurements. The membrane model systems were neutral or negatively charged small unilamellar vesicles or monolayers with a lateral packing density of bilayers (32 mN/m). S19 partitions into charged monolayers/bilayers with a penetration area A(p) = 140 +/- 30 A(2) and a free energy of binding of DeltaG(0) = -5.7 kcal/mol, thereby forming a partially alpha-helical structure. H18 does not interact with lipid monolayers or bilayers. H17 penetrates into neutral and charged monolayers/bilayers with A(p) = 148 +/- 23 A(2) and A(p) = 160 +/- 15 A(2), respectively, forming an alpha-helix in the membrane-bound state. Membrane partitioning is mainly entropy-driven. Under physiological conditions the free energy of binding to negatively charged membranes is DeltaG(0) = -9. 4 kcal/mol with a hydrophobic contribution of DeltaG(h) = -7.8 kcal/mol, comparable to that of post-translationally attached membrane anchors, and an electrostatic contribution of DeltaG(h) = -1.6 kcal/mol. The latter becomes more negative with decreasing pH. We show that H17 provides the binding energy required for a membrane anchor.     

Ion-Induced Defect Permeation of Lipid Membranes

We have explored the mechanisms of uncatalyzed membrane ion permeation using atomistic simulations and electrophysiological recordings. The solubility-diffusion mechanism of membrane charge transport has prevailed since the 1960s, despite inconsistencies in experimental observations and its lack of consideration for the flexible response of lipid bilayers. We show that direct lipid bilayer translocation of alkali metal cations, Cl^–, and a charged arginine side chain analog occurs via an ion-induced defect mechanism. Contrary to some previous suggestions, the arginine analog experiences a large free-energy barrier, very similar to those for Na⁺, K⁺, and Cl^–. Our simulations reveal that membrane perturbations, due to the movement of an ion, are central for explaining the permeation process, leading to both free-energy and diffusion-coefficient profiles that show little dependence on ion chemistry and charge, despite wide-ranging hydration energies and the membrane’s dipole potential. The results yield membrane permeabilities that are in semiquantitative agreement with experiments in terms of both magnitude and selectivity. We conclude that ion-induced defect-mediated permeation may compete with transient pores as the dominant mechanism of uncatalyzed ion permeation, providing new understanding for the actions of a range of membrane-active peptides and proteins.     

Ion-Induced Defect Permeation of Lipid Membranes

We have explored the mechanisms of uncatalyzed membrane ion permeation using atomistic simulations and electrophysiological recordings. The solubility-diffusion mechanism of membrane charge transport has prevailed since the 1960s, despite inconsistencies in experimental observations and its lack of consideration for the flexible response of lipid bilayers. We show that direct lipid bilayer translocation of alkali metal cations, Cl^–, and a charged arginine side chain analog occurs via an ion-induced defect mechanism. Contrary to some previous suggestions, the arginine analog experiences a large free-energy barrier, very similar to those for Na⁺, K⁺, and Cl^–. Our simulations reveal that membrane perturbations, due to the movement of an ion, are central for explaining the permeation process, leading to both free-energy and diffusion-coefficient profiles that show little dependence on ion chemistry and charge, despite wide-ranging hydration energies and the membrane’s dipole potential. The results yield membrane permeabilities that are in semiquantitative agreement with experiments in terms of both magnitude and selectivity. We conclude that ion-induced defect-mediated permeation may compete with transient pores as the dominant mechanism of uncatalyzed ion permeation, providing new understanding for the actions of a range of membrane-active peptides and proteins.     

pH-dependent self-association of influenza hemagglutinin fusion peptides in lipid bilayers

We have recently designed a host-guest peptide system that allows us to quantitatively measure the energetics of interaction of viral fusion peptides with lipid bilayers. Here, we show that fusion peptides of influenza hemagglutinin reversibly associate with one another at membrane surfaces above critical surface concentrations, which range from one to five peptides per 1000 lipids in the systems that we investigated. It is further demonstrated by using circular dichroism and Fourier transform infrared spectroscopy that monomeric peptides insert into the bilayers in a predominantly alpha-helical conformation, whereas self-associated fusion peptides adopt predominantly antiparallel beta-sheet structures at the membrane surface. The two forms are readily interconvertible and the equilibrium between them is determined by the pH and ionic strength of the surrounding solution. Lowering the pH favors the monomeric alpha-helical conformation, whereas increasing the ionic strength shifts the equilibrium towards the membrane-associated beta-aggregates. The binding data are interpreted in terms of a cooperative binding model that yields free energies of insertion and free energies of self-association for each of the peptides studied at pH 7.4 and pH 5. At pH 5 and 35 mM ionic strength, the insertion energy of the 20 residue influenza hemagglutinin fusion peptide is -7.2 kcal/mol and the self-association energy is -1.9 kcal/mol. We propose that self-association of fusion peptides could be a major driving force for recruiting a small number of hemagglutinin trimers into a fusion site.     

Chemical versus Mechanical Perturbations on the Protonation State of Arginine in Complex Lipid Membranes: Insights from Microscopic pKa Calculations

Charged amino acids such as Arginine play important roles in many membrane-mediated biological processes such as voltage gating of ion channels and membrane translocation of cell penetration peptides. It is well established that local membrane deformation and formation of water defects are crucial to the stabilization of charged species in contact with the membrane, which suggests that mechanical properties of the membrane are relevant although a clear connection has not been established. As a quantitative measure, we study how changes in the composition and therefore mechanical properties of a lipid bilayer influence the pK_a of Arg in the membrane center using free energy simulations. Compared to previous studies in a single-component lipid bilayer containing saturated lipids or lipids with a modest degree of unsaturation, substantially larger pK_a shifts are observed in the presence of highly unsaturated lipid tails and cholesterol. Moreover, the underlying molecular mechanisms for the pK_a perturbation are distinct in different systems, with the unsaturated lipid tails mainly destabilizing the charged state of Arg and the cholesterol stabilizing the neutral state of Arg. The observed behaviors in both cases are at odds with predictions based on mechanical considerations at a mesoscopic level—highlighting that, while mechanical considerations are useful for stimulating hypothesis, their applicability to dissecting phenomena at the molecular-length scale is rather limited.     

Arginine-Rich Peptides Destabilize the Plasma Membrane, Consistent with a Pore Formation Translocation Mechanism of Cell-Penetrating Peptides

Recent molecular-dynamics simulations have suggested that the arginine-rich HIV Tat peptides translocate by destabilizing and inducing transient pores in phospholipid bilayers. In this pathway for peptide translocation, Arg residues play a fundamental role not only in the binding of the peptide to the surface of the membrane, but also in the destabilization and nucleation of transient pores across the bilayer. Here we present a molecular-dynamics simulation of a peptide composed of nine Args (Arg-9) that shows that this peptide follows the same translocation pathway previously found for the Tat peptide. We test experimentally the hypothesis that transient pores open by measuring ionic currents across phospholipid bilayers and cell membranes through the pores induced by Arg-9 peptides. We find that Arg-9 peptides, in the presence of an electrostatic potential gradient, induce ionic currents across planar phospholipid bilayers, as well as in cultured osteosarcoma cells and human smooth muscle cells. Our results suggest that the mechanism of action of Arg-9 peptides involves the creation of transient pores in lipid bilayers and cell membranes.     

Outer membrane phospholipase A in phospholipid bilayers: a model system for concerted computational and experimental investigations of amino acid side chain partitioning into lipid bilayers

Understanding the forces that stabilize membrane proteins in their native states is one of the contemporary challenges of biophysics. To date, estimates of side chain partitioning free energies from water to the lipid environment show disparate values between experimental and computational measures. Resolving the disparities is particularly important for understanding the energetic contributions of polar and charged side chains to membrane protein function because of the roles these residue types play in many cellular functions. In general, computational free energy estimates of charged side chain partitioning into bilayers are much larger than experimental measurements. However, the lack of a protein-based experimental system that uses bilayers against which to vet these computational predictions has traditionally been a significant drawback. Moon & Fleming recently published a novel hydrophobicity scale that was derived experimentally by using a host-guest strategy to measure the side chain energetic perturbation due to mutation in the context of a native membrane protein inserted into a phospholipid bilayer. These values are still approximately an order of magnitude smaller than computational estimates derived from molecular dynamics calculations from several independent groups. Here we address this discrepancy by showing that the free energy differences between experiment and computation become much smaller if the appropriate comparisons are drawn, which suggests that the two fields may in fact be converging. In addition, we present an initial computational characterization of the Moon & Fleming experimental system used for the hydrophobicity scale: OmpLA in DLPC bilayers. The hydrophobicity scale used OmpLA position 210 as the guest site, and our preliminary results demonstrate that this position is buried in the center of the DLPC membrane, validating its usage in the experimental studies. We further showed that the introduction of charged Arg at position 210 is well tolerated in OmpLA and that the DLPC bilayers accommodate this perturbation by creating a water dimple that allows the Arg side chain to remain hydrated. Lipid head groups visit the dimple and can hydrogen bond with Arg, but these interactions are transient. Overall, our study demonstrates the unique advantages of this molecular system because it can be interrogated by both computational and experimental practitioners, and it sets the stage for free energy calculations in a system for which there is unambiguous experimental data. This article is part of a Special Issue entitled: Membrane protein structure and function.     

The membrane dipole potential in a total membrane potential model. Applications to hydrophobic ion interactions with membranes.

The total potential energy profile for hydrophobic ion interactions with lipid bilayers can be written as the sum of four terms: the electrical Born, image and dipole contributions, and a neutral energy term. We introduce a specific model for the membrane dipole potential, treating it as a two-dimensional array of point dipoles located near each membrane-water interface. Together with specific theoretical models for the other energy terms, a total potential profile is developed that successfully describes the complete set of thermodynamic parameters for binding and translocation for the two hydrophobic ion structural analogues, tetraphenylphosphonium (TPP+) and tetraphenylboron (TPB-). A reasonable fit to the data is possible if the dipole potential energy has a magnitude of 5.5 + 0.5 kcal/mol (240 + 20 mV), positive inside, and if the neutral energy contribution for TPP+ and TPB- is -7.0 + 1.0 kcal/mol. These results may also have important implications for small ion interactions with membranes and the energetics of charged groups in membrane proteins.     

Influence of cation-pi interactions in different folding types of membrane proteins

Cation-pi interactions play an important role to the stability of protein structures. In this work, we analyze the influence of cation-pi interactions in three-dimensional structures of membrane proteins. We found that transmembrane strand (TMS) proteins have more number of cation-pi interactions than transmembrane helical (TMH) proteins. In TMH proteins, both the positively charged residues Lys and Arg equally experience favorable cation-pi interactions whereas in TMS proteins, Arg is more likely than Lys to be in such interactions. There is no relationship between number of cation-pi interactions and number of residues in TMH proteins whereas a good correlation was observed in TMS proteins. The average cation-pi interaction energy for TMH proteins is -16 kcal/mol and that for TMS proteins is -27 kcal/mol. The pair-wise cation-pi interaction energy between aromatic and positively charged residues showed that Lys-Trp energy is stronger in TMS proteins than TMH proteins; Arg-Phe, Arg-Tyr and Lys-Phe have higher energy in TMH proteins than TMS proteins. The decomposition of energies into electrostatic and van der Waals revealed that the contribution from electrostatic energy is twice as that from van der Waals energy in both TMH and TMS proteins. The results obtained in the present study would be helpful to understand the contribution of cation-pi interactions to the stability of membrane proteins.     

Lipopolysaccharides in bacterial membranes act like cholesterol in eukaryotic plasma membranes in providing protection against melittin-induced bilayer lysis

Melittin is a small, cationic peptide that, like many other antimicrobial peptides, lyses cell membranes by acting on their lipid bilayers. However, the sensitivity to antimicrobial peptides varies among cell types. We have performed direct binding and vesicle leakage experiments to determine the sensitivity to melittin of bilayers composed of various physiologically relevant lipids, in particular, key components of eukaryotic membranes (cholesterol) and bacterial outer membranes (lipopolysaccharide or LPS). Melittin binds to bilayers composed of both zwitterionic and negatively charged phospholipids, as well as to the highly charged LPS bilayers. The magnitude of the free energy of binding (deltaG degrees ) increases with increasing bilayer charge density; deltaG degrees = -7.6 kcal/mol for phosphatidylcholine (PC) bilayers and -8.9 to -11.0 kcal/mol for negatively charged bilayers containing phosphatidylserine (PS), phospholipids with covalently attached polyethylene glycol (PEG-lipids), or LPS. Comparisons of these data show that binding is not markedly affected by the steric barrier produced by the PEG in PEG-lipids or by the polysaccharide core of LPS. The addition of equimolar cholesterol to PC bilayers reduces the level of binding (deltaG degrees = -6.4 kcal/mol) and reduces the extent of melittin-induced leakage by 20-fold. LPS and 1:1 PC/cholesterol bilayers have similar high resistance to melittin-induced leakage, indicating that cholesterol in eukaryotic plasma membranes and LPS in Gram-negative bacteria provide strong protection against the lytic effects of melittin. We argue that this resistance is due at least in part to the similar tight packing of the lipid acyl chains in PC/cholesterol and LPS bilayers. The addition of bacterial phospholipids to LPS bilayers increases their sensitivity to melittin, helping to explain the higher sensitivity of deep rough bacteria compared to smooth phenotypes.     

Binding of small basic peptides to membranes containing acidic lipids: theoretical models and experimental results.

We measured directly the binding of Lys3, Lys5, and Lys7 to vesicles containing acidic phospholipids. When the vesicles contain 33% acidic lipids and the aqueous solution contains 100 mM monovalent salt, the standard Gibbs free energy for the binding of these peptides is 3, 5, and 7 kcal/mol, respectively. The binding energies decrease as the mol% of acidic lipids in the membrane decreases and/or as the salt concentration increases. Several lines of evidence suggest that these hydrophilic peptides do not penetrate the polar headgroup region of the membrane and that the binding is mainly due to electrostatic interactions. To calculate the binding energies from classical electrostatics, we applied the nonlinear Poisson-Boltzmann equation to atomic models of the phospholipid bilayers and the basic peptides in aqueous solution. The electrostatic free energy of interaction, which arises from both a long-range coulombic attraction between the positively charged peptide and the negatively charged lipid bilayer, and a short-range Born or image charge repulsion, is a minimum when approximately 2.5 A (i.e., one layer of water) exists between the van der Waals surfaces of the peptide and the lipid bilayer. The calculated molar association constants, K, agree well with the measured values: K is typically about 10-fold smaller than the experimental value (i.e., a difference of about 1.5 kcal/mol in the free energy of binding). The predicted dependence of K (or the binding free energies) on the ionic strength of the solution, the mol% of acidic lipids in the membrane, and the number of basic residues in the peptide agree very well with the experimental measurements. These calculations are relevant to the membrane binding of a number of important proteins that contain clusters of basic residues.     

Lipid bilayer disruption by oligomeric α-synuclein depends on bilayer charge and accessibility of the hydrophobic core

Soluble oligomeric aggregates of α-synuclein have been implicated to play a central role in the pathogenesis of Parkinson's disease. Disruption and permeabilization of lipid bilayers by α-synuclein oligomers is postulated as a toxic mechanism, but the molecular details controlling the oligomer–membrane interaction are still unknown. Here we show that membrane disruption strongly depends on the accessibility of the hydrophobic membrane core and that charge interactions play an important but complex role. We systematically studied the influence of the physical membrane properties and solution conditions on lipid bilayer disruption by oligomers using a dye release assay. Varying the lipid headgroup composition revealed that membrane disruption only occurs for negatively charged bilayers. Furthermore, the electrostatic repulsion between the negatively charged α-synuclein and the negative surface charge of the bilayer inhibits vesicle disruption at low ionic strength. The disruption of negatively charged vesicles further depends on lipid packing parameters. Bilayer composition changes that result in an increased lipid headgroup spacing make vesicles more prone to disruption, suggesting that the accessibility of the bilayer hydrocarbon core modulates oligomer–membrane interaction. These data shed important new insights into the driving forces governing the highly debated process of oligomer–membrane interactions.     

Thermodynamics of fusion peptide-membrane interactions

The fusion peptides of viral membrane fusion proteins play a key role in the mechanism of viral spike glycoprotein mediated membrane fusion. These peptides insert into the lipid bilayers of cellular target membranes where they adopt mostly helical secondary structures. To better understand how membranes may be converted to high-energy intermediates during fusion, it is of interest to know how much energy, enthalpy and entropy, is provided by the insertion of fusion peptides into lipid bilayers. Here, we describe a detailed thermodynamic analysis of the binding of analogues of the influenza hemagglutinin fusion peptide of different lengths and amino acid compositions. In small unilamellar vesicles, the interaction of these peptides with lipid bilayers is driven by enthalpy (-16.5 kcal/mol) and opposed by entropy (-30 cal mol(-1) K(-1)). Most of the driving force (deltaG = -7.6 kcal/mol) comes from the enthalpy of peptide insertion deep into the lipid bilayer. Enthalpic gains and entropic losses of peptide folding in the lipid bilayer cancel to a large extent and account for only about 40% of the total binding free energy. The major folding event occurs in the N-terminal segment of the fusion peptide. The C-terminal segment mainly serves to drive the N-terminus deep into the membrane. The fusion-defective mutations G1S, which causes hemifusion, and particularly G1V, which blocks fusion, have major structural and thermodynamic consequences on the insertion of fusion peptides into lipid bilayers. The magnitudes of the enthalpies and entropies of binding of these mutant peptides are reduced, their helix contents are reduced, but their energies of self-association at the membrane surface are increased compared to the wild-type fusion peptide.     

Is arginine charged in a membrane?

“Charged” amino acids play countless important roles in protein structure and function. Yet when these side chains come into contact with membranes we do not fully understand their behavior. This is highlighted by a recent model of voltage-gated ion channel activity and translocon-based experiments that suggest small penalties to expose these side chains to lipids, opposing the prevailing view in membrane biophysics. Here we employ a side chain analog as well as a transmembrane helix model to determine the free energy as a function of protonation state and position for a lipid-exposed arginine (Arg) residue across a membrane. We observe high free energy barriers for both the charged and neutral states. Due to the stabilizing influence of membrane deformations for the protonated form, the Arg side chain experiences a pK_a shift of ≤4.5 units and remains mostly protonated. The cost for exposing Arg to lipid hydrocarbon is prohibitively high with implications for many membrane translocating processes and the activation mechanisms of voltage-gated ion channels.     

Membrane thickening by the antimicrobial peptide PGLa

Using x-ray diffraction, solid-state ²H-NMR, differential scanning calorimetry, and dilatometry, we have observed a perturbation of saturated acyl chain phosphatidylglycerol bilayers by the antimicrobial peptide peptidyl-glycylleucine-carboxyamide (PGLa) that is dependent on the length of the hydrocarbon chain. In the gel phase, PGLa induces a quasi-interdigitated phase, previously reported also for other peptides, which is most pronounced for C18 phosphatidylglycerol. In the fluid phase, we found an increase of the membrane thickness and NMR order parameter for C14 and C16 phosphatidylglycerol bilayers, though not for C18. The data is best understood in terms of a close hydrophobic match between the C18 bilayer core and the peptide length when PGLa is inserted with its helical axis normal to the bilayer surface. The C16 acyl chains appear to stretch to accommodate PGLa, whereas tilting within the bilayer seems to be energetically favorable for the peptide when inserted into bilayers of C14 phosphatidylglycerol. In contrast to the commonly accepted membrane thinning effect of antimicrobial peptides, the data demonstrate that pore formation does not necessarily relate to changes in the overall bilayer structure.     

Dye-release assay for investigation of antimicrobial peptide activity in a competitive lipid environment

A dye-release method for investigating the effect of a competitive lipid environment on the activity of two membrane-disrupting antimicrobial peptides (AMP), maculatin 1.1 and aurein 1.2, is presented. The results support the general conclusion that AMP have greater affinity for negatively charged membranes, for example bacterial membranes, than for the neutral membrane surface found in eukaryotic cells, but only within a competitive lipid environment. Indeed, in a single-model membrane environment, both peptides were more potent against neutral vesicles than against charged vesicles. The approach was also used to investigate the effect of pre-incubating the peptides in a neutral lipid environment then introducing charged lipid vesicles. Maculatin was shown to migrate from the neutral lipid bilayers, where pores had already formed, to the charged membrane bilayers. This result was also observed for charged-to-charged bilayers but, interestingly, not for neutral-to-neutral lipid interfaces. Aurein was able to migrate from either lipid environment, indicating weaker binding to lipid membranes, and a different molecular mechanism for lysis of lipid bilayers. Competitive lipid environments could be used to assess other critical conditions that modulate the activity of membrane peptides or proteins.     

Implicit membrane treatment of buried charged groups: Application to peptide translocation across lipid bilayers

The energetic cost of burying charged groups in the hydrophobic core of lipid bilayers has been controversial, with simulations giving higher estimates than certain experiments. Implicit membrane approaches are usually deemed too simplistic for this problem. Here we challenge this view. The free energy of transfer of amino acid side chains from water to the membrane center predicted by IMM1 is reasonably close to all-atom free energy calculations. The shape of the free energy profile, however, for the charged side chains needs to be modified to reflect the all-atom simulation findings (IMM1-LF). Membrane thinning is treated by combining simulations at different membrane widths with an estimate of membrane deformation free energy from elasticity theory. This approach is first tested on the voltage sensor and the isolated S4 helix of potassium channels. The voltage sensor is stably inserted in a transmembrane orientation for both the original and the modified model. The transmembrane orientation of the isolated S4 helix is unstable in the original model, but a stable local minimum in IMM1-LF, slightly higher in energy than the interfacial orientation. Peptide translocation is addressed by mapping the effective energy of the peptide as a function of vertical position and tilt angle, which allows identification of minimum energy pathways and transition states. The barriers computed for the S4 helix and other experimentally studied peptides are low enough for an observable rate. Thus, computational results and experimental studies on the membrane burial of peptide charged groups appear to be consistent. This article is part of a Special Issue entitled: Interfacially Active Peptides and Proteins. Guest Editors: William C. Wimley and Kalina Hristova.