Interactions of alpha-helices with lipid bilayers: a review of simulation studies

Recent Advances of Studies on Cell-Penetrating Peptides Based on Molecular Dynamics Simulations

With the ability to transport cargo molecules across cell membranes with low toxicity, cell-penetrating peptides (CPPs) have become promising candidates for next generation peptide-based drug delivery vectors. Over the past three decades since the first CPP was discovered, a great deal of work has been done on the cellular uptake mechanisms and the applications for the delivery of therapeutic molecules, and significant advances have been made. But so far, we still do not have a precise and unified understanding of the structure-activity relationship of the CPPs. Molecular dynamics (MD) simulations provide a method to reveal peptide-membrane interactions at the atomistic level and have become an effective complement to experiments. In this paper, we review the progress of the MD simulations on CPP-membrane interactions, including the computational methods and technical improvements in the MD simulations, the research achievements in the CPP internalization mechanism, CPP decoration and coupling, and the peptide-induced membrane reactions during the penetration process, as well as the comparison of simulated and experimental results.

Antimicrobial Peptides with Anti-Candida Activity

Mycoses are accountable for millions of infections yearly worldwide. Invasive candidiasis is the most usual, presenting a high morbidity and mortality. Candida albicans remains the prevalent etiologic agent, but the incidence of other species such as Candida parapsilosis, Candida glabrata and Candida auris keeps increasing. These pathogens frequently show a reduced susceptibility to commonly used antifungal drugs, including polyenes, triazoles and echinocandins, and the incidence of emerging multi-drug-resistant strains of these species continues to increase. Therefore, the need to search for new molecules that target these pathogenic species in a different manner is now more urgent than ever. Nature is an almost endless source of interesting new molecules that could meet this need. Among these molecules, antimicrobial peptides, present in different sources in nature, possess some advantages over conventional antifungal agents, even with their own drawbacks, and are considered as a promising pharmacological option against a wide range of microbial infections. In this review, we describe 20 antimicrobial peptides from different origins that possess an activity against Candida.

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.

The influence of lipid membranes on fluorescent probes' optical properties

BACKGROUND

Organic fluorophores embedded in lipid bilayers can nowadays be described by a multiscale computational approach. Combining different length and time scales, a full characterization of the probe localization and optical properties led to novel insight into the effect of the environments.

SCOPE OF REVIEW

Following an introduction on computational advancements, three relevant probes are reviewed that delineate how a multiscale approach can lead to novel insight into the probes' (non) linear optical properties. Attention is paid to the quality of the theoretical description of the optical techniques.

MAJOR CONCLUSIONS

Computation can assess a priori novel probes' optical properties and guide the analysis and interpretation of experimental data in novel studies. The properties can be used to gain information on the phase and condition of the surrounding biological environment.

GENERAL SIGNIFICANCE

Computation showed that a canonical view on some of the probes should be revisited and adapted.

Saddle-curvature instability of lipid bilayer induced by amphipathic peptides: a molecular model

Amphipathic peptides that partition into lipid bilayers affect the curvature elastic properties of their host. Some of these peptides are able to shift the Gaussian modulus to positive values, thus triggering an instability with respect to the formation of saddle curvatures. To characterize the generic aspects of the underlying mechanism, we employ a molecular lipid model that accounts for the interfacial tension between the polar and apolar regions of the membrane, for interactions between the lipid headgroups, and for the energy to stretch or compress the hydrocarbon chains. Peptides are modeled as cylinders that partition into the host membrane in a parallel orientation where they diminish the space available to the lipid headgroups and chains. The penetration depth into the membrane is determined by the angular size of the peptide's hydrophilic region. We demonstrate that only peptides with a small angular size of their hydrophilic region have an intrinsic tendency to render the Gaussian modulus more positive, and we identify conditions at which the Gaussian modulus adopts a positive sign upon increasing the peptide concentration. Our model allows us to also incorporate electrostatic interactions between cationic peptides and anionic lipids on the level of the linear Debye-Hückel model. We show that electrostatic interactions tend to shift the Gaussian modulus toward more positive values. Steric and electrostatic lipid-peptide interactions jointly decrease the effective interaction strength in the headgroup region of the host membrane thus suggesting a generic mechanisms of how certain amphipathic peptides are able to induce the formation of saddle curvatures.

Molecular dynamics study of membrane permeabilization by wild-type and mutant lytic peptides from the non-enveloped Flock House virus

Flock House virus (FHV) serves as a model system for understanding infection mechanisms utilized by non-enveloped viruses to transport across cellular membranes. During the infection cycle of FHV, a fundamental stage involves disruption of the endosomal membrane by membrane active peptides, following externalization of the peptides from the capsid interior. The FHV lytic agents are the 44 C-terminal amino acids residues of the capsid protein, which are auto-catalytically cleaved during the capsid maturation process. The cleaved peptides are termed γ peptides. In this study, we perform multi-scale molecular dynamics simulations including 40 μs all-atom molecular dynamics simulations to study the behavior of pre-inserted transmembrane lytic peptides at a high concentration in a neutral membrane. We study the dynamical organization among peptides to form oligomeric bundles in four systems including the wild-type γ peptide and three mutant forms; namely, a truncation mutant in which the 23 C-terminal residues are deleted (γ₁), a construct where the 8 C-terminal residues of γ are fused to γ₁ (Δ385-399 γ) and a single-point mutant (F402A γ), all of which have been experimentally shown to drastically affect infectivity and lytic activity compared to the wild-type γ. Our results shed light on the actions of varied forms of the FHV lytic peptide including membrane insertion, trans-membrane stability, peptide oligomerization, water permeation activity and dynamic pore formation. Findings from this study provide detailed structural information and rationale for the differences in lytic activity among variants of FHV γ.

