Conformational features of a peptide model Ac-DTVKLMYKGQPMTFR-NH2, corresponding to an early folding beta hairpin region of staphylococcal nuclease.

Recent H-D exchange 1H NMR studies of the refolding of Staphylococcal nuclease (P117G) variant suggest that, a region of the protein corresponding to a beta hairpin in the native structure folded early in the refolding process. In order to investigate whether the formation of beta hairpin is an early folding event, we investigated the conformational features of the beta hairpin peptide model Ac-DTVKLMYKGQPMTFR-NH2 from Staphylococcal nuclease with 1H NMR techniques. It appears that the peptide aggregates even at a low concentration. However, based on the observation of weak dnn(i, i + 1) NOEs between K8-G9, G9-Q10, an upfield shift of Gly9 NH and a low temperature coefficient (-d delta/dT) for Gly9 NH, we suggest that the sequence YKGQP as part of the beta hairpin peptide model samples conformational forms with reduced conformational entropy and turn potential. The presence of aggregation could be restricting the population of folded conformational forms and formation of beta hairpin at detectable concentrations. We suggest that, formation of beta hairpin could be an early event in the folding of Staphylococcal nuclease and this observation correlates with H-D exchange 1H NMR results and also with the prediction of a protein folding model proposed in literature.     

Molecular dynamics simulations of certain mutant peptide models from staphylococcal nuclease reveal that initial hydrophobic collapse associated with turn propensity drives β‐hairpin folding

An important nucleation event during the folding of staphylococcal nuclease involves the formation of a β‐hairpin by the sequence ²¹DTVKLMYKGQPMTFR³⁵. Earlier studies show that the turn sequence ‘YKGQP’ has an important role in the folding of this β‐hairpin. To understand the active or passive nature of the turn sequence ‘YKGQP’ in the folding of the aforementioned β‐hairpin sequence, we studied glycine mutant peptides Ac‐²DTVKLMYGGQPMTFR¹⁶‐NMe (K9G:15), Ac‐²DTVKLMYKGGPMTFR¹⁶‐NMe (Q11G:15), Ac‐²DTVKLMYGGGPMTFR¹⁶‐NMe (K9G/Q11G:15), and Ac‐²DTVKLMGGGGGMTFR¹⁶‐NMe (penta‐G:15) by using molecular dynamics simulations, starting with two different unfolded states, polyproline II and extended conformational forms. Further, 5mer mutant turn peptides Ac‐²YGGQP⁶‐NMe (K3G:5), Ac‐²YKGGP⁶‐NMe (Q5G:5), Ac‐²YGGGP⁶‐NMe (K3G/Q5G:5), and Ac‐²GGGGG⁶‐NMe (penta‐G:5) were also studied individually. Our results show that an initial hydrophobic collapse and loop closure occurs in all 15mer mutants, but only K9G:15 mutant forms a stable native‐like β‐hairpin. In the other 15mer mutants, the hydrophobic collapsed state would not proceed to β‐hairpin formation. Of the different simulations performed for the penta‐G:15 mutant, in only one simulation a nonnative β‐hairpin conformation is sampled with highly flexible loop region (⁸GGGGG¹²), which has no specific conformational preference as a 5mer. While the sequence ‘YGGQP’ in the K3G:5 simulation shows relatively higher β‐turn propensity, the presence of this sequence in K9G:15 peptide seems to be driving the β‐hairpin formation. Thus, these results seem to suggest that for the formation of a stable β‐hairpin, the initial hydrophobic collapse is to be assisted by a turn propensity. Initial hydrophobic collapse alone is not sufficient to guide β‐hairpin formation. Copyright © 2013 European Peptide Society and John Wiley & Sons, Ltd.     

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.     

Stress and strain in staphylococcal nuclease.

