A molecular dynamics study of the 41-56 beta-hairpin from B1 domain of protein G.

The structural and dynamical behavior of the 41-56 beta-hairpin from the protein G B1 domain (GB1) has been studied at different temperatures using molecular dynamics (MD) simulations in an aqueous environment. The purpose of these simulations is to establish the stability of this hairpin in view of its possible role as a nucleation site for protein folding. The conformation of the peptide in the crystallographic structure of the protein GB1 (native conformation) was lost in all simulations. The new equilibrium conformations are stable for several nanoseconds at 300K (>10 ns), 350 K (>6.5 ns), and even at 450 K (up to 2.5 ns). The new structures have very similar hairpin-like conformations with properties in agreement with available experimental nuclear Overhauser effect (NOE) data. The stability of the structure in the hydrophobic core region during the simulations is consistent with the experimental data and provides further evidence for the role played by hydrophobic interactions in hairpin structures. Essential dynamics analysis shows that the dynamics of the peptide at different temperatures spans basically the same essential subspace. The main equilibrium motions in this subspace involve large fluctuations of the residues in the turn and ends regions. Of the six interchain hydrogen bonds, the inner four remain stable during the simulations. The space spanned by the first two eigenvectors, as sampled at 450 K, includes almost all of the 47 different hairpin structures found in the database. Finally, analysis of the hydration of the 300 K average conformations shows that the hydration sites observed in the native conformation are still well hydrated in the equilibrium MD ensemble.     

A structure for the trimeric MHC class II-associated invariant chain transmembrane domain

The major histocompatibility complex (MHC)-associated invariant chain (Ii) contains a single transmembrane domain that forms trimers. Ii is involved in the assembly of the MHC and antigen presentation, and is thus central to the function of the immune system. Here, we show by attenuated total reflectance, Fourier transform infrared (ATR-FTIR) spectroscopy that the transmembrane domain is alpha-helical and we provide a structural model of the transmembrane domain obtained by a combination of site-specific infrared dichroism and molecular dynamics (MD) simulations. This work resolves the backbone structure of a transmembrane peptide by multiple (13)C=(18)O labelling at ten different residues. A second purely computational approach, based on MD simulations of Ii transmembrane homologous sequences, yields a similar structure that is consistent with our experimental results. The structure presented forms a left-handed coiled coil with an average helix tilt of 13(+/-6) degrees; the residue Gln47 implicated in trimer formation forms strong interhelical contacts, Thr50 points to the inside of the trimeric coil and forms a network of hydrogen bonds.     

Conformational properties of a cyclic peptide bradykinin B2 receptor antagonist using experimental and theoretical methods.

The solution conformation of the cyclic peptide J324 (cyclo0,6-[Lys0,Glu6,D-Phe7]BK), an antagonist targeted at the bradykinin (BK) B2 receptor, has been investigated using experimental and theoretical methods. In order to gain insight into the structural requirements essential for BK antagonism, we carried out molecular dynamics (MD) simulations using simulated annealing as the sampling protocol. Following a free MD simulation we performed simulations using nuclear Overhauser enhancement (NOE) distance constraints determined by NMR experiments. The low-energy structures obtained were compared with each other, grouped into families and analyzed with respect to the presence of secondary structural elements in their backbone. We also introduced new ways of plotting structural data for a more comprehensive analysis of large conformational sets. Finally, the relationship between characteristic backbone conformations and the spatial arrangement of specific pharmacophore centers was investigated.     

The investigation of the secondary structure propensities and free-energy landscapes of peptide ligands by replica exchange molecular dynamics simulations

The conformational states of two peptide sequences that bind to staphylococcal enterotoxin B are sampled by replica exchange molecular dynamic (REMD) simulations in explicit water. REMD simulations were treated with 52 replicas in the range of 280–501 K for both peptides. The conformational ensembles of both peptides are dominated by random coil, bend and turn structures with a small amount of helical structures for each temperature. In addition, while an insignificant presence of β-bridge structures were observed for both peptides, the β-sheet structure was observed only for peptide 3. The results obtained from simulations at 300 K are consistent with the experimental results obtained from circular dichroism spectroscopy. From the analysis of REMD results, we also calculated hydrophobic and hydrophilic solvent accessible surface areas for both peptides, and it was observed that the hydrophobic segments of the peptides tend to form bend or turn structures. Moreover, the free-energy landscapes of both peptides were obtained by principal component analysis to understand how the secondary structural properties change according to their complex space. From the free-energy analysis, we have found several minima for both peptides at decreased temperature. For these obvious minima of both peptides, it was observed that the random coil, bend and turn structures are still dominant and the helix, β-bridge or β-sheet structures can appear or disappear with respect to minima. On the other hand, when we compare the results of REMD with conventional MD simulations for these peptides, the configurations of peptide 3 might be trapped in energy minima during the conventional MD simulations. Hence, it can be said that the REMD simulations have provided a sufficiently high sampling efficiency.     

