The formation of the N-terminal β-hairpin of ubiquitin is thought to be an early event in the folding of this small protein. Previously, we have shown that a peptide corresponding to residues 1–17 of ubiquitin folds autonomously and is likely to have a native-like hairpin register. To investigate the causes of the stability of this fold, we have made mutations in the amino acids at the apex of the turn. We find that in a peptide where Thr9 is replaced by Asp, U(1–17)T9D, the native conformation is stabilized with respect to the wild-type sequence, so much so that we are able to characterize the structure of the mutant peptide fully by NMR spectroscopy. The data indicate that U(1–17)T9D peptide does indeed form a hairpin with a native-like register and a type I turn with a G1 β-bulge, as in the full-length protein. The reason for the greater stability of the U(1–17)T9D mutant remains uncertain, but there are nuclear Overhauser effects between the side chains of Asp9 and Lys11, which may indicate that a charge–charge interaction between these residues is responsible.

Previous conformational analysis of 10-residue linear peptides enabled us to identify some cross-strand side-chain interactions that stabilize β-hairpin conformations. The stabilizing influence of these interactions appeared to be greatly reduced when the interaction was located at the N- and C-termini of these 10-residue peptides. To investigate the effect of the position relative to the turn of favorable interactions on β-hairpin formation, we have designed two 15-residue β-hairpin forming peptides with the same residue composition and differing only in the location of two residues within the strand region. The conformational properties of these two peptides in aqueous solution were studied by 1H and 13C NMR. Differences in the conformational behavior of the two designed 15-residue peptides suggest that the influence of stabilizing factors for β-hairpin formation, in particular, cross-strand side-chain interactions, depends on their proximity to the turn. Residues adjacent to the turn are most efficient in that concern. This result agrees with the proposal that the turn region acts as the driving force in β-hairpin folding.

The structural and dynamical behavior of the 41–56 β-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 300 K (>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 N-terminal 17 residues of ubiquitin have been shown by 1H NMR to fold autonomously into a β-hairpin structure in aqueous solution. This structure has a specific, native-like register, though side-chain contacts differ in detail from those observed in the intact protein. An autonomously folding hairpin has previously been identified in the case of streptococcal protein G, which is structurally homologous with ubiquitin, but remarkably, the two are not in topologically equivalent positions in the fold. This suggests that the organization of folding may be quite different for proteins sharing similar tertiary structures. Two smaller peptides have also been studied, corresponding to the isolated arms of the N-terminal hairpin of ubiquitin, and significant differences from simple random coil predictions observed in the spectra of these subfragments, suggestive of significant limitation of the backbone conformational space sampled, presumably as a consequence of the strongly β-structure favoring composition of the sequences. This illustrates the ability of local sequence elements to express a propensity for β-structure even in the absence of actual sheet formation. Attempts were made to estimate the population of the folded state of the hairpin, in terms of a simple two-state folding model. Using published “random coil” values to model the unfolded state, and values derived from native ubiquitin for the putative unique, folded state, it was found that the apparent population varied widely for different residues and with different NMR parameters. Use of the spectra of the subfragment peptides to provide a more realistic model of the unfolded state led to better agreement in the estimates that could be obtained from chemical shift and coupling constant measurements, while making it clear that some other approaches to population estimation could not give meaningful results, because of the tendency to populate the β-region of conformational space even in the absence of the hairpin structure.

The usefulness of molecular dynamics to assess the structural integrity of mutants containing several mutations has been investigated. Our goal was to determine whether molecular dynamics would be able to discriminate mutants of a protein having a close–to–wild-type fold, from those that are not folded under the same conditions. We used as a model the B1 domain of protein G in which we replaced the unique central α-helix by the sequence of the second β-hairpin, which has a strong intrinsic propensity to form this secondary structure in solution. In the resulting protein, one-third of the secondary structure has been replaced by a non-native one. Models of the mutants were built based on the three-dimensional structure of the wild-type GB1 domain. During 2 ns of molecular dynamics simulations on these models, mutants containing up to 10 mutations in the helix retained the native fold, while another mutant with an additional mutation unfolded. This result is in agreement with our circular dichroism and NMR experiments, which indicated that the former mutants fold into a structure similar to the wild-type, as opposed to the latter mutant which is partly unfolded. Additionally, a mutant containing six mutations scattered through the surface of the domain, and which is unfolded, was also detected by the simulation. This study suggests that molecular dynamics calculations could be performed on molecular models of mutants of a protein to evaluate their foldability, prior to a mutagenesis experiment.