The effects of hydrophobic mismatch between phosphatidylcholine bilayers and transmembrane alpha-helical peptides depend on the nature of interfacially exposed aromatic and charged residues

In this study, we investigated the extent to which different aromatic and positively charged side chains, which often flank transmembrane segments of proteins, can influence lipid-peptide interactions. Model systems consisting of phosphatidylcholine and hydrophobic alpha-helical peptides with different flanking residues were investigated. The peptides were incorporated in relatively thick and in relatively thin lipid bilayers to create a peptide-bilayer hydrophobic mismatch, and the compensating effects on lipid structure were analyzed. When relatively long with respect to the thickness of the bilayer, the peptides that are flanked by the aromatic side chains, Trp, Tyr, and Phe, all induce a significant ordering of the lipid acyl chains, while the peptides flanked by the charged residues Lys, Arg, and His do not. However, when the peptides are relatively short with respect to the thickness of the bilayer, their effect on lipid organization does not depend primarily on their aromatic or charged character. Peptides flanked by Trp, Tyr, Lys, or (at low pH) His residues are effective in inducing mismatch-relieving cubic and inverted hexagonal phases, while analogues flanked by Phe, Arg, or (at neutral pH) His residues cannot induce an inverted hexagonal phase. The different responses to mismatch might reflect the different interfacial affinities of the residues that were investigated.     

Distribution of amino acids in a lipid bilayer from computer simulations

We have calculated the distribution in a lipid bilayer of small molecules mimicking 17 natural amino acids in atomistic detail by molecular dynamics simulation. We considered both charged and uncharged forms for Lys, Arg, Glu, and Asp. The results give detailed insight in the molecular basis of the preferred location and orientation of each side chain as well the preferred charge state for ionizable residues. Partitioning of charged and polar side chains is accompanied by water defects connecting the side chains to bulk water. These water defects dominate the energetic of partitioning, rather than simple partitioning between water and a hydrophobic phase. Lys, Glu, and Asp become uncharged well before reaching the center of the membrane, but Arg may be either charged or uncharged at the center of the membrane. Phe has a broad distribution in the membrane but Trp and Tyr localize strongly to the interfacial region. The distributions are useful for the development of coarse-grained and implicit membrane potentials for simulation and structure prediction. We discuss the relationship between the distribution in membranes, bulk partitioning to cyclohexane, and several amino acid hydrophobicity scales.     

Different membrane anchoring positions of tryptophan and lysine in synthetic transmembrane alpha-helical peptides.

Specific interactions of membrane proteins with the membrane interfacial region potentially define protein position with respect to the lipid environment. We investigated the proposed roles of tryptophan and lysine side chains as "anchoring" residues of transmembrane proteins. Model systems were employed, consisting of phosphatidylcholine lipids and hydrophobic alpha-helical peptides, flanked either by tryptophans or lysines. Peptides were incorporated in bilayers of different thickness, and effects on lipid structure were analyzed. Induction of nonbilayer phases and also increases in bilayer thickness were observed that could be explained by a tendency of Trp as well as Lys residues to maintain interactions with the interfacial region. However, effects of the two peptides were remarkably different, indicating affinities of Trp and Lys for different sites at the interface. Our data support a model in which the Trp side chain has a specific affinity for a well defined site near the lipid carbonyl region, while the lysine side chain prefers to be located closer to the aqueous phase, near the lipid phosphate group. The information obtained in this study may further our understanding of the architecture of transmembrane proteins and may prove useful for refining prediction methods for transmembrane segments.     

Distribution of Pentachlorophenol in Phospholipid Bilayers: A Molecular Dynamics Study

Molecular dynamics computer simulations of pentachlorophenol (PCP) in palmitoyl-oleoyl-phosphatidylethanolamine and palmitoyl-oleoyl-phosphatidylcholine lipid bilayers were carried out to investigate the distribution of PCP and the effects of PCP on the phospholipid bilayer structure. Starting from two extreme starting structures, including PCP molecules outside the lipid bilayer, the PCP distribution converges in simulations of up to 50 ns. PCP preferentially occupies the region between the carbonyl groups and the double bonds in the acyl chains of the lipid molecules in the bilayer. In the presence of PCP, the lipid chain order increases somewhat in both chains, and the average tilt angle of the lipid chains decreases. The increase in the lipid chain order in the presence of PCP was more pronounced in the palmitoyl-oleoyl-phosphatidylcholine bilayer compared to the palmitoyl-oleoyl-phosphatidylethanolamine bilayer. The number of trans conformations of lipid chain dihedrals does not change significantly. PCP aligns parallel to the alkyl chains of the lipid to optimize the packing in the dense ordered chain region of the bilayer. The hydroxyl group of PCP forms hydrogen bonds with both water and lipid oxygen atoms in the water/lipid interface region.     

