NMR Characterization of a Peptide Model Provides Evidence for Significant Structure in the Unfolded State of the Villin Headpiece Helical Subdomain

Quantitative Analysis of Protein Unfolded State Energetics: Experimental and Computational Studies Demonstrate That Non-Native Side-Chain Interactions Stabilize Local Native Backbone Structure

Proteins fold on relatively smooth free energy landscapes which are biased toward the native state, but even simple topologies which fold rapidly can experience roughness on their free energy landscape. The details of these interactions are difficult to elucidate experimentally. Closely related to the problem of deciphering the details of the free energy landscape is the problem of defining the interactions in the denatured state ensemble (DSE) which is populated under native conditions, that is, under conditions where the native state is stable. The DSE of many proteins deviates from random coil models, but quantifying and defining the energetics of the transiently populated interactions in this ensemble is extremely challenging. Characterization of the DSE of proteins which fold to compact structures is also relevant to studies of intrinsically disordered proteins (IDPs) since interactions in the dynamic ensemble populated by IDPs can modulate their behavior. Here we show how experimental thermodynamic and pK_a measurements can be combined with computational thermodynamic integration to quantify interactions in the DSE. We show that non-native side chain interactions can stabilize native backbone structure in the DSE and demonstrate that that even rapidly folding proteins can form energetically significant non-native interactions in their DSE. As an example, we characterize a non-native salt bridge that stabilizes local native backbone structure in the DSE of a widely studied fast-folding protein, the villin headpiece helical domain. The combined computational experimental approach is applicable to other protein unfolded states and provides insight that is impossible to obtain with either method alone.

Optimized purification of a fusion protein by reversed-phase high performance liquid chromatography informed by the linear solvent strength model

Fusion protein systems are commonly used for expression of small proteins and peptides. An important criterion for a fusion protein system to be useful is the ability to separate the protein of interest from the tag. Additionally, because no protease cleaves fusion proteins with 100% efficiency, the ability to separate the desired peptide from any remaining uncleaved protein is also necessary. This is likely to be the more difficult task as at least a portion of the sequence of the fusion protein is identical to that of the protein of interest. When a high level of purity is required, gradient elution reversed-phase HPLC is frequently used as a final purification step. Shallow gradients are often advantageous for maximizing both the purity and yield of the final product; however, the relationship between relative retention times at shallow gradients and those at steeper gradients typically used for analytical HPLC are not always straightforward. In this work, we report reversed-phase HPLC results for the fusion protein system consisting of the N-terminal domain of ribosomal protein L9 (NTL9) and the 36-residue villin headpiece subdomain (HP36) linked by a recognition sequence for the protease factor Xa. This system represents an excellent example of the difficulties in purification that may arise from this unexpected elution behavior at shallow gradients. Additionally, we report on the sensitivity of this elution behavior to the concentration of the additive trifluoroacetic acid in the mobile phase and present optimized conditions for separating HP36 from the full fusion protein by reversed-phase HPLC using a shallow gradient. Finally, we suggest that these findings are relevant to the purification of other fusion protein systems, for which similar problems may arise, and support this suggestion using insights from the linear solvent strength model of gradient elution liquid chromatography.

Impact of Turn Propensity on the Folding Rates of Z34C Protein: Implications for the Folding of Helix-Turn-Helix Motif

The rate-limiting step for the folding of the helix-turn-helix (HTH) protein, Z34C, involves β-turn region ²⁰DPNL²³. This reverse turn has been observed to be part of the transition state in the folding process for Z34C, influencing its folding rates. Molecular dynamics simulations were performed on this turn peptide and its two mutants, D20A and P21A, to study turn formation using GROMOS54A7 force field. We find that this region has a turn propensity of its own, and the highest turn propensity is observed for the wild-type, which correlates well with available experimental results. We also find that a slight unfavorable change in ΔG _{turn folding} causes a drastic change in the folding rates of HTH motif and a mechanistic interpretation is given. Implications of these observations for the folding of the HTH protein Z34C are discussed.

Folding stability modulation of the villin headpiece helical subdomain by 4-fluorophenylalanine and 4-methylphenylalanine