Synergism of cationic antimicrobial peptide WLBU2 with antibacterial agents against biofilms of multi-drug resistant Acinetobacter baumannii and Klebsiella pneumoniae

Purpose

The activity of the cationic antimicrobial peptide WLBU2 was evaluated against planktonic cells and biofilms of multi-drug resistant (MDR) Acinetobacter baumannii and Klebsiella pneumoniae, alone and in combination with classical antimicrobial agents.

Methods

Control American Type Culture Collection (ATCC) strains and MDR clinical isolates of A. baumannii and K. pneumoniae were utilized. The minimum inhibitory concentration (MIC) and minimum bactericidal concentration (MBC) of WLBU2 alone and in combination with antimicrobials were determined by classical methods. The Calgary biofilm device was used to determine the minimum biofilm eradication concentration (MBEC). The MTT assay was used to determine the cytotoxicity of agents on eukaryotic cells. The electrophoretic mobility shift assay was used to evaluate the ability of WLBU2 to bind bacterial DNA.

Results

The WLBU2 MIC and MBC values were identical indicating bactericidal activity. The MIC/MBC values ranged from 1.5625 to 12.5 µM. At these concentrations, Vero cells and human skin fibroblasts were viable. The MBEC of WLBU2 ranged from 25 to 200 µM. A significant loss of eukaryotic cell viability was observed at the MBEC range. The combination of sub-inhibitory concentrations of WLBU2 with amoxicillin-clavulanate or ciprofloxacin for K. pneumoniae, and with tobramycin or imipenem for A. baumannii, demonstrated synergism, leading to a significant decrease in MIC and MBEC values for some isolates and ATCC strains. However, all combinations were associated with considerable loss in eukaryotic cells’ viability. WLBU2 did not demonstrate the ability to bind bacterial plasmid DNA.

Conclusion

WLBU2 in combination with antimicrobials holds promise in eradication of MDR pathogens.

Binding, folding and insertion of a β-hairpin peptide at a lipid bilayer surface: Influence of electrostatics and lipid tail packing

Antimicrobial peptides (AMPs) act as host defenses against microbial pathogens. Here we investigate the interactions of SVS-1 (KVKVKVKV^dP^lPTKVKVKVK), an engineered AMP and anti-cancer β-hairpin peptide, with lipid bilayers using spectroscopic studies and atomistic molecular dynamics simulations. In agreement with literature reports, simulation and experiment show preferential binding of SVS-1 peptides to anionic over neutral bilayers. Fluorescence and circular dichroism studies of a Trp-substituted SVS-1 analogue indicate, however, that it will bind to a zwitterionic DPPC bilayer under high-curvature conditions and folds into a hairpin. In bilayers formed from a 1:1 mixture of DPPC and anionic DPPG lipids, curvature and lipid fluidity are also observed to promote deeper insertion of the fluorescent peptide. Simulations using the CHARMM C36m force field offer complementary insight into timescales and mechanisms of folding and insertion. SVS-1 simulated at an anionic mixed POPC/POPG bilayer folded into a hairpin over a microsecond, the final stage in folding coinciding with the establishment of contact between the peptide's valine sidechains and the lipid tails through a "flip and dip" mechanism. Partial, transient folding and superficial bilayer contact are seen in simulation of the peptide at a zwitterionic POPC bilayer. Only when external surface tension is applied does the peptide establish lasting contact with the POPC bilayer. Our findings reveal the influence of disruption to lipid headgroup packing (via curvature or surface tension) on the pathway of binding and insertion, highlighting the collaborative effort of electrostatic and hydrophobic interactions on interaction of SVS-1 with lipid bilayers.

Folding Simulations of the Transmembrane Helix of Virus Protein U in an Implicit Membrane Model

Vpu is an 81-amino-acid auxiliary membrane protein encoded by human immunodeficiency virus type 1 (HIV-1). One of its roles is to amplify viral release by self-assembling in homo-oligomers to form functional water-filled pores enabling the flux of ions across the membrane. Various NMR and CD studies have shown that the transmembrane domain of Vpu has a helical conformation. With a recently developed implicit membrane model and an efficient Monte Carlo (MC) algorithm using concerted backbone rotations, we simulate the folding of the transmembrane domain of Vpu at atomic resolution. The implicit membrane environment is based on the generalized Born theory and enables very long time scale events, such as folding to be observed using detailed all-atom representation of the protein. Such studies are currently computationally unfeasible with fully explicit lipid bilayer molecular dynamics simulations. The correct helical transmembrane structure of Vpu is predicted from extended conformations and remains stably inserted. Tilt and kink angles agree well with experimental estimates from NMR measurements. The experimentally observed change in tilt angle in membranes of varying hydrophobic width is accurately reproduced. The extensive simulation of a pentamer of the Vpu transmembrane domain in the implicit membrane gives results similar to the ones reported previously for fully explicit bilayer simulations.

Modeling the structure of SARS 3a transmembrane protein using a minimum unfavorable contact approach

3a is an accessory protein from SARS coronavirus that is known to play a significant role in the proliferation of the virus by forming tetrameric ion channels. Although the monomeric units are known to consist of three transmembrane (TM) domains, there are no solved structures available for the complete monomer. The present study proposes a structural model for the transmembrane region of the monomer by employing our previously tested approach, which predicts potential orientations of TM α-helices by minimizing the unfavorable contact surfaces between the different TM domains. The best model structure comprising all three α-helices has been subjected to MD simulations to examine its quality. The TM bundle was found to form a compact and stable structure with significant intermolecular interactions. The structural features of the proposed model of 3a account for observations from previous experimental investigations on the activity of the protein. Further analysis indicates that residues from the TM2 and TM3 domains are likely to line the pore of the ion channel, which is in good agreement with a recent experimental study. In the absence of an experimental structure for the protein, the proposed structure can serve as a useful model for inferring structure-function relationships about the protein.

