Prediction of protein folding from amino acid sequence over discrete conformation spaces

Predicting the three-dimensional structure of a protein given only its amino acid sequence is a long-standing goal in computational chemistry. In the thermodynamic approach, one needs a potential function of conformation that resembles the free energy of the real protein to the extent that the global minimum of the potential is attained by the native conformation and no other. In practice, this has never been achieved with certainty because even with greatly simplified representations of the polypeptide chain, there are an astronomical number of local minima to examine. If one chooses instead a protein representation with only a large but manageable number of discrete conformations, then the global preference of the potential for the native can be directly verified. Representing a protein as a walk on a two-dimensional square lattice makes it easy to see that simple functions of the interresidue contacts are sufficient to globally favor a given "native" conformation, as long as it is a compact, globular structure. Explicit representation of the solvent is not required. Another more realistic way to confine the conformational search to a finite set is to draw alternative conformations from fragments of larger proteins having known crystal structure. Then it is possible to construct a simple function of interresidue contacts in three dimensions such that only 8 proteins are required to determine the adjustable parameters, and the native conformations of 37 other proteins are correctly preferred over all alternative conformations. The deduced function favors short-range backbone-backbone contacts regardless of residue type and long-range hydrophobic associations. Interactions over long distances, such as electrostatics, are not required.     

Prelude and Fugue, predicting local protein structure, early folding regions and structural weaknesses

Prelude&Fugue are bioinformatics tools aiming at predicting the local 3D structure of a protein from its amino acid sequence in terms of seven backbone torsion angle domains, using database-derived potentials. Prelude(&Fugue) computes all lowest free energy conformations of a protein or protein region, ranked by increasing energy, and possibly satisfying some interresidue distance constraints specified by the user. (Prelude&)Fugue detects sequence regions whose predicted structure is significantly preferred relative to other conformations in the absence of tertiary interactions. These programs can be used for predicting secondary structure, tertiary structure of short peptides, flickering early folding sequences and peptides that adopt a preferred conformation in solution. They can also be used for detecting structural weaknesses, i.e. sequence regions that are not optimal with respect to the tertiary fold. AVAILABILITY: http://babylone.ulb.ac.be/Prelude_and_Fugue.     

Nuclear Overhauser effects as probes of peptide structure. A diagnostic for left-handed helical conformations

The conformational dependence of the interresidue interproton distances in peptides, C alpha H ... Ni + 1 H and NiH ... Ni + 1 H, have been used to identify zones of sterically allowed phi, psi space, where both distances are less than 3A and expected to yield nuclear Overhauser effects (NOEs). L-residues in left-handed helical conformations are expected to yield both interresidue NOEs and also an appreciable intraresidue NiH----C alpha iH NOE. The effect of cutoff distances has been evaluated. Experimental results on three model peptides illustrate the utility of these NOEs in identifying L-residues at the i + 2 position of Type II and I' beta-turns. Simultaneous observation of both interresidue NOEs may also be indicative of conformational heterogeneity in specific cases, as illustrated for a single residue in a decapeptide.     

Structure of spin-labeled methylmethanethiolsulfonate in solution and bound to TEM-1 beta-lactamase determined by electron nuclear double resonance spectroscopy.

Site-directed spin-labeling of proteins whereby the spin-label methyl 3-(2,2,5,5-tetramethyl-1-oxypyrrolinyl)methanethiolsulfonate (SLMTS) is reacted with the -SH groups of cysteinyl residues incorporated into a protein by mutagenesis has been successfully applied to investigate secondary structure and conformational transitions of proteins. In these studies, it is expected that the spin-label moiety adopts different conformations dependent on its local environment. To determine the conformation of SLMTS in solution reacted with L-cysteine (SLMTCys) and bound in the active site of the Glu240Cys mutant of TEM-1 beta-lactamase, we have synthesized SLMTS both of natural abundance isotope composition and in site-specifically deuterated forms for electron nuclear double resonance (ENDOR) studies. ENDOR-determined electron-proton distances from the unpaired electron of the nitroxyl group of the spin-label to the methylene and methyl protons of SLMTS showed three conformations of the oxypyrrolinyl ring with respect to rotation around the S-S bond dependent on the solvent dielectric constant. For SLMTCys, two conformations of the molecule were compatible with the ENDOR-determined electron-nucleus distances to the side-chain methylene protons and to H(alpha) and H(beta1,2) of cysteine. To determine SLMTS conformation reacted with the Glu240Cys mutant of TEM-1 beta-lactamase, enzyme was overexpressed in both ordinary and perdeuterated minimal medium. Resonance features of H(alpha) and H(beta1,2) of the Cys240 residue of the mutant and of the side-chain methylene protons within the spin-label moiety yielded electron-proton distances that sterically accommodated the two conformations of free SLMTCys in solution.     