Protein molecules generally adopt a tertiary structure in which all backbone and side chain conformations are arranged in local energy minima; however, in several well-refined protein structures examples of locally strained geometries, such as cis peptide bonds, have been observed. Staphylococcal nuclease A contains a single cis peptide bond between residues Lys 116 and Pro 117 within a type VIa beta-turn. Alternative native folded forms of nuclease A have been detected by NMR spectroscopy and attributed to a mixture of cis and trans isomers at the Lys 116-Pro 117 peptide bond. Analyses of nuclease variants K116G and K116A by NMR spectroscopy and X-ray crystallography are reported herein. The structure of K116A is indistinguishable from that of nuclease A, including a cis 116-117 peptide bond (92% populated in solution). The overall fold of K116G is also indistinguishable from nuclease A except in the region of the substitution (residues 112-117), which contains a predominantly trans Gly 116-Pro 117 peptide bond (80% populated in solution). Both Lys and Ala would be prohibited from adopting the backbone conformation of Gly 116 due to steric clashes between the beta-carbon and the surrounding residues. One explanation for these results is that the position of the ends of the residue 112-117 loop only allow trans conformations where the local backbone interactions associated with the phi and psi torsion angles are strained. When the 116-117 peptide bond is cis, less strained backbone conformations are available. Thus the relaxation of the backbone strain intrinsic to the trans conformation compensates for the energetically unfavorable cis X-Pro peptide bond. With the removal of the side chain from residue 116 (K116G), the backbone strain of the trans conformation is reduced to the point that the conformation associated with the cis peptide bond is no longer favorable.     

Understanding the roles of amino acid residues in tertiary structure formation of chignolin by using molecular dynamics simulation

Chignolin is a 10-residue peptide (GYDPETGTWG) that forms a stable beta-hairpin structure in water. However, its design template, GPM12 (GYDDATKTFG), does not have a specific structure. To clarify which amino acids give it the ability to form the beta-hairpin structure, we calculated the folding free-energy landscapes of chignolin, GPM12, and their chimeric peptides using multicanonical molecular dynamics (MD) simulation. Cluster analysis of the conformational ensembles revealed that the native structure of chignolin was the lowest in terms of free energy while shallow local minima were widely distributed in the free energy landscape of GPM12, in agreement with experimental observations. Among the chimeric peptides, GPM12(D4P/K7G) stably formed the same beta-hairpin structure as that of chignolin in the MD simulation. This was confirmed by nuclear magnetic resonance (NMR) spectroscopy. A comparison of the free-energy landscapes showed that the conformational distribution of the Asp3-Pro4 sequence was inherently biased in a way that is advantageous both to forming hydrogen bonds with another beta-strand and to initiating loop structure. In addition, Gly7 helps stabilize the loop structure by having a left-handed alpha-helical conformation. Such a conformation is necessary to complete the loop structure, although it is not preferred by other amino acids. Our results suggest that the consistency between the short-range interactions that determine the local geometries and the long-range interactions that determine the global structure is important for stable tertiary structure formation.     

Loop and stem dynamics during RNA hairpin folding and unfolding

2-Aminopurine (2AP) is a fluorescent adenine analog that probes mainly base stacking in nucleic acids. We labeled the loop or the stem of the RNA hairpin gacUACGguc with 2AP to study folding thermodynamics and kinetics at both loci. Thermal melts and fast laser temperature jumps detected by 2AP fluorescence monitored the stability and folding/unfolding kinetics. The observed thermodynamic and kinetic traces of the stem and loop mutants, though strikingly different at a first glance, can be fitted to the same free-energy landscape. The differences between the two probe locations arise because base stacking decreases upon unfolding in the stem, whereas it increases in the loop. We conclude that 2AP is a conservative adenine substitution for mapping out the contributions of different RNA structural elements to the overall folding process. Molecular dynamics (MD) totaling 0.6 μsec were performed to look at the conformations populated by the RNA at different temperatures. The combined experimental data, and MD simulations lead us to propose a minimal four-state free-energy landscape for the RNA hairpin. Analysis of this landscape shows that a sequential folding model is a good approximation for the full folding dynamics. The frayed state formed initially from the native state is a heterogeneous ensemble of structures whose stem is frayed either from the end or from the loop.     

β‐Hairpin folding and stability: molecular dynamics simulations of designed peptides in aqueous solution

