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.     

Conformation and environment of channel-forming peptides: a simulation study

Ion channel-forming peptides enable us to study the conformational dynamics of a transmembrane helix as a function of sequence and environment. Molecular dynamics simulations are used to study the conformation and dynamics of three 22-residue peptides derived from the second transmembrane domain of the glycine receptor (NK4-M2GlyR-p22). Simulations are performed on the peptide in four different environments: trifluoroethanol/water; SDS micelles; DPC micelles; and a DMPC bilayer. A hierarchy of alpha-helix stabilization between the different environments is observed such that TFE/water < micelles < bilayers. Local clustering of trifluoroethanol molecules around the peptide appears to help stabilize an alpha-helical conformation. Single (S22W) and double (S22W,T19R) substitutions at the C-terminus of NK4-M2GlyR-p22 help to stabilize a helical conformation in the micelle and bilayer environments. This correlates with the ability of the W22 and R19 side chains to form H-bonds with the headgroups of lipid or detergent molecules. This study provides a first atomic resolution comparison of the structure and dynamics of NK4-M2GlyR-p22 peptides in membrane and membrane-mimetic environments, paralleling NMR and functional studies of these peptides.     

Molecular dynamics simulations of model trans-membrane peptides in lipid bilayers: a systematic investigation of hydrophobic mismatch

Hydrophobic mismatch, which is the difference between the hydrophobic length of trans-membrane segments of a protein and the hydrophobic width of the surrounding lipid bilayer, is known to play a role in membrane protein function. We have performed molecular dynamics simulations of trans-membrane KALP peptides (sequence: GKK(LA)nLKKA) in phospholipid bilayers to investigate hydrophobic mismatch alleviation mechanisms. By varying systematically the length of the peptide (KALP15, KALP19, KALP23, KALP27, and KALP31) and the lipid hydrophobic length (DLPC, DMPC, and DPPC), a wide range of mismatch conditions were studied. Simulations of durations of 50-200 ns show that under positive mismatch, the system alleviates the mismatch predominantly by tilting the peptide and to a smaller extent by increased lipid ordering in the immediate vicinity of the peptide. Under negative mismatch, alleviation takes place by a combination of local bilayer bending and the snorkeling of the lysine residues of the peptide. Simulations performed at a higher peptide/lipid molar ratio (1:25) reveal slower dynamics of both the peptide and lipid relative to those at a lower peptide/lipid ratio (1:128). The lysine residues have favorable interactions with specific oxygen atoms of the phospholipid headgroups, indicating the preferred localization of these residues at the lipid/water interface.     

Bicyclic peptide inhibitor reveals large contact interface with a protease target

From a large combinatorial library of chemically constrained bicyclic peptides we isolated a selective and potent (K(i) = 53 nM) inhibitor of human urokinase-type plasminogen activator (uPA) and crystallized the complex. This revealed an extended structure of the peptide with both peptide loops engaging the target to form a large interaction surface of 701 ?(2) with multiple hydrogen bonds and complementary charge interactions, explaining the high affinity and specificity of the inhibitor. The interface resembles that between two proteins and suggests that these constrained peptides have the potential to act as small protein mimics.     

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.     

Conformational aspects and rotational dynamics of synthetic adrenocorticotropin-(1-24) and glucagon in reverse micelles