Conformations of oligonucleotides in solution as determined by sequence-specific antibodies

The conformations of some related oligoribonucleotides in solution were investigated immunochemically with antisera specific for two synthetic oligonucleotide sequences, A-A-U and A-A-U-U. Radioimmunoassay showed differences of as much as three or more orders of magnitude in binding among oligonucleotides with commonly shared sequences. These large differences, which reflect the loss of many points of contact with antibody because of changes in overall conformation, allow the following conclusions: (1) A-A-U and A-A-U-U have conformations distinct from any present in the Un family of oligomers. (2) Conformations of A-A-U and A-A-U-U differ markedly from those of oligomers of A. The dinucleotide A-A, in particular, bears little resemblance in conformation to the A-A sequence in A-A-U and A-A-U-U. (3) The recognizable conformational unit appears to be the triplet A-A-U, which binds as well as A-A-U-U and far better than its component dimers. Interactions between non-adjacent bases may be a factor here, as well as in codon recognition. The immunological data support the conclusion that, in oligonucleotides, as in polypeptides, primary sequence can determine conformation in solution.     

Contact potential that recognizes the correct folding of globular proteins.

We have devised a continuous function of interresidue contacts in globular proteins such that the X-ray crystal structure has a lower function value than that of thousands of protein-like alternative conformations. Although we fit the adjustable parameters of the potential using only 10,000 alternative structures for a selected training set of 37 proteins, a grand total of 530,000 constraints was satisfied, derived from 73 proteins and their numerous alternative conformations. In every case where the native conformation is adequately globular and compact, according to objective criteria we have developed, the potential function always favors the native over all alternatives by a substantial margin. This is true even for an additional three proteins never used in any way in the fitting procedure. Conformations differing only slightly from the native, such as those coming from crystal structures of the same protein complexed with different ligands or from crystal structures of point mutants, have function values very similar to the native's and always less than those of alternatives derived from substantially different crystal structures. This holds for all 95 structures that are homologous to one or another of various proteins we used. Realizing that this potential should be useful for modeling the conformation of new protein sequences from the body of protein crystal structures, we suggest a test for deciding whether a nearly correct approximation to the native conformation has been found.     

A geometrical constraint approach for reproducing the native backbone conformation of a protein

It is known that the backbone conformation of a protein can be reproduced with precision once a correct contact map (two-dimensional representation showing residue pairs in contact) is given as geometrical constraints. There is, however, no way to infer the correct contact map for a protein of unknown structure. We started with one-dimensional constraints using the quantity N14 (the number of neighboring residues within the radius of 14 Å). Since the plot of N14 along a chain shows a good correlation with the corresponding amino acid sequence, the N14 profile obtained from the X-ray structure is predictable from the sequence. Construction of backbone conformations under a given N14 profile was carried out in the following two steps: (1) a contact map from the N14 profile was produced by taking the product of N14 values of every two residues; (2) backbone conformations were generated by applying the distance geometry technique to distance constraints given by the contact map. If present, disulfide bonds in a protein, as well as the secondary structure, were treated as additional constraints, and both cases with or without the additional information were examined. The method was tested for 11 proteins of known structure, and the results indicated that the reproduced conformation was fairly good, using an X-ray structure for comparison, for small proteins of less than 80 residues long. The basic assumption and effectiveness of the present method were compared with those of previous studies employing the geometrical constraint approach. It has become clear that the specific, one-dimensional information (e.g., N14 profile) is more effective than nonspecific, two-dimensional constraints, such as average interresidue distances between particular types of amino acids. © 1993 Wiley-Liss, Inc.     

Abasic frameshift in DNA. Solution conformation determined by proton NMR and molecular mechanics calculations

