The structure of the Alzheimer amyloid beta 10-35 peptide probed through replica-exchange molecular dynamics simulations in explicit solvent

The conformational states sampled by the Alzheimer amyloid beta (10-35) (Abeta 10-35) peptide were probed using replica-exchange molecular dynamics (REMD) simulations in explicit solvent. The Abeta 10-35 peptide is a fragment of the full-length Abeta 40/42 peptide that possesses many of the amyloidogenic properties of its full-length counterpart. Under physiological temperature and pressure, our simulations reveal that the Abeta 10-35 peptide does not possess a single unique folded state. Rather, this peptide exists as a mixture of collapsed globular states that remain in rapid dynamic equilibrium with each other. This conformational ensemble is dominated by random coil and bend structures with insignificant presence of an alpha-helical or beta-sheet structure. The 3D structure of Abeta 10-35 is seen to be defined by a salt bridge formed between the side-chains of K28 and D23. This salt bridge is also observed in Abeta fibrils and our simulations suggest that monomeric conformations of Abeta 10-35 contain pre-folded structural motifs that promote rapid aggregation of this peptide.     

Hydrogen bonding in helical polypeptides from molecular dynamics simulations and amide hydrogen exchange analysis: alamethicin and melittin in methanol.

Molecular dynamics simulations of ion channel peptides alamethicin and melittin, solvated in methanol at 27 degrees C, were run with either regular alpha-helical starting structures (alamethicin, 1 ns; melittin 500 ps either with or without chloride counterions), or with the x-ray crystal coordinates of alamethicin as a starting structure (1 ns). The hydrogen bond patterns and stabilities were characterized by analysis of the dynamics trajectories with specified hydrogen bond angle and distance criteria, and were compared with hydrogen bond patterns and stabilities previously determined from high-resolution NMR structural analysis and amide hydrogen exchange measurements in methanol. The two alamethicin simulations rapidly converged to a persistent hydrogen bond pattern with a high level of 3(10) hydrogen bonding involving the amide NH's of residues 3, 4, 9, 15, and 18. The 3(10) hydrogen bonds stabilizing amide NH's of residues C-terminal to P2 and P14 were previously proposed to explain their high amide exchange stabilities. The absence, or low levels of 3(10) hydrogen bonds at the N-terminus or for A15 NH, respectively, in the melittin simulations, is also consistent with interpretations from amide exchange analysis. Perturbation of helical hydrogen bonding in the residues before P14 (Aib10-P14, alamethicin; T11-P14, melittin) was characterized in both peptides by variable hydrogen bond patterns that included pi and gamma hydrogen bonds. The general agreement in hydrogen bond patterns determined in the simulations and from spectroscopic analysis indicates that with suitable conditions (including solvent composition and counterions where required), local hydrogen-bonded secondary structure in helical peptides may be predicted from dynamics simulations from alpha-helical starting structures. Each peptide, particularly alamethicin, underwent some large amplitude structural fluctuations in which several hydrogen bonds were cooperatively broken. The recovery of the persistent hydrogen bonding patterns after these fluctuations demonstrates the stability of intramolecular hydrogen-bonded secondary structure in methanol (consistent with spectroscopic observations), and is promising for simulations on extended timescales to characterize the nature of the backbone fluctuations that underlie amide exchange from isolated helical polypeptides.     

Characterization of a common intermediate of pea lectin in the folding pathway induced by TFE and HFIP

When pea lectin was exposed to a low pH range, it was found that the secondary structure of the lectin resisted conformational changes to a large extent up to pH 2.4 and below this pH, a sharp transition was observed which could be due to the presence of 27 acidic amino acid residues present in the protein. The effects of 1,1,1,3,3,3 hexafluoro-isopropanol (HFIP) and 2,2,2-Trifluoroethanol (TFE) on the conformation of pea lectin at pH 2.4 were studied using circular dichroism and fluorescence spectroscopy. Analysis varying the TFE concentration showed that up to 80% TFE (v/v) protein retained the residual beta-structure accompanied by a loss in tertiary structure. A similar conformation is presumed to exist at 4% HFIP (v/v), with an increase in HFIP concentration structural rearrangements occurred and a transition from beta-structure to alpha-helical structure started from 12% HFIP which completed at 30% HFIP. Our studies show the occurrence of a common intermediate in the folding pathway of pea lectin induced by two different fluoroalcohols, which differ in their mode of action to stabilize the secondary structure of a given protein. While TFE was not found to induce any alpha-helical structure, HFIP caused the transition of pea lectin, which is predominantly a beta-sheet protein, to a structure rich in alpha-helical contacts. Thus, our results also point out the possibility of a non-hierarchical model of protein folding in lectins.     