The tryptophan (Trp) rotational dynamics and the secondary structure of the peptide hormones adrenocorticotropin-(1-24) [ACTH(1-24)]--the fully active N-terminal fragment of adrenocorticotropin-(1-39)--and glucagon were studied in aqueous solutions and in reverse micelles of sodium bis(2-ethylhexyl) sulfosuccinate (AOT)/water/isooctane, a system selected to mimic the membrane-water interface. In aqueous solutions, the total fluorescence intensity decays of their single Trp residue [Trp-9 and Trp-25 for ACTH(1-24) and glucagon, respectively] are multiexponential. This is also the case for ACTH(5-10), a fragment of the adrenocorticotropin "message" region. Time-resolved fluorescence anisotropy data evidence a high degree of rotational freedom of the single Trp residue. Transfer of these peptides from water to the aqueous core of reverse micelles induces severe restrictions of the Trp internal motion and of its local environment. The results indicate that the Trp-9 residue in ACTH(1-24 is maintained in the close neighborhood of the water-AOT molecular interface where the water molecules are strongly immobilized. By contrast, the Trp residues in ACTH(5-10) and glucagon are likely to be located closer to the center of the micellar aqueous core where the water molecules are in a more mobile state. Furthermore, the above location of Trp can be extended to the peptide chains themselves as evidenced by the overall correlation time values of the peptide-containing micelles. Nevertheless, in all peptides, the indole ring remains susceptible to oxidation by N-bromosuccinimide. Circular dichroism measurements evidence the induction in glucagon of alpha-helices remaining unaffected by the micellar water content. Conversely, beta-sheet structures are favored in ACTH(1-24) at low water-to-surfactant molar ratios (w0) but are disrupted by subsequent additions of water. These results are discussed in terms of the possible role of the micellar interfaces in selecting the preferred peptide dynamical conformation(s)     

Charge pairing of headgroups in phosphatidylcholine membranes: A molecular dynamics simulation study

Molecular dynamics simulation of the hydrated dimyristoylphosphatidylcholine (DMPC) bilayer membrane in the liquid-crystalline phase was carried out for 5 ns to study the interaction among DMPC headgroups in the membrane/water interface region. The phosphatidylcholine headgroup contains a positively charged choline group and negatively charged phosphate and carbonyl groups, although it is a neutral molecule as a whole. Our previous study (Pasenkiewicz-Gierula, M., Y. Takaoka, H. Miyagawa, K. Kitamura, and A. Kusumi. 1997. J. Phys. Chem. 101:3677-3691) showed the formation of water cross-bridges between negatively charged groups in which a water molecule is simultaneously hydrogen bonded to two DMPC molecules. Water bridges link 76% of DMPC molecules in the membrane. In the present study we show that relatively stable charge associations (charge pairs) are formed between the positively and negatively charged groups of two DMPC molecules. Charge pairs link 93% of DMPC molecules in the membrane. Water bridges and charge pairs together form an extended network of interactions among DMPC headgroups linking 98% of all membrane phospholipids. The average lifetimes of DMPC-DMPC associations via charge pairs, water bridges and both, are at least 730, 1400, and over 1500 ps, respectively. However, these associations are dynamic states and they break and re-form several times during their lifetime.     

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.     

Effect of ionic strength on the organization and dynamics of membrane-bound melittin

Melittin, a cationic hemolytic peptide, is intrinsically fluorescent due to the presence of a single functionally important tryptophan residue. We have previously shown that the sole tryptophan of melittin is localized in a motionally restricted environment in the membrane interface. We have monitored the effect of ionic strength on the organization and dynamics of membrane-bound melittin utilizing fluorescence and circular dichroism (CD) spectroscopic approaches. Our results show that red edge excitation shift (REES) of melittin bound to membranes is sensitive to the change in ionic strength of the medium. This could be attributed to a change in the immediate environment around melittin tryptophan with increasing ionic strength due to differential solvation of ions. Interestingly, the rotational mobility of melittin does not appear to be affected with change in ionic strength. In addition, fluorescence parameters such as lifetime and acrylamide quenching of melittin indicate an increase in water penetration in the membrane interface upon increasing ionic strength. Our results suggest that the solvent dynamics and water penetration in the interfacial region of the membranes are significantly affected at physiologically relevant ionic strength. These results assume significance in the overall context of the influence of ionic strength in the organization and dynamics of membrane proteins and membrane-active peptides.     

Modeling of Peptides in Implicit Membrane-Mimetic Media

Abstract

Lipid bilayer plays a crucial role in folding of membrane peptides and their stabilization in the membrane-bound state. Correct treatment of the media effects is thus essential for realistic simulations of peptides in bilayers. Previously (Volynsky et al., 1999), we proposed an efficient solvation model which mimics heterogeneous membrane-water system. The model is based on combined employment of atomic solvation parameters for water and hydrocarbon, which approximate hydrated headgroups and acyl chains of lipids, respectively. In this study, the model is employed in non-restrained Monte Carlo simulations of several peptides: totally apolar 20-residue poly-L-Leu, hydrophobic peptide with polar edges, and strongly amphiphilic pep-tide. The principal goals are: to explore energy landscape of these peptides in membrane; to characterize the structures of low-energy states and their orientations with respect to the bilayer. Simulations were performed starting from different structures (unordered or helical) and orientations. It was found that the membrane environment significantly promotes an α-helical conformation for all the peptides, while their energetically favourable orientations are quite different. Thus, poly-Leu was immobilized inside the membrane, the hydrophobic peptide with polar termini adapted transbilayer orientation, whereas the amphiphilic peptide stayed on the lipid-water interface in peripherial orientation. Energy barriers between different states were characterized. The computational results were compared with the experimental structural data.     

Dynamics of a membrane-bound tryptophan analog in environments of varying hydration: a fluorescence approach

Tryptophan octyl ester (TOE) represents an important model for membrane-bound tryptophan residues. In this article, we have employed a combination of wavelength-selective fluorescence and time-resolved fluorescence spectroscopies to monitor the effect of varying degrees of hydration on the dynamics of TOE in reverse micellar environments formed by sodium bis(2-ethylhexyl) sulfosuccinate (AOT) in isooctane. Our results show that TOE exhibits red edge excitation shift (REES) and other wavelength-selective fluorescence effects when bound to reverse micelles of AOT. Fluorescence parameters such as intensity, emission maximum, anisotropy, and lifetime of TOE in reverse micelles of AOT depend on [water]/[surfactant] molar ratio (w (o)). These results are relevant and potentially useful for analyzing dynamics of proteins or peptides bound to membranes or membrane-mimetic media under conditions of changing hydration.     

How Lipid Headgroups Sense the Membrane Environment: An Application of 14N NMR

The orientation of lipid headgroups may serve as a powerful sensor of electrostatic interactions in membranes. As shown previously by ²H NMR measurements, the headgroup of phosphatidylcholine (PC) behaves like an electrometer and varies its orientation according to the membrane surface charge. Here, we explored the use of solid-state ¹⁴N NMR as a relatively simple and label-free method to study the orientation of the PC headgroup in model membrane systems of varying composition. We found that ¹⁴N NMR is sufficiently sensitive to detect small changes in headgroup orientation upon introduction of positively and negatively charged lipids and we developed an approach to directly convert the ¹⁴N quadrupolar splittings into an average orientation of the PC polar headgroup. Our results show that inclusion of cholesterol or mixing of lipids with different length acyl chains does not significantly affect the orientation of the PC headgroup. In contrast, measurements with cationic (KALP), neutral (Ac-KALP), and pH-sensitive (HALP) transmembrane peptides show very systematic changes in headgroup orientation, depending on the amount of charge in the peptide side chains and on their precise localization at the interface, as modulated by varying the extent of hydrophobic peptide/lipid mismatch. Finally, our measurements suggest an unexpectedly strong preferential enrichment of the anionic lipid phosphatidylglycerol around the cationic KALP peptide in ternary mixtures with PC. We believe that these results are important for understanding protein/lipid interactions and that they may help parametrization of membrane properties in computational studies.     

Towards understanding structure-stability and surface properties of laminin peptide YIGSR and mutants

Properties of laminin peptide YIGSR and its mutated sequences YIGSD, YIGSS, YIGSN and YIGSQ have been investigated using molecular dynamics simulations (MDS) and Langmuir films at air/water interface. Simulation studies on laminin peptide YIGSR were performed in the isothermal-isobaric (N, P, T) ensemble, with run up to 5 ns in water as well as lipid environment at 298 K. From different initial configurations, shape transformations of the peptides on the timescale of nanoseconds were observed. The results showed YIGSR to be the most stable peptide with the order of minimized energy being YIGSR相似文献   

Protein–cofactor binding and ultrafast electron transfer in riboflavin binding protein under the spatial confinement of nanoscopic reverse micelles

In this contribution, we study the effect of confinement on the ultrafast electron transfer (ET) dynamics of riboflavin binding protein (RBP) to the bound cofactor riboflavin (Rf, vitamin B2), an important metabolic process, in anionic sodium bis(2‐ethylhexyl) sulfosuccinate reverse micelles (AOT‐RMs) of various hydration levels. Notably, in addition to excluded volume effect, various nonspecific interactions like ionic charge of the confining surface can influence the biochemical reactions in the confined environment of the cell. To this end, we have also studied the ET dynamics of RBP–Rf complex under the confinement of a cationic hexadecyltrimethylammonium bromide (CTAB) RMs with similar water pool size to the anionic AOT‐RMs towards simulating equal restricted volume effect. It has been found that the spatial confinement of RBP in the AOT‐RM of w₀ = 10 leads to the loss of its tertiary structure and hence vitamin binding capacity. Although, RBP regains its binding capacity and tertiary structure in AOT‐RMs of w₀ ≥20 due to its complete hydration, the ultrafast ET from RBP to Rf merely occurs in such systems. However, to our surprise, the ET process is found to occur in cationic CTAB‐RMs of similar volume restriction. It is found that under the spatial confinement of anionic AOT‐RM, the isoalloxazine ring of Rf is improperly placed in the protein nanospace so that ET between RBP and Rf is not permitted. This anomaly in the binding behaviour of Rf to RBP in AOT‐RMs is believed to be the influence of repulsive potential of the anionic AOT‐RM surface to the protein. Our finding thus suggests that under similar size restriction, both the hydration and surface charge of the confining volume could have major implication in the intraprotein ET dynamics in real cellular environments. Copyright © 2013 John Wiley & Sons, Ltd.     

Membrane interactions of two arginine-rich peptides with different cell internalization capacities

Cell penetrating peptides (CPPs) can cross cell membranes in a receptor independent manner and transport cargo molecules inside cells. These peptides can internalize through two independent routes: energy dependent endocytosis and energy independent translocation across the membrane, but the exact mechanisms are still unknown. The interaction of the CPP with different membrane components is certainly a preliminary key point that triggers internalization, such as the interaction with lipids to lead to the translocation process. In this study, we used two arginine-rich peptides, RW9 (RRWWRRWRR-NH(2)), which is a potent CPP, and RL9 (RRLLRRLRR-NH(2)) that, although binding tightly and accumulating on membranes, does not enter into cells. Using a set of experimental and theoretical techniques, we studied the binding, insertion and orientation of the peptides into different model membranes as well as the subsequent membrane reorganization. Herein we show that although the two peptides had rather similar behavior regarding lipid membrane interaction, subtle differences were found concerning the depth of peptide insertion, effect on the lipid chain ordering and kinetics of peptide insertion in the membrane, which altogether might explain their different cell internalization capacities. Molecular dynamics simulation studies show that some peptide molecules flipped their orientation over the course of the simulation such that the hydrophobic residues penetrated deeper in the lipid core region while Arg-residues maintained H-bonds with the lipid headgroups, serving as a molecular hinge in a conformation that appeared to correspond to the equilibrium one.     

Transmembrane peptides from tyrosine kinase receptor. Mutation-related behavior in a lipid bilayer investigated by molecular dynamics simulations

Polar mutations in transmembrane alpha helices may alter the structural details of the hydrophobic sequences and control intermolecular contacts. We have performed molecular dynamics simulations on the transmembrane domain of the proto-oncogenic and the oncogenic forms of the Neu receptor in a fluid DMPC bilayer to test whether the Glu mutation which replaces the Val residue at position 664 may alter the helical structure and its insertion in the membrane. The simulations show that the wild and the mutant forms of the transmembrane domain have a different behavior in the bilayer. The native transmembrane sequence is found to be more flexible than in the presence of the Glu mutation, characterized by a tendency to pi deformation to accommodate the helix length to the membrane thickness. The mutant form of this domain does not evidence helical deformation in the present simulation. Hydrophobic matching is achieved both by a larger helix tilt and a vertical shift of the helix towards the membrane interface, favoring the accessibility of the Glu side chain to the membrane environment. A rapid exchange of hydrogen bond interactions with the surrounding water molecules and the lipid headgroups is observed. The difference in the behavior between the two peptides in a membrane environment was also observed experimentally. Both simulation and experimental results agree with the hypothesis that water may act as an intermediate for the formation of cross links between the facing Glu side chains stabilizing the dimer.     

The interactions of antimicrobial peptides derived from lysozyme with model membrane systems

Two peptides, RAWVAWR-NH2 and IVSDGNGMNAWVAWR-NH2, derived from human and chicken lysozyme, respectively, exhibit antimicrobial activity. A comparison between the L-RAWVAWR, D-RAWVAWR, and the longer peptide has been carried out in membrane mimetic conditions to better understand how their interaction with lipid and detergent systems relates to the reported higher activity for the all L-peptide. Using CD and 2D 1H NMR spectroscopy, the structures were studied with DPC and SDS micelles. Fluorescence spectroscopy was used to study peptide interactions with POPC and POPG vesicles and DOPC, DOPE, and DOPG mixed vesicle systems. Membrane-peptide interactions were also probed by ITC and DSC. The ability of fluorescein-labeled RAWVAWR to rapidly enter both E. coli and Staphylococcus aureus was visualized using confocal microscopy. Reflecting the bactericidal activity, the long peptide interacted very weakly with the lipids. The RAWVAWR-NH2 peptides preferred lipids with negatively charged headgroups and interacted predominantly in the solvent-lipid interface, causing significant perturbation of membrane mimetics containing PG headgroups. Peptide structures determined by 1H NMR indicated a well-ordered coiled structure for the short peptides and the C-terminus of the longer peptide. Using each technique, the two enantiomers of RAWVAWR-NH2 interacted in an identical fashion with the lipids, indicating that any difference in activity in vivo is limited to interactions not involving the membrane lipids.     

Novel african trypanocidal agents: membrane rigidifying peptides

The bloodstream developmental forms of pathogenic African trypanosomes are uniquely susceptible to killing by small hydrophobic peptides. Trypanocidal activity is conferred by peptide hydrophobicity and charge distribution and results from increased rigidity of the plasma membrane. Structural analysis of lipid-associated peptide suggests a mechanism of phospholipid clamping in which an internal hydrophobic bulge anchors the peptide in the membrane and positively charged moieties at the termini coordinate phosphates of the polar lipid headgroups. This mechanism reveals a necessary phenotype in bloodstream form African trypanosomes, high membrane fluidity, and we suggest that targeting the plasma membrane lipid bilayer as a whole may be a novel strategy for the development of new pharmaceutical agents. Additionally, the peptides we have described may be valuable tools for probing the biosynthetic machinery responsible for the unique composition and characteristics of African trypanosome plasma membranes.