Molecular dynamics simulation of a palmitoyl-oleoyl phosphatidylserine bilayer with Na+ counterions and NaCl

Two 40 ns molecular dynamics simulations of a palmitoyl-oleoyl phosphatidylserine (POPS) lipid bilayer in the liquid crystalline phase with Na(+) counterions and NaCl were carried out to investigate the structure of the negatively charged lipid bilayer and the effect of salt (NaCl) on the lipid bilayer structure. Na(+) ions were found to penetrate deep into the ester region of the water/lipid interface of the bilayer. Interaction of the Na(+) ions with the lipid bilayer is accompanied by a loss of water molecules around the ion and a simultaneous increase in the number of ester carbonyl oxygen atoms binding the ion, which define an octahedral and square pyramidal geometry. The amine group of the lipid molecule is involved in the formation of inter- and intramolecular hydrogen bonds with the carboxylate and the phosphodiester groups of the lipid molecule. The area per lipid of the POPS bilayer is unaffected by the presence of 0.15M NaCl. There is a small increase in the order parameter of carbon atoms in the beginning of the alkyl chain in the presence of NaCl. This is due to a greater number of Na(+) ions being coordinated by the ester carbonyl oxygen atoms in the water/lipid interface region of the POPS bilayer.     

Effect of lysine side chain length on intra-helical glutamate--lysine ion pairing interactions

Ion-pairing interactions are important for protein stabilization. Despite the apparent electrostatic nature of these interactions, natural positively charged amino acids Lys and Arg have multiple methylenes linking the charged functionality to the backbone. Interestingly, the amino acids Lys and Orn have positively charged side chains that differ by only one methylene. However, only Lys is encoded and incorporated into proteins. To investigate the effect of side chain length of Lys on ion-pairing interactions, a series of 12 monomeric alpha-helical peptides containing potential Glu-Xaa (i, i+3), (i, i+4) and (i, i+5) (Xaa = Lys, Orn, Dab, Dap) interactions were studied by circular dichroism (CD) spectroscopy at pH 7 and 2. At pH 7, no Glu-Xaa (i, i+5) interaction was observed, regardless of the Xaa side chain length. Furthermore, only Lys was capable of supporting Glu-Xaa (i, i+3) interactions, whereas any Xaa side chain length supported Glu-Xaa (i, i+4) interactions. Side chain conformational analysis by molecular mechanics calculations showed that the side chain length of Lys enables the Glu-Xaa (i, i+3) interaction with lower energy conformations compared to residues with side chain lengths shorter than that of Lys. Furthermore, these calculated low energy conformers were consistent with conformations of intra-helical Glu-Lys salt bridges in a non-redundant protein structure database. Importantly, the CD spectra for peptides with Glu-Lys interactions did not alter significantly upon changing the pH because of a greater contribution to these interactions by forces other than electrostatics. Incorporating side chains just one methylene shorter (Orn) resulted in significant pH dependence or lack of interaction, suggesting that nature has chosen Lys to form durable interactions with negatively charged functional groups.     

Lipid-protein interactions of integral membrane proteins: a comparative simulation study