Graphical Abstract

Effects of PEGylation on the binding interaction of magainin 2 and tachyplesin I with lipid bilayer surface

Poly(ethylene glycol) (PEG)-grafted magainin 2 and tachyplesin I were simulated with lipid bilayers. In the simulations of PEGylated magainin 2 and tachyplesin I in water, both peptides are wrapped by PEG chains. The α-helical structure of PEGylated magainin 2 is broken in water, while the β-sheet of PEGylated tachyplesin I keeps stable, similar to the structural behavior of unPEGylated peptides, in agreement with experiments. Simulations of PEGylated peptides with lipid bilayers show that PEG chains block the electrostatic interaction between cationic residues of peptides and anionic phosphates of lipids, leading to the less binding of the peptide to the bilayer surface, which is observed more significantly for magainin 2 than for tachyplesin I. Since the random-coiled magainin 2 can be more completely covered by PEGs than does the β-sheet tachyplesin I, the PEGylation effect on the decreased binding is larger for magainin 2, showing the dependence of PEGylation on the peptide structure. These simulation findings qualitatively support the experimental observation of the different extents of the reduced membrane-permeabilizing activity for PEGylated magainin 2 and tachyplesin I.

Peptide-lipid interactions: experiments and applications

The interactions between peptides and lipids are of fundamental importance in the functioning of numerous membrane-mediated cellular processes including antimicrobial peptide action, hormone-receptor interactions, drug bioavailability across the blood-brain barrier and viral fusion processes. Moreover, a major goal of modern biotechnology is obtaining new potent pharmaceutical agents whose biological action is dependent on the binding of peptides to lipid-bilayers. Several issues need to be addressed such as secondary structure, orientation, oligomerization and localization inside the membrane. At the same time, the structural effects which the peptides cause on the lipid bilayer are important for the interactions and need to be elucidated. The structural characterization of membrane active peptides in membranes is a harsh experimental challenge. It is in fact accepted that no single experimental technique can give a complete structural picture of the interaction, but rather a combination of different techniques is necessary.

Plastid genetic engineering in Solanaceae

Plastid genetic engineering has come of age, becoming today an attractive alternative approach for the expression of foreign genes, as it offers several advantages over nuclear transformants. Significant progress has been made in plastid genetic engineering in tobacco and other Solanaceae plants, through the use of improved regeneration procedures and transformation vectors with efficient promoters and untranslated regions. Many genes encoding for industrially important proteins and vaccines, as well as genes conferring important agronomic traits, have been stably integrated and expressed in the plastid genome. Despite these advances, it remains a challenge to achieve marked levels of plastid transgene expression in non-green tissues. In this review, we summarize the basic requirements of plastid genetic engineering and discuss the current status, limitations, and the potential of plastid transformation for expanding future studies relating to Solanaceae plants.

Biophysical properties of membrane-active peptides based on micelle modeling: a case study of cell-penetrating and antimicrobial peptides

We investigated the molecular mechanisms of short peptides interacting with membrane-mimetic systems. Three short peptides were selected for this study: penetratin as a cell-penetrating peptide (CPP), and temporin A and KSL as antimicrobial peptides (AMP). We investigated the detailed interactions of the peptides with dodecylphosphocholine (DPC) and sodium dodecyl sulfate (SDS) micelles, and the subsequent peptide insertion based on free energy calculations by using all-atomistic molecular dynamics simulations with the united atom force field and explicit solvent models. First, we found that the free energy barrier to insertion for the three peptides is dependent on the chemical composition of the micelles. Because of the favorable electrostatic interactions between the peptides and the headgroups of lipids, the insertion barrier into an SDS micelle is less than a DPC micelle. Second, the peptides' secondary structures may play a key role in their binding and insertion ability, particularly for amphiphilic peptides such as penetratin and KSL. The secondary structures with a stronger ability to bind with and insert into micelles are the ones that account for a smaller surface area of hydrophobic core, thus offering a possible criterion for peptide design with specific functionalities.

The thermodynamics of simple biomembrane mimetic systems

Insight into the forces governing a system is essential for understanding its behavior and function. Thermodynamic investigations provide a wealth of information that is not, or is hardly, available from other methods. This article reviews thermodynamic approaches and assays to measure collective properties such as heat adsorption / emission and volume variations. These methods can be successfully applied to the study of lipid vesicles (liposomes) and biological membranes. With respect to instrumentation, differential scanning calorimetry, pressure perturbation calorimetry, isothermal titration calorimetry, dilatometry, and acoustic techniques aimed at measuring the isothermal and adiabatic processes, two- and three-dimensional compressibilities are considered. Applications of these techniques to lipid systems include the measurement of different thermodynamic parameters and a detailed characterization of thermotropic, barotropic, and lyotropic phase behavior. The membrane binding and / or partitioning of solutes (proteins, peptides, drugs, surfactants, ions, etc.) can also be quantified and modeled. Many thermodynamic assays are available for studying the effect of proteins and other additives on membranes, characterizing non-ideal mixing, domain formation, bilayer stability, curvature strain, permeability, solubilization, and fusion. Studies of membrane proteins in lipid environments elucidate lipid-protein interactions in membranes. Finally, a plethora of relaxation phenomena toward equilibrium thermodynamic structures can be also investigated. The systems are described in terms of enthalpic and entropic forces, equilibrium constants, heat capacities, partial volume changes, volume and area compressibility, and so on, also shedding light on the stability of the structures and the molecular origin and mechanism of the structural changes.

