Aggregation of amphipathic peptides at an aqueous–organic interface using coarse-grained simulations

The effect of an aqueous/organic interface on the folding and aggregation of amphipathic peptides is examined by applying discontinuous molecular dynamics (DMD) simulations combined with an intermediate resolution protein model, PRIME20, to a peptide/interface system. The systems contain 48 (KLLK)₄ peptides in random coil or α-helical conformations interacting with both strong and weak interfaces. In the absence of an interface, most of the oligomers form helical bundles, a small fraction of which convert to β-sheets when the temperature is above the folding transition. Adding a weak interface decreases oligomer formation above the folding temperature and increases it below. Little monolayer formation is observed at the weak interface; instead reversible adsorption increases the local peptide concentration near the interface, promoting helical bundle formation in the aqueous phase below the folding temperature and β-sheet formation above the folding temperature. Introducing a strong interface leads to irreversible adsorption, promoting formation of helical monolayers below the folding temperature and mixed β-sheet/amorphous monolayers above the folding temperature. The (KLLK)₄ peptide is more likely to adsorb to the interface when it is in an α-helical conformation, as opposed to a random coil, because of its larger hydrophobic moment.     

Molecular Dynamics Simulations of Homo-oligomeric Bundles Embedded Within a Lipid Bilayer

Using molecular dynamics simulations, we studied the structure, interhelix interactions, and dynamics of transmembrane proteins. Specifically, we investigated homooligomeric helical bundle systems consisting of synthetic α-helices with either the sequence Ac-(LSLLLSL)₃-NH₂ (LS2) or Ac-(LSSLLSL)₃-NH₂ (LS3). The LS2 and LS3 helical peptides are designed to have amphipathic characteristics that form ion channels in membrane. We simulated bundles containing one to six peptides that were embedded in palmitoyl-oleoyl-phosphatidylcholine (POPC) lipid bilayer and placed between two lamellae of water. We aim to provide a fundamental understanding of how amphipathic helical peptides interact with each other and their dynamical behaviors in different homooligomeric states. To understand structural properties, we examined the helix lengths, tilt angles of individual helices and the entire bundle, interhelix distances, interhelix cross-angles, helix hydrophobic-to-hydrophilic vector projections, and the average number of interhelix hydrophilic (serine–serine) contacts lining the pore of the transmembrane channel. To analyze dynamical properties, we calculated the rotational autocorrelation function of each helix and the cross-correlation of the rotational velocity between adjacent helices. The observed structural and dynamical characteristics show that higher order bundles containing four to six peptides are composed of multiple lower order bundles of one to three peptides. For example, the LS2 channel was found to be stable in a tetrameric bundle composed of a “dimer of dimers.” In addition, we observed that there is a minimum of two strong hydrophilic contacts between a pair of adjacent helices in the dimer to tetramer systems and only one strong hydrophilic interhelix contact in helix pairs of the pentamer and hexamer systems. We believe these results are general and can be applied to more complex ion channels, providing insight into ion channel stability and assembly.     

Scrutiny of the mechanism of small molecule inhibitor preventing conformational transition of amyloid-β42 monomer: insights from molecular dynamics simulations

Alzheimer’s disease (AD) is a progressive neurodegenerative disorder that is characterized by loss of intellectual functioning of brain and memory loss. According to amyloid cascade hypothesis, aggregation of amyloid-β₄₂ (Aβ₄₂) peptide can generate toxic oligomers and their accumulation in the brain is responsible for the onset of AD. In spite of carrying out a large number of experimental studies on inhibition of Aβ₄₂ aggregation by small molecules, the detailed inhibitory mechanism remains elusive. In the present study, comparable molecular dynamics (MD) simulations were performed to elucidate the inhibitory mechanism of a sulfonamide inhibitor C1 (2,5-dichloro-N-(4-piperidinophenyl)-3-thiophenesulfonamide), reported for its in vitro and in vivo anti-aggregation activity against Aβ₄₂. MD simulations reveal that C1 stabilizes native α-helix conformation of Aβ₄₂ by interacting with key residues in the central helix region (13–26) with hydrogen bonds and π–π interactions. C1 lowers the solvent-accessible surface area of the central hydrophobic core (CHC), KLVFF (16–20), that confirms burial of hydrophobic residues leading to the dominance of helical conformation in the CHC region. The binding free energy analysis with MM–PBSA demonstrates that Ala2, Phe4, Tyr10, Gln15, Lys16, Leu17, Val18, Phe19, Phe20, Glu22, and Met35 contribute maximum to binding free energy (?43.1 kcal/mol) between C1 and Aβ₄₂ monomer. Overall, MD simulations reveal that C1 inhibits Aβ₄₂ aggregation by stabilizing native helical conformation and inhibiting the formation of aggregation-prone β-sheet conformation. The present results will shed light on the underlying inhibitory mechanism of small molecules that show potential in vitro anti-aggregation activity against Aβ₄₂.     