HP36, the helical subdomain of villin headpiece, contains a hydrophobic core composed of three phenylalanine residues (Phe47, Phe51, and Phe58). Hydrophobic effects and electrostatic interactions were shown to be the critical factors in stabilizing this core and the global structure. To assess the interactions among Phe47, Phe51, and Phe58 residues and investigate how they affect the folding stability, we implanted 4-fluorophenylalanine (Z) and 4-methylphenylalanine (X) into the hydrophobic core of HP36. We chemically synthesized HP36 and its seven variants including four single mutants whose Phe51 or Phe58 was replaced with Z or X, and three double mutants whose Phe51 and Phe58 were both substituted. Circular dichroism and nuclear magnetic resonance measurements show that the variants exhibit a native HP36 like fold, of which F51Z and three double mutants are more stable than the wild type. Molecular modeling provided detailed interaction energy within the phenylalanine residues, revealing that electrostatic interactions dominate the stability modulation upon the introduction of 4-fluorophenylalanine and 4-methylphenylalanine. Our results show that these two non-natural amino acids can successfully tune the interactions in a relatively compact hydrophobic core and the folding stability without inducing dramatic steric effects. Such an approach may be applied to other folded motifs or proteins.

Characterizing a partially ordered miniprotein through folding molecular dynamics simulations: Comparison with the experimental data

The villin headpiece helical subdomain (HP36) is one of the best known model systems for computational studies of fast-folding all-α miniproteins. HP21 is a peptide fragment-derived from HP36-comprising only the first and second helices of the full domain. Experimental studies showed that although HP21 is mostly unfolded in solution, it does maintain some persistent native-like structure as indicated by the analysis of NMR-derived chemical shifts. Here we compare the experimental data for HP21 with the results obtained from a 15-μs long folding molecular dynamics simulation performed in explicit water and with full electrostatics. We find that the simulation is in good agreement with the experiment and faithfully reproduces the major experimental findings, namely that (a) HP21 is disordered in solution with <10% of the trajectory corresponding to transiently stable structures, (b) the most highly populated conformer is a native-like structure with an RMSD from the corresponding portion of the HP36 crystal structure of <1 Å, (c) the simulation-derived chemical shifts-over the whole length of the trajectory-are in reasonable agreement with the experiment giving reduced χ(2) values of 1.6, 1.4, and 0.8 for the Δδ(13) C(α) , Δδ(13) CO, and Δδ(13) C(β) secondary shifts, respectively (becoming 0.8, 0.7, and 0.3 when only the major peptide conformer is considered), and finally, (d) the secondary structure propensity scores are in very good agreement with the experiment and clearly indicate the higher stability of the first helix. We conclude that folding molecular dynamics simulations can be a useful tool for the structural characterization of even marginally stable peptides.

The unusual internal motion of the villin headpiece subdomain

The thermostable 36-residue subdomain of the villin headpiece (HP36) is the smallest known cooperatively folding protein. Although the folding and internal dynamics of HP36 and close variants have been extensively studied, there has not been a comprehensive investigation of side-chain motion in this protein. Here, the fast motion of methyl-bearing amino acid side chains is explored over a range of temperatures using site-resolved solution nuclear magnetic resonance deuterium relaxation. The squared generalized order parameters of methyl groups extensively spatially segregate according to motional classes. This has not been observed before in any protein studied using this methodology. The class segregation is preserved from 275 to 305 K. Motions detected in Helix 3 suggest a fast timescale of conformational heterogeneity that has not been previously observed but is consistent with a range of folding and dynamics studies. Finally, a comparison between the order parameters in solution with previous results based on solid-state nuclear magnetic resonance deuterium line shape analysis of HP36 in partially hydrated powders shows a clear disagreement for half of the sites. This result has significant implications for the interpretation of data derived from a variety of approaches that rely on partially hydrated protein samples.

A compact native 24-residue supersecondary structure derived from the villin headpiece subdomain

Many small proteins fold highly cooperatively in an all-or-none fashion and thus their native states are well protected from thermal fluctuations by an extensive network of interactions across the folded structure. Because protein structures are stabilized by local and nonlocal interactions among distal residues, dissecting individual substructures from the context of folded proteins results in large destabilization and loss of unique three-dimensional structure. Thus, mini-protein substructures can only rarely be derived from natural templates. Here, we describe a compact native 24-residues-long supersecondary structure derived from the hyperstable villin headpiece subdomain consisting of helices 2 and 3 (HP24). Using a combination of experimental techniques, including NMR and small-angle x-ray scattering, as well as all-atom replica exchange molecular-dynamics simulations, we show that a variant with stabilizing substitutions (HP24stab) forms a densely packed and compact conformation. In HP24stab, interactions between helices 2 and 3 are similar to those observed in native folded HP35, and the two helices cooperatively stabilize each other by completing the hydrophobic core lining the central part of HP35. Interestingly, even though the HP24wt fragment shows a more expanded and less structured conformation, NMR and simulations demonstrate a preference for a native-like topology. Thus, the two stabilizing residues in HP24stab shift the energy balance toward the native state, leading to a minimal folding motif.

Rational modification of protein stability by targeting surface sites leads to complicated results

