Structure of the antimicrobial β-hairpin peptide protegrin-1 in a DLPC lipid bilayer investigated by molecular dynamics simulation

All atom molecular dynamics simulations of the 18-residue β-hairpin antimicrobial peptide protegrin-1 (PG-1, RGGRLCYCRRRFCVCVGR-NH₂) in a fully hydrated dilauroylphosphatidylcholine (DLPC) lipid bilayer have been implemented. The goal of the reported work is to investigate the structure of the peptide in a membrane environment (previously solved only in solution [R.L. Fahrner, T. Dieckmann, S.S.L. Harwig, R.I. Lehrer, D. Eisenberg, J. Feigon, Solution structure of protegrin-1, a broad-spectrum antimicrobial peptide from porcine leukocytes. Chemistry and Biology, 3 (1996) 543-550]), and to delineate specific peptide-membrane interactions which are responsible for the peptide's membrane binding properties. A novel, previously unknown, “kick” shaped conformation of the peptide was detected, where a bend at the C-terminal β-strand of the peptide caused the peptide backbone at residues 16-18 to extend perpendicular to the β-hairpin plane. This bend was driven by a highly persistent hydrogen-bond between the polar peptide side-chain of TYR7 and the unshielded backbone carbonyl oxygen atom of GLY17. The H-bond formation relieves the unfavorable free energy of insertion of polar groups into the hydrophobic membrane core. PG-1 was anchored to the membrane by strong electrostatic binding of the protonated N-terminus of the peptide to the lipid head group phosphate anions. The orientation of the peptide in the membrane, and its influence on bilayer structural and dynamic properties are in excellent agreement with solid state NMR measurements [S. Yamaguchi, T. Hong, A. Waring, R.I. Lehrer, M. Hong, Solid-State NMR Investigations of Peptide-Lipid Interaction and Orientation of a b-Sheet Antimicrobial Peptide, Protegrin, Biochemistry, 41 (2002) 9852-9862]. Importantly, two simulations which started from different initial orientations of the peptide converged to the same final equilibrium orientation of the peptide relative to the bilayer. The kick-shaped conformation was observed only in one of the two simulations.     

Ultrafast folding and molecular dynamics of a linear hydrophobic β-hairpin

The first computational study of the folding and dynamics of a hydrophobic β-hairpin containing a central heterochiral diproline segment is reported. Linear hydrophobic sequences containing centrally positioned diproline motifs, heterochiral (DL/LD) and homochiral (LL/DD)), are investigated for their ability to form β-hairpins. Heterochiral diproline motifs (LD/DL) reveal the formation of stable β-hairpins with the backbone adopting β-turn conformation and the formation of backbone hydrogen bonds with antiparallel cross-strand registry, whereas the homochiral diproline (LL/DD) containing sequences tend to adopt PP_II helix conformation. The competition between the β-turn formation and the backbone H-bond ladder of the antiparallel β-strands in heterochiral diproline containing sequences is employed to validate the hypothesis that β-turn formation precedes inter-strand registry in the folding of a β-hairpin (“zipper” mechanism). The observation of noncanonical hydrogen bonds leads to a folded β-hairpin-like conformation and points to the existence of relatively stable transition state intermediates, between the unfolded (extended) and folded (β-hairpin) states. The MD simulations are in excellent agreement with the experimental studies on the model system and constitute the very first computational investigation of the folding and dynamics of a completely hydrophobic synthetic β-hairpin containing heterogeneous residues of mixed chirality.     

Binding of Aβ peptide creates lipid density depression in DMPC bilayer

Using isobaric–isothermal replica exchange molecular dynamics and all-atom explicit water model we study the impact of Aβ monomer binding on the equilibrium properties of DMPC bilayer. We found that partial insertion of Aβ peptide into the bilayer reduces the density of lipids in the binding “footprint” and indents the bilayer thus creating a lipid density depression. Our simulations also reveal thinning of the bilayer and a decrease in the area per lipid in the proximity of Aβ. Although structural analysis of lipid hydrophobic core detects disordering in the orientations of lipid tails, it also shows surprisingly minor structural perturbations in the tail conformations. Finally, partial insertion of Aβ monomer does not enhance water permeation through the DMPC bilayer and even causes considerable dehydration of the lipid–water interface. Therefore, we conclude that Aβ monomer bound to the DMPC bilayer fails to perturb the bilayer structure in both leaflets. Limited scope of structural perturbations in the DMPC bilayer caused by Aβ monomer may constitute the molecular basis of its low cytotoxicity.     

