Influence of the membrane dipole potential on peptide binding to lipid bilayers

Cholesterol Changes Interfacial Water Alignment in Model Cell Membranes

The nanoscopic layer of water that directly hydrates biological membranes plays a critical role in maintaining the cell structure, regulating biochemical processes, and managing intermolecular interactions at the membrane interface. Therefore, comprehending the membrane structure, including its hydration, is essential for understanding the chemistry of life. While cholesterol is a fundamental lipid molecule in mammalian cells, influencing both the structure and dynamics of cell membranes, its impact on the structure of interfacial water has remained unknown. We used surface-specific vibrational sum-frequency generation spectroscopy to study the effect of cholesterol on the structure and hydration of monolayers of the lipids 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC), and egg sphingomyelin (SM). We found that for the unsaturated lipid DOPC, cholesterol intercalates in the membrane without significantly changing the orientation of the lipid tails and the orientation of the water molecules hydrating the headgroups of DOPC. In contrast, for the saturated lipids DPPC and SM, the addition of cholesterol leads to clearly enhanced packing and ordering of the hydrophobic tails. It is also observed that the orientation of the water hydrating the lipid headgroups is enhanced upon the addition of cholesterol. These results are important because the orientation of interfacial water molecules influences the cell membranes' dipole potential and the strength and specificity of interactions between cell membranes and peripheral proteins and other biomolecules. The lipid nature-dependent role of cholesterol in altering the arrangement of interfacial water molecules offers a fresh perspective on domain-selective cellular processes, such as protein binding.

Ebola Virus Glycoprotein Strongly Binds to Membranes in the Absence of Receptor Engagement

Ebola virus (EBOV) is an enveloped virus that must fuse with the host cell membrane in order to release its genome and initiate infection. This process requires the action of the EBOV envelope glycoprotein (GP), encoded by the virus, which resides in the viral envelope and consists of a receptor binding subunit, GP1, and a membrane fusion subunit, GP2. Despite extensive research, a mechanistic understanding of the viral fusion process is incomplete. To investigate GP-membrane association, a key step in the fusion process, we used two approaches: high-throughput measurements of single-particle diffusion and single-molecule measurements with optical tweezers. Using these methods, we show that the presence of the endosomal Niemann-Pick C1 (NPC1) receptor is not required for primed GP-membrane binding. In addition, we demonstrate this binding is very strong, likely attributed to the interaction between the GP fusion loop and the membrane's hydrophobic core. Our results also align with previously reported findings, emphasizing the significance of acidic pH in the protein-membrane interaction. Beyond Ebola virus research, our approach provides a powerful toolkit for studying other protein-membrane interactions, opening new avenues for a better understanding of protein-mediated membrane fusion events.

Peptide Flexibility and the Hydrophobic Moment are Determinants to Evaluate the Clinical Potential of Magainins

Using a flexibility prediction algorithm and in silico structural modeling, we have calculated the intrinsic flexibility of several magainin derivatives. In the case of magainin-2 (Mag-2) and magainin H2 (MAG-H2) we have found that MAG-2 is more flexible than its hydrophobic analog, Mag-H2. This affects the degree of bending of both peptides, with a kink around two central residues (R10, R11), whereas, in Mag-H2, W10 stiffens the peptide. Moreover, this increases the hydrophobic moment of Mag-H2, which could explain its propensity to form pores in POPC model membranes, which exhibit near-to-zero spontaneous curvatures. Likewise, the protective effect described in DOPC membranes for this peptide regarding its facilitation in pore formation would be related to the propensity of this lipid to form membranes with negative spontaneous curvature. The flexibility of another magainin analog (MSI-78) is even greater than that of Mag-2. This facilitates the peptide to present a kind of hinge around the central F12 as well as a C-terminal end prone to be disordered. Such characteristics are key to understanding the broad-spectrum antimicrobial actions exhibited by this peptide. These data reinforce the hypothesis on the determinant role of spontaneous membrane curvature, intrinsic peptide flexibility, and specific hydrophobic moment in assessing the bioactivity of membrane-active antimicrobial peptides.