Proline-induced hinges in transmembrane helices: possible roles in ion channel gating

A number of ion channels contain transmembrane (TM) alpha-helices that contain proline-induced molecular hinges. These TM helices include the channel-forming peptide alamethicin (Alm), the S6 helix from voltage-gated potassium (Kv) channels, and the D5 helix from voltage-gated chloride (CLC) channels. For both Alm and KvS6, experimental data implicate hinge-bending motions of the helix in an aspect of channel gating. We have compared the hinge-bending motions of these TM helices in bilayer-like environments by multi-nanosecond MD simulations in an attempt to describe motions of these helices that may underlie possible modes of channel gating. Alm is an alpha-helical channel-forming peptide, which contains a central kink associated with a Gly-x-x-Pro motif in its sequence. Simulations of Alm in a TM orientation for 10 ns in an octane slab indicate that the Gly-x-x-Pro motif acts as a molecular hinge. The S6 helix from Shaker Kv channels contains a Pro-Val-Pro motif. Modeling studies and recent experimental data suggest that the KvS6 helix may be kinked in the vicinity of this motif. Simulations (10 ns) of an isolated KvS6 helix in an octane slab and in a POPC bilayer reveal hinge-bending motions. A pattern-matching approach was used to search for possible hinge-bending motifs in the TM helices of other ion channel proteins. This uncovered a conserved Gly-x-Pro motif in TM helix D5 of CLC channels. MD simulations of a model of hCLC1-D5 spanning an octane slab suggest that this channel also contains a TM helix that undergoes hinge-bending motion. In conclusion, our simulations suggest a model in which hinge-bending motions of TM helices may play a functional role in the gating mechanisms of several different families of ion channels.     

Conformational interconversion in compstatin probed with molecular dynamics simulations

Compstatin is a 13-residue cyclic peptide that has the potential to become a therapeutic agent against unregulated complement activation. In our effort to understand the structural and dynamic characteristics of compstatin that form the basis for rational and combinatorial optimization of structure and activity, we performed 1-ns molecular dynamics (MD) simulations. We used as input in the MD simulations the ensemble of 21 lowest energy NMR structures, the average minimized structure, and a global optimization structure. At the end of the MD simulations we identified five conformations, with populations ranging between 9% and 44%. These conformations are as follows: 1) coil with alphaR-alphaR beta-turn, as was the conformation of the initial ensemble of NMR structures; 2) beta-hairpin with epsilon-alphaR beta-turn; 3) beta-hairpin with alphaR-alphaR beta-turn; 4) beta-hairpin with alphaR-beta beta-turn; and 5) alpha-helical. Conformational switch was possible with small amplitude backbone motions of the order of 0.1-0.4 A and free energy barrier crossing of 2-11 kcal/mol. All of the 21 MD structures corresponding to the NMR ensemble possessed a beta-turn, with 14 structures retaining the alphaR-alphaR beta-turn type, but the average minimized structure and the global optimization structures were converted to alpha-helical conformations. Overall, the MD simulations have aided to gain insight into the conformational space sampled by compstatin and have provided a measure of conformational interconversion. The calculated conformers will be useful as structural and possibly dynamic templates for optimization in the design of compstatin using structure-activity relations (SAR) or dynamics-activity relations (DAR).     

Equilibrium unfolding of the poly(glutamic acid)20 helix

The equilibrium structural ensemble of a 20-residue polyglutamic acid peptide (E(20)) was studied with FRET, circular dichroism, and molecular dynamics (MD) simulations. A FRET donor, o-aminobenzamide, and acceptor, 3-nitrotyrosine, were introduced at the N- and C-termini, respectively. Circular dichroism, steady state FRET, and time-resolved FRET measurements were employed to characterize the fraction helix and end-to-end distance under different pH conditions: pH 4 (60% alpha-helix), pH 6 (0% alpha-helix), and pH 9 (0% alpha-helix). At pH 4, the end-to-end distance was measured at 24 A and determined to be considerably less than the 31 A predicted for an alpha-helix of the same length. At pH 6 and 9, the end-to-end distance was measured at > 31 and 39 A respectively, both which are determined to be considerably greater than the 27 A predicted for a freely jointed random coil of the same length. To better understand the physical forces underlying the unusual helix-coil transition in this peptide, three theoretical MD models of E(20) were constructed: (1) a pure alpha-helix, (2) an alpha-helix with equivalent attractive intramolecular contacts, and (3) a weak alpha-helix with termini-weighted intramolecular contacts ("sticky ends"). Using MD simulations, the bent helix structure calculated from Model 3 was found to be the closest in agreement with the experimental data.     