The rational modification of protein stability is an important goal of protein design. Protein surface electrostatic interactions are not evolutionarily optimized for stability and are an attractive target for the rational redesign of proteins. We show that surface charge mutants can exert stabilizing effects in distinct and unanticipated ways, including ones that are not predicted by existing methods, even when only solvent-exposed sites are targeted. Individual mutation of three solvent-exposed lysines in the villin headpiece subdomain significantly stabilizes the protein, but the mechanism of stabilization is very different in each case. One mutation destabilizes native-state electrostatic interactions but has a larger destabilizing effect on the denatured state, a second removes the desolvation penalty paid by the charged residue, whereas the third introduces unanticipated native-state interactions but does not alter electrostatics. Our results show that even seemingly intuitive mutations can exert their effects through unforeseen and complex interactions.

Glassy dynamics of protein methyl groups revealed by deuteron NMR

We investigated site-specific dynamics of key methyl groups in the hydrophobic core of chicken villin headpiece subdomain (HP36) over the temperature range between 298 and 140 K using deuteron solid-state NMR longitudinal relaxation measurements. The relaxation of the longitudinal magnetization is weakly nonexponential (glassy) at high temperatures and exhibits a stronger degree of nonexponentiality below about 175 K. In addition, the characteristic relaxation times deviate from the simple Arrhenius law. We interpret this behavior via the existence of distribution of activation energy barriers for the three-site methyl jumps, which originates from somewhat different methyl environments within the local energy landscape. The width of the distribution of the activation barriers for methyl jumps is rather significant, about 1.4 kJ/mol. Our experimental results and modeling allow for the description of the apparent change at about 175 K without invoking a specific transition temperature. For most residues in the core, the relaxation behavior at high temperatures points to the existence of conformational exchange between the substates of the landscape, and our model takes into account the kinetics of this process. The observed dynamics are the same for dry and hydrated protein. We also looked at the effect of F58L mutation inside the hydrophobic core on the dynamics of one of the residues and observed a significant increase in its conformational exchange rate constant at high temperatures.

Low folding cooperativity of HP35 revealed by single-molecule force spectroscopy and molecular dynamics simulation

Some small proteins, such as HP35, fold at submicrosecond timescale with low folding cooperativity. Although these proteins have been extensively investigated, still relatively little is known about their folding mechanism. Here, using single-molecule force spectroscopy and steered molecule dynamics simulation, we study the unfolding of HP35 under external force. Our results show that HP35 unfolds at extremely low forces without a well-defined unfolding transition state. Subsequently, we probe the structure of unfolded HP35 using the persistence length obtained in the force spectroscopy. We found that the persistence length of unfolded HP35 is around 0.72 nm, >40% longer than typical unstructured proteins, suggesting that there are a significant amount of residual secondary structures in the unfolded HP35. Molecular dynamics simulation further confirmed this finding and revealed that many native contacts are preserved in HP35, even its two ends have been extended up to 8 nm. Our results therefore suggest that retaining a significant amount of secondary structures in the unfolded state of HP35 may be an efficient way to reduce the entropic cost for the formation of tertiary structure and increase the folding speed, although the folding cooperativity is compromised. Moreover, we anticipate that the methods we used in this work can be extended to the study of other proteins with complex folding behaviors and even intrinsically disordered ones.

Fast dynamics of HP35 for folded and urea-unfolded conditions

The changes in fast dynamics of HP35 with a double CN vibrational dynamics label (HP35-P(2)) as a function of the extent of denaturation by urea were investigated with two-dimensional infrared (2D IR) vibrational echo spectroscopy. Cyanophenylalanine (PheCN) replaces the native phenylalanine at two residues in the hydrophobic core of HP35, providing vibrational probes. NMR data show that HP35-P(2) maintains the native folded structure similar to wild type and that both PheCN residues share essentially the same environment within the peptide. A series of time-dependent 2D IR vibrational echo spectra were obtained for the folded peptide and the increasingly unfolded peptide. Analysis of the time dependence of the 2D spectra yields the system's spectral diffusion, which is caused by the sampling of accessible structures of the peptide under thermal equilibrium conditions. The structural dynamics become faster as the degree of unfolding is increased.

Thermally induced protein unfolding probed by isotope-edited IR spectroscopy

Infrared (IR) spectroscopy has been widely utilized for the study of protein folding, unfolding, and misfolding processes. We have previously developed a theoretical method for calculating IR spectra of proteins in the amide I region. In this work, we apply this method, in combination with replica-exchange molecular dynamics simulations, to study the equilibrium thermal unfolding transition of the villin headpiece subdomain (HP36). Temperature-dependent IR spectra and spectral densities are calculated. The spectral densities correctly reflect the unfolding conformational changes in the simulation. With the help of isotope labeling, we are able to capture the feature that helix 2 of HP36 loses its secondary structure before global unfolding occurs, in agreement with experiment.