A physical perspective to understand myelin II: The physical origin of myelin development

The physical principle of myelin development is obtained from our previous study by explaining Peter’s quadrant mystery: an externally applied negative and positive E-field can promote and inhibit the growth of the inner tongue of the myelin sheath, respectively. In this study, this principle is considered as a fundamental hypothesis, named Hypothesis-E, to explain more phenomena about myelin development systematically. Specifically, the g-ratio and the fate of the Schwann cell’s differentiation are explained in terms of the E-field. Moreover, an experiment is proposed to validate this theory.

COGRIMEN: Coarse-Grained Method for Modeling of Membrane Proteins in Implicit Environments

The presented methodology is based on coarse-grained representation of biomolecules in implicit environments and is designed for the molecular dynamics simulations of membrane proteins and their complexes. The membrane proteins are not only found in the cell membrane but also in all membranous compartments of the cell: Golgi apparatus, mitochondria, endosomes and lysosomes, and they usually form large complexes. To investigate such systems the methodology is proposed based on two independent approaches combining the coarse-grained MARTINI model for proteins and the effective energy function to mimic the water/membrane environments. The latter is based on the implicit environment developed for all-atom simulations in the IMM1 method. The force field solvation parameters for COGRIMEN were initially calculated from IMM1 all-atom parameters and then optimized using Genetic Algorithms. The new methodology was tested on membrane proteins, their complexes and oligomers. COGRIMEN method is implemented as a patch for NAMD program and can be useful for fast and brief studies of large membrane protein complexes.

Structure-activity relationships of mitochondria-targeted tetrapeptide pharmacological compounds

Mitochondria play a central role in metabolic homeostasis, and dysfunction of this organelle underpins the etiology of many heritable and aging-related diseases. Tetrapeptides with alternating cationic and aromatic residues such as SS-31 (elamipretide) show promise as therapeutic compounds for mitochondrial disorders. In this study, we conducted a quantitative structure-activity analysis of three alternative tetrapeptide analogs, benchmarked against SS-31, that differ with respect to aromatic side chain composition and sequence register. We present the first structural models for this class of compounds, obtained with Nuclear Magnetic Resonance (NMR) and molecular dynamics approaches, showing that all analogs except for SS-31 form compact reverse turn conformations in the membrane-bound state. All peptide analogs bound cardiolipin-containing membranes, yet they had significant differences in equilibrium binding behavior and membrane interactions. Notably, analogs had markedly different effects on membrane surface charge, supporting a mechanism in which modulation of membrane electrostatics is a key feature of their mechanism of action. The peptides had no strict requirement for side chain composition or sequence register to permeate cells and target mitochondria in mammalian cell culture assays. All four peptides were pharmacologically active in serum withdrawal cell stress models yet showed significant differences in their abilities to restore mitochondrial membrane potential, preserve ATP content, and promote cell survival. Within our peptide set, the analog containing tryptophan side chains, SPN10, had the strongest impact on most membrane properties and showed greatest efficacy in cell culture studies. Taken together, these results show that side chain composition and register influence the activity of these mitochondria-targeted peptides, helping provide a framework for the rational design of next-generation therapeutics with enhanced potency.

Membrane Dipole Potential: An Emerging Approach to Explore Membrane Organization and Function

Biological membranes are complex organized molecular assemblies of lipids and proteins that provide cells and membrane-bound intracellular organelles their individual identities by morphological compartmentalization. Membrane dipole potential originates from the electrostatic potential difference within the membrane due to the nonrandom arrangement (orientation) of amphiphile and solvent (water) dipoles at the membrane interface. In this Feature Article, we will focus on the measurement of dipole potential using electrochromic fluorescent probes and highlight interesting applications. In addition, we will focus on ratiometric fluorescence microscopic imaging technique to measure dipole potential in cellular membranes, a technique that can be used to address novel problems in cell biology which are otherwise difficult to address using available approaches. We envision that membrane dipole potential could turn out to be a convenient tool in exploring the complex interplay between membrane lipids and proteins and could provide novel insights in membrane organization and function.