Molecular Dynamics Simulations of Homo-oligomeric Bundles Embedded Within a Lipid Bilayer

Using molecular dynamics simulations, we studied the structure, interhelix interactions, and dynamics of transmembrane proteins. Specifically, we investigated homooligomeric helical bundle systems consisting of synthetic α-helices with either the sequence Ac-(LSLLLSL)₃-NH₂ (LS2) or Ac-(LSSLLSL)₃-NH₂ (LS3). The LS2 and LS3 helical peptides are designed to have amphipathic characteristics that form ion channels in membrane. We simulated bundles containing one to six peptides that were embedded in palmitoyl-oleoyl-phosphatidylcholine (POPC) lipid bilayer and placed between two lamellae of water. We aim to provide a fundamental understanding of how amphipathic helical peptides interact with each other and their dynamical behaviors in different homooligomeric states. To understand structural properties, we examined the helix lengths, tilt angles of individual helices and the entire bundle, interhelix distances, interhelix cross-angles, helix hydrophobic-to-hydrophilic vector projections, and the average number of interhelix hydrophilic (serine–serine) contacts lining the pore of the transmembrane channel. To analyze dynamical properties, we calculated the rotational autocorrelation function of each helix and the cross-correlation of the rotational velocity between adjacent helices. The observed structural and dynamical characteristics show that higher order bundles containing four to six peptides are composed of multiple lower order bundles of one to three peptides. For example, the LS2 channel was found to be stable in a tetrameric bundle composed of a “dimer of dimers.” In addition, we observed that there is a minimum of two strong hydrophilic contacts between a pair of adjacent helices in the dimer to tetramer systems and only one strong hydrophilic interhelix contact in helix pairs of the pentamer and hexamer systems. We believe these results are general and can be applied to more complex ion channels, providing insight into ion channel stability and assembly.     

Membrane partitioning of peptide aggregates: coarse-grained molecular dynamics simulations

Abstract

Coarse-grained molecular dynamics (CGMD) simulation technique (MARTINI force field) is applied to monitor the aggregation of helical peptides representing the transmembrane sequence and its extension of bone marrow stromal cell antigen 2 (BST-2). One of the peptides is coupled with a protein transducing domain (PTD) of nine arginine residues (R9) at its N-terminal side as well as a peptide, pep11**, which has been shown to bind to human papilloma virus 16 (HPV16) E6 oncoprotein. A short hydrophobic stretch of the transmembrane domain (TMD) of BST-2 aggregates the fastest and inserts into a lipid membrane. An aggregate of R9-pep11** attaches to the membrane via simultaneous contact of many arginine residues. Monomers from the aggregates of the shortest of the hydrophobic TMDs dissolve into the opposing leaflet when the aggregate spans the bilayer. A ‘flipping’ of the individual monomeric peptides is not observed.

Communicated by Ramaswamy H. Sarma     

Effect of one D‐Leu residue on right‐handed helical ‐L‐Leu‐Aib‐ peptides in the crystal state