Slow motions in the hydrophobic core of chicken villin headpiece subdomain and their contributions to configurational entropy and heat capacity from solid-state deuteron NMR measurements

We have investigated microsecond to millisecond time scale dynamics in several key hydrophobic core methyl groups of chicken villin headpiece subdomain protein (HP36) using a combination of single-site labeling, deuteron solid-state NMR line shape analysis, and computational modeling. Deuteron line shapes of hydrated powder samples are dominated by rotameric jumps and show a large variability of rate constants, activation energies, and rotameric populations. Site-specific activation energies vary from 6 to 38 kJ/mol. An additional mode of diffusion on a restricted arc is significant for some sites. In dry samples, the dynamics is quenched. Parameters of the motional models allow for calculations of configurational entropy and heat capacity, which, together with the rate constants, allow for observation of interplay between thermodynamic and kinetic picture of the landscape. Mutations at key phenylalanine residues at both distal (F47L&F51L) and proximal (F58L) locations to a relatively rigid side chain of L69 have a pronounced effect on alleviating the rigidity of this side chain at room temperature and demonstrate the sensitivity of the hydrophobic core environment to such perturbations.

Bridging the Gap between Folding Simulations and Experiments: The Case of the Villin Headpiece

The increasing accuracy of molecular dynamics force fields parameters and the increasing resolution of experimental results allow one to carefully compare and complement in silico data with experimental observations. Here, we study the human villin headpiece C-terminal helical subdomain (HP35) with the recent highly optimized Amber99SB*-ILDN force field and compare the results with recent high resolution triplet-triplet energy transfer (TTET) experiments. The correct reproduction of the main structural features reveals a good agreement between experimental data and simulations.

Hydrophobic core formation and dehydration in protein folding studied by generalized-ensemble simulations

Despite its small size, chicken villin headpiece subdomain HP36 folds into the native structure with a stable hydrophobic core within several microseconds. How such a small protein keeps up its conformational stability and fast folding in solution is an important issue for understanding molecular mechanisms of protein folding. In this study, we performed multicanonical replica-exchange simulations of HP36 in explicit water, starting from a fully extended conformation. We observed at least five events of HP36 folding into nativelike conformations. The smallest backbone root mean-square deviation from the crystal structure was 1.1 A. In the nativelike conformations, the stably formed hydrophobic core was fully dehydrated. Statistical analyses of the simulation trajectories show the following sequential events in folding of HP36: 1), Helix 3 is formed at the earliest stage; 2), the backbone and the side chains near the loop between Helices 2 and 3 take nativelike conformations; and 3), the side-chain packing at the hydrophobic core and the dehydration of the core side chains take place simultaneously at the later stage of folding. This sequence suggests that the initial folding nucleus is not necessarily the same as the hydrophobic core, consistent with a recent experimental phi-value analysis.

Dynamics of protein folding: probing the kinetic network of folding-unfolding transitions with experiment and theory

The problem of spontaneous folding of amino acid chains into highly organized, biologically functional three-dimensional protein structures continues to challenge the modern science. Understanding how proteins fold requires characterization of the underlying energy landscapes as well as the dynamics of the polypeptide chains in all stages of the folding process. In recent years, important advances toward these goals have been achieved owing to the rapidly growing interdisciplinary interest and significant progress in both experimental techniques and theoretical methods. Improvements in the experimental time resolution led to determination of the timescales of the important elementary events in folding, such as formation of secondary structure and tertiary contacts. Sensitive single molecule methods made possible probing the distributions of the unfolded and folded states and following the folding reaction of individual protein molecules. Discovery of proteins that fold in microseconds opened the possibility of atomic-level theoretical simulations of folding and their direct comparisons with experimental data, as well as of direct experimental observation of the barrier-less folding transition. The ultra-fast folding also brought new questions, concerning the intrinsic limits of the folding rates and experimental signatures of barrier-less "downhill" folding. These problems will require novel approaches for even more detailed experimental investigations of the folding dynamics as well as for the analysis of the folding kinetic data. For theoretical simulations of folding, a main challenge is how to extract the relevant information from overwhelmingly detailed atomistic trajectories. New theoretical methods have been devised to allow a systematic approach towards a quantitative analysis of the kinetic network of folding-unfolding transitions between various configuration states of a protein, revealing the transition states and the associated folding pathways at multiple levels, from atomistic to coarse-grained representations. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Thermodynamic consequences of incorporating 4-substituted proline derivatives into a small helical protein