Extension of the CAVS model to the simulation of helical peptides in a membrane environment

Considering the effect of peptide insertion on the dipole potential of the lipid membrane, we extend the CAVS coarse-grained (CG) model to the simulation of helical peptides in a membrane environment. In this approach, the CG scheme for a peptide backbone is similar to the treatment in the united-atom model, while we treated the side chain of an amino acid by grouping 1-3 heavy atoms into a CG unit. The CAVS CG force field for peptides is optimized by reproducing the experimental results for the backbone (φ, ψ) distribution and predicting the PMF profiles of transferring organic molecules in a lipid bilayer membrane obtained from all-atom simulations. The CAVS simulation of a helical peptide in a phosphatidylcholine (PC) lipid bilayer revealed that the insertion of a peptide increases the dipole potential of the PC lipid bilayer, in which the peptide and its neutralized ions make a significant contribution. Finally, we carried out the CAVS simulation for five different helical peptides in the PC lipid bilayer to explore the behavior of peptide tilt, showing excellent agreement with the all-atom simulations. Our work suggests that the peptide tilt should relieve the deformation stress from the lipid bilayer, and the peptide aggregation could reduce the peptide tilt by resisting the deformation stress from the surrounding lipids.

The effect of acidic pH on the adsorption and lytic activity of the peptides Polybia-MP1 and its histidine-containing analog in anionic lipid membrane: a biophysical study by molecular dynamics and spectroscopy

Antimicrobial peptides (AMPs) are part of the innate immune system of many species. AMPs are short sequences rich in charged and non-polar residues. They act on the lipid phase of the plasma membrane without requiring membrane receptors. Polybia-MP1 (MP1), extracted from a native wasp, is a broad-spectrum bactericide, an inhibitor of cancer cell proliferation being non-hemolytic and non-cytotoxic. MP1 mechanism of action and its adsorption mode is not yet completely known. Its adsorption to lipid bilayer and lytic activity is most likely dependent on the ionization state of its two acidic and three basic residues and consequently on the bulk pH. Here we investigated the effect of bulk acidic (pH 5.5) and neutral pH (7.4) solution on the adsorption, insertion, and lytic activity of MP1 and its analog H-MP1 to anionic (7POPC:3POPG) model membrane. H-MP1 is a synthetic analog of MP1 with lysines replaced by histidines. Bulk pH changes could modulate this peptide efficiency. The combination of different experimental techniques and molecular dynamics (MD) simulations showed that the adsorption, insertion, and lytic activity of H-MP1 are highly sensitive to bulk pH in opposition to MP1. The atomistic details, provided by MD simulations, showed peptides contact their N-termini to the bilayer before the insertion and then lay parallel to the bilayer. Their hydrophobic faces inserted into the acyl chain phase disturb the lipid-packing.

An ω-3, but Not an ω-6 Polyunsaturated Fatty Acid Decreases Membrane Dipole Potential and Stimulates Endo-Lysosomal Escape of Penetratin

Although the largely positive intramembrane dipole potential (DP) may substantially influence the function of transmembrane proteins, its investigation is deeply hampered by the lack of measurement techniques suitable for high-throughput examination of living cells. Here, we describe a novel emission ratiometric flow cytometry method based on F66, a 3-hydroxiflavon derivative, and demonstrate that 6-ketocholestanol, cholesterol and 7-dehydrocholesterol, saturated stearic acid (SA) and ω-6 γ-linolenic acid (GLA) increase, while ω-3 α-linolenic acid (ALA) decreases the DP. These changes do not correlate with alterations in cell viability or membrane fluidity. Pretreatment with ALA counteracts, while SA or GLA enhances cholesterol-induced DP elevations. Furthermore, ALA (but not SA or GLA) increases endo-lysosomal escape of penetratin, a cell-penetrating peptide. In summary, we have developed a novel method to measure DP in large quantities of individual living cells and propose ALA as a physiological DP lowering agent facilitating cytoplasmic entry of penetratin.

