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.     

UV resonance Raman investigation of a 3(10)-helical peptide reveals a rough energy landscape

We used UVRRS at 194 and 204 nm excitation to examine the backbone conformation of a 13-residue polypeptide (gp41(659-671)) that has been shown by NMR to predominantly fold into a 3(10)-helix. Examination of the conformation sensitive AmIII(3) region indicates the peptide has significant populations of beta-turn, PPII, 3(10)-helix, and pi-helix-like conformations but little alpha-helix. We estimate that at 1 degree C on average six of the 12 peptide bonds are in folded conformations (predominantly 3(10)- and pi-helix), while the other six are in unfolded (beta-turn/PPII) conformations. The folded and unfolded populations do not change significantly as the temperature is increased from 1 to 60 degrees C, suggesting a unique energy landscape where the folded and unfolded conformations are essentially degenerate in energy and exhibit identical temperature dependences.     

The role of alpha-, 3(10)-, and pi-helix in helix-->coil transitions

The conformational equilibrium between 3(10)- and alpha-helical structure has been studied via high-resolution NMR spectroscopy by Millhauser and coworkers using the MW peptide Ac-AMAAKAWAAKA AAARA-NH2. Their 750-MHz nuclear Overhauser effect spectroscopy (NOESY) spectra were interpreted to reflect appreciable populations of 3(10)-helix throughout the peptide, with the greatest contribution at the N and C termini. The presence of simultaneous alphaN(i,i + 2) and alphaN(i,i + 4) NOE cross-peaks was proposed to represent conformational averaging between 3(10)- and alpha-helical structures. In this study, we describe 25-nsec molecular dynamics simulations of the MW peptide at 298 K, using both an 8 A and a 10 A force-shifted nonbonded cutoff. The ensemble averages of both simulations are in reasonable agreement with the experimental helical content from circular dichroism (CD), the (3)J(HNalpha) coupling constants, and the 57 observed NOEs. Analysis of the structures from both simulations revealed very little formation of contiguous i --> i + 3 hydrogen bonds (3(10)-helix); however, there was a large population of bifurcated i --> i + 3 and i --> i + 4 alpha-helical hydrogen bonds. In addition, both simulations contained considerable populations of pi-helix (i --> i + 5 hydrogen bonds). Individual turns formed over residues 1-9, which we predict contribute to the intensities of the experimentally observed alphaN(i,i + 2) NOEs. Here we show how sampling of both folded and unfolded structures can provide a structural framework for deconvolution of the conformational contributions to experimental ensemble averages.     

A conformational analysis of leucine enkephalin as a function of pH

The conformations of Leu enkephalin in aqueous solution have been investigated as a function of pH using molecular dynamics simulations. The simulations suggest the peptide backbone exists as a mixture of folded and unfolded forms (approximately 50% each) at neutral pH, but is always unfolded at low or high pH. The folded form at neutral pH possesses a 2 --> 5 hydrogen bond and a close head to tail separation. No significant intramolecular hydrogen bonding of the carbonyl oxygens was observed in either the folded or unfolded forms of the peptide. Analysis of the Gly carbonyl oxygens and terminal groups indicated that, while the conformational population distribution of Leu enkephalin did vary noticeably as a function of pH, their hydration was essentially independent of pH and in agreement with the available NMR data. Further study indicated that the unfolded state of the peptide was not random in nature and consisted of one major unfolded backbone arrangement stabilized by a persistent hydrophobic interaction between the side chains of Tyr and Leu.     

Folding-unfolding thermodynamics of a beta-heptapeptide from equilibrium simulations

The thermodynamics of folding and unfolding of a beta-heptapeptide in methanol solution has been studied at four different temperatures, 298 K, 340 K, 350 K, and 360 K, by molecular dynamics simulation. At each of these temperatures, the 50-ns simulations were sufficient to generate an equilibrium distribution between a relatively small number of conformations (approximately 10(2)), showing that, even above the melting temperature (approximately 340 K), the peptide does not randomly sample conformational space. The free energy of folding and the free energy difference between pairs of conformations have been calculated from their relative populations. The experimentally determined folded conformation at 298 K, a left-handed 3(1)-helix, is at each of the four temperatures the predominant conformation, with its probability and average lifetime decreasing with increasing temperature. The most common intermediates of folding and unfolding are also the same at the four temperatures. Paths and rates of interconversion between different conformations have been determined. It has been found that folding can occur through multiple pathways, not necessarily downhill in free energy, although the final step involves a reduced number of intermediates.     