NMR structure of a viral peptide inserted in artificial membranes: a view on the early steps of the birnavirus entry process

Nonenveloped virus must penetrate the cellular membrane to access the cytoplasm without the benefit of membrane fusion. For birnavirus, one of the peptides present in the virus capsid, pep46 for infectious bursal disease virus, is able to induce pores into membranes as an intermediate step of the birnavirus-penetration pathway. Using osmotic protection experiments, we demonstrate here that pep46 and its pore-forming N-terminal moiety (pep22) form pores of different diameters, 5-8 and 2-4 nm, respectively, showing that both pep46 moieties participate to pore formation. The solution structures of pep46, pep22, and pep24 (the pep46 C-terminal moiety) in different hydrophobic environments and micelles determined by (1)H NMR studies provide structural insights of the pep46 domain interaction. In CDCl(3)/CD(3)OH mixture and in dodecylphosphocholine micelles, the N-terminal domain of pep46 is structured in a long kinked helix, although the C terminus is structured in one or two helices depending upon the solvents used. We also show that the folding and the proline isomerization status of pep46 depend on the type of hydrophobic environment. NMR spectroscopy with labeled phospholipid micelles, differential scanning calorimetry, and plasmon waveguide resonance studies show the peptides lie parallel to the lipid-water interface, perturbing the fatty acid chain packing. All these data lead to a model in which the two domains of pep46 interact with the membrane to form pores.

Expression of a synthesized gene encoding cationic peptide cecropin B in transgenic tomato plants protects against bacterial diseases

The cationic lytic peptide cecropin B (CB), isolated from the giant silk moth (Hyalophora cecropia), has been shown to effectively eliminate Gram-negative and some Gram-positive bacteria. In this study, the effects of chemically synthesized CB on plant pathogens were investigated. The S(50)s (the peptide concentrations causing 50% survival of a pathogenic bacterium) of CB against two major pathogens of the tomato, Ralstonia solanacearum and Xanthomonas campestris pv. vesicatoria, were 529.6 microg/ml and 0.29 microg/ml, respectively. The CB gene was then fused to the secretory signal peptide (sp) sequence from the barley alpha-amylase gene, and the new construct, pBI121-spCB, was used for the transformation of tomato plants. Integration of the CB gene into the tomato genome was confirmed by PCR, and its expression was confirmed by Western blot analyses. In vivo studies of the transgenic tomato plant demonstrated significant resistance to bacterial wilt and bacterial spot. The levels of CB expressed in transgenic tomato plants ( approximately 0.05 microg in 50 mg of leaves) were far lower than the S(50) determined in vitro. CB transgenic tomatoes could therefore be a new mode of bioprotection against these two plant diseases with significant agricultural applications.

Dermaseptins from Phyllomedusa oreades and Phyllomedusa distincta: Secondary structure, antimicrobial activity, and mammalian cell toxicity

The present study reports the structural characteristics, the biological activities, and preliminary clinical investigations of three synthetic members of the dermaseptin family of antimicrobial peptides. The three peptides showed similar tendencies to form alpha-helical structures in non-polar media. The antimicrobial activity towards bacteria and fungi was determined in the micromolar concentration and the peptides did not influenced peritoneal cells viability. One of the peptides was intravenously administered in mice at concentrations similar to those of antibiotics employed in bacterial/fungal infections and it did not cause any detectable changes in cells and tissues.

Characterization of an ideal amphipathic peptide as a procoagulant agent

On the basis of previous evidence that amphipathic helical peptides accelerate Factor IXa activation of Factor X [Blostein, Rigby, Furie, Furie and Gilbert (2000) Biochemistry 39, 12000–12006], the present study was designed to assess the procoagulant activity of an IAP (ideal amphipathic peptide) of Lys7Leu15 composition. The results show that IAP accelerates Factor X activation by Factor IXa in a concentration-dependent manner and accelerates thrombin generation by Factor Xa with a comparable peptide- and substrate-concentration-dependence. A scrambled helical peptide with the same amino acid composition as IAP, but with its amphipathicity abolished, eliminated most of the aforementioned effects. The Gla (γ-carboxyglutamic acid)-rich domain of Factor X is required for IAP activity, suggesting that this peptide behaves as a phospholipid membrane. This hypothesis was confirmed, using fluorescence spectroscopy, by demonstrating direct binding between IAP and the Gla-rich domain of Factor X. In addition, the catalytic efficiencies of the tenase and prothrombinase enzymatic complexes, containing cofactors Factor VIIIa and Factor Va respectively, are enhanced by IAP. Finally, we show that IAP delays clot lysis in vitro. In summary, these observations demonstrate that IAP not only enhances essential procoagulant reactions required for fibrin generation, but also inhibits fibrinolysis, suggesting a potential role for IAP as a haemostatic agent.

The impact of peptides on lipid membranes

We review the fundamental strategies used by small peptides to associate with lipid membranes and how the different strategies impact on the structure and dynamics of the lipids. In particular we focus on the binding of amphiphilic peptides by electrostatic and hydrophobic forces, on the anchoring of peptides to the bilayer by acylation and prenylation, and on the incorporation of small peptides that form well-defined channels. The effect of lipid-peptide interactions on the lipids is characterized in terms of lipid acyl-chain order, membrane thickness, membrane elasticity, permeability, lipid-domain and annulus formation, as well as acyl-chain dynamics. The different situations are illustrated by specific cases for which experimental observations can be interpreted and supplemented by theoretical modeling and simulations. A comparison is made with the effect on lipids of trans-membrane proteins. The various cases are discussed in the context of the possible roles played by lipid-peptide interactions for the biological, physiological, and pharmacological function of peptides.