Lipid Membrane Interaction of Peptide/DNA Complexes Designed for Gene Delivery

Among gene delivery systems, peptide-based gene carriers have received significant attention because of their selectivity, biocompatibility, and biodegradability. Since cellular membranes function as a barrier toward exogenous molecules, cell-penetrating peptides (CPPs), which are usually cationic and/or amphiphilic, can serve as efficient carriers to deliver cargo into the cytosol. Here, we examined the interactions of carrier peptides and their DNA complexes with lipid membranes using a quartz crystal microbalance (QCM) and high-speed atomic force microscopy (HS-AFM). The carrier peptides are a 12-residue partial presequence of yeast cytochrome c oxidase subunit IV (Cytcox) and BP100, which are a mitochondria-targeting signal peptide and a CPP, respectively. QCM data showed that BP100 has a higher binding affinity than Cytcox to both plasma membrane- and mitochondrial membrane-mimicking lipid bilayers. The DNA complexes with either Cytcox or BP100 exhibited the same tendency. Furthermore, HS-AFM data demonstrated that the DNA complexes of either peptide can disrupt the lipid membranes, forming larger pores in the case of Cytcox. Our results suggest that the binding affinity of the peptide/DNA complex to the plasma membrane is more critical than its membrane disruption ability in enhancing the cellular uptake of DNA.

Computational studies of peptide-induced membrane pore formation

A variety of peptides induce pores in biological membranes; the most common ones are naturally produced antimicrobial peptides (AMPs), which are small, usually cationic, and defend diverse organisms against biological threats. Because it is not possible to observe these pores directly on a molecular scale, the structure of AMP-induced pores and the exact sequence of steps leading to their formation remain uncertain. Hence, these questions have been investigated via molecular modelling. In this article, we review computational studies of AMP pore formation using all-atom, coarse-grained, and implicit solvent models; evaluate the results obtained and suggest future research directions to further elucidate the pore formation mechanism of AMPs.This article is part of the themed issue 'Membrane pores: from structure and assembly, to medicine and technology'.

Effect of N-terminal acetylation on lytic activity and lipid-packing perturbation induced in model membranes by a mastoparan-like peptide

L1A (IDGLKAIWKKVADLLKNT-NH2) is a peptide that displays a selective antibacterial activity to Gram-negative bacteria without being hemolytic. Its lytic activity in anionic lipid vesicles was strongly enhanced when its N-terminus was acetylated (ac-L1A). This modification seems to favor the perturbation of the lipid core of the bilayer by the peptide, resulting in higher membrane lysis. In the present study, we used lipid monolayers and bilayers as membrane model systems to explore the impact of acetylation on the L1A lytic activity and its correlation with lipid-packing perturbation. The lytic activity investigated in giant unilamellar vesicles (GUVs) revealed that the acetylated peptide permeated the membrane at higher rates compared with L1A, and modified the membrane's mechanical properties, promoting shape changes. The peptide secondary structure and the changes in the environment of the tryptophan upon adsorption to large unilamellar vesicles (LUVs) were monitored by circular dichroism (CD) and red-edge excitation shift experiments (REES), respectively. These experiments showed that the N-terminus acetylation has an important effect on both, peptide secondary structure and peptide insertion into the bilayer. This was also confirmed by experiments of insertion into lipid monolayers. Compression isotherms for peptide/lipid mixed films revealed that ac-L1A dragged lipid molecules to the more disordered phase, generating a more favorable environment and preventing the lipid molecules from forming stiff films. Enthalpy changes in the main phase transition of the lipid membrane upon peptide insertion suggested that the acetylated peptide induced higher impact than the non-acetylated one on the thermotropic behavior of anionic vesicles.