Angular variances for internal bond rotations of side chains in GXG-based tripeptides derived from (13)C-NMR relaxation measurements: Implications to protein folding

The study of backbone and side-chain internal motions in proteins and peptides is crucial to having a better understanding of protein/peptide "structure" and to characterizing unfolded and partially folded states of proteins and peptides. To achieve this, however, requires establishing a baseline for internal motions and motional restrictions for all residues in the fully, solvent-exposed "unfolded state." GXG-based tripeptides are the simpliest peptides where residue X is fully solvent exposed in the context of an actual peptide. In this study, a series of GXG-based tripeptides has been synthesized with X being varied to include all twenty common amino acid residues. Proton-coupled and -decoupled (13)C-nmr relaxation measurements have been performed on these twenty tripeptides and various motional models (Lipari-Szabo model free approach, rotational anisotropic diffusion, rotational fluctuations within a potential well, rotational jump model) have been used to analyze relaxation data for derivation of angular variances and motional correlation times for backbone and side-chain chi(1) and chi(2) bonds and methyl group rotations. At 298 K, backbone motional correlation times range from about 50 to 85 ps, whereas side-chain motional correlation times show a much broader spread from about 18 to 80 ps. Angular variances for backbone phi,psi bond rotations range from 11 degrees to 23 degrees and those for side chains vary from 5 degrees to 24 degrees for chi(1) bond rotations and from 5 degrees to 27 degrees for chi(2) bond rotations. Even in these peptide models of the "unfolded state," side-chain angular variances can be as restricted as those for backbone and beta-branched (valine, threonine, and isoleucine) and aromatic side chains display the most restricted motions probably due to steric hinderence with backbone atoms. Comparison with motional data on residues in partially folded, beta-sheet-forming peptides indicates that side-chain motions of at least hydrophobic residues are less restricted in the partially folded state, suggesting that an increase in side-chain conformational entropy may help drive early-stage protein folding. Copyright 1999 John Wiley & Sons, Inc.     

A shorter peptide model from staphylococcal nuclease for the folding-unfolding equilibrium of a beta-hairpin shows that unfolded state has significant contribution from compact conformational states

It is important to understand the conformational features of the unfolded state in equilibrium with folded state under physiological conditions. In this paper, we consider a short peptide model LMYKGQPM from staphylococcal nuclease to model the conformational equilibrium between a hairpin conformation and its unfolded state using molecular dynamics simulation under NVT conditions at 300K using GROMOS96 force field. The free energy landscape has overall funnel-like shape with hairpin conformations sampling the minima. The "unfolded" state has a higher free energy of approximately 12kJ/mol with respect to native hairpin minimum and occupies a plateau region. We find that the unfolded state has significant contributions from compact conformations. Many of these conformations have hairpin-like topology. Further, these compact conformational forms are stabilized by hydrophobic interactions. Conversion between native and non-native hairpins occurs via unfolded states. Frequent conversions between folded and unfolded hairpins are observed with single exponential kinetics. We compare our results with the emerging picture of unfolded state from both experimental and theoretical studies.     

Basis of substrate binding by the chaperonin GroEL.

The molecular chaperonins are essential proteins involved in protein folding, complex assembly, and polypeptide translocation. While there is abundant structural information about the machinery and the mechanistic details of its action are well studied, it is yet unresolved how chaperonins recognize a large number of structurally unrelated polypeptides in their unfolded or partially folded forms. To determine the nature of chaperonin-substrate recognition, we have characterized by NMR methods the interactions of GroEL with synthetic peptides that mimic segments of unfolded proteins. In previous work, we found using transferred nuclear Overhauser effect (trNOE) analysis that two 13 amino acid peptides bound GroEL in an amphipathic alpha-helical conformation. By extending the study to a variety of peptides with differing sequence motifs, we have observed that peptides can adopt conformations other than alpha-helix when bound to GroEL. Furthermore, peptides of the same composition exhibited significantly different affinities for GroEL as manifested by the magnitude of trNOEs. Binding to GroEL correlates well with the ability of the peptide to cluster hydrophobic residues on one face of the peptide, as determined by the retention time on reversed-phase (RP) HPLC. We conclude that the molecular basis of GroEL-substrate recognition is the presentation of a hydrophobic surface by an incompletely folded polypeptide and that many backbone conformations can be accommodated.     

