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.     

Exploring the energy landscape of protein folding using replica-exchange and conventional molecular dynamics simulations

Two independent replica-exchange molecular dynamics (REMD) simulations with an explicit water model were performed of the Trp-cage mini-protein. In the first REMD simulation, the replicas started from the native conformation, while in the second they started from a nonnative conformation. Initially, the first simulation yielded results qualitatively similar to those of two previously published REMD simulations: the protein appeared to be over-stabilized, with the predicted melting temperature 50-150K higher than the experimental value of 315K. However, as the first REMD simulation progressed, the protein unfolded at all temperatures. In our second REMD simulation, which starts from a nonnative conformation, there was no evidence of significant folding. Transitions from the unfolded to the folded state did not occur on the timescale of these simulations, despite the expected improvement in sampling of REMD over conventional molecular dynamics (MD) simulations. The combined 1.42 micros of simulation time was insufficient for REMD simulations with different starting structures to converge. Conventional MD simulations at a range of temperatures were also performed. In contrast to REMD, the conventional MD simulations provide an estimate of Tm in good agreement with experiment. Furthermore, the conventional MD is a fraction of the cost of REMD and continuous, realistic pathways of the unfolding process at atomic resolution are obtained.     

Reproducible polypeptide folding and structure prediction using molecular dynamics simulations

The folding of a polypeptide from an extended state to a well-defined conformation is studied using microsecond classical molecular dynamics (MD) simulations and replica exchange molecular dynamics (REMD) simulations in explicit solvent and in vacuo. It is shown that the solvated peptide folds many times in the REMD simulations but only a few times in the conventional simulations. From the folding events in the classical simulations we estimate an approximate folding time of 1-2 micros. The REMD simulations allow enough sampling to deduce a detailed Gibbs free energy landscape in three dimensions. The global minimum of the energy landscape corresponds to the native state of the peptide as determined previously by nuclear magnetic resonance (NMR) experiments. Starting from an extended state it takes about 50 ns before the native structure appears in the REMD simulations, about an order of magnitude faster than conventional MD. The calculated melting curve is in good qualitative agreement with experiment. In vacuo, the peptide collapses rapidly to a conformation that is substantially different from the native state in solvent.     

Computer simulations of protein folding by targeted molecular dynamics

We have performed 128 folding and 45 unfolding molecular dynamics runs of chymotrypsin inhibitor 2 (CI2) with an implicit solvation model for a total simulation time of 0.4 microseconds. Folding requires that the three-dimensional structure of the native state is known. It was simulated at 300 K by supplementing the force field with a harmonic restraint which acts on the root-mean-square deviation and allows to decrease the distance to the target conformation. High temperature and/or the harmonic restraint were used to induce unfolding. Of the 62 folding simulations started from random conformations, 31 reached the native structure, while the success rate was 83% for the 66 trajectories which began from conformations unfolded by high-temperature dynamics. A funnel-like energy landscape is observed for unfolding at 475 K, while the unfolding runs at 300 K and 375 K as well as most of the folding trajectories have an almost flat energy landscape for conformations with less than about 50% of native contacts formed. The sequence of events, i.e., secondary and tertiary structure formation, is similar in all folding and unfolding simulations, despite the diversity of the pathways. Previous unfolding simulations of CI2 performed with different force fields showed a similar sequence of events. These results suggest that the topology of the native state plays an important role in the folding process.     

Microscopic reversibility of protein folding in molecular dynamics simulations of the engrailed homeodomain

The principle of microscopic reversibility states that at equilibrium the number of molecules entering a state by a given path must equal those exiting the state via the same path under identical conditions or, in structural terms, that the conformations along the two pathways are the same. There has been some indirect evidence indicating that protein folding is such a process, but there have been few conclusive findings. In this study, we performed molecular dynamics simulations of an ultrafast unfolding and folding protein at its melting temperature to observe, on an atom-by-atom basis, the pathways the protein followed as it unfolded and folded within a continuous trajectory. In a total of 0.67 micros of simulation in water, we found six transient denaturing events near the melting temperature (323 and 330 K) and an additional refolding event following a previously identified unfolding event at a high temperature (373 K). In each case, unfolding and refolding transition state ensembles were identified, and they agreed well with experiment on the basis of a comparison of S and Phi values. On the basis of several structural properties, these 13 transition state ensembles agreed very well with each other and with four previously identified transition states from high-temperature denaturing simulations. Thus, not only were the unfolding and refolding transition states part of the same ensemble, but in five of the seven cases, the pathway the protein took as it unfolded was nearly identical to the subsequent refolding pathway. These events provide compelling evidence that protein folding is a microscopically reversible process. In the other two cases, the folding and unfolding transition states were remarkably similar to each other but the paths deviated.     