The interactions between membrane proteins and their lipid bilayer environment play important roles in the stability and function of such proteins. Extended (15-20 ns) molecular dynamics simulations have been used to explore the interactions of two membrane proteins with phosphatidylcholine bilayers. One protein (KcsA) is an alpha-helix bundle and embedded in a palmitoyl oleoyl phosphatidylcholine bilayer; the other (OmpA) is a beta-barrel outer-membrane protein and is in a dimyristoyl phosphatidylcholine bilayer. The simulations enable analysis in detail of a number of aspects of lipid-protein interactions. In particular, the interactions of aromatic amphipathic side chains (i.e., Trp, Tyr) with lipid headgroups, and "snorkeling" interactions of basic side chains (i.e., Lys, Arg) with phosphate groups are explored. Analysis of the number of contacts and of H-bonds reveal fluctuations on an approximately 1- to 5-ns timescale. There are two clear bands of interacting residues on the surface of KcsA, whereas there are three such bands on OmpA. A large number of Arg-phosphate interactions are seen for KcsA; for OmpA, the number of basic-phosphate interactions is smaller and shows more marked fluctuations with respect to time. Both classes of interaction occur in clearly defined interfacial regions of width approximately 1 nm. Analysis of lateral diffusion of lipid molecules reveals that "boundary" lipid molecules diffuse at about half the rate of bulk lipid. Overall, these simulations present a dynamic picture of lipid-protein interactions: there are a number of more specific interactions but even these fluctuate on an approximately 1- to 5-ns timescale.     

Bilayer conformation of fusion peptide of influenza virus hemagglutinin: a molecular dynamics simulation study

Unraveling the conformation of membrane-bound viral fusion peptides is essential for understanding how those peptides destabilize the bilayer topology of lipids that is important for virus-cell membrane fusion. Here, molecular dynamics (MD) simulations were performed to investigate the conformation of the 20 amino acids long fusion peptide of influenza hemagglutinin of strain X31 bound to a dimyristoyl phosphatidylcholine (DMPC) bilayer. The simulations revealed that the peptide adopts a kinked conformation, in agreement with the NMR structures of a related peptide in detergent micelles. The peptide is located at the amphipathic interface between the headgroups and hydrocarbon chains of the lipid by an energetically favorable arrangement: The hydrophobic side chains of the peptides are embedded into the hydrophobic region and the hydrophilic side chains are in the headgroup region. The N-terminus of the peptide is localized close to the amphipathic interface. The molecular dynamics simulations also revealed that the peptide affects the surrounding bilayer structure. The average hydrophobic thickness of the lipid phase close to the N-terminus is reduced in comparison with the average hydrophobic thickness of a pure dimyristoyl phosphatidylcholine bilayer.     

Molecular dynamics simulation of melittin in a dimyristoylphosphatidylcholine bilayer membrane.

Molecular dynamics trajectories of melittin in an explicit dimyristoyl phosphatidylcholine (DMPC) bilayer are generated to study the details of lipid-protein interactions at the microscopic level. Melittin, a small amphipathic peptide found in bee venom, is known to have a pronounced effect on the lysis of membranes. The peptide is initially set parallel to the membrane-solution interfacial region in an alpha-helical conformation with unprotonated N-terminus. Solid-state nuclear magnetic resonance (NMR) and polarized attenuated total internal reflectance Fourier transform infrared (PATIR-FTIR) properties of melittin are calculated from the trajectory to characterize the orientation of the peptide relative to the bilayer. The residue Lys7 located in the hydrophobic moiety of the helix and residues Lys23, Arg24, Gln25, and Gln26 at the C-terminus hydrophilic form hydrogen bonds with water molecules and with the ester carbonyl groups of the lipids, suggesting their important contribution to the stability of the helix in the bilayer. Lipid acyl chains are closely packed around melittin, contributing to the stable association with the membrane. Calculated density profiles and order parameters of the lipid acyl chains averaged over the molecular dynamics trajectory indicate that melittin has effects on both layers of the membrane. The presence of melittin in the upper layer causes a local thinning of the bilayer that favors the penetration of water through the lower layer. The energetic factors involved in the association of melittin at the membrane surface are characterized using an implicit mean-field model in which the membrane and the surrounding solvent are represented as structureless continuum dielectric material. The results obtained by solving the Poisson-Bolztmann equation numerically are in qualitative agreement with the detailed dynamics. The influence of the protonation state of the N-terminus of melittin is examined. After 600 ps, the N-terminus of melittin is protonated and the trajectory is continued for 400 ps, which leads to an important penetration of water molecules into the bilayer. These observations provide insights into how melittin interacts with membranes and the mechanism by which it enhances their lysis.     

Relative reactivity of lysine and other peptide-bound amino acids to oxidation by hypochlorite