Real-time structural investigation of a lipid bilayer during its interaction with melittin using sum frequency generation vibrational spectroscopy

Interactions between membrane bilayers and peptides/proteins are ubiquitous throughout a cell. To determine the structure of membrane bilayers and the associated peptides/proteins, model systems such as supported lipid bilayers are often used. It has been difficult to directly investigate the interactions between a single membrane bilayer and peptides/proteins without exogenous labeling. In this work we demonstrate that sum frequency generation vibrational spectroscopy can be employed to study the interactions between peptides/proteins and a single lipid bilayer in real time, in situ, and without exogenous labeling. Using melittin and a dipalmitoyl phosphatidylglycerol bilayer as a model system, we monitored the C-H and C-D stretching signals from isotopically symmetric or asymmetric dipalmitoyl phosphatidylglycerol bilayers during their interaction with melittin. It has been found that the extent and kinetics of bilayer perturbation induced by melittin are very sensitive to melittin concentration. Such concentration dependence is correlated to melittin's mode of action. Melittin is found to function via the early and late stage of the carpet model at low and high concentrations, respectively, whereas the toroidal model is probable at intermediate concentrations. This research illustrates the potential of sum frequency generation as a biophysical technique to monitor individual leaflet structure of lipid bilayers in real time during their interactions with biomolecules.

A generalized born implicit-membrane representation compared to experimental insertion free energies

An implicit-membrane representation based on the generalized Born theory of solvation has been developed. The method was parameterized against the water-to-cyclohexane insertion free energies of hydrophobic side-chain analogs. Subsequently, the membrane was compared with experimental data from translocon inserted polypeptides and validated by comparison with an independent dataset of six membrane-associated peptides and eight integral membrane proteins of known structure and orientation. Comparison of the insertion energy of alpha-helical model peptides with the experimental values from the biological hydrophobicity scale of Hessa et al. gave a correlation of 93% with a mean unsigned error of 0.64 kcal/mol, when charged residues were ignored. The membrane insertion energy was found to be dependent on residue position. This effect is particularly pronounced for charged and polar residues, which strongly prefer interfacial locations. All integral membrane proteins investigated orient and insert correctly into the implicit-membrane model. Remarkably, the membrane model correctly predicts a partially inserted configuration for the monotopic membrane protein cyclooxygenase, matching experimental and theoretical predictions. To test the applicability and usefulness of the implicit-membrane method, molecular simulations of influenza A M2 as well as the glycophorin A dimer were performed. Both systems remain structurally stable and integrated into the membrane.

Solution structure of Ser14Gly-humanin, a potent rescue factor against neuronal cell death in Alzheimer's disease

The NMR solution study of Ser14Gly-humanin (S14G-HN), a 1000-fold more potent derivative of humanin (HN), is reported. HN is 24-residue peptide that selectively suppresses neuronal cell death caused by Alzheimer's disease (AD)-specific insults and offers hope for the development of a cure against AD. In aqueous solution the NMR data show that S14G-HN is a flexible peptide with turn-like structures in its conformational ensemble distributed over an extensive part of its sequence from Pro3 to Glu15. In the more lipophilic environment of 30% TFE, an alpha-helical structure spanning residues Phe6 to Thr13 is identified. Comparison of these findings to the NMR structure of the parent HN and to existing structure-function relationship literature data outlines the important for activity structural features for this class of neuroprotective peptides, and brings forth flexibility as an important characteristic that may facilitate interactions with functional counterparts of the neuroprotection pathway.

Structure of magainin and alamethicin in model membranes studied by x-ray reflectivity

We have investigated the structure of lipid bilayers containing varied molar ratios of different lipids and the antimicrobial peptides magainin and alamethicin. For this structural study, we have used x-ray reflectivity on highly aligned solid-supported multilamellar lipid membranes. The reflectivity curves have been analyzed by semi-kinematical reflectivity theory modeling the bilayer density profile rho(z). Model simulations of the reflectivity curves cover a large range of vertical momentum transfer q(z), and yield excellent agreement between data and theory. The structural changes observed as a function of the molar peptide/lipid concentration P/L are discussed in a comparative way.

Association of transmembrane helices: what determines assembling of a dimer?

Self-association of two hydrophobic alpha-helices is studied via unrestrained Monte Carlo (MC) simulations in a hydrophobic slab described by an effective potential. The system under study represents two transmembrane (TM) segments of human glycophorin A (GpA), which form homo-dimers in membranes. The influence of TM electrostatic potential, thickness and hydrophobicity degree of lipid bilayer is investigated. It is shown that the membrane environment stabilizes alpha-helical conformation of GpA monomers, induces their TM insertion and facilitates inter-helical contacts. Head-to-head orientation of the helices is promoted by the voltage difference across the membrane. Subsequent "fine-tuned" assembling of the dimer is mediated by van der Waals interactions. Only the models of dimer, calculated in a hydrophobic slab with applied voltage agree with experimental data, while simulations in vacuo or without TM voltage fail to give reasonable results. The moderate structural heterogeneity of GpA dimers (existence of several groups of states with close energies) is proposed to reflect their equilibrium dynamics in membrane-mimics. The calculations performed for GpA mutants G83A and G86L permit rationalization of mutagenesis data for them. The results of Monte Carlo simulations critically depend on the parameters of the membrane model: adequate description of helix association is achieved in the water-cyclohexane-water system with the membrane thickness 30-34 A, while in membranes with different hydrophobicities and thickness unrealistic conformations of the dimer are found. The computational approach permits efficient prediction of TM helical oligomers based solely on the sequences of interacting peptides.