Phi-analysis at the experimental limits: mechanism of beta-hairpin formation

The 37-residue Formin-binding protein, FBP28, is a canonical three-stranded beta-sheet WW domain. Because of its small size, it is so insensitive to chemical denaturation that it is barely possible to determine accurately a denaturation curve, as the transition spans 0-7 M guanidinium hydrochloride (GdmCl). It is also only marginally stable, with a free energy of denaturation of just 2.3 kcal/mol at 10 degrees Celsius so only small changes in energy upon mutation can be tolerated. But these properties and relaxation times for folding of 25 micros-400 micros conspire to allow the rapid acquisition of accurate and reproducible kinetic data for Phi-analysis using classical temperature-jump methods. The transition state for folding is highly polarized with some regions having Phi-values of 0 and others 1, as readily seen in chevron plots, with Phi-values of 0 having the refolding arms overlaying and those of 1 the unfolding arms superimposable. Good agreement is seen with transition state structures identified from independent molecular dynamics (MD) simulations at 60, 75, and 100 degrees Celsius, which allows us to explore further the details of the folding and unfolding pathway of FBP28. The first beta-turn is near native-like in the transition state for folding (experimental) and unfolding (MD and experiment). The simulations show that there are transient contacts between the aromatic side-chains of the beta-strands in the denatured state and that these interactions provide the driving force for folding of the first beta-hairpin of this three-stranded sheet. Only after the backbone hydrogen bonds are formed between beta1 and beta2 does a hydrogen bond form to stabilize the intervening turn, or the first beta-turn.     

Equilibrium structure and folding of a helix-forming peptide: circular dichroism measurements and replica-exchange molecular dynamics simulations

We have performed experimental measurements and computer simulations of the equilibrium structure and folding of a 21-residue alpha-helical heteropeptide. Far ultraviolet circular dichroism spectroscopy is used to identify the presence of helical structure and to measure the thermal unfolding curve. The observed melting temperature is 296 K, with a folding enthalpy of -11.6 kcal/mol and entropy of -39.6 cal/(mol K). Our simulations involve 45 ns of replica-exchange molecular dynamics of the peptide, using eight replicas at temperatures between 280 and 450 K, and the program CHARMM with a continuum solvent model. In a 30-ns simulation started from a helical structure, conformational equilibrium at all temperatures was reached after 15 ns. This simulation was used to calculate the peptide melting curve, predicting a folding transition with a melting temperature in the 330-350 K range, enthalpy change of -10 kcal/mol, and entropy change of -30 cal/(mol K). The simulation results were also used to analyze the peptide structural fluctuations and the free-energy surface of helix unfolding. In a separate 15-ns replica-exchange molecular dynamics simulation started from the extended structure, the helical conformation was first attained after approximately 2.8 ns, and equilibrium was reached after 10 ns of simulation. These results showed a sequential folding process with a systematic increase in the number of hydrogen bonds until the helical state is reached, and confirmed that the alpha-helical state is the global free-energy minimum for the peptide at low temperatures.     

Tuning lambda6-85 towards downhill folding at its melting temperature