Antibacterial and inflammatory responses of neutrophils and macrophages produce hypochlorite as a major oxidant. Numerous side chains of amino acids found in extracellular proteins can be modified by hypochlorite, including His, Arg, Tyr, Lys, Trp, and Met. We studied the relative reactivity of each of these amino acid residues in short N-blocked peptides, where other residues in the peptide were highly resistant to hypochlorite attack. Hypochlorite treatment led to modified peptides in each case, which were detected by changes in retention on reversed-phase HPLC. A distinct single product, consuming two equivalents of hypochlorite per equivalent of peptide, was obtained from the Lys-containing peptides. UV spectroscopy, nuclear magnetic resonance (NMR), and electrospray/mass spectroscopy identified this product as the dichloramine at the epsilon-amino group of the Lys side chain. The dichloramine at Lys did not decompose to form a detectable amount of carbonyl reactive with dinitrophenylhydrazine. The dichloramine at Lys did however quantitatively revert back to Lys during HCl digestion of the tetrapeptide for amino acid analysis, with simultaneous modification of the adjacent Phe residue. The formation of the dichloramine at Lys was not blocked by peptides or acetylated amino acids that contained Tyr, His, or Arg. In contrast, the presence of equimolar Met-containing peptide, or N-Acetyl-Trp, both inhibited the formation of the dichloramine at Lys. Thus, Met and Trp side chains of proteins might be able to protect Lys from chloramine formation under some circumstances, but this interpretation must consider that Met and Trp are typically found in relatively inaccessible hydrophobic sites, whereas lysine is typically exposed on the protein surface. The hierarchy of amino acid reactivities examined here will aid in the prediction of residues in biological samples most likely to be modified by hypochlorite.     

Investigating the cationic side chains of the antimicrobial peptide tritrpticin: hydrogen bonding properties govern its membrane-disruptive activities

The positively charged side chains of cationic antimicrobial peptides are generally thought to provide the initial long-range electrostatic attractive forces that guide them towards the negatively charged bacterial membranes. Peptide analogs were designed to examine the role of the four Arg side chains in the cathelicidin peptide tritrpticin (VRRFPWWWPFLRR). The analogs include several noncoded Arg and Lys derivatives that offer small variations in side chain length and methylation state. The peptides were tested for bactericidal and hemolytic activities, and their membrane insertion and permeabilization properties were characterized by leakage assays and fluorescence spectroscopy. A net charge of +5 for most of the analogs maintains their high antimicrobial activity and directs them towards preferential insertion into model bacterial membrane systems with a similar extent of burial of the Trp side chains. However the peptides exhibit significant functional differences. Analogs with methylated cationic side chains cause lower levels of membrane leakage and are associated with lower hemolytic activities, making them potentially attractive pharmaceutical candidates. Analogs containing the Arg guanidinium groups cause more membrane disruption than those containing the Lys amino groups. Peptides in the latter group with shorter side chains have increased membrane activity and conversely, elongating the Arg residue causes slightly higher membrane activity. Altogether, the potential for strong hydrogen bonding between the four positive Arg side chains with the phospholipid head groups seems to be a determinant for the membrane disruptive properties of tritrpticin and many related cationic antimicrobial peptides.     

Helicity of short E-R/K peptides

Understanding the secondary structure of peptides is important in protein folding, enzyme function, and peptide‐based drug design. Previous studies of synthetic Ala‐based peptides (>12 a.a.) have demonstrated the role for charged side chain interactions involving Glu/Lys or Glu/Arg spaced three (i, i + 3) or four (i, i + 4) residues apart. The secondary structure of short peptides (<9 a.a.), however, has not been investigated. In this study, the effect of repetitive Glu/Lys or Glu/Arg side chain interactions, giving rise to E‐R/K helices, on the helicity of short peptides was examined using circular dichroism. Short E‐R/K–based peptides show significant helix content. Peptides containing one or more E‐R interactions display greater helicity than those with similar E‐K interactions. Significant helicity is achieved in Arg‐based E‐R/K peptides eight, six, and five amino acids long. In these short peptides, each additional i + 3 and i + 4 salt bridge has substantial contribution to fractional helix content. The E‐R/K peptides exhibit a strongly linear melt curve indicative of noncooperative folding. The significant helicity of these short peptides with predictable dependence on number, position, and type of side chain interactions makes them an important consideration in peptide design.     