2H-NMR and MD Simulations Reveal Membrane-Bound Conformation of Magainin 2 and Its Synergy with PGLa

Magainin 2 (MAG2) and PGLa are two α-helical antimicrobial peptides found in the skin of the African frog Xenopus laevis. They act by permeabilizing bacterial membranes and exhibit an exemplary synergism. Here, we determined the detailed molecular alignment and dynamical behavior of MAG2 in oriented lipid bilayers by using ²H-NMR on Ala-d₃-labeled peptides, which yielded orientation-dependent quadrupolar splittings of the labels. The amphiphilic MAG2 helix was found to lie flat on the membrane surface in 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC)/1,2-dimyristoyl-sn-glycero-3-phosphatidylglycerol (DMPG) and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC)/1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylglycerol (POPG), as expected, with a tilt angle close to 90°. This orientation fits well with all-atom molecular-dynamics simulations of MAG2 performed in DMPC and DMPC/DMPG. In the presence of an equimolar amount of PGLa, the NMR analysis showed that MAG2 becames tilted at an angle of 120°, and its azimuthal rotation angle also changes. Since this interaction was found to occur in a concentration range where the peptides per se do not interact with their own type, we propose that MAG2 forms a stable heterodimer with PGLa. Given that the PGLa molecules in the complex are known to be flipped into a fully upright orientation, with a helix tilt close to 180°, they must make up the actual transmembrane pore. We thus suggest that the two negative charges on the C-terminus of the obliquely tilted MAG2 peptides neutralize some of the cationic groups on the upright PGLa helices. This would stabilize the assembly of PGLa into a toroidal pore with an overall reduced charge density, which could explain the mechanism of synergy.

Molecular dynamics simulations of biological membranes and membrane proteins using enhanced conformational sampling algorithms

This paper reviews various enhanced conformational sampling methods and explicit/implicit solvent/membrane models, as well as their recent applications to the exploration of the structure and dynamics of membranes and membrane proteins. Molecular dynamics simulations have become an essential tool to investigate biological problems, and their success relies on proper molecular models together with efficient conformational sampling methods. The implicit representation of solvent/membrane environments is reasonable approximation to the explicit all-atom models, considering the balance between computational cost and simulation accuracy. Implicit models can be easily combined with replica-exchange molecular dynamics methods to explore a wider conformational space of a protein. Other molecular models and enhanced conformational sampling methods are also briefly discussed. As application examples, we introduce recent simulation studies of glycophorin A, phospholamban, amyloid precursor protein, and mixed lipid bilayers and discuss the accuracy and efficiency of each simulation model and method. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Free Energy Calculations for the Peripheral Binding of Proteins/Peptides to an Anionic Membrane. 1. Implicit Membrane Models

The binding of peptides and proteins to the surface of complex lipid membranes is important in many biological processes such as cell signaling and membrane remodeling. Computational studies can aid experiments by identifying physical interactions and structural motifs that determine the binding affinity and specificity. However, previous studies focused on either qualitative behaviors of protein/membrane interactions or the binding affinity of small peptides. Motivated by this observation, we set out to develop computational protocols for bimolecular binding to charged membranes that are applicable to both peptides and large proteins. In this work, we explore a method based on an implicit membrane/solvent model (generalized Born with a simple switching in combination with the Gouy-Chapman-Stern model for a charged interface), which we expect to lead to useful results when the binding does not implicate significant membrane deformation and local demixing of lipids. We show that the binding free energy can be efficiently computed following a thermodynamic cycle similar to protein-ligand binding calculations, especially when a Bennett acceptance ratio based protocol is used to consider both the membrane bound and solution conformational ensembles. Test calculations on a series of peptides show that our computational approach leads to binding affinities in encouraging agreement with experimental data, including for the challenging example of the bringing of flexible MARCKS-ED peptides to membranes. The calculations highlight that for a membrane with a significant fraction of anionic lipids, it is essential to include the effect of ion adsorption using the Stern model, which significantly modifies the effective surface charge. This implicit membrane model based computational protocol helps lay the groundwork for more systematic analysis of protein/peptide binding to membranes of complex shape and composition.