The five-helix bundle lambda6-85* is a fast two-state folder. Several stabilized mutants have been reported to fold kinetically near-downhill or downhill. These mutants undergo a transition to two-state folding kinetics when heated. It has been suggested that this transition is caused by increased hydrophobicity at higher temperature. Here we investigate two histidine-containing mutants of lambda6-85* to see if a weaker hydrophobic core can extend the temperature range of downhill folding. The very stable lambdaHA is the fastest-folding lambda repressor to date (k(f)(-1) approximately k(obs)(-1)=2.3 micros at 44 degrees C). It folds downhill at low temperature, but transits back to two-state folding at its unfolding midpoint. lambdaHG has a weakened hydrophobic core. It is less stable than some slower folding mutants of lambda6-85*, and it has more exposed hydrophobic surface area in the folded state. This mutant nonetheless folds very rapidly, and has the non-exponential folding kinetics of an incipient downhill folder even at the unfolding midpoint (k(m)(-1) approximately 2 micros, k(a)(-1)=15 micros at 56 degrees C). We also compare the thermodynamic melting transition of lambdaHG with the nominal two-state folding mutant lambdaQG, which has a similar melting temperature. Unlike lambdaQG, lambdaHG yields fluorescence wavelength-dependent cooperativities and probe-dependent melting temperatures. This result combined with previous work shows that the energy landscapes of lambda repressor mutants support all standard folding mechanisms.     

Temperature-dependent downhill unfolding of ubiquitin. I. Nanosecond-to-millisecond resolved nonlinear infrared spectroscopy

Transient thermal unfolding of ubiquitin is investigated using nonlinear infrared spectroscopy after a nanosecond laser temperature jump (T-jump). The abrupt change in the unfolding free energy surface and the ns time resolution allow us to observe a fast response on ns to micros time-scales, which we attribute to downhill unfolding, before a cross-over to ms kinetics. The downhill unfolding by a sub-population of folded proteins is induced through a shift of the barrier toward the native state. By adjusting the T-jump width, the effect of the initial (T(i)) and final (T(f)) temperature on the unfolding dynamics can be separated. From the amplitude of the fast downhill unfolding, the fractional population prepared at the unfolding transition state is obtained. This population increases with both T(i) and with T(f). A two-state kinetic analysis of the ms refolding provides thermodynamic information about the barrier height. By a combination of the fast and slow unfolding and folding parameters, a quasi-two-state kinetic analysis is performed to calculate the time-dependent population changes of the folded state. This calculation coincides with the experimentally obtained population changes at low temperature but deviations are found in the T-jump from 67 to 78 degrees C. Using temperature-dependent barrier height changes, a temperature Phi value analysis is performed. The result shows a decreasing trend of Phi(T) with temperature, which indicates an increase of the heterogeneity of the transition state. We conclude that ubiquitin unfolds along a well-defined pathway at low temperature which expands with increasing temperature to include multiple routes.     

Temperature dependence of the free energy landscape of the src-SH3 protein domain

The temperature dependence of the free energy landscape of the src-SH3 protein domain is investigated through fully atomic simulations in explicit solvent. Simulations are performed above and below the folding transition temperature, enabling an analysis of both protein folding and unfolding. The transition state for folding and unfolding, identified from the free energy surfaces, is found to be very similar, with structure in the central hydrophobic sheet and little structure throughout the rest of the protein. This is a result of a polarized folding (unfolding) mechanism involving early formation (late loss) of the central hydrophobic sheet at the transition state. Unfolding simulations map qualitatively well onto low-temperature free energy surfaces but appear, however, to miss important features observed in folding simulations. In particular, details of the folding mechanism involving the opening and closing of the hydrophobic core are not captured by unfolding simulations performed under strongly denaturing conditions. In addition, free energy surfaces at high temperatures do not display a desolvation barrier found at lower temperatures, involving the expulsion of water molecules from the hydrophobic core.     

Structural and dynamic effects of alpha-helix deletion in Sso7d: implications for protein thermal stability

Sso7d is a 62-residue protein from the hyperthemophilic archaeon Sulfolobus solfataricus with a denaturation temperature close to 100 degrees C around neutral pH. An engineered form of Sso7d truncated at leucine 54 (L54Delta) is significantly less stable, with a denaturation temperature of 53 degrees C. Molecular dynamics (MD) studies of Sso7d and its truncated form at two different temperatures have been performed. The results of the MD simulations at 300 K indicate that: (1) the flexibility of Sso7d chain at 300 K agrees with that detected from X-ray and NMR structural studies; (2) L54Delta remains stable in the native folded conformation and possesses an overall dynamic behavior similar to that of the parent protein. MD simulations performed at 500 K, 10 ns long, indicate that, while Sso7d is in-silico resistant to high temperature, the truncated variant partially unfolds, revealing the early phases of the thermal unfolding pathway of the protein. Analysis of the trajectories of L54Delta suggests that the unzipping of the N-terminal and C-terminal beta-strands should be the first event of the unfolding pathway, and points out the regions more resistant to thermal unfolding. These findings allow one to understand the role played by specific interactions connecting the two ends of the chain for the high thermal stability of Sso7d, and support recent hypotheses on its folding mechanism emerged from site-directed mutagenesis studies.     