N‐terminal diproline and charge group effects on the stabilization of helical conformation in alanine‐based short peptides: CD studies with water and methanol as solvent

Protein folding problem remains a formidable challenge as main chain, side chain and solvent interactions remain entangled and have been difficult to resolve. Alanine‐based short peptides are promising models to dissect protein folding initiation and propagation structurally as well as energetically. The effect of N‐terminal diproline and charged side chains is assessed on the stabilization of helical conformation in alanine‐based short peptides using circular dichroism (CD) with water and methanol as solvent. A1 (Ac–Pro–Pro–Ala–Lys–Ala–Lys–Ala–Lys–Ala–NH₂) is designed to assess the effect of N‐terminal homochiral diproline and lysine side chains to induce helical conformation. A2 (Ac–Pro–Pro–Glu–Glu–Ala–Ala–Lys–Lys–Ala–NH₂) and A3 (Ac–d Pro–Pro–Glu–Glu–Ala–Ala–Lys–Lys–Ala–NH₂) with N‐terminal homochiral and heterochiral diproline, respectively, are designed to assess the effect of Glu...Lys (i , i  + 4) salt bridge interactions on the stabilization of helical conformation. The CD spectra of A1 , A2 and A3 in water manifest different amplitudes of the observed polyproline II (PPII) signals, which indicate different conformational distributions of the polypeptide structure. The strong effect of solvent substitution from water to methanol is observed for the peptides, and CD spectra in methanol evidence A2 and A3 as helical folds. Temperature‐dependent CD spectra of A1 and A2 in water depict an isodichroic point reflecting coexistence of two conformations, PPII and β‐strand conformation, which is consistent with the previous studies. The results illuminate the effect of N‐terminal diproline and charged side chains in dictating the preferences for extended‐β, semi‐extended PPII and helical conformation in alanine‐based short peptides. The results of the present study will enhance our understanding on stabilization of helical conformation in short peptides and hence aid in the design of novel peptides with helical structures. Copyright © 2017 European Peptide Society and John Wiley & Sons, Ltd.     

Crystal structures of human saposins C andD: implications for lipid recognition and membrane interactions

Human saposins are essential proteins required for degradation of sphingolipids and lipid antigen presentation. Despite the conserved structural organization of saposins, their distinct modes of interaction with biological membranes are not fully understood. We describe two crystal structures of human saposin C in an "open" configuration with unusual domain swapped homodimers. This form of SapC dimer supports the "clip-on" model for SapC-induced vesicle fusion. In addition, we present the crystal structure of SapD in two crystal forms. They reveal the monomer-monomer interface for the SapD dimer, which was confirmed in solution by analytical ultracentrifugation. The crystal structure of SapD suggests that side chains of Lys10 and Arg17 are involved in initial association with the preferred anionic biological membranes by forming salt bridges with sulfate or phosphate lipid headgroups.     

Molecular dynamics simulations predict a tilted orientation for the helical region of dynorphin A(1-17) in dimyristoylphosphatidylcholine bilayers

The structural properties of the endogenous opioid peptide dynorphin A(1-17) (DynA), a potential analgesic, were studied with molecular dynamics simulations in dimyristoylphosphatidylcholine bilayers. Starting with the known NMR structure of the peptide in dodecylphosphocholine micelles, the N-terminal helical segment of DynA (encompassing residues 1-10) was initially inserted in the bilayer in a perpendicular orientation with respect to the membrane plane. Parallel simulations were carried out from two starting structures, systems A and B, that differ by 4 A in the vertical positioning of the peptide helix. The complex consisted of approximately 26,400 atoms (dynorphin + 86 lipids + approximately 5300 waters). After >2 ns of simulation, which included >1 ns of equilibration, the orientation of the helical segment of DynA had undergone a transition from parallel to tilted with respect to the bilayer normal in both the A and B systems. When the helix axis achieved a approximately 50 degrees angle with the bilayer normal, it remained stable for the next 1 ns of simulation. The two simulations with different starting points converged to the same final structure, with the helix inserted in the bilayer throughout the simulations. Analysis shows that the tilted orientation adopted by the N-terminal helix is due to specific interactions of residues in the DynA sequence with phospholipid headgroups, water, and the hydrocarbon chains. Key elements are the "snorkel model"-type interactions of arginine side chains, the stabilization of the N-terminal hydrophobic sequence in the lipid environment, and the specific interactions of the first residue, Tyr. Water penetration within the bilayer is facilitated by the immersed DynA, but it is not uniform around the surface of the helix. Many water molecules surround the arginine side chains, while water penetration near the helical surface formed by hydrophobic residues is negligible. A mechanism of receptor interaction is proposed for DynA, involving the tilted orientation observed from these simulations of the peptide in the lipid bilayer.     