Application of NMR,molecular simulation,and hydrodynamics to conformational analysis of trisaccharides

The preferred conformations and conformational flexibilities of the trisaccharides alpha-D-Glcp-(1-->2)-beta-D-Glcp-(1-->3)-alpha-D-Glcp-OMe (I) and alpha-D-Glcp-(1-->3)[beta-D-Glcp-(1-->4)]-alpha-D-Glcp-OMe (II) in aqueous solution were determined using nuclear magnetic resonance (NMR) spectroscopy, molecular dynamics (MD) and Langevin dynamics (LD) simulations, and hydrodynamics calculations. Both trisaccharides have a vicinal substitution pattern in which long range (nonsequential) interactions may play an important role. LD simulation at 600 K indicated that the all-syn conformation predominated, though other conformations were apparent. NOE data and MD and LD simulations at 298 K all indicated that trisaccharide I is a single all-syn conformer in solution. Given that previous studies showed evidence of anti-conformers in beta-D-Glcp-(1-->2)-beta-D-Glcp-(1-->3)-alpha-D-Glcp-OMe, this result provides an example of how changing the anomeric configuration of one residue from beta to alpha can make an oligosaccharide more rigid. Discrepancies in inter-ring distances obtained by experiment and by simulation of the all-syn conformer suggest the presence of an anti-psi conformation at the beta-(1-->4)-linkage for II. A combined analysis of measured and calculated translational diffusion constants and (13)C T(1) relaxation times yield order parameters of 0.9 for each trisaccharide. This implies that any interconversion among conformations is significantly slower than tumbling. Anisotropies of approximately 1.6 and 1.3 calculated for I and II, respectively, are consistent with the observed relatively flat T(1) profiles because the tumbling is not in the motional narrowing regime.     

Folding and stability of the three-stranded beta-sheet peptide Betanova: insights from molecular dynamics simulations

