A Combined Pulse EPR and Monte Carlo Simulation Study Provides Molecular Insight on Peptide−Membrane Interactions

Intracellular Protein-Lipid Interactions Studied by Rapid-Scan Electron Paramagnetic Resonance Spectroscopy

Protein-membrane interactions play key roles in essential cellular processes; studying these interactions in the cell is a challenging task of modern biophysical chemistry. A prominent example is the interaction of human α-synuclein (αS) with negatively charged membranes. It has been well-studied in vitro, but in spite of the huge amount of lipid membranes in the crowded environment of biological cells, to date, no interactions have been detected in cells. Here, we use rapid-scan (RS) electron paramagnetic resonance (EPR) spectroscopy to study αS interactions with negatively charged vesicles in vitro and upon transfection of the protein and lipid vesicles into model cells, i.e., oocytes of Xenopus laevis. We show that protein-vesicle interactions are reflected in RS spectra in vitro and in cells, which enables time-resolved monitoring of protein-membrane interaction upon transfection into cells. Our data suggest binding of a small fraction of αS to endogenous membranes.

Perspective of Use of Antiviral Peptides against Influenza Virus

The threat of a worldwide influenza pandemic has greatly increased over the past decade with the emergence of highly virulent avian influenza strains. The increased frequency of drug-resistant influenza strains against currently available antiviral drugs requires urgent development of new strategies for antiviral therapy, too. The research in the field of therapeutic peptides began to develop extensively in the second half of the 20(th) century. Since then, the mechanisms of action for several peptides and their antiviral prospect received large attention due to the global threat posed by viruses. Here, we discussed the therapeutic properties of peptides used in influenza treatment. Peptides with antiviral activity against influenza can be divided into three main groups. First, entry blocker peptides such as a Flupep that interact with influenza hemagglutinin, block its binding to host cells and prevent viral fusion. Second, several peptides display virucidal activity, disrupting viral envelopes, e.g., Melittin. Finally, a third set of peptides interacts with the viral polymerase complex and act as viral replication inhibitors such as PB1 derived peptides. Here, we present a review of the current literature describing the antiviral activity, mechanism and future therapeutic potential of these influenza antiviral peptides.

Multistep Molecular Dynamics Simulations Identify the Highly Cooperative Activity of Melittin in Recognizing and Stabilizing Membrane Pores

The prototypical antimicrobial peptide, melittin, is well-known for its ability to induce pores in zwitterionic model lipid membranes. However, the mechanism by which melittin accomplishes this is not fully understood. We have conducted all-atom and coarse-grained molecular dynamics simulations which suggest that melittin employs a highly cooperative mechanism for the induction of both small and large membrane pores. The process by which this peptide induces membrane pores appears to be driven by its affinity to membrane defects via its N-terminus region. In our simulations, a membrane defect was deliberately created through either lipid flip-flop or the reorientation of one adsorbed melittin peptide. In a cooperative response, other melittin molecules also inserted their N-termini into the created defect, thus lowering the overall free energy. The insertion of these peptide molecules ultimately allowed the defect to develop into a small transmembrane pore, with an estimated diameter of ∼1.5 nm and a lifetime of the order of tens of milliseconds. In the presence of a finite membrane tension, we show that this small pore can act as a nucleation site for the stochastic rupture of the lipid bilayer, so as to create a much larger pore. We found that a threshold membrane tension of 25 mN/m was needed to create a ruptured pore. Furthermore, by actively accumulating at its edge, adsorbed peptides are able to cooperatively stabilize this larger pore. The defect-mediated pore formation mechanism revealed in this work may also apply to other amphipathic membrane-active peptides.

Probing water density and dynamics in the chaperonin GroEL cavity