Cooperative alpha-helix formation of beta-lactoglobulin and melittin induced by hexafluoroisopropanol.

Alcohols denature the native state of proteins, and also stabilize the alpha-helical conformation in unfolded proteins and peptides. Among various alcohols, trifluoroethanol (TFE) and hexafluoroisopropanol (HFIP) are often used because of their high potential to induce such effects. However, the reason why TFE and HFIP are more effective than other alcohols is unknown. Using CD, we studied the effects of TFE and HFIP as well as reference alcohols, i.e., methanol, ethanol, and isopropanol, on the conformation of bovine beta-lactoglobulin and the bee venom melittin at pH 2. Upon addition of alcohols, beta-lactoglobulin exhibited a transformation from the native state, consisting of beta-sheets, to the alpha-helical state, whereas melittin folded from the unfolded state to the alpha-helical state. In both cases, the order of effectiveness of alcohols was shown to be: HFIP > TFE > isopropanol > ethanol > methanol. The alcohol-induced transitions were analyzed assuming a two-state mechanism to obtain the m value, a measure of the dependence of the free energy change on alcohol concentration. Comparison of the m values indicates that the high potential of TFE can be explained by the additive contribution of constituent groups, i.e., F atoms and alkyl group. On the other hand, the high potential of HFIP is more than that expected from the additive effects, suggesting that the cooperative formation of micelle-like clusters of HFIP is important.     

The interaction of bee melittin with lipid bilayer membranes

The influence of melittin and the related 8-26 peptide on the stability and electrical properties of bilayer lipid membranes is reported. Melittin, unlike the 8-26 peptide, has a dramatic influence on lipid membranes, causing rupture at dilute concentrations. The circular dichroism of melittin demonstrated that under physiological conditions, in water, melittin is in extended conformation, which is enhanced in aqueous ethanol. However in 'membrane-like' conditions it is essentially alpha-helical. Secondary structure predictions were used to locate possible alpha-helical nucleation centres and a model of melittin was built according to these predictions. It is postulated that melittin causes a wedge effect in membranes.     

Conformational sampling of peptides in cellular environments

Biological systems provide a complex environment that can be understood in terms of its dielectric properties. High concentrations of macromolecules and cosolvents effectively reduce the dielectric constant of cellular environments, thereby affecting the conformational sampling of biomolecules. To examine this effect in more detail, the conformational preference of alanine dipeptide, poly-alanine, and melittin in different dielectric environments is studied with computer simulations based on recently developed generalized Born methodology. Results from these simulations suggest that extended conformations are favored over alpha-helical conformations at the dipeptide level at and below dielectric constants of 5-10. Furthermore, lower-dielectric environments begin to significantly stabilize helical structures in poly-alanine at epsilon = 20. In the more complex peptide melittin, different dielectric environments shift the equilibrium between two main conformations: a nearly fully extended helix that is most stable in low dielectrics and a compact, V-shaped conformation consisting of two helices that is preferred in higher dielectric environments. An additional conformation is only found to be significantly populated at intermediate dielectric constants. Good agreement with previous studies of different peptides in specific, less-polar solvent environments, suggest that helix stabilization and shifts in conformational preferences in such environments are primarily due to a reduced dielectric environment rather than specific molecular details. The findings presented here make predictions of how peptide sampling may be altered in dense cellular environments with reduced dielectric response.     

Effects of solvent on the structure of the Alzheimer amyloid-beta(25-35) peptide

The free energy landscape for folding of the Alzheimer's amyloid-beta(25-35) peptide is explored using replica exchange molecular dynamics in both pure water and in HFIP/water cosolvent. This amphiphilic peptide is a natural by-product of the Alzheimer's amyloid-beta(1-40) peptide and retains the toxicity of its full-length counterpart as well as the ability to aggregate into beta-sheet-rich fibrils. Our simulations reveal that the peptide preferentially populates a helical structure in apolar organic solvent, while in pure water, the peptide adopts collapsed coil conformations and to a lesser extent beta-hairpin conformations. The beta-hairpin is characterized by a type II' beta-turn involving residues G29 and A30 and two short beta-strands involving residues N27, K28, I31, and I32. The hairpin is stabilized by backbone hydrogen-bonding interactions between residues K28 and I31; S26 and G33; and by side-chain-to-side-chain interactions between N27 and I32. Implications regarding the mechanism of aggregation of this peptide into fibrils and the role of the environment in modulating secondary structure are discussed.     

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.     