Solvent-free simulations of fluid membrane bilayers

A molecular level model for lipid bilayers is presented. Lipids are represented by rigid, asymmetric, soft spherocylinders in implicit solvent. A simple three parameter potential between pairs of lipids gives rise to a rich assortment of phases including (but not limited to) micelles, fluid bilayers, and gel-like bilayers. Monte Carlo simulations have been carried out to verify self-assembly, characterize the phases corresponding to different potential parametrizations, and to quantify the physical properties associated with those parameter sets corresponding to fluid bilayer behavior. The studied fluid bilayers have compressibility moduli in agreement with experimental systems, but display bending moduli at least three times larger than typical biological membranes without cholesterol.

High-Field NMR Studies of Molecular Recognition and Structure−Function Relationships in Antimicrobial Piscidins at the Water−Lipid Bilayer Interface

High magnetic field solid-state NMR was performed on amphipathic cationic antimicrobial peptides from fish to characterize their secondary structure and orientation in hydrated phospholipid bilayers. High-resolution distance and orientational restraints on 13C- and 15N-labeled amidated piscidins 1 and 3 provided site-specific information establishing alpha-helicity and an orientation parallel to the membrane surface. Few membrane-bound natural peptides with this topology have been structurally studied at high resolution in the presence of hydrated lipid bilayers. This orientation was foreseen since the partitioning of amphipathic cationic antimicrobial peptides at the water-bilayer interface allows for favorable peptide-lipid interactions, and it may be related to the mechanism of action. The enhanced resolution obtained at 900 MHz evidences a determinant advantage of ultra-high-field NMR for the structural determination of multiple-labeled peptides and proteins.

Properties of integral membrane protein structures: derivation of an implicit membrane potential

Distributions of each amino acid in the trans-membrane domain were calculated as a function of the membrane normal using all currently available alpha-helical membrane protein structures with resolutions better than 4 A. The results were compared with previous sequence- and structure-based analyses. Calculation of the average hydrophobicity along the membrane normal demonstrated that the protein surface in the membrane domain is in fact much more hydrophobic than the protein core. While hydrophobic residues dominate the membrane domain, the interfacial regions of membrane proteins were found to be abundant in the small residues glycine, alanine, and serine, consistent with previous studies on membrane protein packing. Charged residues displayed nonsymmetric distributions with a preference for the intracellular interface. This effect was more prominent for Arg and Lys resulting in a direct confirmation of the positive inside rule. Potentials of mean force along the membrane normal were derived for each amino acid by fitting Gaussian functions to the residue distributions. The individual potentials agree well with experimental and theoretical considerations. The resulting implicit membrane potential was tested on various membrane proteins as well as single trans-membrane alpha-helices. All membrane proteins were found to be at an energy minimum when correctly inserted into the membrane. For alpha-helices both interfacial (i.e. surface bound) and inserted configurations were found to correspond to energy minima. The results demonstrate that the use of trans-membrane amino acid distributions to derive an implicit membrane representation yields meaningful residue potentials.

Amphipathic alpha-helices in proteins: results from analysis of protein structures

Amphipathic alpha-helices play a crucial role in mediating the interaction of peptides and proteins with membranes. We have analyzed protein structures for the occurrence of 18-residue amphipathic helices. We find several of these alpha-helices having average hydrophobic moments and average hydrophobicities that would favor their interaction with membranes. We have analyzed the distribution of net charge, helix length, normalized frequency of occurrence, and propensities of the 20 amino acids in the delineated 18-residue helices. We have observed distinct differences in the frequencies of occurrence of polar and hydrophobic amino acids at positions 1-18 in amphipathic and nonamphipathic helices. There are also differences in propensities of the 20 amino acids to occur at positions 1-18 of amphipathic and nonamphipathic helices. Synthetic peptides corresponding to some of these surface-seeking helices do possess antibacterial and/or hemolytic activities. Knowledge of the distribution of charges in 18-residue surface-seeking amphipathic alpha-helices, as well as propensity of occurrence of amino acids at various positions, would be useful inputs in the de novo design of amphipathic peptides.

Evaluation of Transgenic ‘Chardonnay’ (Vitis vinifera) Containing Magainin Genes for Resistance to Crown Gall and Powdery Mildew

Magainins, short peptides with broad-spectrum antimicrobial activity in vitro, were assayed for their ability to confer resistance to pathogens in transgenic grapevines. Embryogenic cell suspensions of 'Chardonnay' (Vitis vinifera L.) were co-transformed by microprojectile bombardment with a plasmid carrying the npt-II gene and a second plasmid harboring either a natural magainin-2 (mag2) or a synthetic derivative (MSI99) gene. Magainin genes and the marker gene were driven by Arabidopsis ubiquitin-3 and ubiquitin-11 promoters, respectively. A total of 10 mag2 and 9 MSI99 regenerated lines were studied by Southern blot hybridization, which showed 1-6 transgene integration events into the plant genome. Semi-quantitative reverse transcription-polymerase chain reaction (RT-PCR) revealed a variable range in transcription levels among mag2 and MSI99 lines. A positive correlation between number of integration events and transcription level was observed (p<0.05). Plants were acclimated and challenged in the greenhouse with either Agrobacterium vitis strains (bacterial crown gall pathogen) at 10(8) cfu/ml or Uncinula necator (fungal powdery mildew pathogen) at 10(5) conidia/ml for evaluation of disease resistance. A total of 6 mag2 and 5 MSI99 lines expressing the antimicrobial genes exhibited significant reductions of crown gall symptoms as compared to non-transformed controls. However, only two mag2 lines showed measurable symptom reductions in response to U. necator, but not strong resistance. Our results suggest that the expression of magainin-type genes in grapevines may be more effective against bacteria than fungi. Additional strategies to enhance transgene expression and the spectrum of resistance to grape diseases are suggested.