Folding of β-barrel membrane proteins in lipid bilayers — Unassisted and assisted folding and insertion

In cells, β-barrel membrane proteins are transported in unfolded form to an outer membrane into which they fold and insert. Model systems have been established to investigate the mechanisms of insertion and folding of these versatile proteins into detergent micelles, lipid bilayers and even synthetic amphipathic polymers. In these experiments, insertion into lipid membranes is initiated from unfolded forms that do not display residual β-sheet secondary structure. These studies therefore have allowed the investigation of membrane protein folding and insertion in great detail. Folding of β-barrel membrane proteins into lipid bilayers has been monitored from unfolded forms by dilution of chaotropic denaturants that keep the protein unfolded as well as from unfolded forms present in complexes with molecular chaperones from cells. This review is aimed to provide an overview of the principles and mechanisms observed for the folding of β-barrel transmembrane proteins into lipid bilayers, the importance of lipid–protein interactions and the function of molecular chaperones and folding assistants. This article is part of a Special Issue entitled: Lipid–protein interactions.     

Role of hydrophobic interactions and salt-bridges in β-hairpin folding

β-Hairpins are the simplest form of β-sheets which, due to the presence of long-range interactions, can be considered as tertiary structures. Molecular dynamics simulation is a powerful tool that can unravel whole pathways of protein folding/unfolding at atomic resolution. We have performed several molecular dynamics simulations, to a total of over 250 ns, of a β-hairpin peptide in water using GROMACS. We show that hydrophobic interactions are necessary for initiating the folding of the peptide. Once formed, the peptide is stabilized by hydrogen bonds and disruption of hydrophobic interactions in the folded peptide does not denature the structure. In the absence of hydrophobic interactions, the peptide fails to fold. However, the introduction of a salt-bridge compensates for the loss of hydrophobic interactions to a certain extent. Figure Model of b-hairpin folding: Folding is initiated by hydrophobic interactions (Brown circles). The folded structure, once formed, is stabilized by hydrogen bonds (red lines) and is unaffected by loss of hydrophobic contacts     

Effect of heterochiral inversions on the structure of a β-hairpin peptide

We study computationally a family of β-hairpin peptides with systematically introduced chiral inversions, in explicit water, and we investigate the extent to which the backbone structure is able to fold in the presence of heterochiral perturbations. In contrast to the recently investigated case of a helical peptide, we do not find a monotonic change in secondary structure content as a function of the number of L- to D-inversions. The effects of L- to D-inversions are instead found to be highly position-specific. Additionally, in contrast to the helical peptide, some inversions increase the stability of the folded peptide: in such cases, we compute an increase in β-sheet content in the aqueous solution equilibrium ensemble. However, the tertiary structures of the stable (folded) configurations for peptides for which inversions cause an increase in β-sheet content show differences from one another, as well as from the native fold of the nonchirally perturbed β-hairpin. Our results suggest that although some chiral perturbations can increase folding stability, chirally perturbed proteins may still underperform functionally, given the relationship between structure and function.     

Molecular mechanism of action of β-hairpin antimicrobial peptide arenicin: oligomeric structure in dodecylphosphocholine micelles and pore formation in planar lipid bilayers

The membrane-active, cationic, β-hairpin peptide, arenicin, isolated from marine polychaeta Arenicola marina exhibits a broad spectrum of antimicrobial activity. The peptide in aqueous solution adopts the significantly twisted β-hairpin conformation without pronounced amphipathicity. To assess the mechanism of arenicin action, the spatial structure and backbone dynamics of the peptide in membrane-mimicking media and its pore-forming activity in planar lipid bilayers were studied. The spatial structure of the asymmetric arenicin dimer stabilized by parallel association of N-terminal strands of two β-hairpins was determined using triple-resonance nuclear magnetic resonance (NMR) spectroscopy in dodecylphosphocholine (DPC) micelles. Interaction of arenicin with micelles and its oligomerization significantly decreased the right-handed twist of the β-hairpin, increased its amphipathicity, and led to stabilization of the peptide backbone on a picosecond to nanosecond time scale. Relaxation enhancement induced by water-soluble (Mn(2+)) and lipid-soluble (16-doxylstearate) paramagnetic probes pointed to the dimer transmembrane arrangement. Qualitative NMR and circular dichroism study of arenicin-2 in mixed DPC/1,2-dioleoyl-sn-glycero-3-phosphoglycerol bicelles, sodium dodecyl sulfate micelles, and lipid vesicles confirmed that a similar dimeric assembly of the peptide was retained in membrane-mimicking systems containing negatively charged lipids and detergents. Arenicin-induced conductance was dependent on the lipid composition of the membrane. Arenicin low-conductivity pores were detected in the phosphatidylethanolamine-containing lipid mixture, whereas the high-conductivity pores were observed in an exclusively anionic lipid system. The measured conductivity levels agreed with the model in which arenicin antimicrobial activity was mediated by the formation of toroidal pores assembled of two, three, or four β-structural peptide dimers and lipid molecules. The structural transitions involved in arenicin membrane-disruptive action are discussed.     