Although proline residues are incompatible with an alpha-helix conformation, they fit well into the N-terminal end of alpha-helices. Proline can form either a C(gamma)-exo ring pucker or a C(gamma)-endo ring pucker. An electron-withdrawing substituent on the 4R position of proline favors an exo ring pucker while an endo ring pucker is preferred if the substituent is on the 4S position due to stereoelectronic effects. The villin headpiece subdomain (HP36) is a small helical protein composed of three alpha-helices and contains a proline residue (Pro62) at the N-terminus of its C-terminal alpha-helix. Pro62 has a C(gamma)-exo ring pucker and forms an aromatic-proline interaction, with Trp64 in the native structure. This work reports the use of 4-substituted proline derivatives, including (2S,4R)-4-hydroxyproline (Hyp), (2S,4R)-4-fluoroproline (Flp), (2S,4R)-4-methoxyproline (Mop), (2S,4S)-4-hydroxyproline (hyp), (2S,4S)-4-fluoroproline (flp), and (2S,4S)-4-methoxyproline (mop), to replace Pro62 and study how the pucker conformation affects the proline-aromatic interaction and the stability of HP36. CD and NMR measurements indicate that all of the HP36 variants incorporated with proline derivatives maintain a structure similar to that of the wild type. Thermal unfolding and urea-induced denaturation measurements have shown that all of the mutants with the exception of the one with the flp substitution are less stable than the wild type. Our results reveal that, upon the replacement of Pro62 to proline derivatives, not only do stereoelectronic effects influence the aromatic-proline interaction but the steric and hydrophobic effects induced by the substituents also play an important role in modulating the stability of HP36.

A critical assessment of putative gatekeeper interactions in the villin headpiece helical subdomain

The helical subdomain of the villin headpiece (HP36) is one of the smallest naturally occurring proteins that folds cooperatively. Its small size, rapid folding, and simple three-helix topology have made it an extraordinary popular model system for computational, theoretical, and experimental studies of protein folding. Aromatic-proline interactions involving Trp64 and Pro62 have been proposed to play a critical role in specifying the subdomain fold by acting as gatekeeper residues. Note that the numbering corresponds to full-length headpiece. Mutation of Pro62 has been shown to lead to a protein that does not fold, but this may arise for two different reasons: The residue may make interactions that are critical for the specificity of the fold or the mutation may simply destabilize the domain. In the first case, the protein cannot fold, while in the second, the small fraction of molecules that do fold adopt the correct structure. The modest stability of the wild type prevents a critical analysis of these interactions because even moderately destabilizing mutations lead to a very small folded state population. Using a hyperstable variant of HP36, denoted DM HP36, as our new wild type, we characterized a set of mutants designed to assess the role of the putative gatekeeper interactions. Four single mutants, DM Pro62Ala, DM Trp64Leu, DM Trp64Lys, and DM Trp64Ala, and a double mutant, DM Pro62Ala Trp64Leu, were prepared. All mutants are less stable than DM HP36, but all are well folded as judged by CD and (1)H NMR. All of the mutants display sigmoidal thermal unfolding and urea-induced unfolding curves. Double-mutant cycle analysis shows that the interactions between Pro62 and Trp64 are weak but favorable. Interactions involving Pro62 and proline-aromatic interactions are, thus, not required for specifying the subdomain fold. The implications for the design and thermodynamics of miniature proteins are discussed.

An unlocking/relocking barrier in conformational fluctuations of villin headpiece subdomain

A reversible structural unlocking reaction, in which the close-packed van der Waals interactions break cooperatively, has been found for the villin headpiece subdomain (HP35) using triplet-triplet-energy transfer to monitor conformational fluctuations from equilibrium. Unlocking is associated with an unfavorable enthalpy change (DeltaH(0) = 35 +/- 4 kJ/mol) which is nearly compensated in free energy by the entropy change (DeltaS(0) = 112 +/- 20 Jxmol(-1)xK(-1)). The unlocking reaction has a time constant of about 1 mus at 5 degrees C and is enthalpy-limited with an activation energy of 32 +/- 1 kJ/mol and a large Arrhenius preexponential factor of A = 7.5 x 10(11) s(-1). In the unlocked state a fast local conformational fluctuation with a time constant of 170 ns and a low activation barrier of 17 +/- 1 kJ/mol leads to unfolding of the C-terminal helix and to its undocking from the core. On a much slower time scale, global unfolding occurs from the unlocked state. These results suggest that native protein structures are locked into conformations with low amplitude motions. Large scale motions and global unfolding require an initial structural unlocking step leading to a state with properties of a dry molten globule. The experiments additionally yielded information on the dynamics of loop formation between different positions in unfolded HP35. Comparison of the results with dynamics in unstructured model peptides indicates slightly decelerated kinetics of local loop formation in the C-terminal region which points at residual, nonrandom structure. Dynamics of long-range loop formation, in contrast, are not influenced by residual structure, which argues against unfolded state properties as molecular origin for ultrafast folding of HP35.