Calcium-regulated interactions of human alpha-lactalbumin with bee venom melittin

Affinity chromatography, fluorescence and circular dichroism spectroscopy methods have been used to study the interaction of melittin, a 26-residue peptide from bee venom, with Ca2(+)-binding alpha-lactalbumin from human milk. It has been revealed that melittin binds to the apo- and acidic states of alpha-lactalbumin while the presence of Ca2+ makes the interaction essentially weaker. The association constant for the complex of melittin with apo-alpha-lactalbumin determined from spectropolarimetric melittin-titration data is 2 X 10(7) M-1. The complexation of alpha-lactalbumin with melittin decreases its affinity to Ca2+ by three orders of magnitude. The interaction of apo-alpha-lactalbumin with melittin causes some changes in the environment of its aromatic amino acid residues and drastically alters the conformation of melittin, increasing its alpha-helical content but leaving its single tryptophan residue accessible to water. In the case of the acidic state of alpha-lactalbumin the interaction does not induce an increase in alpha-helical content of melittin.     

Thermal unfolding of tetrameric melittin: comparison with the molten globule state of cytochrome c.

Whereas melittin at micromolar concentrations is unfolded under conditions of low salt at neutral pH, it transforms to a tetrameric alpha-helical structure under several conditions, such as high peptide concentration, high anion concentration, or alkaline pH. The anion- and pH-dependent stabilization of the tetrameric structure is similar to that of the molten globule state of several acid-denatured proteins, suggesting that tetrameric melittin might be a state similar to the molten globule state. To test this possibility, we studied the thermal unfolding of tetrameric melittin using far-UV CD and differential scanning calorimetry. The latter technique revealed a broad but distinct heat absorption peak. The heat absorption curves were consistent with the unfolding transition observed by CD and were explainable by a 2-state transition mechanism between the tetrameric alpha-helical state and the monomeric unfolded state. From the peptide or salt-concentration dependence of unfolding, the heat capacity change upon unfolding was estimated to be 5 kJ (mol of tetramer)-1 K-1 at 30 degrees C and decreased with increasing temperature. The observed change in heat capacity was much smaller than that predicted from the crystallographic structure (9.2 kJ (mol of tetramer)-1 K-1), suggesting that the hydrophobic residues of tetrameric melittin in solution are exposed in comparison with the crystallographic structure. However, the results also indicate that the structure is more ordered than that of a typical molten globule state. We consider that the conformation is intermediate between the molten globule state and the native state of globular proteins.     

Binding and reorientation of melittin in a POPC bilayer: Computer simulations

We performed, using an all-atom force field, molecular dynamics computer simulations to study the binding of melittin to the POPC bilayer and its subsequent reorientation in this bilayer. The binding process involves a simultaneous folding and adsorption of the peptide to the bilayer, followed by the creation of a "U shaped" conformation. The reorientation of melittin from the parallel to the perpendicular conformation requires charged residues to cross the hydrophobic core of the bilayer. This is accomplished by a creation of defects in the bilayer that are filled out with water. The defects are caused by peptide charged residues dragging the lipid headgroup atoms along with them, as they reorient. With increased concentration of melittin water defects form stable pores; this makes it easier for the peptide N-terminus to reorient. Our results complement experimental and computational observations of the melittin/lipid bilayer interaction.     

In situ study by polarization modulated Fourier transform infrared spectroscopy of the structure and orientation of lipids and amphipathic peptides at the air-water interface.

Free amphipathic peptides and peptides bound to dimyristoylphosphatidylcholine (DMPC) were studied directly at the air/water interface using polarization modulation infrared reflection absorption spectroscopy (PMIRRAS). Such differential reflectivity measurements proved to be a sensitive and efficient technique to investigate in situ the respective conformations and orientations of lipid and peptide molecules in pure and mixed films. Data obtained for melittin, a natural hemolytic peptide, are compared to those of L15K7, an ideally amphipathic synthetic peptide constituted by only apolar Leu and polar Lys residues. For pure peptidic films, the intensity, shape, and position of the amide I and II bands indicate that the L15K7 peptide adopts a totally alpha-helical structure, whereas the structure of melittin is mainly alpha-helical and presents some unordered domains. The L15K7 alpha-helix axis is oriented essentially parallel to the air-water interface plane; it differs for melittin. When injected into the subphase, L15K7 and melittin insert into preformed expanded DMPC monolayers and can be detected by PMIRRAS, even at low peptide content (> 50 DMPC molecules per peptide). In such conditions, peptides have the same secondary structure and orientation as in pure peptidic films.     

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 ³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.     