Molecular dynamics simulations of a highly charged peptide from an SH3 domain: possible sequence-function relationship

A seven-residue peptide that is highly conserved in SH3 domains despite being far from the active site has been shown by NMR to be stable in solution. This peptide, biologically important because it is a likely folding nucleus for SH3 domains, provides a challenging subject for molecular dynamics because it is highly charged. We present stable, 10-ns simulations of both the native-like diverging turn structure and a helical model. Free energies of these two conformations, estimated through MM-PBSA analysis using several force fields, suggest a comparable free energy (DeltaDeltaG < or =6 kcal/mol) for native and helix conformations. NOE intensities calculated from the native trajectory reproduce experimental data quite well, suggesting that the conformations sampled by the trajectory reasonably represent those observed in the NMR experiment. The molecular dynamics results, as well as sequence analysis of a diverse 690-member family of SH3 domain proteins, suggest that the presence of two elements is essential for formation of the diverging turn structure: a pair of residues with low helical propensity in the turn region and, as previously recognized, two hydrophobic residues to close the end of the diverging turn. Thus, these two sequence features may form the structural basis for the function of this peptide as a folding nucleus in this family of proteins.     

Folding study of an Aib-rich peptide in DMSO by molecular dynamics simulations.

To evaluate the ability of molecular dynamics (MD) simulations using atomic force-fields to correctly predict stable folded conformations of a peptide in solution, we show results from MD simulations of the reversible folding of an octapeptide rich in alpha-aminoisobutyric acid (2-amino-2-methyl-propanoic acid, Aib) solvated in di-methyl-sulfoxide (DMSO). This solvent generally prevents the formation of secondary structure, whereas Aib-rich peptides show a high propensity to form secondary structural elements, in particular 3(10)- and alpha-helical structures. Aib is, moreover, achiral, so that Aib-rich peptides can form left- or right-handed helices depending on the overall composition of the peptide, the temperature, and the solvation conditions. This makes the system an interesting case to study the ensembles of peptide conformations as a function of temperature by MD simulation. Simulations involving the folding and unfolding of the peptide were performed starting from two initial structures, a right-handed alpha-helical structure and an extended structure, at three temperatures, 298 K, 340 K, and 380 K, and the results are compared with experimental nuclear magnetic resonance (NMR) data measured at 298 K and 340 K. The simulations generally reproduce the available experimental nuclear Overhauser effect (NOE) data, even when a wide range of conformations is sampled at each temperature. The importance of adequate statistical sampling in order to reliably interpret the experimental data is discussed.     

A comparison of the dynamic behavior of monomeric and dimeric insulin shows structural rearrangements in the active monomer

Molecular dynamics (MD) simulations (5-10ns in length) and normal mode analyses were performed for the monomer and dimer of native porcine insulin in aqueous solution; both starting structures were obtained from an insulin hexamer. Several simulations were done to confirm that the results obtained are meaningful. The insulin dimer is very stable during the simulation and remains very close to the starting X-ray structure; the RMS fluctuations calculated from the MD simulation agree with the experimental B-factors. Correlated motions were found within each of the two monomers; they can be explained by persistent non-bonded interactions and disulfide bridges. The correlated motions between residues B24 and B26 of the two monomers are due to non-bonded interactions between the side-chains and backbone atoms. For the isolated monomer in solution, the A chain and the helix of the B chain are found to be stable during 5ns and 10ns MD simulations. However, the N-terminal and the C-terminal parts of the B chain are very flexible. The C-terminal part of the B chain moves away from the X-ray conformation after 0.5-2.5ns and exposes the N-terminal residues of the A chain that are thought to be important for the binding of insulin to its receptor. Our results thus support the hypothesis that, when monomeric insulin is released from the hexamer (or the dimer in our study), the C-terminal end of the monomer (residues B25-B30) is rearranged to allow binding to the insulin receptor. The greater flexibility of the C-terminal part of the beta chain in the B24 (Phe-->Gly) mutant is in accord with the NMR results. The details of the backbone and side-chain motions are presented. The transition between the starting conformation and the more dynamic structure of the monomers is characterized by displacements of the backbone of Phe B25 and Tyr B26; of these, Phe B25 has been implicated in insulin activation.     