Probing the dynamics of a protein hydrophobic core by deuteron solid-state nuclear magnetic resonance spectroscopy

With the goal of investigating dynamical features of hydrophobic cores of proteins over a wide range of temperatures, the chicken villin headpiece subdomain protein (HP36) was labeled at a "single" site corresponding to any one of the two C(delta)D(3) groups of leucine-69, which is located in a key position of the core. The main techniques employed are deuteron NMR quadrupolar echo line shape analysis, and T(1Z) (Zeeman) and T(1Q) (quadrupolar order) relaxation experiments performed at 11.7 and 17.6 T over the temperature range of 112 to 298 K. The experimental data are compared with computer simulations. The deuteron line shapes give an excellent fit to a three-mode motional model that consists of (a) fast three-site rotational jumps about the pseudo C(3) methyl spinning axis, (b) slower reorientation of the spinning axis, described by diffusion along a restricted arc, and (c) large angle jumps between traces of rotameric conformers. Relaxation behavior is described by a phenomenological distribution of activation energies for three-site hops at high temperatures that collapses to a single, distinctly smaller value for lower temperatures.

Native like structure in the unfolded state of the villin headpiece helical subdomain, an ultrafast folding protein

The villin headpiece subdomain, HP36, is the smallest naturally occurring protein that folds cooperatively. Its small size, rapid folding, and simple three-helix topology have made it an extremely popular system for computational studies of protein folding. The role of unfolded state structure in rapid folding is an area of active investigation, but relatively little is known about the properties of unfolded states under native conditions. A peptide fragment, HP21, which contains the first and second helices of HP36 has been shown to be a good model for structure in the unfolded state of the intact domain but a detailed description of the conformational propensities of HP21 is lacking and the balance between native and nonnative interactions is not known. A series of three-dimensional NMR experiments were performed on (13)C, (15)N-labeled HP21 to investigate in detail its conformational propensities. Analysis of (13)C(alpha), (13)C(beta), (13)CO chemical shifts, Deltadelta(13)C(alpha) - Deltadelta(13)C(beta) secondary shifts, the secondary structure propensity scores, NOEs, (15)N R(2) values and comparison of experimental chemical shifts with those of HP36 and with chemical shifts calculated using the SHIFTS and SHIFTX programs all indicate that there is significant native like structure in the HP21 ensemble, and thus by implication in the unfolded state of HP36.

Probing the folding transition state structure of the villin headpiece subdomain via side chain and backbone mutagenesis

Backbone-backbone hydrogen bonds are a common feature of native protein structures, yet their thermodynamic and kinetic influence on folding has long been debated. This is reflected by the disparity between current protein folding models, which place hydrogen bond formation at different stages along the folding trajectory. For example, previous studies have suggested that the denatured state of the villin headpiece subdomain contains a residual helical structure that may provide a bias toward the folded state by confining the conformational search associated with its folding. Although helical hydrogen bonds clearly stabilize the folded state, here we show, using an amide-to-ester mutation strategy, that the formation of backbone hydrogen bonds within helices is not rate-limiting in the folding of the subdomain, thereby suggesting that such hydrogen bonds are unlikely to be formed en route from the denatured to the transition state. On the other hand, elimination of hydrogen bonds within the turn region elicits a slower folding rate, consistent with the hypothesis that these residues are involved in the formation of a folding nucleus. While illustrating a potentially conserved aspect of helix-turn-helix folding, our results further underscore the inherent importance of turns in protein supersecondary structure formation.

Quantitative determination of site-specific conformational distributions in an unfolded protein by solid-state nuclear magnetic resonance