Vectorial insertion of a β-helical peptide into membrane: a theoretical study on polytheonamide B

Spontaneous unidirectional, or vectorial, insertion of transmembrane peptides is a fundamental biophysical process for toxin and viral actions. Polytheonamide B (pTB) is a potent cytotoxic peptide with a β^6.3-helical structure. Previous experimental studies revealed that the pTB inserts into the membrane in a vectorial fashion and forms a channel with its single molecular length long enough to span the membrane. Also, molecular dynamics simulation studies demonstrated that the pTB is prefolded in aqueous solution. These are unique features of pTB because most of the peptide toxins form channels through oligomerization of transmembrane helices. Here, we performed all-atom molecular dynamics simulations to examine the dynamic mechanism of the vectorial insertion of pTB, providing underlying elementary processes of the membrane insertion of a prefolded single transmembrane peptide. We find that the insertion of pTB proceeds with only the local lateral compression of the membrane in three successive phases: “landing,” “penetration,” and “equilibration” phases. The free energy calculations using the replica-exchange umbrella sampling simulations present an energy cost of 4.3 kcal/mol at the membrane surface for the membrane insertion of pTB from bulk water. The trajectories of membrane insertion revealed that the insertion process can occur in two possible pathways, namely “trapped” and “untrapped” insertions; in some cases, pTB is trapped in the upper leaflet during the penetration phase. Our simulations demonstrated the importance of membrane anchoring by the hydrophobic N-terminal blocking group in the landing phase, leading to subsequent vectorial insertion.     

Internal and environmental effects on folding and dimerisation of Alzheimer's β-amyloid peptide

Amyloid deposits are a hallmark of many diseases. In the case of Alzheimer's disease, a turn between 21Ala and 30Ala, stabilised by a salt bridge between 22Glu/23Asp and 28Lys, may nucleate folding and aggregation of the amyloid β (Aβ) peptide. In the present paper, we test this hypothesis by studying how salt bridge and turn formation vary with intrinsic and environmental changes, and how these changes affect folding and aggregation of the Aβ-peptide.     

The role of surface electrostatics on the stability,function and regulation of human cystathionine β-synthase,a complex multidomain and oligomeric protein

Human cystathionine β-synthase (hCBS) is a key enzyme of sulfur amino acid metabolism, controlling the commitment of homocysteine to the transsulfuration pathway and antioxidant defense. Mutations in hCBS cause inherited homocystinuria (HCU), a rare inborn error of metabolism characterized by accumulation of toxic homocysteine in blood and urine. hCBS is a complex multidomain and oligomeric protein whose activity and stability are independently regulated by the binding of S-adenosyl-methionine (SAM) to two different types of sites at its C-terminal regulatory domain. Here we study the role of surface electrostatics on the complex regulation and stability of hCBS using biophysical and biochemical procedures. We show that the kinetic stability of the catalytic and regulatory domains is significantly affected by the modulation of surface electrostatics through noticeable structural and energetic changes along their denaturation pathways. We also show that surface electrostatics strongly affect SAM binding properties to those sites responsible for either enzyme activation or kinetic stabilization. Our results provide new insight into the regulation of hCBS activity and stability in vivo with implications for understanding HCU as a conformational disease. We also lend experimental support to the role of electrostatic interactions in the recently proposed binding modes of SAM leading to hCBS activation and kinetic stabilization.     

Binding of phospholipids to β-Lactoglobulin and their transfer to lipid bilayers