We have determined the three-dimensional structure of a non-self-complementary oligodeoxynucleotide duplex that contains a model abasic site. The duplex contains six GC base pairs plus the abasic site at the center of one strand and corresponds to an abasic frameshift. Two-dimensional NMR studies on the nonexchangeable protons show that the guanine bases on either side of the abasic site are stacked over each other and that the abasic site is rotated out of the helix. Close proton-proton interactions are observed between the H4' proton of the abasic site and sugar protons of the guanosine in the 5' direction, which allows the position of the free sugar to be well-defined. NOE buildup curves from NOESY spectra recorded at very short mixing times were used to calculate a set of interproton distances. This data set was incorporated into the refinement of the oligonucleotide structure by molecular mechanics calculations. Two conformations that differ in the sugar conformation of the guanosine next to the abasic site in the 3' direction were necessary to fit all the NMR data. One of these two conformations could only be stabilized by addition of counterions at specific sites.     

Effects of distance constraints on macromolecular conformation. II. Simulation of experimental results and theoretical predictions

By generating classes of random structures for trypsin inhibitor and carp myogen, each consistent with a given set of experimental or theoretical information, we have assessed the relative utility of various experiments and theories in deducing the conformation of macromolecules. We compare the calculated structures with known x-ray coordinates and compute for each class an average error. Small errors mean that the experimental or theoretical constraints limit the structures to the vicinity of the crystal structure, whereas large errors show that the constraints permit a wide variety of tertiary conformations. We find the following points to hold true: (1) Qualitative information on all the distances, as might be obtained from the correct prediction of interresidue contacts, effectively determines the structure (error approximately 1 Å). (2) Quantitative information on a limited number of distances, as might be obtained from nmr or crosslinking experiments, significantly restricts the range of possible structures only when the number of distances given is comparable to the number of residues (error approximately 3 Å). (3) Quantitative information on the distances of each residue to the center of mass of the molecule, as might in part be obtained from solvent accessibility and solution x-ray studies, is not particularly restrictive by itself (error approximately 5 Å). (4) Complete qualitative local distance information, as might be obtained from secondary prediction and CD/ORD studies, is clearly consistent with a wide variety of tertiary structures (error approximately 7 Å).     

A distance geometry program for determining the structures of small proteins and other macromolecules from nuclear magnetic resonance measurements of intramolecular1H−1H proximities in solution

DISGEO is a new implementation of a distance geometry algorithm which has been specialized for the calculation of macromolecular conformation from distance measurements obtained by two-dimensional nuclear Overhauser enhancement spectroscopy. The improvements include (1) a decomposition of the complete embedding process into two successive, more tractable calculations by the use of “substructures”, (2) a compact data structure for storing incomplete distance information on a structure, (3) a more efficient shortest-path algorithm for computing the triangle inequality limits on all distances from this information, (4) a new algorithm for selecting random metric spaces from within these limits, (5) the use of chirality constraints to obtain good covalent geometry without the use ofad hoc weights or excessive optimization. The utility of the resultant program with nuclear magnetic resonance data is demonstrated by embedding complete spatial structures for the protein basic pancreatic trypsin inhibitor vs all 508 intramolecular, interresidue proton-proton contacts shorter than 4.0 Å that were present in its crystal structure. The crystal structure could be reproduced from this data set to within 1.3 Å minimum root mean square coordinate difference between the backbone atoms. We conclude that the information potentially available from nuclear magnetic resonance experiments in solution is sufficient to define the spatial structure of small proteins.     

An approach to the multiple-minima problem in protein folding by relaxing dimensionality. Tests on enkephalin

An algorithm for locating the region in conformational space containing the global energy minimum of a polypeptide is described. Distances are used as the primary variables in the minimization of an objective function that incorporates both energetic and distance-geometric terms. The latter are obtained from geometry and energy functions, rather than nuclear magnetic resonance experiments, although the algorithm can incorporate distances from nuclear magnetic resonance data if desired. The polypeptide is generated originally in a space of high dimensionality. This has two important consequences. First, all interatomic distances are initially at their energetically most favorable values; i.e. the polypeptide is initially at a global minimum-energy conformation, albeit a high-dimensional one. Second, the relaxation of dimensionality constraints in the early stages of the minimization removes many potential energy barriers that exist in three dimensions, thereby allowing a means of escaping from three-dimensional local minima. These features are used in an algorithm that produces short trajectories of three-dimensional minimum-energy conformations. A conformation in the trajectory is generated by allowing the previous conformation in the trajectory to evolve in a high-dimensional space before returning to three dimensions. The resulting three-dimensional structure is taken to be the next conformation in the trajectory, and the process is iterated. This sequence of conformations results in a limited but efficient sampling of conformational space. Results for test calculations on Met-enkephalin, a pentapeptide with the amino acid sequence H-Tyr-Gly-Gly-Phe-Met-OH, are presented. A tight cluster of conformations (in three-dimensional space) is found with ECEPP energies (Empirical Conformational Energy Program for Peptides) lower than any previously reported. This cluster of conformations defines a region in conformational space in which the global-minimum-energy conformation of enkephalin appears to lie.     