Solid-state nuclear magnetic resonance (NMR) techniques are used to investigate the structure of the 35-residue villin headpiece subdomain (HP35) in folded, partially denatured, and fully denatured states. Experiments are carried out in frozen glycerol/water solutions, with chemical denaturation by guanidine hydrochloride (GdnHCl). Without GdnHCl, two-dimensional solid-state (13)C NMR spectra of samples prepared with uniform (13)C labeling of selected residues show relatively sharp cross-peaks at chemical shifts that are consistent with the known three-helix bundle structure of folded HP35. At high GdnHCl concentrations, most cross-peaks broaden and shift, qualitatively indicating disruption of the folded structure and development of static conformational disorder in the frozen denatured state. Conformational distributions at one residue in each helical segment are probed quantitatively with three solid-state NMR techniques that provide independent constraints on backbone varphi and psi torsion angles in samples with sequential pairs of carbonyl (13)C labels. Without GdnHCl, the combined data are well fit by alpha-helical conformations. At [GdnHCl]=4.5 M, corresponding to the approximate denaturation midpoint, the combined data are well fit by a combination of alpha-helical and partially extended conformations at each site, but with a site-dependent population ratio. At [GdnHCl]=7.0 M, corresponding to the fully denatured state, the combined data are well fit by a combination of partially extended and polyproline II conformations, again with a site-dependent population ratio. Two entirely different models for conformational distributions lead to nearly the same best-fit distributions, demonstrating the robustness of these conclusions. This work represents the first quantitative investigation of site-specific conformational distributions in partially folded and unfolded states of a protein by solid-state NMR.

Probing the lower size limit for protein-like fold stability: ten-residue microproteins with specific, rigid structures in water

Mutational optimization of two long-range interactions first observed in Ac-WINGKWT-NH2, (a) bifurcated H-bonding involving the threonine amide H(N) and side chain OH and the N-terminal acetyl carbonyl and (b) an H-bond between the entgegen-H(N) of the C-terminal amide and the indole ring of Trp6 that stabilizes a face-to-edge indole/indole interaction between Trp1 and Trp6, has afforded < or = 10 residue systems that yield a remarkably stable fold in water. Optimization was achieved by designing a hydrophobic cluster that sequesters these H-bonds from solvent exposure. The structures and extent of amide H/D exchange protection for CH3CH2CO-WI pGXWTGPS (p = D-Pro, X = Leu or Ile) were determined. These two systems are greater than 94% folded at 298 K (97.5% at 280 K) with melting temperatures > 75 degrees C. The fold appears to display minimal fluxionality; a well-converged NMR structure rationalizes all of the large structuring shifts observed, and we suggest that these designed constructs can be viewed as microproteins.

Slow motions in chicken villin headpiece subdomain probed by cross-correlated NMR relaxation of amide NH bonds in successive residues

The villin headpiece subdomain (HP36) is a widely used system for protein-folding studies. Nuclear magnetic resonance cross-correlated relaxation rates arising from correlated fluctuations of two N-H(N) dipole-dipole interactions involving successive residues were measured at two temperatures at which HP36 is at least 99% folded. The experiment revealed the presence of motions slower than overall tumbling of the molecule. Based on the theoretical analysis of the spectral densities we show that the structural and dynamic contributions to the experimental cross-correlated relaxation rate can be separated under certain conditions. As a result, dynamic cross-correlated order parameters describing slow microsecond-to-millisecond motions of N-H bonds in neighboring residues can be introduced for any extent of correlations in the fluctuations of the two bond vectors. These dynamic cross-correlated order parameters have been extracted for HP36. The comparison of their values at two different temperatures indicates that when the temperature is raised, slow motions increase in amplitude. The increased amplitude of these fluctuations may reflect the presence of processes directly preceding the unfolding of the protein.

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.

Conformational analyses of a partially-folded bioactive prodomain of human furin

The 81-residue multifunctional prodomain of human furin adopts only a partially-folded conformational state under near physiological conditions. By use of NMR spectroscopy, we demonstrate that the N-terminal residues 1-46 of the prodomain in 50% trifluoroethanol (TFE) populates backbone conformations containing a short helix, a beta-strand and a helix-loop-helix super-secondary structure with elements of tertiary interactions. (15)N NMR relaxation measurements indicate that the helix-loop-helix region has similar motional characteristics in the fast picosecond to nanosecond timescales. On the other hand, the intervening segment (residues 47-65) is predominantly unstructured with a long and highly flexible region surrounding the protease 'activation loop' followed by a partially helical segment in the C-terminal end. Interestingly, the helix-loop-helix "fold" was found to be populated even when excised out of the full-length prodomain, since a peptide fragment derived from residues Pro16-Arg49 can also form the helix-loop-helix structure in aqueous solution in the absence of TFE. Structure analyses reveal that two helices orient in an antiparallel fashion directed by the sharing of hydrophobic residues involved in helix-capping interactions. Very importantly, a positively-charged Lys residue replacing His43 in the 16-49 fragment imparts stability to the super-secondary structure at both acidic and neutral pH, while a hydrophobic residue Leu at position 43 appears to destabilize the helical conformation in the 31-44 region. As such, this study provides valuable insights into the structural properties of the furin prodomain in relation to its role in the folding of the furin zymogen and its inhibitory action toward furin.

Protein folding under confinement: a role for solvent