The dynamics of the three-stranded beta-sheet peptide Betanova has been studied at four different temperatures (280, 300, 350, and 450 K by molecular dynamics simulation techniques, in explicit water. Two 20-ns simulations at 280 K indicate that the peptide remains very flexible under "folding" conditions sampling a range of conformations that together satisfy the nuclear magnetic resonance (NMR)-derived experimental constraints. Two simulations at 300 K (above the experimental folding temperature) of 20 ns each show partial formation of "native"-like structure, which also satisfies most of the NOE constraints at 280 K. At higher temperature, the presence of compact states, in which a series of hydrophobic contacts remain present, are observed. This is consistent with experimental observations regarding the role of hydrophobic contacts in determining the peptide's stability and in initiating the formation of turns and loops. A set of different structures is shown to satisfy NMR-derived distance restraints and a possible mechanism for the folding of the peptide into the NMR-determined structure is proposed.     

Membrane adsorption,folding, insertion and translocation of synthetic trans-membrane peptides

Spontaneous membrane adsorption, folding and insertion of the synthetic WALP16 and KALP16 peptides was studied by computer simulations starting from completely extended conformations. The peptides were simulated using an unmodified all-atom force field in combination with an efficient Monte Carlo sampling algorithm. The membrane is represented implicitly as a hydrophobic zone inside a continuum solvent modelled using the generalized Born theory of solvation. The method was previously parameterized to match insertion energies of hydrophobic side chain analogs into cyclohexane and no parameters were optimized for the present simulations. Both peptides rapidly precipitate out of bulk solution and adsorb to the membrane surface. Interfacial folding into a helical conformation is followed by membrane insertion. Both the peptide conformations and their location in the membrane are strongly temperature dependent. The temperature dependent behaviour can be summarized by fitting to a four-state model, separating the system into folded and unfolded conformers, which are either inserted into the membrane or located at the interfaces. As the temperature is lowered the dominant peptide conformation of the system changes from unfolded surface bound configurations to folded surface bound states. Folded trans-membrane conformers represent the dominant configuration at low temperatures. The analysis allows direct estimates of the free energies of peptide folding and membrane insertion. In the case of WALP the quality of the fit is excellent and the thermodynamic behaviour is in good agreement with expected theoretical consideration. For KALP the fit is more problematic due to the large solvation energies of the charged lysine residues.     

Membrane adsorption, folding, insertion and translocation of synthetic trans-membrane peptides

Spontaneous membrane adsorption, folding and insertion of the synthetic WALP16 and KALP16 peptides was studied by computer simulations starting from completely extended conformations. The peptides were simulated using an unmodified all-atom force field in combination with an efficient Monte Carlo sampling algorithm. The membrane is represented implicitly as a hydrophobic zone inside a continuum solvent modelled using the generalized Born theory of solvation. The method was previously parameterized to match insertion energies of hydrophobic side chain analogs into cyclohexane and no parameters were optimized for the present simulations. Both peptides rapidly precipitate out of bulk solution and adsorb to the membrane surface. Interfacial folding into a helical conformation is followed by membrane insertion. Both the peptide conformations and their location in the membrane are strongly temperature dependent. The temperature dependent behaviour can be summarized by fitting to a four-state model, separating the system into folded and unfolded conformers, which are either inserted into the membrane or located at the interfaces. As the temperature is lowered the dominant peptide conformation of the system changes from unfolded surface bound configurations to folded surface bound states. Folded trans-membrane conformers represent the dominant configuration at low temperatures. The analysis allows direct estimates of the free energies of peptide folding and membrane insertion. In the case of WALP the quality of the fit is excellent and the thermodynamic behaviour is in good agreement with expected theoretical consideration. For KALP the fit is more problematic due to the large solvation energies of the charged lysine residues.     

Transient intermediary states with high and low folding probabilities in the apparent two-state folding equilibrium of ACBP at low pH

Measurements of the stability as a function of pH for the acyl-coenzyme A binding protein (ACBP) has shown a significant difference in the pH transition midpoint measured by NMR spectroscopy at pH 3.12 and the transition midpoint measured at pH 2.92 and 2.97 by circular dichroism and by fluorescence spectroscopy, respectively. A similar behavior has not been observed in other proteins. It is suggested that these differences arise because the population of the unfolded molecules still contains significant amounts of native like secondary and tertiary structure. NMR spectroscopy measures the concentration of the two components of the folding unfolding equilibrium individually, whereas circular dichroism and fluorescence measure the concentration of the conformations of the light-absorbing chromophores present in both the folded and the unfolded molecules. In the narrow pH range, nascent structure can be detected as the average amount of secondary structure per unfolded molecule and hydrophobic interactions in the population of unfolded molecules. These structures are not observable immediately by NMR spectroscopy; however, a chemical shift analysis of the peptide backbone (13)C chemical shift indicates strongly the existence of short-lived and transient helical structures at pH 2.3. Magnetization transfer studies have been applied to study the equilibrium between folded and unfolded ACBP near the pH transition point measured by NMR. This study has shown that there are two categories of subpopulations in the population of unfolded ACBP. One for which magnetization can be transferred to the folded form during the folding process, and one for which transfer is not observed. The molecules of the latter population of unfolded protein apparently, do not fold within the time-frame of the magnetization transfer experiment. This result suggests the existence of a subpopulation of the acid-unfolded protein molecules with a high propensity for folding. It is suggested that in this subpopulation, a particular set of native like interactions in the peptide backbone and between side-chains in the peptide chain have to be formed.     

Molecular dynamics simulations of beta-hairpin folding

Molecular dynamics simulations of beta-hairpin folding have been carried out with a solvent-referenced potential at 274 K. The model peptide V4DPGV4 formed stable beta-hairpin conformations and the beta-hairpin ratio calculated by the DSSP algorithm was about 56% in the 50-ns simulation. Folding into beta-hairpin conformations is independent of the initial conformations. The simulations provided insights into the folding mechanism. The hydrogen bond often formed in a beta-turn first, and then propagated by forming more hydrogen bonds along the strands. Unfolding and refolding occurred repeatedly during the simulations. Both the hydrogen bonding and the hydrophobic interaction played important roles in forming the ordered structure. Without the hydrophobic effect, stable beta-hairpin conformations did not form in the simulations. With the same energy functions, the alanine-based peptide (AAQAA)3Y folded into helical conformations, in agreement with experiments. Folding into an alpha-helix or a beta-hairpin is amino acid sequence-dependent.     

A small tripeptide AFA undergoes two state cooperative conformational transitions: implications for conformational biases in unfolded states

It is important to understand the conformational biases that are present in unfolded states to understand protein folding. In this context, it is surprising that even a short tripeptide like AFA samples folded/ordered conformation as demonstrated recently by NMR experiments of the peptide in aqueous solution at 280 K. In this paper, we present molecular dynamics simulation of the peptide in explicit water using OPLS-AA/L all-atom force field. The results are in overall agreement with NMR results and provide some further insights. The peptide samples turn and extended conformational forms corresponding to minima in free energy landscape. Frequent transitions between the minima are observed due to modest free energy barriers. The turn conformation seems to be stabilized by hydrophobic interactions and possibly by bridging water molecules between backbone donors and acceptors. Thus the peptide does not sample conformations randomly, but samples well defined conformations. The peptide served as a model for folding-unfolding equilibrium in the context of peptide folding. Further, implications for drug design are also discussed.     

Functional and structural characterization of a synthetic peptide representing the N-terminal domain of prokaryotic pyruvate dehydrogenase

A synthetic peptide (Nterm-E1p) is used to characterize the structure and function of the N-terminal region (amino acid residues 4-45) of the pyruvate dehydrogenase component (E1p) from the pyruvate dehydrogenase multienzyme complex (PDHC) from Azotobacter vinelandii. Activity and binding studies established that Nterm-E1p specifically competes with E1p for binding to the dihydrolipoyl transacetylase component (E2p) of PDHC. Moreover, the experiments show that the N-terminal region of E1p forms an independent folding domain that functions as a binding domain. CD measurements, two-dimensional (2D) (1)H NMR analysis, and secondary structure prediction all indicate that Nterm-E1p has a high alpha-helical content. Here a structural model of the N-terminal domain is proposed. The peptide is present in two conformations, the population of which depends on the sample conditions. The conformations are designated "unfolded" at pH > or =6 and "folded" at pH <5. The 2D (1)H TOCSY spectrum of a mixture of folded and unfolded Nterm-E1p shows exchange cross-peaks that "link" the folded and unfolded state of Nterm-E1p. The rate of exchange between the two species is in the range of 0.5-5 s(-1). Sharp resonances in the NMR spectra of wild-type E1p demonstrate that this 200 kDa enzyme contains highly flexible regions. The observed dynamic character of E1p and of Nterm-E1p is likely required for the binding of the E1p dimer to the two different binding sites on E2p. Moreover, the flexibility might be essential in sustaining the allosteric properties of the enzyme bound in the complex.     

Thermal unfolding in a GCN4-like leucine zipper: 13C alpha NMR chemical shifts and local unfolding curves.

13C alpha chemical shifts and site-specific unfolding curves are reported for 12 sites on a 33-residue, GCN4-like leucine zipper peptide (GCN4-lzK), ranging over most of the chain and sampling most heptad positions. Data were derived from NMR spectra of nine synthetic, isosequential peptides bearing 99% 13C alpha at sites selected to avoid spectral overlap in each peptide. At each site, separate resonances appear for unfolded and folded forms, and most sites show resonances for two folded forms near room temperature. The observed chemical shifts suggest that 1) urea-unfolded GCN4-lzK chains are randomly coiled; 2) thermally unfolded chains include significant transient structure, except at the ends; 3) the coiled-coli structure in the folded chains is atypical near the C-terminus; 4) only those interior sites surrounded by canonical interchain salt bridges fail to show two folded forms. Local unfolding curves, obtained from integrated resonance intensities, show that 1) sites differ in structure content and in melting temperature, so the equilibrium population must comprise more than two molecular conformations; 2) there is significant end-fraying, even at the lowest temperatures, but thermal unfolding is not a progressive unwinding from the ends; 3) residues 9-16 are in the lowest melting region; 4) heptad position does not dictate stability; 5) significant unfolding occurs below room temperature, so the shallow, linear decline in backbone CD seen there has conformational significance. It seems that only a relatively complex array of conformational states could underlie these findings.     