Four diastereomeric‐Leu‐Leu‐Aib‐Leu‐Leu‐Aib‐peptides, Boc‐D ‐Leu‐L ‐Leu‐Aib‐L ‐Leu‐L ‐Leu‐Aib‐OMe (1), Boc‐L ‐Leu‐D ‐Leu‐Aib‐L ‐Leu‐L ‐Leu‐Aib‐OMe (2), Boc‐L ‐Leu‐L ‐Leu‐Aib‐D ‐Leu‐L ‐Leu‐Aib‐OMe (3), and Boc‐L ‐Leu‐L ‐Leu‐Aib‐L ‐Leu‐D ‐Leu‐Aib‐OMe (4), were synthesized. The crystals of the four hexapeptides were characterized by X‐ray crystallographic analysis. Two diastereomeric hexapeptides 1 and 2 having D ‐Leu(1) or D ‐Leu(2) were folded into right‐handed (P) 3 ₁₀ ‐helical structures, while peptide 3 having D ‐Leu(4) was folded into a turn structure nucleated by type III′ and I$' \bf{\beta}$ ‐turns, and peptide 4 having D ‐Leu(5) was folded into a left‐handed (M) 3 ₁₀ ‐helical structure. Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

An Evaluation of the Crystal Structure of C-terminal Truncated Apolipoprotein A-I in Solution Reveals Structural Dynamics Related to Lipid Binding

Apolipoprotein (apo) A-I mediates many of the anti-atherogenic functions attributed to high density lipoprotein. Unfortunately, efforts toward a high resolution structure of full-length apoA-I have not been fruitful, although there have been successes with deletion mutants. Recently, a C-terminal truncation (apoA-I^Δ185–243) was crystallized as a dimer. The structure showed two helical bundles connected by a long, curved pair of swapped helical domains. To compare this structure to that existing under solution conditions, we applied small angle x-ray scattering and isotope-assisted chemical cross-linking to apoA-I^Δ185–243 in its dimeric and monomeric forms. For the dimer, we found evidence for the shared domains and aspects of the N-terminal bundles, but not the molecular curvature seen in the crystal. We also found that the N-terminal bundles equilibrate between open and closed states. Interestingly, this movement is one of the transitions proposed during lipid binding. The monomer was consistent with a model in which the long shared helix doubles back onto the helical bundle. Combined with the crystal structure, these data offer an important starting point to understand the molecular details of high density lipoprotein biogenesis.     

Nascent β-Hairpin Formation of a Natively Unfolded Peptide Reveals the Role of Hydrophobic Contacts

Despite the important role of the unfolded states in protein stability, folding, and aggregation, they remain poorly understood due to the lack of residue-specific experimental data. Here, we explore features of the unfolded state of the NTL9 protein by applying all-atom replica-exchange simulations to the two fragment peptides NTL9(1–22) and NTL9(6–17). We found that while NTL9(6–17) is unstructured, NTL9(1–22) transiently folds as various β-hairpins, a fraction of which contain a native β-sheet. Interestingly, despite a large number of charged residues, the formation of backbone hydrogen bonds is concomitant with hydrophobic but not electrostatic contacts. Although the fragment peptides lack a proposed specific contact between Asp⁸ and Lys¹², the individually weak, nonspecific interactions with lysines together stabilize the charged Asp⁸, leading to a pK_a shift of nearly 0.5 units, in agreement with the NMR data. Taken together, our data suggest that the unfolded state of NTL9 likely contains a β-hairpin in segment 1–22 with sequence-distant hydrophobic contacts, thus lending support to a long-standing hypothesis that the unfolded states of proteins exhibit native-like topology with hydrophobic clusters.     

Role of strong anchor residues in the effective binding of 10-mer and 11-mer peptides to HLA-A*2402 molecules

The binding capacity of one-hundred-and-seventy-two 8-mer to 11-mer peptides carrying HLA-A24 anchor residues to HLA-A^*2402 molecules was analyzed by using a HLA class I stabilization assay. Most (76.2%) of these peptides bound to HLA-A^*2402 molecules. These results confirmed previous findings that Tyr and Phe at P2 as well as Phe, Trp, Ile, and Leu at the C-terminus were main anchor residues for HLA-A^*2402. Tyr at P2 was a stronger anchor residue than Phe, while bulky aromatic hydrophobic residues Phe and Trp at the C-terminus are stronger anchors than aliphatic hydrophobic residues Ile and Leu. These results were also supported by an analysis using a panel of mutated 9-mer peptides at P2 and P9. Taken together, these results suggest that HLA-A^*2402 molecules have deep B- and F-pockets because they favor peptides carrying bulky aromatic hydrophobic residues at P2 and the C-terminus. The affinity of 8-mer peptides was significantly lower than that of 9-mer to 11-mer peptides, while there was no difference in affinity between 9-mer, 10-mer, and 11-mer peptides. The affinity of peptides carrying bulky aromatic hydrophobic residues at the C-terminus was higher than that of peptides carrying aliphatic hydrophobic residues in each of the 8-mer to 11-mer peptides, though the greatest difference in affinity was observed in 11-mer peptides. The strong interaction of side chains of these anchor residues with the corresponding pockets may permit the effective binding of 10-mer and 11-mer peptides to HLA-A^*2402 molecules.     