Effects of Dipole Potential Modifiers on Heterogenic Lipid Bilayers

In this work, we examine the ability of dipole modifiers, flavonoids, and RH dyes to affect the dipole potential (φ d) and phase separation in membranes composed of ternary mixtures of POPC with different sphingolipids and sterols. Changes in the steady-state conductance induced by cation-ionophore complexes have been measured to evaluate the changes in dipole potential of planar lipid bilayers. Confocal fluorescence microscopy has been employed to investigate lipid segregation in giant unilamellar vesicles. The effects of flavonoids on φ d depend on lipid composition and dipole modifier type. The effectiveness of RH dyes to increase φ d depends on sphingolipid type but is not influenced by sterol content. Tested modifiers lead to partial or complete disruption of gel domains in bilayers composed of POPC, sphingomyelin, and cholesterol. Substitution of cholesterol to ergosterol or 7-dehydrocholesterol leads to a loss of fluidizing effects of modifiers except phloretin. This may be due to various compositions of gel domains. The lack of influence of modifiers on phase scenario in vesicles composed of ternary mixtures of POPC, cholesterol, and phytosphingosine or sphinganine is related to an absence of gel-like phase. It was concluded that the membrane lateral heterogeneity affects the dipole-modifying abilities of the agents that influence the magnitude of φ d by intercalation into the bilayer and orientation of its own large dipole moments (phloretin and RH dyes). The efficacy of modifiers that do not penetrate deeply and affect φ d through water adsorption (phlorizin, quercetin, and myricetin) is not influenced by lateral heterogeneity of membrane.

Modifiers of membrane dipole potentials as tools for investigating ion channel formation and functioning

Electrostatic fields generated on and within biological membranes play a fundamental role in key processes in cell functions. The role of the membrane dipole potential is of particular interest because of its powerful impact on membrane permeability and lipid-protein interactions, including protein insertion, oligomerization, and function. The membrane dipole potential is defined by the orientation of electric dipoles of lipid headgroups, fatty acid carbonyl groups, and membrane-adsorbed water. As a result, the membrane interior is several hundred millivolts more positive than the external aqueous phase. This potential decrease depends on the lipid, and especially sterol, composition of the membrane. The adsorption of certain electroneutral molecules known as dipole modifiers may also lead to significant changes in the magnitude of the potential decrease. These agents are widely used to study the effects of the dipole potential on membrane transport. This review presents a critical analysis of a variety of data from studies dedicated to ion channel formation and functioning in membranes with different dipole potentials. The types of ion channels found in cellular membranes and pores formed by antimicrobial agents and toxins in artificial lipid membranes are summarized. The mechanisms underlying the influence of the membrane dipole potential on ion channel activity, including dipole-dipole and charge-dipole interactions in the pores and in membranes, are discussed. A hypothesis, in which lipid rafts in both model and cellular membranes also modulate ion channel activity by virtue of an increased or decreased dipole potential, is also considered.

Implicit Membrane Investigation of the Stability of Antimicrobial Peptide β-Barrels and Arcs