'Wave-type' structure of a synthetic hexaglycosylated decapeptide: A part of the extracellular domain of human glycophorin A

The three-dimensional structure of a glycopeptide, His-Thr*-Ser*-Thr*-Ser*-Ser*-Ser*-Val-Thr-Lys, with 2-acetamido-2-deoxy--D-galactose (GalNAc) residues linked to six adjacent amino acids from Thr-10 to Ser-15, was studied by NMR spectroscopy and molecular dynamics (MD) simulations. The hexaglycosylated decapeptide is part of the extracellular domain of human glycophorin A and shows an extended structure of the peptide backbone due to O-glycosylation. Furthermore, each GalNAc residue exhibits one and only one NOE contact from the NHAc proton to the backbone amide proton of the amino acid that the sugar is directly bound to. This indicates a strong preference for the orientation of all GalNAc residues towards the N-terminus. NOE build-up curves were used to determine 42 inter-proton distances that, in connection with angles of the peptide backbone obtained from 3J-coupling constants, resulted in constraints for a MD simulation in water. The NMR data and the MD simulations show a preference for an extended backbone structure. The GalNAc residues are located alternatingly on opposite sides of the backbone and reduce the flexibility of the peptide backbone. The conformation of the molecule is relatively rigid and shows a 'wave-type' 3D structure of the peptide backbone within the glycosylation cluster. This new structural element is also supported by the unusual CD spectrum of the glycopeptide.     

Metropolis Monte Carlo calculations of DNA structure using internal coordinates and NMR distance restraints: An alternative method for generating a high-resolution solution structure

Summary A new method, a restrained Monte Carlo (rMC) calculation, is demonstrated for generating high-resolution structures of DNA oligonucleotides in solution from interproton distance restraints and bounds derived from complete relaxation matrix analysis of two-dimensional nuclear Overhauser effect (NOE) spectral peak intensities. As in the case of restrained molecular dynamics (rMD) refinement of structures, the experimental distance restraints and bounds are incorporated as a pseudo-energy term (or penalty function) into the mathematical expression for the molecular energy. However, the use of generalized helical parameters, rather than Cartesian coordinates, to define DNA conformation increases efficiency by decreasing by an order of magnitude the number of parameters needed to describe a conformation and by simplifying the potential energy profile. The Metropolis Monte Carlo method is employed to simulate an annealing process. The rMC method was applied to experimental 2D NOE data from the octamer duplex d(GTA-TAATG)·d(CATTATAC). Using starting structures from different locations in conformational space (e.g. A-DNA and B-DNA), the rMC calculations readily converged, with a root-mean-square deviation (RMSD) of <0.3 Å between structures generated using different protocols and starting structures. Theoretical 2D NOE peak intensities were calculated for the rMC-generated structures using the complete relaxation matrix program CORMA, enabling a comparison with experimental intensities via residual indices. Simulation of the vicinal proton coupling constants was carried out for the structures generated, enabling a comparison with the experimental deoxyribose ring coupling constants, which were not utilized in the structure determination in the case of the rMC simulations. Agreement with experimental 2D NOE and scalar coupling data was good in all cases. The rMC structures are quite similar to that refined by a traditional restrained MD approach (RMSD<0.5 Å) despite the different force fields used and despite the fact that MD refinement was conducted with additional restraints imposed on the endocyclic torsion angles of deoxyriboses. The computational time required for the rMC and rMD calculations is about the same. A comparison of structural parameters is made and some limitations of both methods are discussed with regard to the average nature of the experimental restraints used in the refinement.Abbreviations MC Monte Carlo - rMC restrained Monte Carlo - MD molecular dynamics - rMD restrained molecular dynamics - DG distance geometry - EM energy minimization - 2D NOE two-dimensional nuclear Overhauser effect - DQF-COSY double-quantum-filtered correlation spectroscopy - RMSD root-mean-square deviation To whom correspondence should be addressed.     

The structure of Escherichia coli heat-stable enterotoxin b by nuclear magnetic resonance and circular dichroism.