ATP-dependent binding of the chaperonin GroEL to its cofactor GroES forms a cavity in which encapsulated substrate proteins can fold in isolation from bulk solution. It has been suggested that folding in the cavity may differ from that in bulk solution owing to steric confinement, interactions with the cavity walls, and differences between the properties of cavity-confined and bulk water. However, experimental data regarding the cavity-confined water are lacking. Here, we report measurements of water density and diffusion dynamics in the vicinity of a spin label attached to a cysteine in the Tyr71 → Cys GroES mutant obtained using two magnetic resonance techniques: electron-spin echo envelope modulation and Overhauser dynamic nuclear polarization. Residue 71 in GroES is fully exposed to bulk water in free GroES and to confined water within the cavity of the GroEL–GroES complex. Our data show that water density and translational dynamics in the vicinity of the label do not change upon complex formation, thus indicating that bulk water-exposed and cavity-confined GroES surface water share similar properties. Interestingly, the diffusion dynamics of water near the GroES surface are found to be unusually fast relative to other protein surfaces studied. The implications of these findings for chaperonin-assisted folding mechanisms are discussed.

Dual-fluorescence L-amino acid reports insertion and orientation of melittin peptide in cell membranes

Monitoring insertion and orientation of peptides in situ on cell membranes remains a challenge. To this end, we synthesized an l-amino acid (AFaa) containing a dual-fluorescence dye of the 3-hydroxyflavone family, as a side chain. In contrast to other labeling approaches using a flexible linker, the AFaa fluorophore, introduced by solid phase synthesis into desired position of a peptide, is attached closely to its backbone with well-defined orientation, and, therefore, could reflect its localization in the membrane. This concept was validated by replacing the leucine-9 (L9) and tryptophan-19 (W19) residues by AFaa in melittin, a well-studied membrane-active peptide. Due to high sensitivity of AFaa dual emission to the environment polarity, we detected a much deeper insertion of L9 peptide position into the bilayer, compared to the W19 position. Moreover, using fluorescence microscopy with a polarized light excitation, we found different orientation of AFaa at L9 and W19 positions of melittin in the bilayers of giant vesicles and cellular membranes. These results suggested that in the natural membranes, similarly to the model lipid bilayers, melittin is preferentially oriented parallel to the membrane surface. The developed amino acid and the proposed methodology will be of interest to study other membrane peptides.

Melittin creates transient pores in a lipid bilayer: results from computer simulations

To study the interaction between melittin peptides and lipid bilayer, we performed coarse-grained simulations on systems containing melittin interacting with a bilayer containing zwitterionic dipalmitoylphosphatidylcholine (DPPC) and anionic palmitoyloleoylphosphatidylglycerol (POPG) phospholipids in a 7:3 ratio. Eight different systems were considered: four at low and four at high peptide to lipid (P/L) ratios. In case of low P/L ratio we did not observe any pore creation in the bilayer. In two out of four of the simulations with the high P/L ratio, appearance of transient pores in the bilayer was observed. These pores were created due to an assembly of 3-5 melittin peptides. Not all of the peptides in the pores were in a transmembrane conformation; many of them had their termini residues anchored to the same leaflet, and these peptides assumed bent, U-shaped, conformations. We propose that when an assembly of melittin peptides creates pores, such an assembly acts as a "wedge" that splits the bilayer. To get a more detailed description of melittin on the bilayer surface and in transient pores, we performed coarse-grained to united-atom scale transformations and after that performed 50 ns molecular dynamics simulations using the united atom description of the systems. While these simulations did not show much of the change in the pore structure during the 50 ns time interval, they clearly showed the presence of water in the transient pores. The appearance of transient pores together with the translocation of peptides across the membranes is consistent with the mechanism proposed to explain graded dye leakage from large vesicles in the presence of melittin.

Topology of the trans-membrane peptide WALP23 in model membranes under negative mismatch conditions