Although most experimental and theoretical studies of protein folding involve proteins in vitro, the effects of spatial confinement may complicate protein folding in vivo. In this study, we examine the folding dynamics of villin (a small fast folding protein) with explicit solvent confined to an inert nanopore. We have calculated the probability of folding before unfolding (P(fold)) under various confinement regimes. Using P(fold) correlation techniques, we observed two competing effects. Confining protein alone promotes folding by destabilizing the unfolded state. In contrast, confining both protein and solvent gives rise to a solvent-mediated effect that destabilizes the native state. When both protein and solvent are confined we see unfolding to a compact unfolded state different from the unfolded state seen in bulk. Thus, we demonstrate that the confinement of solvent has a significant impact on protein kinetics and thermodynamics. We conclude with a discussion of the implications of these results for folding in confined environments such as the chaperonin cavity in vivo.

Reconciling the solution and X-ray structures of the villin headpiece helical subdomain: molecular dynamics simulations and double mutant cycles reveal a stabilizing cation-pi interaction

The 36-residue helical subdomain of the villin headpiece, HP36, is one of the smallest cooperatively folded proteins, folding on the microsecond time scale. The domain is an extraordinarily popular model system for both experimental and computational studies of protein folding. The structure of HP36 has been determined using X-ray crystallography and NMR spectroscopy, with the resulting structures exhibiting differences in helix packing, van der Waals contacts, and hydrogen bonding. It is important to determine the solution structure of HP36 with as much accuracy as possible since this structure is widely used as a reference for simulations and experiments. We complement the existing data by using all-atom molecular dynamics simulations with explicit solvent to evaluate which of the experimental models is the better representation of HP36 in solution. After simulation for 50 ns initiated with the NMR structure, we observed that the protein spontaneously adopts structures with a backbone conformation, core packing, and C-capping motif on the third helix that are more consistent with the crystal structure. We also examined hydrogen bonding and side chain packing interactions between D44 and R55 and between F47 and R55, respectively, which were observed in the crystal structure but not in the NMR-based solution structure. Simulations showed large fluctuations in the distance between D44 and R55, while the distance between F47 and R55 remained stable, suggesting the formation of a cation-pi interaction between those residues. Experimental double mutant cycles confirmed that the F47-R55 pair has a larger energetic coupling than the D44-R55 interaction. Overall, these combined experimental and computational studies show that the X-ray crystal structure is the better reference structure for HP36 in solution at neutral pH. Our analysis also shows how detailed molecular dynamics simulations combined with experimental validation can help bridge the gap between NMR and crystallographic methods.

Atomic-level characterization of disordered protein ensembles

The roles of unfolded states of proteins in normal folding and in diseases involving aggregation, as well as the prevalence and regulatory functions of intrinsically disordered proteins, have become increasingly recognized. The structural representation of these disordered states as ensembles of interconverting conformers can therefore provide critical insights. Experimental methods can be used to probe ensemble-averaged structural properties of disordered states and computational approaches generate representative ensembles of conformers using experimental restraints. In particular, NMR and small-angle X-ray scattering provide quantitative data that can readily be incorporated into calculations. These techniques have gleaned structural information about denatured, unfolded and intrinsically disordered proteins. The use of experimental data in different computational approaches, including ensemble molecular dynamics simulations and algorithms that assign populations to pregenerated conformers, has highlighted the presence of both local and long-range structure, and the occurrence of native-like and non-native interactions in unfolded and denatured states. Analysis of the resulting ensembles has suggested important implications of this fluctuating structure for folding, aggregation and binding.

Fluctuations of primary ubiquitin folding intermediates in a force clamp

Folding experiments of single ubiquitin molecules under force clamp using an atomic force microscope revealed a dynamic long-lived intermediate with nanometer scale end-to-end distance fluctuations along an unexpectedly complex folding pathway. To examine the nature of this intermediate at the atomic level as well as the driving forces that give rise to the observed fluctuations, we performed molecular dynamics refolding simulations of unfolded ubiquitin under constant force. After an initial fast collapse, we find a highly dynamic, broad ensemble of conformations with partial and continuously changing secondary structure and side chain interactions. This ensemble resembles a molten-globule-like state, similar in nature to the previously described non-native state of ubiquitin in solution, but stretched by the external force. The scale of the end-to-end distance fluctuations derived from the simulations compares well with experiment. Transient formation of unspecific and metastable hydrophobic clusters along the chain are found to give rise to the observed end-to-end distance fluctuations. These distinct collapses, interpreted as folding attempts, imply an upper limit for the folding attempt frequency of approximately 10 ns. Our results suggest possible relations between force-induced unfolding and temperature or chemically induced denaturation.