The Cecropin Superfamily of Toxic Peptides

All sequenced peptide toxins of the cecropin, pleurocidin and dermaceptin/ceratotoxin families in the National Center for Biotechnology Information (NCBI) database as of May 2005 were identified and shown to comprise a single superfamily. The peptide sequences were multiply aligned, revealing conserved residues that may play roles in structure and function. Signature sequences were derived for each of the 3 constituent families. Phylogenetic analyses revealed the relationships of these peptides to each other, and average hydropathy/amphipathicity studies provided structural information. This study serves to characterize a large superfamily of toxic peptides that perform host defense functions in a range of animal kingdoms.

Evaluating tilt angles of membrane-associated helices: comparison of computational and NMR techniques

A computational method to calculate the orientation of membrane-associated alpha-helices with respect to a lipid bilayer has been developed. It is based on a previously derived implicit membrane representation, which was parameterized using the structures of 46 alpha-helical membrane proteins. The method is validated by comparison with an independent data set of six transmembrane and nine antimicrobial peptides of known structure and orientation. The minimum energy orientations of the transmembrane helices were found to be in good agreement with tilt and rotation angles known from solid-state NMR experiments. Analysis of the free-energy landscape found two types of minima for transmembrane peptides: i), Surface-bound configurations with the helix long axis parallel to the membrane, and ii), inserted configurations with the helix spanning the membrane in a perpendicular orientation. In all cases the inserted configuration also contained the global energy minimum. Repeating the calculations with a set of solution NMR structures showed that the membrane model correctly distinguishes native transmembrane from nonnative conformers. All antimicrobial peptides investigated were found to orient parallel and bind to the membrane surface, in agreement with experimental data. In all cases insertion into the membrane entailed a significant free-energy penalty. An analysis of the contributions of the individual residue types confirmed that hydrophobic residues are the main driving force behind membrane protein insertion, whereas polar, charged, and aromatic residues were found to be important for the correct orientation of the helix inside the membrane.

A molecular dynamics study of the bee venom melittin in aqueous solution, in methanol, and inserted in a phospholipid bilayer

The structural properties of melittin, a small amphipathic peptide found in the bee venom, are investigated in three different environments by molecular dynamics simulation. Long simulations have been performed for monomeric melittin solvated in water, in methanol, and shorter ones for melittin inserted in a dimyristoylphosphatidylcholine bilayer. The resulting trajectories were analysed in terms of structural properties of the peptide and compared to the available NMR data. While in water and methanol solution melittin is observed to partly unfold, the peptide retains its structure when embedded in a lipid bilayer. The latter simulation shows good agreement with the experimentally derived (3)J-coupling constants. Generally, it appears that higher the stability of the helical conformation of melittin, lower is the dielectric permittivity of the environment. In addition, peptide-lipid interactions were investigated showing that the C-terminus of the peptide provides an anchor to the lipid bilayer by forming hydrogen bonds with the lipid head groups.

Interaction of the antimicrobial peptide cyclo(RRWWRF) with membranes by molecular dynamics simulations

Antimicrobial peptides have gained a lot of interest in recent years due to their potential use as a new generation of antibiotics. It is believed that this type of relatively short, amphipathic, cationic peptide targets the bacterial membrane, and destroys the chemical gradients over the membrane via formation of stable or transient pores. Here we use the NMR structure of cyclo(RRWWRF) in a series of molecular dynamics simulations in membranes at various peptide/lipid ratios. We observe that the NMR structure of the peptide is still stable after 100 ns simulation. At a peptide/lipid ratio of 2:128, the membrane is only a little affected compared to a pure dipalmitoylphosphatidylcholine lipid membrane, but at a ratio of 12:128, the water-lipid interface becomes more fuzzy, the water molecules can reach deeper into the hydrophobic core, and the water penetration free-energy barrier changes. Moreover, we observe that the area per lipid decreases and the deuterium order parameters increase in the presence of the peptide. We suggest that the changes in the hydrophobic core, together with the changes in the headgroups, result in an imbalance of the membrane and that it is thus not an efficient hydrophobic barrier in the presence of the peptides, independent of pore formation.

De novo generation of cationic antimicrobial peptides: influence of length and tryptophan substitution on antimicrobial activity