The Alzheimer's peptides Abeta40 and 42 adopt distinct conformations in water: a combined MD / NMR study

The role of peptides Abeta40 and Abeta42 in the early pathogenesis of Alzheimer's disease (AD) is frequently emphasized in the literature. It is known that Abeta42 is more prone to aggregation than Abeta40, even though they differ in only two (IA) amino acid residues at the C-terminal end. A direct comparison of the ensembles of conformations adopted by the monomers in solution has been limited by the inherent flexibility of the unfolded peptides. Here, we characterize the conformations of Abeta40 and Abeta42 in water by using a combination of molecular dynamics (MD) and measured scalar (3)J(HNHalpha) data from NMR experiments. We perform replica exchange MD (REMD) simulations and find that classical forcefields reproduce the NMR data quantitatively when the sampling is extended to the microseconds time-scale. Using the quantitative agreement of the NMR data as a validation of the model, we proceed to compare the conformational ensembles of the Abeta40 and Abeta42 peptide monomers. Our analysis confirms the existence of structured regions within the otherwise flexible Abeta peptides. We find that the C terminus of Abeta42 is more structured than that of Abeta40. The formation of a beta-hairpin in the sequence (31)IIGLMVGGVVIA involving short strands at residues 31-34 and 38-41 (in bold) reduces the C-terminal flexibility of the Abeta42 peptide and may be responsible for the higher propensity of this peptide to form amyloids.     