The structural properties of a 10‐residue and a 15‐residue peptide in aqueous solution were investigated by molecular dynamics simulation. The two designed peptides, SYINSDGTWT and SESYINSDGTWTVTE, had been studied previously by NMR at 278 K and the resulting model structures were classified as 3:5 β‐hairpins with a type I + G1 β‐bulge turn. In simulations at 278 K, starting from the NMR model structure, the 3:5 β‐hairpin conformers proved to be stable over the time period evaluated (30 ns). Starting from an extended conformation, simulations of the decapeptide at 278 K, 323 K and 353 K were also performed to study folding. Over the relatively short time scales explored (30 ns at 278 K and 323 K, 56 ns at 353 K), folding to the 3:5 β‐hairpin could only be observed at 353 K. At this temperature, the collapse to β‐hairpin‐like conformations is very fast. The conformational space accessible to the peptide is entirely dominated by loop structures with different degrees of β‐hairpin character. The transitions between different types of ordered loops and β‐hairpins occur through two unstructured loop conformations stabilized by a single side‐chain interaction between Tyr2 and Trp9, which facilitates the changes of the hydrogen‐bond register. In agreement with previous experimental results, β‐hairpin formation is initially driven by the bending propensity of the turn segment. Nevertheless, the fine organization of the turn region appears to be a late event in the folding process. Copyright © 2004 European Peptide Society and John Wiley & Sons, Ltd.     

Dissecting the stability of a beta-hairpin peptide that folds in water: NMR and molecular dynamics analysis of the beta-turn and beta-strand contributions to folding.

NMR studies of the folding and conformational properties of a beta-hairpin peptide, several peptide fragments of the hairpin, and sequence-modified analogues, have enabled the various contributions to beta-hairpin stability in water to be dissected. Temperature and pH-induced unfolding studies indicate that the folding-unfolding equilibrium approximates to a two-state model. The hairpin is highly resistant to denaturation and is still significantly folded in 7 M urea at 298 K. Thermodynamic analysis shows the hairpin to fold in water with a significant change in heat capacity, however, DeltaCp degrees in 7 M urea is reduced. V/Y-->A mutations on one strand of the hairpin reduce folding to <10 %, consistent with a hydrophobic stabilisation model. We show that in a truncated peptide (residues 6-16) lacking the hydrophobic residues on one beta-strand, the type I' Asn-Gly turn in the sequence SINGKK is significantly populated in water in the absence of interstrand hydrophobic contacts. Unrestrained molecular dynamics simulations of unfolding, using an explicit solvation model, show that the conformation of the NG turn persists for longer than the AG analogue, which has a much lower propensity for type I' turn formation from a data base analysis of preferred turns. The origin of the high stability of the Asn-Gly turn is not entirely clear; data base analysis of 66 NG turns, together with molecular dynamics simulations, reveals no participation of the Asn side-chain in turn-stabilising interactions with the peptide backbone. However, hydration analysis of the molecular dynamics simulations reveals a pocket of "high density" water bridging between the Asn side-chain and peptide main-chain that suggests solvent-mediated interactions may play an important role in modulating phi,psi propensities in the NG turn region.     

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.     

Folding landscapes of the Alzheimer amyloid-beta(12-28) peptide

The energy landscape for folding of the 12-28 fragment of the Alzheimer amyloid beta (Abeta) peptide is characterized using replica-exchange molecular dynamics simulations with an all-atom peptide model and explicit solvent. At physiological temperatures, the peptide exists mostly as a collapsed random coil, populating a small fraction (less than 10%) of hairpins with a beta-turn at position V18F19, with another 10% of hairpin-like conformations possessing a bend rather than a turn in the central VFFA positions. A small fraction of the populated states, approximately 14%, adopt polyproline II (PPII) conformations. Folding of the structured hairpin states proceeds through the assembly of two locally stable segments, VFFAE and EDVGS. The interactions stabilizing these locally folded structural motifs are in conflict with those stabilizing the global fold of A12-28, a signature of underlying residual frustration in this peptide. At increased temperature, the population of both beta-strand and PPII conformations diminishes in favor of beta-turn and random-coil states. On the basis of the conformational preferences of Abeta 12-28 monomers, two models for the molecular structure of amyloid fibrils formed by this peptide are proposed.     

Conformations of alanine-based peptides in water probed by FTIR, Raman, vibrational circular dichroism, electronic circular dichroism, and NMR spectroscopy