Kinetic analysis of molecular dynamics simulations reveals changes in the denatured state and switch of folding pathways upon single-point mutation of a beta-sheet miniprotein

The effects of a single-point mutation on folding thermodynamics and kinetics are usually interpreted by focusing on the native structure and the transition state. Here, the entire conformational spaces of a 20-residue three-stranded antiparallel beta-sheet peptide (double hairpin) and of its single-point mutant W10V are sampled close to the melting temperature by equilibrium folding-unfolding molecular dynamics simulations for a total of 40 micros. The folded state as well as the most populated free energy basins in the denatured state are isolated by grouping conformations according to fast relaxation at equilibrium. Such kinetic analysis provides more detailed and useful information than a simple projection of the free energy. The W10V mutant has the same native structure as the wild type peptide, and similar folding rate and stability. In the denatured state, the N-terminal hairpin is about 20% more structured in W10V than the wild type mainly because of van der Waals interactions. Notably, the W10V mutation influences also the van der Waals energy at the transition state ensemble causing a shift in the ratio of fluxes between two different transition state regions on parallel folding pathways corresponding to nucleation at either of the two beta-hairpins. Previous experimental studies have focused on the effects of denaturant-dependent or temperature-dependent changes in the structure of the denatured state. The atomistic simulations show that a single-point mutation in the central strand of a beta-sheet peptide results in remarkable changes in the topography of the denatured state ensemble. These changes modulate the relative accessibility of parallel folding pathways because of kinetic partitioning of the denatured state. Therefore, the observed dependence of the folding process on the starting ensemble raises questions on the biological significance of in vitro folding studies under strongly denaturing conditions.     

Conformational transition states of a beta-hairpin peptide between the ordered and disordered conformations in explicit water

The conformational transition states of a beta-hairpin peptide in explicit water were identified from the free energy landscapes obtained from the multicanonical ensemble, using an enhanced conformational sampling calculation. The beta-hairpin conformations were significant at 300 K in the landscape, and the typical nuclear Overhauser effect signals were reproduced, consistent with the previously reported experiment. In contrast, the disordered conformations were predominant at higher temperatures. Among the stable conformations at 300 K, there were several free energy barriers, which were not visible in the landscapes formed with the conventional parameters. We identified the transition states around the saddle points along the putative folding and unfolding paths between the beta-hairpin and the disordered conformations in the landscape. The characteristic features of these transition states are the predominant hydrophobic contacts and the several hydrogen bonds among the side-chains, as well as some of the backbone hydrogen bonds. The unfolding simulations at high temperatures, 400 K and 500 K, and their principal component analyses also provided estimates for the transition state conformations, which agreed well with those at 400 K and 500 K deduced from the current free energy landscapes at 400 K and 500 K, respectively. However, the transition states at high temperatures were much more widely distributed on the landscape than those at 300 K, and their conformations were different.     

Role of native topology investigated by multiple unfolding simulations of four SH3 domains

The relative importance of amino acid sequence and native topology in the unfolding process of two SH3 domains and two circular permutants was investigated by 120 molecular dynamics runs at 375 K for a total simulation time of 0.72 micros. The alpha-spectrin (aSH3) and src SH3 (sSH3) domains, which have the same topology and a sequence identity of only 34%, show similar unfolding pathways. The disappearance of the three-stranded antiparallel beta-sheet is the last unfolding event, in agreement with a large repertoire of kinetic data derived from point mutations as well as glycine insertions and disulfide crosslinks. Two alternative routes of beta-sheet unfolding have emerged from the analysis of the trajectories. One is statistically preferred in aSH3 (n-src loop breaks before distal hair-pin) and the inverse in sSH3. An elongation of the beta2-beta3 hairpin was observed during the unfolding of sSH3 at 375 K and in 300 K simulations started from the putative transition state of sSH3 in accord with unusual kinetic data for point mutations at the n-src loop. The change of connectivity in the permutants influenced the sequence of unfolding events mainly at the permutation site. Regions where the connectivity remained unaffected showed the same chronology of contact disappearance. Taken together with previous folding simulations of two designed three-stranded antiparallel beta-sheet peptides, these results indicate that, at least for small beta-sheet proteins, the folding mechanism is primarily defined by the native state topology, whilst specific interactions determine the statistically predominant folding route.     