Beta-hairpin folding mechanism of a nine-residue peptide revealed from molecular dynamics simulations in explicit water

The beta-hairpin fold mechanism of a nine-residue peptide, which is modified from the beta-hairpin of alpha-amylase inhibitor tendamistat (residues 15-23), is studied through direct folding simulations in explicit water at native folding conditions. Three 300-nanosecond self-guided molecular dynamics (SGMD) simulations have revealed a series of beta-hairpin folding events. During these simulations, the peptide folds repeatedly into a major cluster of beta-hairpin structures, which agree well with nuclear magnetic resonance experimental observations. This major cluster is found to have the minimum conformational free energy among all sampled conformations. This peptide also folds into many other beta-hairpin structures, which represent some local free energy minimum states. In the unfolded state, the N-terminal residues of the peptide, Tyr-1, Gln-2, and Asn-3, have a confined conformational distribution. This confinement makes beta-hairpin the only energetically favored structure to fold. The unfolded state of this peptide is populated with conformations with non-native intrapeptide interactions. This peptide goes through fully hydrated conformations to eliminate non-native interactions before folding into a beta-hairpin. The folding of a beta-hairpin starts with side-chain interactions, which bring two strands together to form interstrand hydrogen bonds. The unfolding of the beta-hairpin is not simply the reverse of the folding process. Comparing unfolding simulations using MD and SGMD methods demonstrate that SGMD simulations can qualitatively reproduce the kinetics of the peptide system.     

Conformational characteristics of unstructured peptides: alpha-synuclein

We have performed replica-exchange molecular dynamics simulations on 41 residue peptides containing NAC region of alpha-synuclein in various force fields and solvent conditions. Alpha-synuclein is known to be the major cause of Parkinson's disease by amyloid-like aggregation, and one of the natively unfolded proteins. To investigate conformational characteristics of intrinsically unstructured peptides, we carried out structural analysis by introducing 'representative structure' for ensemble of structures occurring during the overall trajectory. Representative structures may be defined by using either coordinate averaging or distance averaging. When applied to the natively folded proteins such as villin headpiece, structural analysis based on representative structure was found to yield consistent results with those obtained from conventional analysis. Individual conformations obtained from the simulations of NAC peptide for various conditions show flexible structures close to random coil. Secondary structure contents and free energy surfaces showed dependency on solvent conditions, which may be interpreted as another manifestation of structural diversity. It is found that representative structures can provide useful information about structural characteristics of intrinsically unstructured proteins.