Previous simulations showed that the β-hairpin antimicrobial peptide (AMP) protegrin-1 can form stable octameric β-barrels and tetrameric arcs (half barrels) in both implicit and explicit membranes. Here, we extend this investigation to several AMPs of similar structure: tachyplesin, androctonin, polyphemusin, gomesin, and the retrocyclin θ-defensin. These peptides form short β-hairpins stabilized by 2-3 disulfide bonds. We also examine synthetic β-sheet peptides selected from a combinatorial library for their ability or inability to form pores in lipid membranes. When heptameric, octameric, and decameric β-barrels and tetrameric arcs of these peptides were embedded in pre-formed neutral or anionic lipid pores (i.e., pores in neutral or anionic membranes, respectively), a variety of behaviors and membrane binding energies were observed. Due to the cationic charge of the peptides, more favorable transfer energies and more stable binding were observed in anionic than neutral pores. The synthetic peptides bound very strongly and formed stable barrels and arcs in both neutral and anionic pores. The natural AMPs exhibited unfavorable or marginally favorable binding energy and kinetic stability in neutral pores, consistent with the lower hemolytic activity of some of them compared with protegrin-1. Binding to anionic pores was more favorable, but significant distortions of the barrel or arc structures were sometimes noted. These results are discussed in light of the available experimental data. The diversity of behaviors obtained makes it unlikely that the barrel and arc mechanisms are valid for the entire family of β-hairpin AMPs.

Protein arcs may form stable pores in lipid membranes

Electron microscopy and atomic force microscopy images of cholesterol-dependent cytolysins and related proteins that form large pores in lipid membranes have revealed the presence of incomplete rings, or arcs. Some evidence indicates that these arcs are inserted into the membrane and induce membrane leakage, but other experiments seem to refute that. Could such pores, only partially lined by protein, be kinetically and thermodynamically stable? How would the lipids be structured in such a pore? Using the antimicrobial peptide protegrin-1 as a model, we test the stability of pores only partially lined by peptide using all-atom molecular dynamics simulations in POPC and POPE/POPG membranes. The data show that, whereas pure lipid pores close rapidly, pores partially lined by protegrin arcs are stable for at least 300 ns. Estimates of the thermodynamic stability of these arcs using line tension data and implicit solvent calculations show that these arcs can be marginally stable in both zwitterionic and anionic membranes. Arcs provide an explanation for the observed ion selectivity in protegrin electrophysiology experiments and could possibly be involved in other membrane permeabilization processes where lipids are thought to participate, such as those induced by antimicrobial peptides and colicins, as well as the Bax apoptotic pore.

A thermodynamic approach to alamethicin pore formation

The structure and energetics of alamethicin Rf30 monomer to nonamer in cylindrical pores of 5 to 11Å radius are investigated using molecular dynamics simulations in an implicit membrane model that includes the free energy cost of acyl chain hydrophobic area exposure. Stable, low energy pores are obtained for certain combinations of radius and oligomeric number. The trimer and the tetramer formed 6Å pores that appear closed while the larger oligomers formed open pores at their optimal radius. The hexamer in an 8Å pore and the octamer in an 11Å pore give the lowest effective energy per monomer. However, all oligomers beyond the pentamer have comparable energies, consistent with the observation of multiple conductance levels. The results are consistent with the widely accepted "barrel-stave" model. The N terminal portion of the molecule exhibits smaller tilt with respect to the membrane normal than the C terminal portion, resulting in a pore shape that is a hybrid between a funnel and an hourglass. Transmembrane voltage has little effect on the structure of the oligomers but enhances or decreases their stability depending on its orientation. Antiparallel bundles are lower in energy than the commonly accepted parallel ones and could be present under certain experimental conditions. Dry aggregates (without an aqueous pore) have lower average effective energy than the corresponding aggregates in a pore, suggesting that alamethicin pores may be excited states that are stabilized in part by voltage and in part by the ion flow itself.

Inclusion of lateral pressure/curvature stress effects in implicit membrane models