The heat-stable enterotoxin b (STb) is secreted by enterotoxigenic Escherichia coli that cause secretory diarrhea in animals and humans. It is a 48-amino acid peptide containing two disulfide bridges, between residues 10 and 48 and 21 and 36, which are crucial for its biological activity. Here, we report the solution structure of STb determined by two- and three-dimensional NMR methods. Approximate interproton distances derived from NOE data were used to construct structures of STb using distance-geometry and simulated annealing procedures. The NMR-derived structure shows that STb is helical between residues 10 and 22 and residues 38 and 44. The helical structure in the region 10-22 is amphipathic and exposes several polar residues to the solvent, some of which have been shown to be important in determining the toxicity of STb. The hydrophobic residues on the opposite face of this helix make contacts with the hydrophobic residues of the C-terminal helix. The loop region between residues 21 and 36 has another cluster of hydrophobic residues and exposes Arg 29 and Asp 30, which have been shown to be important for intestinal secretory activity. CD studies show that reduction of disulfide bridges results in a dramatic loss of structure, which correlates with loss of function. Reduced STb adopts a predominantly random-coil conformation. Chromatographic measurements of concentrations of native, fully reduced, and single-disulfide species in equilibrium mixtures of STb in redox buffers indicate that the formation of the two disulfide bonds in STb is only moderately cooperative. Similar measurements in the presence of 8 M urea suggest that the native secondary structure significantly stabilizes the disulfide bonds.     

Bending of the calmodulin central helix: a theoretical study.

The crystal structure of calcium-calmodulin (CaM) reveals a protein with a typical dumbbell structure. Various spectroscopic studies have suggested that the central linker region of CaM, which is alpha-helical in the crystal structure, is flexible in solution. In particular, NMR studies have indicated the presence of a flexible backbone between residues Lys 77 and Asp 80. This flexibility is related directly to the function of the protein because it enables the N- and C-terminal domains of the protein to move toward each other and bind to the CaM-binding domain of a target protein. We have investigated the flexibility of the CaM central helix by a variety of computational techniques: molecular dynamics (MD) simulations, normal mode analysis (NMA), and essential dynamics (ED) analysis. Our MD results reproduce the experimentally determined location of the bend in a simulation of only the CaM central helix, indicating that the bending point is an intrinsic property of the alpha-helix, for which the remainder of the protein is not important. Interestingly, the modes found by the ED analysis of the MD trajectory are very similar to the lowest frequency modes from the NM analysis and to modes found by an ED analysis of different structures in a set of NMR structures. Electrostatic interactions involving residues Arg 74 and Asp 80 seem to be important for these bending motions and unfolding, which is in line with pH-dependent NMR and CD studies.     

Structural similarities of micelle-bound peptide YY (PYY) and neuropeptide Y (NPY) are related to their affinity profiles at the Y receptors

Here, we investigate the structure of porcine peptide YY (pPYY) both when unligated in solution at pH 4.2 and when bound to dodecylphosphocholine (DPC) micelles at pH 5.5. pPYY in solution displays the PP-fold, with the N-terminal segment being back-folded onto the C-terminal alpha-helix, which extends from residue 17 to 31. In contrast to the solution structure of Keire et al. published in the year 2000 the C-terminal helix does not display a kink around residue 23-25. The root mean square deviation (RMSD) for backbone atoms of the NMR ensemble of conformers to the mean structure is 0.99(+/-0.35) Angstrom for residues 14-31. The back-fold is supported by values of 0.60+/-0.1 for the (15)N(1)H-NOE and by generalized order parameters S(2) of 0.74+/-0.1 for residues 5-31 which indicate that the peptide is folded in that segment. We have additionally used DPC micelles as a membrane model and determined the structure of pPYY when bound to it. Therein, an alpha-helix occurs in the segment comprising residues 17-31 and the N terminus freely diffuses in solution. The hydrophobic side of the amphipathic helix forms the micelle-binding interface and hydrophobic side-chains extend into the micelle interior. A significant stabilization of helical conformation occurs in the C-terminal pentapeptide, which is important for receptor binding. The latter is supported by positive values of the heteronuclear NOE in that segment (0.52+/-0.1 compared to 0.08+/-0.4 for the unligated form) and by values of S(2) of 0.6+/-0.2 (versus 0.38+/-0.2 for the unligated form). The structures of micelle-bound pPYY and pNPY are much more similar than those of pPYY and bPP with pairwise RMSDs of 1.23(+/-0.21)A or 3.21(+/-0.39) Angstrom, respectively. In contrast to the conformational similarities in the DPC-bound state their structures in solution are very different. In fact pPYY is more similar to bPP, which with its strong preference for the Y(4) receptor displays a completely different binding profile. Considering the high degree of sequence homology of pNPY and pPYY (>80%) and the fact, that their binding affinities at all receptor subtypes are high and, more importantly, rather similar, it is much more likely that PYY and NPY are recognized by the Y receptors from the membrane-bound state. As a consequence of the latter the PP-fold is not important for recognition of PYY or NPY at the Y receptors. To our knowledge this work provides for the first time strong arguments derived from structural data that support a membrane-bound receptor recognition pathway.     