Comparison of human immunodeficiency virus lentiviral lytic peptide 1 with other host-derived peptides indicates that antimicrobial properties of membrane-active peptides are markedly influenced by their cationic, hydrophobic, and amphipathic properties. Many common themes, such as Arg composition of the cationic face of an amphipathic helix and the importance of maintaining the hydrophobic face, have been deduced from these observations. These studies suggest that a peptide with these structural properties can be derived de novo by using only a few strategically positioned amino acids. However, the effects of length and helicity on antimicrobial activity and selectivity have not been objectively evaluated in the context of this motif. To address these structure-function issues, multimers of a 12-residue lytic base unit (LBU) peptide composed only of Arg and Val residues aligned to form idealized amphipathic helices were designed. Bacterial killing assays and circular dichroism analyses reveal a strong correlation between antibacterial activity, peptide length, and propensity to form a helix in solvent mimicking the environment of a membrane. Increasing peptide length beyond two LBUs (24-residue peptides) resulted in no appreciable increase in antimicrobial activity. Derivatives (WLBU) of the LBU series were further engineered by substituting Trp residues in the hydrophobic domains. The 24-residue WLBU2 peptide was active at physiologic NaCl concentrations against Staphylococcus aureus and mucoid and nonmucoid strains of Pseudomonas aeruginosa. Further, WLBU2 displayed the highest antibacterial selectivity of all peptides evaluated in the present study by using a coculture model of P. aeruginosa and primary human skin fibroblasts. These findings provide fundamental information toward the de novo design of an antimicrobial peptide useful for the management of infectious diseases.

Conformation of peptides in lipid membranes studied by x-ray grazing incidence scattering

Although the antimicrobial, fungal peptide alamethicin has been extensively studied, the conformation of the peptide and the interaction with lipid bilayers as well as the mechanism of channel gating are still not completely clear. As opposed to studies of the crystalline state, the polypeptide structures in the environment of fluid bilayers are difficult to probe. We have investigated the conformation of alamethicin in highly aligned stacks of model lipid membranes by synchrotron-based x-ray scattering. The (wide-angle) scattering distribution has been measured by reciprocal space mappings. A pronounced scattering signal is observed in samples of high molar peptide/lipid ratio which is distinctly different from the scattering distribution of an ideal helix in the transmembrane state. Beyond simple models of ideal helices, the data is analyzed in terms of models based on atomic coordinates from the Brookhaven Protein Data Bank, as well as from published molecular dynamics simulations. The results can be explained by assuming a wide distribution of helix tilt angles with respect to the membrane normal and a partial insertion of the N-terminus into the membrane.

Computer simulations of membrane proteins

Computer simulations are rapidly becoming a standard tool to study the structure and dynamics of lipids and membrane proteins. Increasing computer capacity allows unbiased simulations of lipid and membrane-active peptides. With the increasing number of high-resolution structures of membrane proteins, which also enables homology modelling of more structures, a wide range of membrane proteins can now be simulated over time spans that capture essential biological processes. Longer time scales are accessible by special computational methods. We review recent progress in simulations of membrane proteins.

Electrostatic sequestration of PIP2 on phospholipid membranes by basic/aromatic regions of proteins

The basic effector domain of myristoylated alanine-rich C kinase substrate (MARCKS), a major protein kinase C substrate, binds electrostatically to acidic lipids on the inner leaflet of the plasma membrane; interaction with Ca2+/calmodulin or protein kinase C phosphorylation reverses this binding. Our working hypothesis is that the effector domain of MARCKS reversibly sequesters a significant fraction of the L-alpha-phosphatidyl-D-myo-inositol 4,5-bisphosphate (PIP2) on the plasma membrane. To test this, we utilize three techniques that measure the ability of a peptide corresponding to its effector domain, MARCKS(151-175), to sequester PIP2 in model membranes containing physiologically relevant fractions (15-30%) of the monovalent acidic lipid phosphatidylserine. First, we measure fluorescence resonance energy transfer from Bodipy-TMR-PIP2 to Texas Red MARCKS(151-175) adsorbed to large unilamellar vesicles. Second, we detect quenching of Bodipy-TMR-PIP2 in large unilamellar vesicles when unlabeled MARCKS(151-175) binds to vesicles. Third, we identify line broadening in the electron paramagnetic resonance spectra of spin-labeled PIP2 as unlabeled MARCKS(151-175) adsorbs to vesicles. Theoretical calculations (applying the Poisson-Boltzmann relation to atomic models of the peptide and bilayer) and experimental results (fluorescence resonance energy transfer and quenching at different salt concentrations) suggest that nonspecific electrostatic interactions produce this sequestration. Finally, we show that the PLC-delta1-catalyzed hydrolysis of PIP2, but not binding of its PH domain to PIP2, decreases markedly as MARCKS(151-175) sequesters most of the PIP2.

Membrane anchoring of the AgrD N-terminal amphipathic region is required for its processing to produce a quorum-sensing pheromone in Staphylococcus aureus

Quorum-sensing pheromones are signal molecules that are secreted from Gram-positive bacteria and utilized by these bacteria to communicate among individual cells to regulate their activities as a group through a cell density-sensing mechanism. Typically, these pheromones are processed from precursor polypeptides. The mechanisms of trafficking, processing, and modification of the precursor to generate a mature pheromone are unclear. In Staphylococcus aureus, AgrD is the propeptide for an autoinducing peptide (AIP) pheromone that triggers the Agr cell density-sensing system upon reaching a threshold and subsequently regulates expression of virulence factor genes. The transmembrane protein AgrB, encoded in the agr locus, is necessary for the processing of AgrD to produce mature AIP; however, it is not clear how AgrD interacts with AgrB and how this interaction results in the generation of mature AIP. In this study, we found that the AgrD propeptide was integrated into the cytoplasmic membrane by a conserved alpha-helical amphipathic motif in its N-terminal region. We demonstrated that membrane targeting of AgrD by this motif was required for the stabilization of AgrD and the production of mature AIP, although this region was not specifically involved in the interaction with AgrB. An artificial amphipathic peptide replacing the N-terminal amphipathic motif of AgrD directed the protein to the cytoplasmic membrane and enabled the production of AIP. Analysis of Bacillus ComX precursor protein sequences suggested that the amphipathic membrane-targeting motif might also exist in pheromone precursors of other Gram-positive bacteria.