Stereochemistry of Schellman motifs in peptides: Crystal structure of a hexapeptide with a C‐terminus 6 → 1 hydrogen bond

The Schellman motif is a widely observed helix terminating structural motif in proteins, which is generated when the C‐terminus residue adopts a left‐handed helical (α_L) conformation. The resulting hydrogen‐bonding pattern involves the formation of an intramolecular 6 → 1 interaction. This helix terminating motif is readily mimicked in synthetic helical peptides by placing an achiral residue at the penultimate position of the sequence. Thus far, the Schellman motif has been characterized crystallographically only in peptide helices of length 7 residues or greater. The structure of the hexapeptide Boc–Pro–Aib–Gly–Leu–Aib–Leu–OMe in crystals reveal a short helical stretch terminated by a Schellman motif, with the formation of 6 → 1 C‐terminus hydrogen bond. The crystals are in the space group P2₁2₁2₁ with a = 18.155(3) Å, b = 18.864(8) Å, c = 11.834(4) Å, and Z = 4 . The final R₁ and wR₂ values are 7.68 and 14.6%, respectively , for 1524 observed reflections [F_o ≥ 3ς(F_o)]. A 6 → 1 hydrogen bond between Pro(1)CO · · · Leu(6)NH and a 5 → 2 hydrogen bond between Aib(2)CO · · · Aib(5)NH are observed. An analysis of the available oligopeptides having an achiral Aib residue at the penultimate position suggests that chain length and sequence effects may be the other determining factors in formation of Schellman motifs. © 1999 John Wiley & Sons, Inc. Biopoly 50: 13–22, 1999     

Asymmetric inclusion by de novo designed proteins: fluorescence probe studies on amphiphilic alpha-helix bundles

The inclusion feature and supersecondary structure of the de novo designed proteins which are constructed with several amphiphilic alpha-helices and flexible linkage parts were investigated with fluorescence probes. Five types of small proteins (or peptides) have been designed, which are composed of 2, 3, 4, 4, and 6 helices, respectively, and are linked with only linear junctions except for one of 4-helix proteins. All of these proteins have inclusion ability for hydrophobic fluorophores. Further, by the analysis of fluorescence polarization anisotropy, it was suggested that these proteins include guest molecules in compact helix bundles constructed with about 4 helices. Asymmetric inclusion of both monomer and stacked dimer of acridine orange derivatives was found by means of induced circular dichroism except for the 4-helix protein with cross-junction. The chirality of the included dimer proved to be in accordance with the chiral sense of alpha-helical coiled-coil. The 6-helix protein has especially high efficiency in inclusion for any fluorophores examined in this study and brings about a significant blue-shift of maximal emission for 8-anilino-1-naphthalenesulfonate.     

Computational structure–activity study directs synthesis of novel antitumor enkephalin analogs

The capability of a Support Vector Machines QSAR model to predict the antiproliferative ability of small peptides was evaluated by screening a virtual library of enkephalin-like analogs modified by incorporation of the (R,S)-(1-adamantyl)glycine (Aaa) residue. From an initial set of 390 compounds, the peptides, Tyr-Aaa-Gly-Phe-Met (2), Tyr-Aaa-Gly-Phe-Phe (3), Phe-Aaa-Gly-Phe-Phe (4) and Phe-Aaa-Gly-Phe-Met (5) were selected, synthesized and their antitumor activity was tested and compared to that of Met-enkephalin (1). The antiproliferative activity correlated with the computational prediction and with the foldamer-forming ability of the studied peptides. The most active compounds were the hydrophobic peptides, Phe-Aaa-Gly-Phe-Phe (4) and Phe-Aaa-Gly-Phe-Met (5), having a greater propensity to adopt folded structures than the other peptides.     