Nuclear magnetic resonance studies and molecular dynamics simulations of the solution conformation of a 'designed', alpha-helical peptide.

The complete three-dimensional structure in methanol of an amphipathic alpha-helical peptide, that has been designed by taking into account the three-dimensional structures of small haemolytic peptides, secondary structure prediction algorithms and the well documented literature on alpha-helix stabilizing factors, has been elucidated by two-dimensional NMR spectroscopy. Initially various two-dimensional spectra (COSY, TOCSY, and NOESY) allowed the complete sequence specific assignment of all signals in the 1H spectrum. Consequently trial structures were generated which were then subjected to molecular dynamics simulations using 121 NOE-derived distances and 25 vicinal coupling constant values as structural restraints to give a final set of calculated structures. These structures are in complete agreement with the results of a circular dichroism study and reveal that the peptide adopted a highly ordered alpha-helical conformation. Details of the structure which throw light on future peptide/protein design are discussed.     

Cooperative antimicrobial action of melittin on lipid membranes: A coarse-grained molecular dynamics study

We conducted a series of coarse-grained molecular dynamics (CG-MD) simulations to investigate the complicated actions of melittin, which is an antimicrobial peptide (AMP) derived from honey bee venom, on a lipid membrane. To accurately simulate the AMP action, we developed and used a protein CG model as an extension of the pSPICA force field (FF), which was designed to reproduce several thermodynamic quantities and structural properties. At a low peptide-to-lipid (P/L) ratio (1/102), no defect was detected. At P/L = 1/51, toroidal pore formation was observed due to collective insertion of multiple melittin peptides from the N-termini. The pore formation was initiated by a local increase in membrane curvature in the vicinity of the peptide aggregate. At a higher P/L ratio (1/26), two more modes were detected, seemingly not controlled by the P/L ratio but by a local arrangement of melittin peptides: 1. Pore formation accompanied by lipid extraction by melittin peptides:a detergent-like mechanism. 2. A rapidly formed large pore in a significantly curved membrane: bursting. Thus, we observed three pore formation modes (toroidal pore formation, lipid extraction, and bursting) depending on the peptide concentration and local arrangement. These observations were consistent with experimental observations and hypothesized melittin modes. Through this study, we found that the local arrangements and population of melittin peptides and the area expansion rate by membrane deformation were key to the initiation of and competition among the multiple pore formation mechanisms.     

Fluoroalcohols induced unfolding of Succinylated Con A: native like beta-structure in partially folded intermediate and alpha-helix in molten globule like state

Concanavalin A (Con A) exists in dimeric state at pH 5. In concentration range 20-60% (v/v) 2,2,2-trifluoroethanol (TFE) and 2-40% (v/v) 1,1,1,3,3,3-hexafluoroisopropanol (HFIP), Con A at pH 5.0 shows visible aggregation. However, when succinyl Con A was used, no aggregation was observed in the entire concentration range of fluoroalcohols (0-90% v/v TFE and HFIP) and resulted in stable alpha-helix formation. Temperature-induced concentration-dependent aggregation in Con A was also found to be prevented/reduced in succinylated form. Possible role of electrostatic repulsion among residues in the prevention of hydrophobically driven aggregation has been discussed. Results indicate that succinylation of a protein resulted in greater stability (in both beta-sheet and alpha-helical forms) against alcohol-induced and temperature-induced concentration-dependent aggregation and this observation may play significant role in amyloid-forming proteins. Effect of TFE and HFIP on the conformation of a dimeric protein, Succinylated Con A, has been investigated by circular dichroism (CD), fluorescence emission spectroscopy, binding of hydrophobic dye ANS (8-anilinonaphthalene-1-sulfonic acid). Far UV-CD, a probe for secondary structure shows loss of native secondary structure in the presence of low concentration of both the alcohols, TFE (10% v/v) and HFIP (4% v/v). Upon addition of higher concentration of these alcohols, Succinylated Con A exhibited transformation from beta-sheet to alpha-helical structure. Intrinsic tryptophan fluorescence studies, ANS binding and near UV-CD experiments indicate the protein is more expanded, have more exposed hydrophobic surfaces and highly disrupted tertiary structure at 60% (v/v) TFE and 30% (v/v) HFIP concentrations. Taken together, these results it might be concluded that TFE and HFIP induce two intermediate states at their low and high concentrations in Succinyl Con A.     