Single-molecule nanopositioning: structural transitions of a helicase-DNA complex during ATP hydrolysis

The conformational states of Escherichia coli Rep helicase undergoing ATP hydrolysis while bound to a partial-duplex DNA (pdDNA) were studied using single-molecule FRET. Crystallographic studies showed that Rep bound to single-stranded DNA can exist in open and closed conformations that differ in the orientation of the 2B subdomain. FRET measurements between eight Rep mutants donor-labeled at different residues and pdDNA acceptor-labeled at the junction were conducted at each of the four nucleotide states. The positions of donor-labeled residues, based on crystal structure, and FRET measurements between these donor molecules and the acceptor fluorophore at the DNA junction were used to predict the most likely position for the DNA junction using a triangulation algorithm. These predicted junction positions are compared with the crystal structure to determine whether the open or closed conformation is more consistent with the FRET data. Our data revealed that there are two distinct Rep-pdDNA conformations in the ATPγS and ADP states, an unexpected finding. The primary conformation is similar to that observed in nucleotide-free and ADP.Pi states, and the secondary conformation is a novel conformation where the duplex DNA and 2B subdomain moved as a unit by 13 Å relative to the rest of the protein. The primary conformation found in all nucleotide states is consistent with the closed conformation of the crystal structure however; the secondary conformation is a new conformation that has not been observed before. We discuss the possible implications of this newly observed conformation.     

Empirical solvation models in the context of conformational energy searches: application to bovine pancreatic trypsin inhibitor.

Continuum solvation models that estimate free energies of solvation as a function of solvent accessible surface area are computationally simple enough to be useful for predicting protein conformation. The behavior of three such solvation models has been examined by applying them to the minimization of the conformational energy of bovine pancreatic trypsin inhibitor. The models differ only with regard to how the constants of proportionality between free energy and surface area were derived. Each model was derived by fitting to experimentally measured equilibrium solution properties. For two models, the solution property was free energy of hydration. For the third, the property was NMR coupling constants. The purpose of this study is to determine the effect of applying these solvation models to the nonequilibrium conformations of a protein arising in the course of global searches for conformational energy minima. Two approaches were used: (1) local energy minimization of an ensemble of conformations similar to the equilibrium conformation and (2) global search trajectories using Monte Carlo plus minimization starting from a single conformation similar to the equilibrium conformation. For the two models derived from free energy measurements, it was found that both the global searches and local minimizations yielded conformations more similar to the X-ray crystallographic structures than did searches or local minimizations carried out in the absence of a solvation component of the conformational energy. The model derived from NMR coupling constants behaved similarly to the other models in the context of a global search trajectory. For one of the models derived from measured free energies of hydration, it was found that minimization of an ensemble of near-equilibrium conformations yielded a new ensemble in which the conformation most similar to the X-ray determined structure PTI4 had the lowest total free energy. Despite the simplicity of the continuum solvation models, the final conformation generated in the trajectories for each of the models exhibited some of the characteristics that have been reported for conformations obtained from molecular dynamics simulations in the presence of a bath of explicit water molecules. They have smaller root mean square (rms) deviations from the experimentally determined conformation, fewer incorrect hydrogen bonds, and slightly larger radii of gyration than do conformations derived from search trajectories carried out in the absence of solvent.     

Probing the folding capacity and residual structures in 1-79 residues fragment of staphylococcal nuclease by biophysical and NMR methods

1-79 residues SNase fragment (SNase79) has chain length containing a sequence for helix alpha(1), omega-loop, beta(I)-sheet, and partial beta(II)-sheet of native SNase. The incomplete "beta-barrel" structural region of SNase79 makes this fragment to be interested in investigation of its conformation. For this study, we use CD, fluorescence, and NMR spectroscopy to probe the folding capacity and the residual structures in SNase79. The optical spectra obtained for SNase79 and its mutants reveal the presence of retained capacity for folding of the fragment. The NMR derived (13)C(alpha) secondary chemical shifts, (3)J(NH-Halpha) coupling constants, amide-proton temperature coefficients, interresidue NOEs, and (15)N relaxation data determine the intrinsic propensities for helix- and turn- or beta-sheet-like conformations of SNase79, which is not the result of stabilizing inter-molecular interactions by oligomerization effects. The residual turn- and helix-like structures may serve as potential local nucleation sites, whereas the residual beta(I)-sheet-like structure can be regarded as a potential non-local nucleation site in the folding of SNase79. The intrinsic local and non-local interactions in these potential initiation sites are insufficient to stabilize the folding of SNase79 due to the shortage of relevant long-range interactions from other part of the fragment. The conformational ensemble of SNase79 is a highly heterogeneous collection of interconverting conformations having transiently populated helix- and beta-sheet- or turn-like structures.     