The bovine milk lipocalin, β-Lactoglobulin (β-LG), has been associated with the binding and transport of small hydrophobic and amphiphilic compounds, whereby it is proposed to increase their bioavailability. We have studied the binding of the fluorescent phospholipid-derivative, NBD-didecanoylphosphatidylethanolamine (NBD-diC₁₀PE) to β-LG by following the increase in amphiphile fluorescence upon binding to the protein using established methods. The equilibrium association constant, K_B, was (1.2 ± 0.2) × 10⁶ M^− 1 at 25 °C, pH 7.4 and I = 0.15 M. Dependence of K_B on pH and on the monomer-dimer equilibrium of β-LG gave insight on the nature of the binding site which is proposed to be the hydrophobic calyx formed by the β-barrel in the protein. The monomer-dimer equilibrium of β-LG was re-assessed using fluorescence anisotropy of Tryptophan. The equilibrium constant for dimerization, K_D, was (7.0 ± 1.5) × 10⁵ M^− 1 at 25 °C, pH 7.4, and 0.15 M ionic strength. The exchange of NBD-diC₁₀PE between β-LG and POPC lipid bilayers was followed by the change in NBD fluorescence. β-LG was shown to be a catalyst of phospholipid exchange between lipid bilayers, the mechanism possibly involving adsorption of the protein at the bilayer surface.     

Binding of small alcohols to a lipid bilayer membrane: does the partitioning coefficient express the net affinity?

The total vapor pressures at 26 degreesC of binary (water-alcohol) and ternary (water-alcohol-vesicle) systems were measured for six short chain alcohols. The vesicles were unilamellar dipalmitoyl phosphatidylcholine (DMPC). The data was used to evaluate the effect of vesicles on the chemical potential of alcohols expressed as the preferential binding parameter of the alcohol-lipid interaction, gamma23. This quantity is a thermodynamic (model-free) measure of the net strength of membrane-alcohol interactions. For the smaller investigated alcohols (methanol, ethanol and 1-propanol) gamma23 was negative. This is indicative of so-called preferential hydration, a condition where the affinity of the membrane for water is higher than the affinity for the alcohol. For the longer alcohols (1-butanol, 1-pentanol, 1-hexanol) gamma23 was positive and increasing with increasing chain length. This demonstrates preferential binding, i.e. enrichment of alcohol in the membrane and a concomitant depletion of the solute in the aqueous bulk. The measured values of gamma23 were compared to the number of alcohol-membrane contacts specified by partitioning coefficients from the literature. It was found that for the small alcohols the number of alcohol-membrane contacts is much larger than the number of preferentially bound solutes. This discrepancy, which is theoretically expected in cases of very weak binding, becomes less pronounced with increasing alcohol chain length, and when the partitioning coefficient exceeds approximately 3 on the molal scale (10(2) in mole fraction units) it vanishes. Based on this, relationships between structural and thermodynamic interpretations of membrane partitioning are discussed.     

Binding of benzo(α)pyrene, ellipticine, and cis-parinaric acid to β-lactoglobulin: Influence of protein modifications

The binding of benzo()pyrene, ellipticine, and cis-parinaric acid to native, esterified, and alkylated -lactoglobulin was followed by enhancement of the ligand fluorescence. Three studied ligands bind to native or modified -lactoglobulin in apparent molar ratios varying between 1/8 and 2/1, with apparent dissociation constants in the range of 10^–8 M for ligand/-lactoglobulin complexes. The studied, chemically modified -lactoglobulin derivatives display higher binding affinities for all studied ligands, cis-parinaric acid excluded. The reductive alkylation of -NH₂ lysyl residues of -lactoglobulin increases the apparent molar ratios of benzo()pyrene and cis-parinaric acid, and decreases it for ellipticine. The esterified and native -lactoglobulin complexed to the investigated ligands display similar stoichiometries. Dynamic light scattering study of ligand--lactoglobulin complexes in solution shows the formation of aggregates: the apparent hydrodynamic radius value of -lactoglobulin dimer (3.4 nm) reaches 49, 46, and 74 nm upon addition and binding of benzo()pyrene, ellipticine, and cis-parinaric acid, respectively.     

AMPs and OMPs: Is the folding and bilayer insertion of β-stranded outer membrane proteins governed by the same biophysical principles as for α-helical antimicrobial peptides?