Molten-globule like partially folded states of human serum albumin induced by fluoro and alkyl alcohols at low pH

Human serum albumin (HSA) exists in a molten-globule like state at low pH (pH 2.0). We studied the effects of trifluoroethanol (TFE) and hexafluoroisopropanol (HFIP) on the acid-denatured state of HSA by far-UV circular dichroism (CD), near-UV CD, tryptophan fluorescence, and 1-anilinonaphthalene-8-sulfonic acid (ANS) binding. At pH 2.0, these alcohols induced the formation of alpha-helical structure as evident from the increase in mean residue ellipticity (MRE) value at 222 nm. On addition of different alcohols, HSA exhibited a transition from the acid-denatured state to the alpha-helical state and loss of ANS-binding sites reflected by the decrease in ANS fluorescence at 480 nm. However, the concentration range required to bring about the transition varied greatly among different alcohols. HFIP was found to have highest potential whereas methanol was least effective in inducing the transition. The order of effectiveness of alcohols was shown to be: HFIP > TFE > 2-chloroethanol > tert-butanol > isopropanol > ethanol > methanol as evident from the Cm values. The near-UV CD spectra and tryptophan fluorescence showed the differential effects of halogenated alcohols with those of alkanols. A comparison of the m values, showing the dependence of Delta GH on alcohol concentration, suggests that the helix stabilizing potential of different alcohols is due to the additive effect of different constituent groups present whereas remarkably higher potential of HFIP involves some other factor in addition to the contribution of constituent groups.     

Structure and solvation of melittin in 1,1,1,3,3,3-hexafluoro-2-propanol/water

Fluorinated alcohols can induce peptides and proteins to take up helical conformations. Nuclear Overhauser effect (NOE) spectroscopy experiments and analysis of C(alpha)H proton chemical shifts show that the conformation of melittin in 35% hexafluoro-2-propanol/water is alpha-helical from residues Ile-2 to Val-8 and from Leu-13 to Gln-25. As has been found in other solvent systems, the two helical regions are not colinear; the interhelix angle (73 +/- 15 degrees ) in 35% 1,1,1,3,3,3-hexafluoro-2-propanol/water is smaller than the angle found in other fluoroalcohol-water mixtures or in the crystal. Intermolecular (1)H(19)F and (1)H(1)H nuclear Overhauser effects were used to explore interaction of solvent components with melittin dissolved in this solvent mixture. The NOEs observed indicate that fluoroalcohol and water molecules are both tightly bound to the peptide in the vicinity of the interhelix bend. For the remainder of the molecule, solute-solvent NOEs are consistent with preferential solvation of the peptide by the fluoroalcohol component of the solvent mixture.     

Effect of urea on peptide conformation in water: molecular dynamics and experimental characterization

Molecular dynamics simulations of a ribonuclease A C-peptide analog and a sequence variant were performed in water at 277 and 300 K and in 8 M urea to clarify the molecular denaturation mechanism induced by urea and the early events in protein unfolding. Spectroscopic characterization of the peptides showed that the C-peptide analog had a high alpha-helical content, which was not the case for the variant. In the simulations, interdependent side-chain interactions were responsible for the high stability of the alpha-helical C-peptide analog in the different solvents. The other peptide displayed alpha-helical unwinding that propagated cooperatively toward the N-terminal. The conformations sampled by the peptides depended on their sequence and on the solvent. The ability of water molecules to form hydrogen bonds to the peptide as well as the hydrogen bond lifetimes increased in the presence of urea, whereas water mobility was reduced near the peptide. Urea accumulated in excess around the peptide, to which it formed long-lived hydrogen bonds. The unfolding mechanisms induced by thermal denaturation and by urea are of a different nature, with urea-aqueous solutions providing a better peptide solvation than pure water. Our results suggest that the effect of urea on the chemical denaturation process involves both the direct and indirect mechanisms.     

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.