The organization and orientation of membrane-inserted helices is important for better understanding the mode of action of membrane-active peptides and of protein-membrane interactions. Here we report on the application of ESEEM (electron spin-echo envelope modulation) and DEER (double electron-electron resonance) techniques to probe the orientation and oligomeric state of an α-helical trans-membrane model peptide, WALP23, under conditions of negative mismatch between the hydrophobic cores of the model membrane and the peptide. Using ESEEM, we measured weak dipolar interactions between spin-labeled WALP23 and (2)H nuclei of either the solvent (D2O) or of lipids specifically deuterated at the choline group. The ESEEM data obtained from the deuterated lipids were fitted using a model that provided the spin label average distance from a layer of (2)H nuclei in the hydrophilic region of the membrane and the density of the (2)H nuclei in the layer. DEER was used to probe oligomerization through the dipolar interaction between two spin-labels on different peptides. We observed that the center of WALP23 does not coincide with the bilayer midplane and its N-terminus is more buried than the C-terminus. In addition, the ESEEM data fitting yielded a (2)H layer density that was much lower than expected. The DEER experiments revealed the presence of oligomers, the presence of which was attributable to the negative mismatch and the electrostatic dipole of the peptide. A discussion of a possible arrangement of the individual helices in the oligomers that is consistent with the ESEEM and DEER data is presented.

Monte Carlo simulations of peptide-membrane interactions with the MCPep web server