Conformational Analysis of the β-amyloid Peptide Fragment, β(12–28)

Abstract

NMR and CD spectroscopy have been used to examine the conformation of the peptide, β(12–28), (VHHQKLVFFAEDVGSNK) in aqueous and 60% TFE/40% H₂0 solution at pH 2.4. In 60% TFE solution, the peptide is helical as confirmed by the CD spectrum and by the pattern of the NOE cross peaks detected in the NOESY spectrum of the peptide. In aqueous solution, the peptide adopts a more extended and flexible conformation. Broadening of resonances at low temperature, temperature-dependent changes in the chemical shifts of several of the CH_α resonances and the observation of a number of NOE contacts between the hydrophobic side-chain protons of the peptide are indicative of aggregation in aqueous solution. The behavior of β(12–28) in 60% TFE and in aqueous solution are consistent with the overall conformation and aggregation behavior reported for the larger peptide fragment, β(1–28) and the parent β-amyloid peptide.     

Studies on the metabolic interactions between 4-methylpyrazole and methanol using the monkey as an animal model

Bovine adrenal chromaffin granules have been shown to contain [Met]enkephalin and [Leu]enkephalin and at least seven other small peptides that exhibit specific binding to opiate receptors. Six of these peptides have been characterized and their structures established as (O)-[Met]enkephalin, [Met]enkephalin-Arg⁶-Arg⁷, [Met]enkephalin-Lys⁶, [Met]enkephalin-Arg⁶, [Leu]enkephalin-Arg⁶, and [Met]enkephalin-Arg⁶-Phe⁷. Many of these hexa- and heptapeptides are also present in bovine and human brain. It is suggested that the presence of these peptides in a tissue is evidence of a common biosynthetic pathway of the enkephalins from a large precursor protein.     

A nonhelical,multiple β-turn conformation in a glycine-rich heptapeptide fragment of trichogin A IV containing a single central α-aminoisobutyric acid residue

The conformational properties of the protected seven-residue C-terminal fragment the lipopeptaibol antibiotic Trichogin A IV (Boc-Gly-Gly-Leu-Aib-Gly-Ile-Leu-OMe) has been examined in CDCl₃ and (CD₃)₂SO by ¹H-nmr. Evidence for a multiple β-turn conformation [type I′ at Gly(1)-Gly(2), type II at Leu(3)-Aib(4), and a type I′ at Aib(4)-Gly(5)] suggests that Leu(3) has preferred an extended or semiextended conformation over a helical conformation in CDCl₃. This structure is thus in contrast to earlier observations of seven-residue peptides containing a single central Aib preferring helical conformations in both solution and crystalline slates. A structural transition to a frayed right-handed helix is absented in (CD₃)₂SO. These results suggest that nonhelical conformations may be important in Gly-rich peptides containing Aib. Further, the presence of amino acids with contradictory influences on backbone conformational freedom can lead to well-defined conformational transitions even in small peptides. © 1995 John Wiley & Sons, Inc.     

Studies of the minimum hydrophobicity of alpha-helical peptides required to maintain a stable transmembrane association with phospholipid bilayer membranes

The effects of the hydrophobicity and the distribution of hydrophobic residues on the surfaces of some designed alpha-helical transmembrane peptides (acetyl-K2-L(m)-A(n)-K2-amide, where m + n = 24) on their solution behavior and interactions with phospholipids were examined. We find that although these peptides exhibit strong alpha-helix forming propensities in water, membrane-mimetic media, and lipid model membranes, the stability of the helices decreases as the Leu content decreases. Also, their binding to reversed phase high-performance liquid chromatography columns is largely determined by their hydrophobicity and generally decreases with decreases in the Leu/Ala ratio. However, the retention of these peptides by such columns is also affected by the distribution of hydrophobic residues on their helical surfaces, being further enhanced when peptide helical hydrophobic moments are increased by clustering hydrophobic residues on one side of the helix. This clustering of hydrophobic residues also increases peptide propensity for self-aggregation in aqueous media and enhances partitioning of the peptide into lipid bilayer membranes. We also find that the peptides LA3LA2 [acetyl-K2-(LAAALAA)3LAA-K2-amide] and particularly LA6 [acetyl-K2-(LAAAAAA)3LAA-K2-amide] associate less strongly with and perturb the thermotropic phase behavior of phosphatidylcholine bilayers much less than peptides with higher L/A ratios. These results are consistent with free energies calculated for the partitioning of these peptides between water and phospholipid bilayers, which suggest that LA3LA2 has an equal tendency to partition into water and into the hydrophobic core of phospholipid model membranes, whereas LA6 should strongly prefer the aqueous phase. We conclude that for alpha-helical peptides of this type, Leu/Ala ratios of greater than 7/17 are required for stable transmembrane associations with phospholipid bilayers.     