Solution structure of monomeric peptide YY supports the functional significance of the PP-fold

Peptide YY (PYY) belongs to a family of peptides including neuropeptide Y (NPY) and pancreatic peptide (PP) that regulate numerous functions through both central and peripheral receptors. The solution structure of these peptides is hypothesized to be critically important in receptor selectivity and activation, based on prior demonstration of a stable tertiary conformation of PP called the "PP-fold". Circular dichroism (CD) spectra show a pH-dependent structural transition in the pH range 3-4. Thus we describe the tertiary structure of porcine PYY in water at pH 5.5, 25 degrees C, and 150 mM NaCl, as determined from 2D (1)H NMR data recorded at 500 MHz. A constraint set consisting of 396 interproton distances from NOE data was used as input for distance geometry, simulated annealing, and restrained energy minimization calculations in X-PLOR. The RMSDs of the 20 X-PLOR-generated structures were 0.71 +/- 0.14 and 1.16 +/- 0.17 A, respectively, for backbone and heavy atom overlays of residues 1-34. The resulting structure consists of two C-terminal helical segments from residues 17 to 22 and 25 to 33 separated by a kink at residues 23, 24, and 25, a turn centered around residues 12-14, and the N-terminus folded near residues 30 and 31. The well-defined portions of the PYY structure reported here bear a marked similarity to the structure of PP. Our findings strongly support the importance of the stable folded structure of this family of peptides for binding and activation of Y receptor subtypes.     

3D-Structure of the interior fusion peptide of HGV/GBV-C by 1H NMR, CD and molecular dynamics studies

In this work, we present a structural characterization of the putative fusion peptide E2(279-298) corresponding to the E2 envelope protein of the HGV/GBV-C virus by (1)H NMR, CD and MD studies performed in H(2)O/TFE and in lipid model membranes. The peptide is largely unstructured in water, whereas in H(2)O/TFE and in model membranes it adopts an helical structure (approximately 65-70%). The partitioning free energy DeltaG ranges from -6 to -7.5 kcal mol(-1). OCD measurements on peptide-containing hydrated and oriented lipid multilayers showed that the peptide adopts a predominantly surface orientation. The (1)H NMR data (observed NOEs, deuterium exchange rates, Halpha chemical shift index and vicinal coupling constants) and the molecular dynamics calculations support the conclusions that the peptide adopts a stable helix in the C-terminal 9-18 residues slightly inserted into the lipid bilayer and a major mobility in the amino terminus of the sequence (1-8 residues).     

Conformational study on glycosylated asparagine-oligopeptides by NMR spectroscopy and molecular dynamics calculations.

The conformational properties of the homo oligomers of increasing chain length Boc-(Asn)(n)-NHMe (n = 2, 4, 5), (GlcNAc-beta-Asn)(n)-NHMe (n = 2, 4, 5, 8) and Boc-[GlcNAc(Ac)(3)-beta-Asn](n)-NHMe (n = 2, 4, 5) were studied by using NOE experiments and molecular dynamic calculations (MD). Sequential NOEs and medium range NOEs, including (i,i+2) interactions, were detected by ROESY experiments and quantified. The calculated inter-proton distances are longer than those characteristic of beta-turn secondary structures. Owing to the large conformational motions expected for linear peptides, MD simulations were performed without NMR constraints, with explicit water and by applying different treatments of the electrostatic interactions. In agreement with the NOE results, the simulations showed, for all peptides, the presence of both folded and unfolded structures. The existence of significant populations of beta-turn structures can be excluded for all the examined compounds, but two families of structures were more often recognized. The first one with sinusoidal or S-shaped forms, and another family of large turns together with some more extended conformations. Only the glycosylated pentapeptide shows in vacuo a large amount of structures with helical shaped form. The results achieved in water and in DMSO are compared and discussed, together with the effect of the glycosylation.