Binding and insertion of alpha-helical anti-microbial peptides in POPC bilayers studied by molecular dynamics simulations

We have performed molecular dynamics simulations of the interactions of two alpha-helical anti-microbial peptides, magainin2 and its synthetic analog of MSI-78, with palmitoyl-oleoyl-phosphatidylcholine (POPC) lipid bilayers. We used various initial positions and orientations of the peptide with respect to the lipid bilayer, including a surface-bound state parallel to the interface, a trans-membrane state, and a partially inserted state. Our 20 ns long simulations show that both magainin2 and MSI-78 are most stable in the lipid environment, with the peptide destabilized to different extents in both aqueous and lipid/water interfacial environments. We found that there are strong specific interactions between the lysine residues of the peptides and the lipid head-group regions. MSI-78, owing to its large number of lysines, shows better binding characteristics and overall stability when compared to magainin2. We also find that both peptides destabilize the bilayer environment, as observed by the increase in lipid tail disorder and the induction of local curvature on the lipid head-groups by the peptides. From all the simulations, we conclude that the hydrogen bonding interactions between the lysines of the peptides and the oxygens of the polar lipid head-groups are the strongest and determine the overall peptide binding characteristics to the lipids.     

Electrostatic couplings in OmpA ion-channel gating suggest a mechanism for pore opening

The molecular forces that drive structural transitions between the open and closed states of channels and transporters are not well understood. The gate of the OmpA channel is formed by the central Glu52-Arg138 salt bridge, which can open to form alternate ion pairs with Lys82 and Glu128. To gain deeper insight into the channel-opening mechanism, we measured interaction energies between the relevant side chains by double-mutant cycle analysis and correlated these with the channel activities of corresponding point mutants. The closed central salt bridge has a strong interaction energy of -5.6 kcal mol(-1), which can be broken by forming the open-state salt bridge Glu52-Lys82 (DeltaDeltaG(Inter) = -3.5 kcal mol(-1)) and a weak interaction between Arg138 and Glu128 (DeltaDeltaG(Inter) = -0.6 kcal mol(-1)). A covalent disulfide bond in place of the central salt bridge completely blocks the channel. Growth assays indicate that this gating mechanism could physiologically contribute to the osmoprotection of Escherichia coli cells from environmental stress.     

Solution structure of BmKK4, the first member of subfamily alpha-KTx 17 of scorpion toxins

BmKK4 is a 30 amino acid peptide purified from the venom of the Chinese scorpion Buthus martensi Karsch. It has been classified as the first member of scorpion toxin subfamily alpha-KTx 17. The 3D structure of BmKK4 in solution has been determined by 2D NMR spectroscopy. This toxin adopts a common alpha/beta-motif, but shows a distinctive local conformation. The most novel feature is that the regular arrangements of the side chains of the residues involved in the beta-sheet of BmKK4 are distorted by a classic beta-bulge structure, which involves two residues (Asp18 and Arg19) in the first strand opposite a single residue (Tyr26) in the second strand. The bulge produces two main changes in the structure of the antiparallel beta-sheet: (1) It disrupts the normal alteration of the side chain direction; the side chain of Asp18 turns over to form a salt bridge with that of Arg19. (2) It accentuates the twist of the sheet, and alters the direction of the antiparallel beta-sheet. The unusual structural feature of the toxin is attributed to the shorter peptide segment (Leu15-Arg19) between the third and fourth Cys residues and two unique residues (Asp18 and Arg19) at the position preceding the fourth Cys. In addition, the lower affinity of the peptide for the Kv channel is correlated to the structural features: residue Arg19 instead of a Lys residue at the critical position for binding and the salt bridge formed between residues Arg19 and Asp18.