Side-chain entropy effects on protein secondary structure formation

Loss of conformational entropy is one of the primary factors opposing protein folding. Both the backbone and side-chain of each residue in a protein will have their freedom of motion restricted in the final folded structure. The type of secondary structure of which a residue is part will have a significant impact on how much side-chain entropy is lost. Side-chain conformational entropies have previously been determined for folded proteins, simple models of unfolded proteins, alpha-helices, and a dipeptide model for beta-strands, but not for polyproline II (PII) helices. In this work, we present side-chain conformational estimates for the three regular secondary structure types: alpha-helices, beta-strands, and PII helices. Entropies are estimated from Monte Carlo computer simulations. Beta-strands are modeled as two structures, parallel and antiparallel beta-strands. Our data indicate that restraining a residue to the PII helix or antiparallel beta-strand conformations results in side-chain entropies equal to or higher than those obtained by restraining residues to the parallel beta-strand conformation. Side-chains in the alpha-helix conformation have the lowest side-chain entropies. The observation that extended structures retain the most side-chain entropy suggests that such structures would be entropically favored in unfolded proteins under folding conditions. Our data indicate that the PII helix conformation would be somewhat favored over beta-strand conformations, with antiparallel beta-strand favored over parallel. Notably, our data imply that, under some circumstances, residues may gain side-chain entropy upon folding. Implications of our findings for protein folding and unfolded states are discussed.     

Effects of limited input distance constraints upon the distance geometry algorithm.

In this paper we examine the distance geometry (DG) algorithm in the form used to determine the structure of proteins. We focus on three aspects of the algorithm: bound smoothing with the triangle inequality, the random selection of distances within the bounds, and the number of distances needed to specify a structure. Computational experiments are performed using simulated and real data for basic pancreatic trypsin inhibitor (BPTI) from nmr and crystallographic measurements. We find that the upper bounds determined by bound smoothing to be a linear function of the true crystal distance. A simple model that describes the results obtained with randomly selected trial distances is proposed. Using this representation of the trial distances, we show that BPTI DG structures are more compact than the true crystal structure. We also show that the DG-generated structures no longer resemble test structures when the number of these interresidue distance constraints is less than the number of degrees of freedom of the protein backbone. While the actual model will be sensitive the way distances are chosen, our conclusions are likely to apply to other versions of the DG algorithm.     

Refinement of the spatial structure of the gramicidin A ion channel]

The spatial structure of the gramicidin A (GA) transmembrane ion-channel was refined on the base of cross-peak volumes measured in NOESY spectra (mixing time tau m = 100 and 200 ms). The refinement methods included the comparison of experimental cross-peak volumes with those calculated for low-energy GA conformations, dynamic averaging of the low-energy conformation set and restrained energy minimization. Accuracy of the spatial structure determination was estimated by the penalty function Fr defined as a root mean square deviation of interproton distances corresponding to the calculated and experimental cross-peak volumes. As the initial conformation we used the right-handed pi 6,3 LD pi 6,3 LD helix established on the base of NMR data regardless of the cross-peak volumes. The conformation is in a good agreement with NOE cross-peak volumes (Fr 0.2 to 0.5 A depending on NOESY spectrum). For a number of NOEs formed by the side chain protons, distances errors were found as much as 0.5-2.0 A. Restrained energy minimization procedure had little further success. However some of these errors were eliminated by the change in torsional angle chi 2 of D-Leu12 and dynamic averaging of the Val7 side chain conformations. Apparently, majority of deviations of the calculated and experimental cross-peak volumes are due to the intramolecular mobility of GA and cannot be eliminated within the framework of rigid globule model. In summary the spatial structure of GA ion-channel can be thought as a set of low-energy conformations, differing by the side chain torsion angles chi 1 Val7 and chi 2 D-Leu4 and D-Leu10 and the orientation of the C-terminal ethanolamine group. Root mean square differences between the atomic coordinates of conformations are in the range of 0.3-0.8 A.     