We have used a combination of FTIR, VCD, ECD, Raman, and NMR spectroscopies to probe the solution conformations sampled by H-(AAKA)-OH by utilizing an excitonic coupling model and constraints imposed by the 3JCalphaHNH coupling constants of the central residues to simulate the amide I' profile of the IR, isotropic Raman, anisotropic Raman, and VCD spectra in terms of a mixture of three conformations, i.e., polyproline II, beta-strand and right-handed helical. The representative coordinates of the three conformations were obtained from published coil libraries. Alanine was found to exhibit PPII fractions of 0.60 or greater, mixed with smaller fractions of helices and beta-strand conformations. Lysine showed no clear conformational propensity in that it samples polyproline II, beta-strand, and helical conformations with comparable probability. This is at variance with results obtained earlier for ionized polylysine, which suggest a high polyproline II propensity. We reanalyzed previously investigated tetra- and trialanine by combining published vibrational spectroscopy data with 3JCalphaHNH coupling constants and obtained again blends dominated by PPII with smaller admixtures of beta-strand and right-handed helical conformations. The polyproline II propensity of alanine was found to be higher in tetraalanine than in trialanine. For all peptides investigated, our results rule out a substantial population of turn-like conformations. Our results are in excellent agreement with MD simulations on short alanine peptides by Gnanakaran and Garcia [(2003) J. Phys. Chem. B 107, 12555-12557] but at variance with multiple MD simulations particularly for the alanine dipeptide.     

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

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

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

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

Cooperative folding mechanism of a beta-hairpin peptide studied by a multicanonical replica-exchange molecular dynamics simulation

G-peptide is a 16-residue peptide of the C-terminal end of streptococcal protein G B1 domain, which is known to fold into a specific beta-hairpin within 6 micros. Here, we study molecular mechanism on the stability and folding of G-peptide by performing a multicanonical replica-exchange (MUCAREM) molecular dynamics simulation with explicit solvent. Unlike the preceding simulations of the same peptide, the simulation was started from an unfolded conformation without any experimental information on the native conformation. In the 278-ns trajectory, we observed three independent folding events. Thus MUCAREM can be estimated to accelerate the folding reaction more than 60 times than the conventional molecular dynamics simulations. The free-energy landscape of the peptide at room temperature shows that there are three essential subevents in the folding pathway to construct the native-like beta-hairpin conformation: (i) a hydrophobic collapse of the peptide occurs with the side-chain contacts between Tyr45 and Phe52, (ii) then, the native-like turn is formed accompanying with the hydrogen-bonded network around the turn region, and (iii) finally, the rest of the backbone hydrogen bonds are formed. A number of stable native hydrogen bonds are formed cooperatively during the second stage, suggesting the importance of the formation of the specific turn structure. This is also supported by the accumulation of the nonnative conformations only with the hydrophobic cluster around Tyr45 and Phe52. These simulation results are consistent with high phi-values of the turn region observed by experiment.     

Searching for folding initiation sites of staphylococcal nuclease: a study of N-terminal short fragments

The N-terminal short fragments of staphylococcal nuclease (SNase), SNase20, SNase28, and SNase36, corresponding to the sequence regions, Ala1-Gly20, Ala1-Lys28, and Ala1-Leu36, respectively, as well as an 8-residue peptide (Ala17-Ile18-Asp19-Gly20-Asp21-Thr22-Val23-Lys24) have been synthesized. The conformational states of these fragments were investigated using CD and NMR spectroscopy in aqueous solution and in trifluoroethanol (TFE)-H(2)O mixture. SNase20 containing a sequence corresponding to a bent peptide in native SNase shows a transient population of bend-like conformation around Ala12-Thr13-Leu14 in TFE-H(2)O mixture. The sequence region of Ala17-Thr22 of SNase28 displays a localized propensity for turn-like conformation in both aqueous solution and TFE-H(2)O mixture. The conformational ensemble of SNase36 in aqueous solution includes populated turn-like conformations localized in sequence regions Ala17-Thr22 and Tyr27-Gln30. The analysis suggests that these sequence regions, which form the regular secondary structures in native protein, may serve as the folding nucleation sites of SNase fragments of different chain lengths starting from the N-terminal end. Thus, the formation of bend- and turn-like conformations of these sequence regions may be involved in the early folding events of the SNase polypeptide chain in vitro.     

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.     