The MCPep server (http://bental.tau.ac.il/MCPep/) is designed for non-experts wishing to perform Monte Carlo (MC) simulations of helical peptides in association with lipid membranes. MCPep is a web implementation of a previously developed MC simulation model. The model has been tested on a variety of peptides and protein fragments. The simulations successfully reproduced available empirical data and provided new molecular insights, such as the preferred locations of peptides in the membrane and the contribution of individual amino acids to membrane association. MCPep simulates the peptide in the aqueous phase and membrane environments, both described implicitly. In the former, the peptide is subjected solely to internal conformational changes, and in the latter, each MC cycle includes additional external rigid body rotational and translational motions to allow the peptide to change its location in the membrane. The server can explore the interaction of helical peptides of any amino-acid composition with membranes of various lipid compositions. Given the peptide’s sequence or structure and the natural width and surface charge of the membrane, MCPep reports the main determinants of peptide–membrane interactions, e.g. average location and orientation in the membrane, free energy of membrane association and the peptide’s helical content. Snapshots of example simulations are also provided.

The Transmembrane Helix Tilt May Be Determined by the Balance between Precession Entropy and Lipid Perturbation

Hydrophobic helical peptides interact with lipid bilayers in various modes, determined by the match between the length of the helix’s hydrophobic core and the thickness of the hydrocarbon region of the bilayer. For example, long helices may tilt with respect to the membrane normal to bury their hydrophobic cores in the membrane, and the lipid bilayer may stretch to match the helix length. Recent molecular dynamics simulations and potential of mean force calculations have shown that some TM helices whose lengths are equal to, or even shorter than, the bilayer thickness may also tilt. The tilt is driven by a gain in the helix precession entropy, which compensates for the free energy penalty resulting from membrane deformation. Using this free energy balance, we derived theoretically an equation of state, describing the dependence of the tilt on the helix length and membrane thickness. To this end, we conducted coarse-grained Monte Carlo simulations of the interaction of helices of various lengths with lipid bilayers of various thicknesses, reproducing and expanding the previous molecular dynamics simulations. Insight from the simulations facilitated the derivation of the theoretical model. The tilt angles calculated using the theoretical model agree well with our simulations and with previous calculations and measurements.

DEER distance measurements on proteins

Distance distributions between paramagnetic centers in the range of 1.8 to 6 nm in membrane proteins and up to 10 nm in deuterated soluble proteins can be measured by the DEER technique. The number of paramagnetic centers and their relative orientation can be characterized. DEER does not require crystallization and is not limited with respect to the size of the protein or protein complex. Diamagnetic proteins are accessible by site-directed spin labeling. To characterize structure or structural changes, experimental protocols were optimized and techniques for artifact suppression were introduced. Data analysis programs were developed, and it was realized that interpretation of the distance distributions must take into account the conformational distribution of spin labels. First methods have appeared for deriving structural models from a small number of distance constraints. The present scope and limitations of the technique are illustrated.

The mode of α-synuclein binding to membranes depends on lipid composition and lipid to protein ratio

Interactions of the presynaptic protein α-synuclein with membranes are involved in its physiological action as well as in the pathological misfolding and aggregation related to Parkinsons's disease. We studied the conformation and orientation of α-synuclein bound to model vesicular membranes using multiparametric response polarity-sensitive fluorescent probes together with CD and EPR measurements. At low lipid to α-synuclein ratio the protein binds membranes through its N-terminal domain. When lipids are in excess, the α-helical content and the role of the C-terminus in binding increase. Highly rigid membranes also induce a greater α-helical content and a lower polarity of the protein microenvironment.

Membrane interactions of novicidin, a novel antimicrobial peptide: phosphatidylglycerol promotes bilayer insertion

Novicidin is an antimicrobial peptide derived from ovispirin, a cationic peptide which originated from the ovine cathelicidin SMAP-29. Novicidin, however, has been designed to minimize the cytotoxic properties of SMAP-29 and ovisipirin toward achieving potential therapeutic applications. We present an analysis of membrane interactions and lipid bilayer penetration of novicidin, using an array of biophysical techniques and biomimetic membrane assemblies, complemented by Monte Carlo (MC) simulations. The data indicate that novicidin interacts minimally with zwitterionic bilayers, accounting for its low hemolytic activity. Negatively charged phosphatidylglycerol, on the other hand, plays a significant role in initiating membrane binding of novicidin, and promotes peptide insertion into the interface between the lipid headgroups and the acyl chains. The significant insertion into bilayers containing negative phospholipids might explain the enhanced antibacterial properties of novicidin. Overall, this study highlights two distinct outcomes for membrane interactions of novicidin, and points to a combination between electrostatic attraction to the lipid/water interface and penetration into the subsurface lipid headgroups region as important determinants for the biological activity of novicidin.

Interaction of an antimicrobial peptide with membranes: experiments and simulations with NKCS

We used Monte Carlo simulations and biophysical measurements to study the interaction of NKCS, a derivative of the antimicrobial peptide NK-2, with a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) membrane. The simulations showed that NKCS adsorbed on the membrane surface and the dominant conformation featured two amphipathic helices connected by a hinge region. We designed two mutants in the hinge to investigate the interplay between helicity and membrane affinity. Simulations with a Leu-to-Pro substitution showed that the helicity and membrane affinity of the mutant (NKCS-[LP]) decreased. Two Ala residues were added to NKCS to produce a sequence that is compatible with a continuous amphipathic helix structure (NKCS-[AA]), and the simulations showed that the mutant adsorbed on the membrane surface with a particularly high affinity. The circular dichroism spectra of the three peptides also showed that NKCS-[LP] is the least helical and NKCS-[AA] is the most. However, the activity of the peptides, determined in terms of their antimicrobial potency and influence on the temperature of the transition of the lipid to hexagonal phase, displayed a complex behavior: NKCS-[LP] was the least potent and had the smallest influence on the transition temperature, and NKCS was the most potent and had the largest effect on the temperature.

Electron spin-echo envelope modulation (ESEEM) reveals water and phosphate interactions with the KcsA potassium channel

Electron spin-echo envelope modulation (ESEEM) spectroscopy is a well-established technique for the study of naturally occurring paramagnetic metal centers. The technique has been used to study copper complexes, hemes, enzyme mechanisms, micellar water content, and water permeation profiles in membranes, among other applications. In the present study, we combine ESEEM spectroscopy with site-directed spin labeling (SDSL) and X-ray crystallography in order to evaluate the technique's potential as a structural tool to describe the native environment of membrane proteins. Using the KcsA potassium channel as a model system, we demonstrate that deuterium ESEEM can detect water permeation along the lipid-exposed surface of the KcsA outer helix. We further demonstrate that (31)P ESEEM is able to identify channel residues that interact with the phosphate headgroup of the lipid bilayer. In combination with X-ray crystallography, the (31)P data may be used to define the phosphate interaction surface of the protein. The results presented here establish ESEEM as a highly informative technique for SDSL studies of membrane proteins.