Detection of native-like models for amino acid sequences of unknown three-dimensional structure in a data base of known protein conformations.

We present an approach which can be used to identify native-like folds in a data base of protein conformations in the absence of any sequence homology to proteins in the data base. The method is based on a knowledge-based force field derived from a set of known protein conformations. A given sequence is mounted on all conformations in the data base and the associated energies are calculated. Using several conformations and sequences from the globin family we show that the native conformation is identified correctly. In fact the resolution of the force field is high enough to discriminate between a native fold and several closely related conformations. We then apply the procedure to several globins of known sequence but unknown three dimensional structure. The homology of these sequences to globins of known structures in the data base ranges from 49 to 17%. With one exception we find that for all globin sequences one of the known globin folds is identified as the most favorable conformation. These results are obtained using a force field derived from a data base devoid of globins of known structure. We briefly discuss useful applications in protein structural research and future development of our approach.     

Stabilization of a strained protein loop conformation through protein engineering.

Staphylococcal nuclease is found in two folded conformations that differ in the isomerization of the Lys 116-Pro 117 peptide bond, resulting in two different conformations of the residue 112-117 loop. The cis form is favored over the trans with an occupancy of 90%. Previous mutagenesis studies have shown that when Lys 116 is replaced by glycine, a trans conformation is stabilized relative to the cis conformation by the release of steric strain in the trans form. However, when Lys 116 is replaced with alanine, the resulting variant protein is identical to the wild-type protein in its structure and in the dominance of the cis configuration. The results of these studies suggested that any nuclease variant with a non-glycine residue at position 116 should also favor the cis form because of steric requirements of the beta-carbon at this position. In this report, we present a structural analysis of four nuclease variants with substitutions at position 116. Two variants, K116E and K116M, follow the "beta-carbon" hypothesis by favoring the cis form. Furthermore, the crystal structure of K116E is nearly identical to that of the wild-type protein. Two additional variants, K116D and K116N, provide exceptions to this simple "beta-carbon" rule in that the trans conformation is stabilized relative to the cis configuration by these substitutions. Crystallographic data indicate that this stabilization is effected through the addition of tertiary interactions between the side chain of position 116 with the surrounding protein and water structure. The detailed trans conformation of the K116D variant appears to be similar to the trans conformation observed in the K116G variant, suggesting that these two mutations stabilize the same conformation but through different mechanisms.     

Conformational and structural analysis of the equilibrium between single- and double-strand beta-helix of a D,L-alternating oligonorleucine

Alternating sequences of D and L residues in peptides are directly related to the formation of several kinds of regular helical conformations usually called beta-helices. The major feature of these structures is that they can be associated with the transmembrane ion-conducting channel activity in some natural antibacterial peptides. The study of alternating D,L synthetic peptides is critical to understand how factors such as surrounding media, main chain length, type of side chain and terminal groups, among others, can determine the adoption of a specific kind of beta-helix. Early studies pointed out that the peptides Boc-(D-NLeu-L-NLeu)(6)-D-MeNLe-L-Nl-D-Nl-L-Nl-OMe (Boc: tert-butyloxycarbonyl) and Boc-L-Nle-(D-Nle-L-Nle)(5)-D-MeNle-L-Nle-D-Nle-L-Nle-OMe adopt in chloroform a unique detectable conformation single beta(4.4)- and double beta(5.6) upward arrow downward arrow -helix, respectively. The influence of terminal groups on the final stable conformation of N-formylated peptides has been studied in this work. The initial basic NMR data analysis of a synthetic alternating D,L-oligopeptide with ten norleucines, N-methylated on the residue 7 and having HCO- and -OMe as terminal groups clearly indicates the coexistence of two different conformations in equilibrium. NMR data and molecular dynamics calculations point to a dimeric antiparallel beta-helix structure beta(5.6) upward arrow downward arrow for the main conformation. On the other hand, NMR data suggest a single beta-helix structure beta(4.4) for the second conformation. Finally, a thermodynamic analysis of the equilibrium between both conformations has been carried out by one-dimensional NMR measurements at ten different temperatures. The temperature at which 50% of dimer conformation is dissociated is 319 K. In addition, the dimer-monomer equilibrium curve obtained shows a DeltaG>0 for the whole range of studied temperatures, and its behavior can be considered similar to the thermodynamic denaturation protein processes.