Insight into the stability of cross-β amyloid fibril from VEALYL short peptide with molecular dynamics simulation

Amyloid fibrils are found in many fatal neurodegenerative diseases such as Alzheimer's disease, Parkinson's disease, type II diabetes, and prion disease. The VEALYL short peptide from insulin has been confirmed to aggregate amyloid-like fibrils. However, the aggregation mechanism of amyloid fibril is poorly understood. Here, we utilized molecular dynamics simulation to analyse the stability of VEALYL hexamer. The statistical results indicate that hydrophobic residues play key roles in stabilizing VEALYL hexamer. Single point and two linkage mutants confirmed that Val1, Leu4, and Tyr5 of VEALYL are key residues. The consistency of the results for the VEALYL oligomer suggests that the intermediate states might be trimer (3-0) and pentamer(3-2). These results can help us to obtain an insight into the aggregation mechanism of amyloid fibril. These methods can be used to study the stability of amyloid fibril from other short peptides.     

Terminal effect of chiral residue on helical screw sense in achiral peptides

To understand the terminal effect of chiral residue for determining a helical screw sense, we adopted five kinds of peptides I – V containing N‐ and/or C‐terminal chiral Leu residue(s): Boc–L ‐Leu–(Aib–ΔPhe)₂–Aib–OMe ( I ), Boc–(Aib–ΔPhe)₂–L ‐Leu–OMe ( II ), Boc–L ‐Leu–(Aib–ΔPhe)₂–L ‐Leu–OMe ( III ), Boc–D ‐Leu–(Aib–ΔPhe)₂–L ‐Leu–OMe ( IV ), and Boc–D ‐Leu–(Aib–ΔPhe)₂–Aib–OMe ( V ). The segment –(Aib–ΔPhe)₂– was used for a backbone composed of two “enantiomeric” (left‐/right‐handed) helices. Actually, this could be confirmed by ¹H‐nmr [nuclear Overhauser effect (NOE) and solvent accessibility of NH resonances] and CD spectroscopy on Boc–(Aib–ΔPhe)₂–Aib–OMe, which took a left‐/right‐handed 3₁₀‐helix. Peptides I – V were also found to take 3₁₀‐type helical conformations in CDCl₃, from difference NOE measurement and solvent accessibility of NH resonances. Chloroform, acetonitrile, methanol, and tetrahydrofuran were used for CD measurement. The CD spectra of peptides I – III in all solvents showed marked exciton couplets with a positive peak at longer wavelengths, indicating that their main chains prefer a left‐handed screw sense over a right‐handed one. Peptide V in all solvents showed exciton couplets with a negative peak at longer wavelengths, indicating it prefers a right‐handed screw sense. Peptide IV in chloroform showed a nonsplit type CD pattern having only a small negative signal around 280 nm, meaning that left‐ and right‐handed helices should exist with almost the same content. In the other solvents, peptide IV showed exciton couplets with a negative peak at longer wavelengths, corresponding to a right‐handed screw sense. From conformational energy calculation and the above ¹H‐nmr studies, an N‐ or C‐terminal L ‐Leu residue in the lowest energy left‐handed 3₁₀‐helical conformation was found to take an irregular conformation that deviates from a left‐handed helix. The positional effect of the L ‐residue on helical screw sense was discussed based on CD data of peptides I – V and of Boc–(L ‐Leu–ΔPhe)_n–L ‐Leu–OMe (n = 2 and 3). © 1999 John Wiley & Sons, Inc. Biopoly 49: 551–564, 1999     

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.