Temperature-dependent downhill unfolding of ubiquitin. II. Modeling the free energy surface

To provide evidence for the interpretation of temperature‐dependent unfolding kinetics and the downhill unfolding scenario presented in the accompanying experimental article (Part I), the free energy surface of ubiquitin unfolding is calculated using statistical mechanical models of the Muñoz‐Eaton (ME) form. The models allow only two states for each amino acid residue, folded or unfolded, and permutations of these states generate an ensemble of microstates. One‐dimensional free energy curves are calculated using the number of folded residues as a reaction coordinate. The proposed sequential unfolding of ubiquitin's β‐sheet is tested by mapping the free energy onto two reaction coordinates inspired by the experiment as follows: the number of folded residues in ubiquitin's stable β‐strands I and II and those of the less stable strands III–V. Although the original ME model successfully captures folding features of zipper‐like one‐dimensional folders, it misses important tertiary interactions between residues that are far from each other in primary sequence. To take tertiary contacts into account, partially folded microstates based on a spherical growth model are included in the calculation and compared with the original model. By calculating the folding probability of each residue for a given point on the free energy surface, the unfolding pathway of ubiquitin is visualized. At low temperature, thermal unfolding occurs along a sequential unfolding pathway as follows: disruption of the β‐strands III–V followed by unfolding of the strands I and II. At high temperature, multiple unfolding routes are formed. The heterogeneity of the transition state explains the global nonexponential unfolding observed in the T‐jump experiment at high temperature. The calculation also reports a high stability for the α‐helix of ubiquitin. Proteins 2008. © 2008 Wiley‐Liss, Inc.     

A breakdown of symmetry in the folding transition state of protein L

The 62 residue IgG binding domain of protein L consists of a central alpha-helix packed on a four-stranded beta-sheet formed by N and C-terminal beta-hairpins. The overall topology of the protein is quite symmetric: the beta-hairpins have similar lengths and make very similar interactions with the central helix. Characterization of the effects of 70 point mutations distributed throughout the protein on the kinetics of folding and unfolding reveals that this symmetry is completely broken during folding; the first beta-hairpin is largely structured while the second beta-hairpin and helix are largely disrupted in the folding transition state ensemble. The results are not consistent with a "hydrophobic core first" picture of protein folding; the first beta-hairpin appears to be at least as ordered at the rate limiting step in folding as the hydrophobic core.     

Folding and aggregation kinetics of a beta-hairpin

We have investigated the solution structure, equilibrium properties, and folding kinetics of a 17-residue beta-hairpin-forming peptide derived from the protein ubiquitin. NMR experiments show that at 4 degrees C the peptide has a highly populated beta-hairpin conformation. At protein concentrations higher than 0.35 mM, the peptide aggregates. Sedimentation equilibrium measurements show that the aggregate is a trimer, while NMR indicates that the beta-hairpin conformation is maintained in the trimer. The relaxation kinetics in nanosecond laser temperature-jump experiments reveal a concentration-independent microsecond phase, corresponding to beta-hairpin unfolding-refolding, and a concentration-dependent millisecond phase due to oligomerization. Kinetic modeling of the relaxation rates and amplitudes yields the folding and unfolding rates for the monomeric beta-hairpin, as well as assembly and disassembly rates for trimer formation consistent with the equilibrium constant determined by sedimentation equilibrium. When the net charge on the peptides and ionic strength were taken into account, the rate of trimer assembly approaches the Debye-Smoluchowski diffusion limit. At 300 K, the rate of formation of the monomeric hairpin is (17 micros)(-1), compared to rates of (0.8 micros)(-1) to (52 micros)(-1) found for other peptides. After using Kramers theory to correct for the temperature dependence of the pre-exponential factor, the activation energy for hairpin formation is near zero, indicating that the barrier to folding is purely entropic. Comparisons with previously measured rates for a series of hairpins are made to distinguish between zipper and hydrophobic collapse mechanisms. Overall, the experimental data are most consistent with the zipper mechanism in which structure formation is initiated at the turn, the mechanism predicted by the Ising-like statistical mechanical model that was developed to explain the equilibrium and kinetic data for the beta-hairpin from protein GB1. In contrast, the majority of simulation studies favor a hydrophobic collapse mechanism. However, with few exceptions, there is little or no quantitative comparison of the simulation results with experimental data.     