Implicit membrane models usually treat the membrane as a hydrophobic slab and neglect lateral pressure/curvature stress effects. As a result, they cannot distinguish, for example, PE from PC lipids. Here, the implicit membrane model IMM1 is extended to include these effects using a combination of classical thermodynamics and membrane elasticity theory. The proposed model is tested by molecular dynamics simulation of the peptides alamethicin, melittin, cyclotide kalata B1, 18A, and KKpL15. The lateral pressure term stabilizes interfacial binding due to the negative pressure at the hydrocarbon-water interface. In agreement with experiment, increase in the peptide/lipid molar ratio shifts the equilibrium from the interfacial to the transmembrane orientation. Simulations of mixed DOPC/DOPE bilayers show that increase of the DOPE mole fraction in general stabilizes interfacial orientations and destabilizes transmembrane orientations. The extent of the stabilization or destabilization varies depending on the exact position of the peptides. The computational results are in good agreement with experiments.

Membrane interactions and pore formation by the antimicrobial peptide protegrin

Protegrin is an antimicrobial peptide with a β-hairpin structure stabilized by a pair of disulfide bonds. It has been extensively studied by solid-state NMR and computational methods. Here we use implicit membrane models to examine the binding of monomers on the surface and in the interior of the membrane, the energetics of dimerization, the binding to membrane pores, and the stability of different membrane barrel structures in pores. Our results challenge a number of conclusions based on previous experimental and theoretical work. The burial of monomers into the membrane interior is found to be unfavorable for any membrane thickness. Because of its imperfect amphipathicity, protegrin binds weakly, at most, on the surface of zwitterionic membranes. However, it binds more favorably onto toroidal pores. Anionic charge on the membrane facilitates the binding due to electrostatic interactions. Solid-state NMR results have suggested a parallel NCCN association of monomers in dimers and association of dimers to form octameric or decameric β-barrels. We find that this structure is not energetically plausible for binding to bilayers, because in this configuration the hydrophobic sides of two monomers point in opposite directions. In contrast, the antiparallel NCCN and especially the parallel NCNC octamers are stable and exhibit a favorable binding energy to the pore. The results of 100-ns simulations in explicit bilayers corroborate the higher stability of the parallel NCNC barrel compared with the parallel NCCN barrel. The ability to form pores in zwitterionic membranes provides a rationalization for the peptide's cytotoxicity. The discrepancies between our results and experiment are discussed, and new experiments are proposed to resolve them and to test the validity of the models.

Modeling peptide binding to anionic membrane pores

Peptide-induced pore formation in membranes can be dissected into two steps: pore formation and peptide binding to the pore. A computational method is proposed to study the second step in anionic membranes. The electrostatic potential is obtained from numerical solutions to the Poisson-Boltzmann equation and is then used in conjunction with IMM1 (implicit membrane model 1). A double charge layer model is used to incorporate the effects of the membrane dipole potential. Inhomogeneity of the charge density in the pore, characterized by explicit membrane simulations of toroidal pores, is included in the model. This approach was applied to two extensively studied peptides, magainin and melittin. In agreement with previous work, binding to toroidal pores is more favorable than binding to the flat membrane. The dependence of binding energy on anionic content exhibits different patterns for the two peptides, in correlation with the different lipid selectivity that has been observed experimentally.

Membrane protein native state discrimination by implicit membrane models

Four implicit membrane models [IMM1, generalized Born (GB)-surface area-implicit membrane (GBSAIM), GB with a simple switching (GBSW), and heterogeneous dielectric GB (HDGB)] were tested for their ability to discriminate the native conformation of five membrane proteins from 450 decoys generated by the Rosetta-Membrane program. The energy ranking of the native state and Z-scores were used to assess the performance of the models. The effect of membrane thickness was examined and was found to be substantial. Quite satisfactory discrimination was achieved with the all-atom IMM1 and GBSW models at 25.4 Å thickness and with the HDGB model at 28.5 Å thickness. The energy components by themselves were not discriminative. Both van der Waals and electrostatic interactions contributed to native state discrimination, to a different extent in each model. Computational efficiency of the models decreased in the order: extended-atom IMM1 > all-atom IMM1 > GBSAIM > GBSW > HDGB. These results encourage the further development and use of implicit membrane models for membrane protein structure prediction.