A pentapeptide model for an early folding step in the refolding of staphylococcal nuclease: The role of its turn propensity

Recently the folding of a staphylococcal nuclease (P117G) variant was examined with the hydrogen-deuterium (H-D) exchange technique. Many of the residues that showed significant protection are located in protection are located in β-sheet regions. About half the residues protected belong to an antiparallel β-hairpin structure (residues 21–35) in the native structure. The β-hairpin structure is formed by strands 2 and 3 of sheet 2 connected by the sequence²⁷ Y KGQP³¹ in a type I′ reverse turn conformation with a 4 → 1 hydrogen bonding between Q30 NH and Y27 C=O. We have targeted the conformational characterization of the peptide model Ac-YKGQP-NH₂ with ¹II two-dimensional nmr techniques in aqueous solution with a view to assessing its propensity to sample turn conformational forms and thus initiate the formation of β-hairpin structure. Based upon the observed dαn (i, i + 1), dαn (i, i + 3), and dnn (i, i + 1) nuclear Overhauser effect connectivities, temperature coefficients for amide protons and conformational analysis with quantum mechanical perturbative configuration interaction over localized orbitals method, we conclude that the model peptide samples turn conformational forms with reduced conformational entropy. We suggest that the turn can nucleate the formation of the β-hairpin structure in the refolding of nuclease. Observation of turn propensity for this sequence is consistent with the folding mechanism of the Greek key motif (present in Staphylococcal nuclease) proposed in the literature. © 1997 John Wiley & Sons, Inc.     

Folding of a DNA hairpin loop structure in explicit solvent using replica-exchange molecular dynamics simulations

Hairpin loop structures are common motifs in folded nucleic acids. The 5'-GCGCAGC sequence in DNA forms a characteristic and stable trinucleotide hairpin loop flanked by a two basepair stem helix. To better understand the structure formation of this hairpin loop motif in atomic detail, we employed replica-exchange molecular dynamics (RexMD) simulations starting from a single-stranded DNA conformation. In two independent 36 ns RexMD simulations, conformations in very close agreement with the experimental hairpin structure were sampled as dominant conformations (lowest free energy state) during the final phase of the RexMDs ( approximately 35% at the lowest temperature replica). Simultaneous compaction and accumulation of folded structures were observed. Comparison of the GCA trinucleotides from early stages of the simulations with the folded topology indicated a variety of central loop conformations, but arrangements close to experiment that are sampled before the fully folded structure also appeared. Most of these intermediates included a stacking of the C(2) and G(3) bases, which was further stabilized by hydrogen bonding to the A(5) base and a strongly bound water molecule bridging the C(2) and A(5) in the DNA minor groove. The simulations suggest a folding mechanism where these intermediates can rapidly proceed toward the fully folded hairpin and emphasize the importance of loop and stem nucleotide interactions for hairpin folding. In one simulation, a loop motif with G(3) in syn conformation (dihedral flip at N-glycosidic bond) accumulated, resulting in a misfolded hairpin. Such conformations may correspond to long-lived trapped states that have been postulated to account for the folding kinetics of nucleic acid hairpins that are slower than expected for a semiflexible polymer of the same size.     

The disordered mobile loop of GroES folds into a defined beta-hairpin upon binding GroEL

The GroES mobile loop is a stretch of approximately 16 amino acids that exhibits a high degree of flexible disorder in the free protein. This loop is responsible for the interaction between GroES and GroEL, and it undergoes a folding transition upon binding to GroEL. Results derived from a combination of transferred nuclear Overhauser effect NMR experiments and molecular dynamics simulations indicate that the mobile loop adopts a beta-hairpin structure with a Type I, G1 Bulge turn. This structure is distinct from the conformation of the loop in the co-crystal of GroES with GroEL-ADP but identical to the conformation of the bacteriophage-panned "strongly binding peptide" in the co-crystal with GroEL. Analysis of sequence conservation suggests that sequences of the mobile loop and strongly binding peptide were selected for the ability to adopt this hairpin conformation.