Peroxidase activity as a tool for studying the folding of c-type cytochromes

The peroxidase activity of c-type cytochromes increases substantially by unfolding. This phenomenon was used to study the equilibrium unfolding of ferricytochrome c. The peroxidase activity is already enhanced at low denaturant concentrations. The lowest free energy folding intermediate is easily detected by this method, while it is invisible using fluorescence or optical spectroscopy. The free energy difference between this folding intermediate and the native state depends on the strength of the sixth ligand of the heme-iron and the increase in peroxidase activity upon unfolding is shown to be a sensitive indicator of the strength of this ligand. Under fully denaturing conditions, the peroxidase activity is inhibited by protein-based ligands. It is shown that at least three different ligand groups can be responsible for this inhibition, and that at neutral or alkaline pH, the predominant ligand is not histidine. The use of peroxidase activity assays as a method to study the unfolding of cytochrome c is evaluated.     

Unfolding of hen egg lysozyme by molecular dynamics simulations at 300K: insight into the role of the interdomain interface

We present the results of two 1.2 ns molecular dynamics (MD) unfolding simulations on hen egg lysozyme in water at 300K, performed using a new procedure called PEDC (Path Exploration With Distance Constraints). This procedure allows exploration of low energy structures as a function of increasing RMSD from the native structure, and offers especially the possibility of extensive exploration of the conformational space during the initial unfolding stages. The two independent MD simulations gave similar chronology of unfolding events: disruption of the active site, kinking of helix C, partial unfolding of the three-stranded beta-sheet to a two-stranded sheet (during which the helices A, B, and D remain to a great extent native), and finally unfolding of the beta-domain and partial unfolding of the alpha-domain in which hydrophobic clusters persist. We show particularly that the loss of hydrophobic contacts between the beta-sheet turn residues Leu55 and Ile56 and the hydrobic patch of the alpha-domain destabilizes the beta-domain and leads to its unfolding, suggesting that the correct embedding of these residues in the alpha-beta interface may constitute the rate limiting step in folding. These results are in accord with experimental observations on the folding/unfolding behavior of hen egg lysozyme at room temperature. They would also explain the loss of stability and the tendency to aggregation observed for the mutant Leu55Thr, and the slow refolding kinetics observed in the analogous amyloidogenic variant of human lysozyme.     

Molecular dynamics simulations of synthetic peptide folding

Because the time scale of protein folding is much greater than that of the widely used simulations of native structures, a detailed report of molecular dynamics simulations of folding has not been available. In this study, we Included the average solvent effect in the potential functions to simplify the calculation of the solvent effect and carried out long molecular dynamics simulations of the alanine-based synthetic peptides at 274 K. From either an extended or a randomly generated conformation, the simulations approached a helix-coil equilibrium in about 3 ns. The multiple minima problem did not prevent helix folding. The calculated helical ratio of Ac-AAQAAAAQAAAAQAAY-NH₂ was 47%, in good agreement with the circular dichroism measurement (about 50%). A helical segment with frayed ends was the most stable conformation, but the hydrophobic interaction favored the compact, distorted helix-turn-helix conformations. The transition between the two types of conformations occurred in a much larger time scale than helix propagation. The transient hydrogen bonds between the glutamine side chain and the backbone carbonyl group could reduce the free energy barrier of helix folding and unfolding. The substitution of a single alanine residue in the middle of the peptide with valine or glycine decreased the average helical ratio significantly, in agreement with experimental observations. © 1996 Wiley-Liss, Inc.