The folding and function of membrane proteins is controlled not only by specific but also by unspecific interactions with the constituent lipids. In this review, we focus on the influence of the spontaneous lipid curvature on the folding and insertion of peptides and proteins in membranes. Amphiphilic α-helical peptides, as represented by various antimicrobial sequences, are compared with β-barrel proteins, which are found in the outer membrane of Gram-negative bacteria. It has been shown that cationic amphiphilic peptides are always surface-bound in lipids with a negative spontaneous curvature like POPC, i.e. they are oriented parallel to the membrane plane. On the other hand, in lipids like DMPC with a positive curvature, these peptides can get tilted or completely inserted in a transmembrane state. Remarkably, the folding and spontaneous membrane insertion of β-barrel outer membrane proteins also proceeds more easily in lipids with a positive intrinsic curvature, while it is hampered by negative curvature. We therefore propose that a positive spontaneous curvature of the lipids promotes the ability of a surface-bound molecule to insert more deeply into the bilayer core, irrespective of the conformation, size, or shape of the peptide, protein, or folding intermediate. This article is part of a Special Issue entitled: Lipid-protein interactions.     

A knowledge-based potential highlights unique features of membrane α-helical and β-barrel protein insertion and folding

Outer membrane β-barrel proteins differ from α-helical inner membrane proteins in lipid environment, secondary structure, and the proposed processes of folding and insertion. It is reasonable to expect that outer membrane proteins may contain primary sequence information specific for their folding and insertion behavior. In previous work, a depth-dependent insertion potential, E(z) , was derived for α-helical inner membrane proteins. We have generated an equivalent potential for TM β-barrel proteins. The similarities and differences between these two potentials provide insight into unique aspects of the folding and insertion of β-barrel membrane proteins. This potential can predict orientation within the membrane and identify functional residues involved in intermolecular interactions.     

Peptide-lipid interactions of the β-hairpin antimicrobial peptide tachyplesin and its linear derivatives from solid-state NMR

The peptide-lipid interaction of a β-hairpin antimicrobial peptide tachyplesin-1 (TP-1) and its linear derivatives are investigated to gain insight into the mechanism of antimicrobial activity. ³¹P and ²H NMR spectra of uniaxially aligned lipid bilayers of varying compositions and peptide concentrations are measured to determine the peptide-induced orientational disorder and the selectivity of membrane disruption by tachyplesin. The disulfide-linked TP-1 does not cause any disorder to the neutral POPC and POPC/cholesterol membranes but induces both micellization and random orientation distribution to the anionic POPE/POPG membranes above a peptide concentration of 2%. In comparison, the anionic POPC/POPG bilayer is completely unaffected by TP-1 binding, suggesting that TP-1 induces negative curvature strain to the membrane as a mechanism of its action. Removal of the disulfide bonds by substitution of Cys residues with Tyr and Ala abolishes the micellization of POPE/POPG bilayers but retains the orientation randomization of both POPC/POPG and POPE/POPG bilayers. Thus, linear tachyplesin derivatives have membrane disruptive abilities but use different mechanisms from the wild-type peptide. The different lipid-peptide interactions between TP-1 and other β-hairpin antimicrobial peptides are discussed in terms of their molecular structure.     

Adsorption of glycinin and β-conglycinin on silica and cellulose: surface interactions as a function of denaturation, pH, and electrolytes

Soybean proteins have found uses in different nonfood applications due to their interesting properties. We report on the kinetics and extent of adsorption on silica and cellulose surfaces of glycinin and β-conglycinin, the main proteins present in soy. Quartz crystal microgravimetry (QCM) experiments indicate that soy protein adsorption is strongly affected by changes in the physicochemical environment. The affinity of glycinin and the mass adsorbed on silica and cellulose increases (by ca. 13 and 89%, respectively) with solution ionic strength (as it increases from 0 to 100 mM NaCl) due to screening of electrostatic interactions. In contrast, β-conglycinin adsorbs on the same substrates to a lower extent and the addition of electrolyte reduces adsorption (by 25 and 57%, respectively). The addition of 10 mM 2-mercaptoethanol, a denaturing agent, reduces the adsorption of both proteins with a significant effect for glycinin. This observation is explained by the cleavage of disulfide bonds which allows unfolding of the molecules and promotes dissociation into subunits that favors more compact adsorbed layer structures. In addition, adsorption of glycinin onto cellulose decreases with lowering the pH from neutral to pH 3 due to dissociation of the macromolecules, resulting in flatter adsorbed layers. The respective adsorption isotherms fit a Langmuir model and QCM shifts in energy dissipation and frequency reveal multiple-step kinetic processes indicative of changes in adlayer structure.