Non-native Contacts Research Articles

Structural Changes to Monomeric CuZn Superoxide Dismutase Caused by the Familial Amyotrophic Lateral Sclerosis-Associated Mutation A4V

Amyotrophic lateral sclerosis (ALS) is a progressive motor neuron degenerative disease, and the inherited form, familial ALS (fALS), has been linked to over 100 different point mutations scattered throughout the Cu-Zn superoxide dismutase protein (SOD1). The disease is likely due to a toxic gain of function caused by the misfolding, oligomerization, and eventual aggregation of mutant SOD1, but it is not yet understood how the structurally diverse mutations result in a common disease phenotype. The behavior of the apo-monomer fALS-associated mutant protein A4V was explored using molecular-dynamics simulations to elucidate characteristic structural changes to the protein that may allow the mutant form to improperly associate with other monomer subunits. Simulations showed that the mutant protein is less stable than the WT protein overall, with shifts in residue-residue contacts that lead to destabilization of the dimer and metal-binding sites, and stabilization of nonnative contacts that leads to a misfolded state. These findings provide a unifying explanation for disparate experimental observations, allow a better understanding of alterations of residue contacts that accompany loss of SOD1 structural integrity, and suggest sites where compensatory changes may stabilize the mutant structure.

The Limited Role of Nonnative Contacts in the Folding Pathways of a Lattice Protein

Models of protein energetics that neglect interactions between amino acids that are not adjacent in the native state, such as the Gō model, encode or underlie many influential ideas on protein folding. Implicit in this simplification is a crucial assumption that has never been critically evaluated in a broad context: Detailed mechanisms of protein folding are not biased by nonnative contacts, typically argued to be a consequence of sequence design and/or topology. Here we present, using computer simulations of a well-studied lattice heteropolymer model, the first systematic test of this oft-assumed correspondence over the statistically significant range of hundreds of thousands of amino acid sequences that fold to the same native structure. Contrary to previous conjectures, we find a multiplicity of folding mechanisms, suggesting that Gō-like models cannot be justified by considerations of topology alone. Instead, we find that the crucial factor in discriminating among topological pathways is the heterogeneity of native contact energies: The order in which native contacts accumulate is profoundly insensitive to omission of nonnative interactions, provided that native contact heterogeneity is retained. This robustness holds over a surprisingly wide range of folding rates for our designed sequences. Mirroring predictions based on the principle of minimum frustration, fast-folding sequences match their Gō-like counterparts in both topological mechanism and transit times. Less optimized sequences dwell much longer in the unfolded state and/or off-pathway intermediates than do Gō-like models. For dynamics that bridge unfolded and unfolded states, however, even slow folders exhibit topological mechanisms and transit times nearly identical with those of their Gō-like counterparts. Our results do not imply a direct correspondence between folding trajectories of Gō-like models and those of real proteins, but they do help to clarify key topological and energetic assumptions that are commonly used to justify such caricatures.

Is the Formation of the Correct Nucleating Loop the Rate-Limiting Step in Hairpin Formation in ss-Polynucleotides?

Numerous kinetics measurements on the formation of single-stranded (ss) DNA and RNA hairpin structures with ∼4-20 nucleotides (nt) in the loop and ∼5-8 base-pairs in the stem, indicate that the time required to form hairpins is ∼10-500 microseconds. If ss-polynucleotide chain is treated as an ideal semiflexible polymer with a statistical segment length of ∼4 nt, the theoretical estimate for the end-to-end contact time for an ∼10-nt long chain is expected to be tens-of-nanoseconds. To explain this discrepancy in time-scale, we proposed that the formation of the nucleating loop, prior to the zipping step, is slowed down as a result of transient trapping in misfolded conformations, with mis-paired base-pairs, non-native hydrogen bonding, or intrastrand stacking interactions in the unfolded state. Experimental measurements of end-to-end contact formation indicate that loop closure times for 4-nt poly(dT) loops are ∼400 ns, and for 4-nt poly(dA) loops are ∼8 microseconds, thus confirming that intrachain interactions slow down the configurational diffusion of the chain (Wang and Nau, J. Am. Chem. Soc. 2004, 126, 808). Interestingly, despite this evidence for intrachain interactions slowing down diffusion, the hairpin closing times for both ssDNA and RNA hairpins are found to scale with the length of the loop as L2.2-2.6, in reasonable agreement with the scaling behavior expected for loop-closure of a semiflexible polymer.Here, we present a kinetic zipper model that explicitly includes all misfolded microstates with non-native contacts, to describe the hairpin relaxation rates. The temperature and loop-size dependence for the relaxation rates is described in terms of two free parameters, the configurational diffusion coefficient that is relevant for the single-strand chain dynamics, and one parameter that characterizes the strength of non-native interactions prior to the formation of the nucleating loop.

Analysing the origin of long-range interactions in proteins using lattice models

BackgroundLong-range communication is very common in proteins but the physical basis of this phenomenon remains unclear. In order to gain insight into this problem, we decided to explore whether long-range interactions exist in lattice models of proteins. Lattice models of proteins have proven to capture some of the basic properties of real proteins and, thus, can be used for elucidating general principles of protein stability and folding.ResultsUsing a computational version of double-mutant cycle analysis, we show that long-range interactions emerge in lattice models even though they are not an input feature of them. The coupling energy of both short- and long-range pairwise interactions is found to become more positive (destabilizing) in a linear fashion with increasing 'contact-frequency', an entropic term that corresponds to the fraction of states in the conformational ensemble of the sequence in which the pair of residues is in contact. A mathematical derivation of the linear dependence of the coupling energy on 'contact-frequency' is provided.ConclusionOur work shows how 'contact-frequency' should be taken into account in attempts to stabilize proteins by introducing (or stabilizing) contacts in the native state and/or through 'negative design' of non-native contacts.

Aromatic Residues Engineered into the β-Turn Nucleation Site of Ubiquitin Lead to a Complex Folding Landscape, Non-Native Side-Chain Interactions, and Kinetic Traps

The fast folding of small proteins is likely to be the product of evolutionary pressures that balance the search for native-like contacts in the transition state with the minimum number of stable non-native interactions that could lead to partially folded states prone to aggregation and amyloid formation. We have investigated the effects of non-native interactions on the folding landscape of yeast ubiquitin by introducing aromatic substitutions into the beta-turn region of the N-terminal beta-hairpin, using both the native G-bulged type I turn sequence (TXTGK) as well as an engineered 2:2 XNGK type I' turn sequence. The N-terminal beta-hairpin is a recognized folding nucleation site in ubiquitin. The folding kinetics for wt-Ub (TLTGK) and the type I' turn mutant (TNGK) reveal only a weakly populated intermediate, however, substitution with X = Phe or Trp in either context results in a high propensity to form a stable compact intermediate where the initial U-->I collapse is visible as a distinct kinetic phase. The introduction of Trp into either of the two host turn sequences results in either complex multiphase kinetics with the possibility of parallel folding pathways, or formation of a highly compact I-state stabilized by non-native interactions that must unfold before refolding. Sequence substitutions with aromatic residues within a localized beta-turn capable of forming non-native hydrophobic contacts in both the native state and partially folded states has the undesirable consequence that folding is frustrated by the formation of stable compact intermediates that evolutionary pressures at the sequence level may have largely eliminated.

The Denatured State of N-PGK Is Compact and Predominantly Disordered

The organisation of the structure present in the chemically denatured N-terminal domain of phosphoglycerate kinase (N-PGK) has been determined by paramagnetic relaxation enhancements (PREs) to define the conformational landscape accessible to the domain. Below 2.0 M guanidine hydrochloride (GuHCl), a species of N-PGK (denoted I b) is detected, distinct from those previously characterised by kinetic experiments [folded (F), kinetic intermediate (I k) and denatured (D)]. The transition to I b is never completed at equilibrium, because F predominates below 1.0 M GuHCl. Therefore, the ability of PREs to report on transient or low population species has been exploited to characterise I b. Five single cysteine variants of N-PGK were labelled with the nitroxide electron spin-label MTSL [(1-oxyl-2,2,5,5-tetramethyl-3-pyrroline-3-methyl)methanesulfonate] and the denaturant dependences of the relaxation properties of the amide NMR signals between 1.2 and 3.6 M GuHCl were determined. Significant PREs for I b were obtained, but these were distributed almost uniformly throughout the sequence. Furthermore, the PREs indicate that no specific short tertiary contacts persist. The data indicate a collapsed state with no coherent three-dimensional structure, but with a restricted radius beyond which the protein chain rarely reaches. The NMR characteristics of I b indicate that it forms from the fully denatured state within 100 μs, and therefore a rapid collapse is the initial stage of folding of N-PGK from its chemically denatured state. By extrapolation, I b is the predominant form of the denatured state under native conditions, and the non-specifically collapsed structure implies that many non-native contacts and chain reversals form early in protein folding and must be broken prior to attaining the native state topology.

Theoretical and experimental demonstration of the importance of specific nonnative interactions in protein folding

Many experimental and theoretical studies have suggested a significant role for nonnative interactions in protein folding pathways, but the energetic contributions of these interactions are not well understood. We have addressed the energetics and the position specificity of nonnative hydrophobic interactions by developing a continuum coarse-grained chain model with a native-centric potential augmented by sequence-dependent hydrophobic interactions. By modeling the effect of different hydrophobicity values at various positions in the Fyn SH3 domain, we predicted energetically significant nonnative interactions that led to acceleration or deceleration of the folding rate depending on whether they were more populated in the transition state or unfolded state. These nonnative contacts were centered on position 53 in the Fyn SH3 domain, which lies in an exposed position in a 3(10)-helix. The energetic importance of the predicted nonnative interactions was confirmed experimentally by folding kinetics studies combined with double mutant thermodynamic cycles. By attaining agreement of theoretical and experimental investigations, this study provides a compelling demonstration that specific nonnative interactions can significantly influence folding energetics. Moreover, we show that a coarse-grained model with a simple consideration of hydrophobicity is sufficient for the accurate prediction of kinetically important nonnative interactions.

Stability of the Glycerol Facilitator in Detergent Solutions

Understanding membrane protein folding and stability is required for a molecular explanation of function and for the development of interventions in membrane protein folding diseases. Stable aqueous detergent solutions of the Escherichia coli glycerol facilitator in its native oligomeric state have been difficult to prepare as the protein readily unfolds and forms nonspecific aggregates. Here, we report a study of the structure and stability of the glycerol facilitator in several detergent solutions by Blue Native and sodium dodecyl sulfate (SDS)-polyacrylamide gel electrophoresis (PAGE), circular dichroism (CD), and fluorescence. Protein tetramers were prepared in neutral dodecyl maltoside (DDM) and in zwitterionic lysomyristoylphosphatidylcholine (LMPC) detergent solutions that are stable during SDS-PAGE. Thermal unfolding experiments show that the protein is more stable in LMPC than in DDM. Tertiary structure unfolds before quaternary and some secondary structure in LMPC, whereas unfolding is more cooperative in DDM. The high stability of the protein in DDM is evident from the unfolding half-life of 8 days in 8 M urea, suggesting that hydrophobic interactions contribute to the stability. The protein unfolds readily in LMPC below pH 6, whereas the tetramer remains intact at pH 4 in DDM. At pH 4 in DDM, the protein is more sensitive than at neutral pH to unfolding by SDS and the effect is reversible. At pH 3 in DDM, the tetramer unfolds, losing its tertiary structure but retaining native helical structure which melts at significantly lower temperatures than in the native tetramer. The glycerol facilitator prepared in SDS is mainly monomeric and has ~10% less alpha-helix than the native protein. CD suggests that it forms a condensed structure with non-native tertiary contacts highly similar to the state observed in LMPC at low pH. The implications of the results for in vitro and in vivo folding of the protein are discussed.

Special issue on protein folding: experimental and theoretical approaches

Special issue on protein folding: experimental and theoretical approaches

RNA misfolding and the action of chaperones

RNA folds to a myriad of three-dimensional structures and performs an equally diverse set of functions. The ability of RNA to fold and function in vivo is all the more remarkable because, in vitro, RNA has been shown to have a strong propensity to adopt misfolded, non-functional conformations. A principal factor underlying the dominance of RNA misfolding is that local RNA structure can be quite stable even in the absence of enforcing global tertiary structure. This property allows non-native structure to persist, and it also allows native structure to form and stabilize non-native contacts or non-native topology. In recent years it has become clear that one of the central reasons for the apparent disconnect between the capabilities of RNA in vivo and its in vitro folding properties is the presence of RNA chaperones, which facilitate conformational transitions of RNA and therefore mitigate the deleterious effects of RNA misfolding. Over the past two decades, it has been demonstrated that several classes of non-specific RNA binding proteins possess profound RNA chaperone activity in vitro and when overexpressed in vivo, and at least some of these proteins appear to function as chaperones in vivo. More recently, it has been shown that certain DExD/H-box proteins function as general chaperones to facilitate folding of group I and group II introns. These proteins are RNA-dependent ATPases and have RNA helicase activity, and are proposed to function by using energy from ATP binding and hydrolysis to disrupt RNA structure and/or to displace proteins from RNA-protein complexes. This review outlines experimental studies that have led to our current understanding of the range of misfolded RNA structures, the physical origins of RNA misfolding, and the functions and mechanisms of putative RNA chaperone proteins.

Conformational Changes below the Tm: Molecular Dynamics Studies of the Thermal Pretransition of Ribonuclease A

Recent work suggests that some native conformations of proteins can vary with temperature. To obtain an atomic-level description of this structural and conformational variation, we have performed all-atom, explicit-solvent molecular dynamics simulations of bovine pancreatic ribonuclease A (RNase A) up to its melting temperature (Tm approximately 337 K). RNase A has a thermal pretransition near 320 K [Stelea, S. D., Pancoska, P., Benight, A. S., and Keiderling, T. A. (2001) Protein Sci. 10, 970-978]. Our simulations identify a conformational change that coincides with this pretransition. Between 310 and 320 K, there is a small but significant decrease in the number of native contacts, beta-sheet hydrogen bonding, and deviation of backbone conformation from the starting structure, and an increase in the number of nonnative contacts. Native contacts are lost in beta-sheet regions and in alpha1, partially due to movement of alpha1 away from the beta-sheet core. At 330 and 340 K, a nonnative helical segment of residues 15-20 forms, corresponding to a helix observed in the N-terminal domain-swapped dimer [Liu, Y. S., Hart, P. J., Schulnegger, M. P., and Eisenberg, D. (1998) Proc. Natl. Acad. Sci. U.S.A. 95, 3437-3432]. The conformations observed at the higher temperatures possess nativelike topology and overall conformation, with many native contacts, but they have a disrupted active site. We propose that these conformations may represent the native state at elevated temperature, or the N' state. These simulations show that subtle, functionally important changes in protein conformation can occur below the Tm.

Phi-value analysis of a three-state protein folding pathway by NMR relaxation dispersion spectroscopy.

Experimental studies of protein folding frequently are consistent with two-state folding kinetics. However, recent NMR relaxation dispersion studies of several fast-folding mutants of the Fyn Src homology 3 (SH3) domain have established that folding proceeds through a low-populated on-pathway intermediate, which could not be detected with stopped-flow experiments. The dispersion experiments provide precise kinetic and thermodynamic parameters that describe the folding pathway, along with a detailed site-specific structural characterization of both the intermediate and unfolded states from the NMR chemical shifts that are extracted. Here we describe NMR relaxation dispersion Phi-value analysis of the A39V/N53P/V55L Fyn SH3 domain, where the effects of suitable point mutations on the energy landscape are quantified, providing additional insight into the structure of the folding intermediate along with per-residue structural information of both rate-limiting transition states that was not available from previous studies. In addition to the advantage of delineating the full three-state folding pathway, the use of NMR relaxation dispersion as opposed to stopped-flow kinetics to quantify Phi values facilitates their interpretation because the obtained chemical shifts monitor any potential structural changes along the folding pathway that might be introduced by mutation, a significant concern in their analysis. Phi-Value analysis of several point mutations of A39V/N53P/V55L Fyn SH3 establishes that the beta(3)-beta(4)-hairpin already is formed in the first transition state, whereas strand beta(1), which forms nonnative interactions in the intermediate, does not fully adopt its native conformation until after the final transition state. The results further support the notion that on-pathway intermediates can be stabilized by nonnative contacts.

Implementation of a k/k0 Method to Identify Long-Range Structure in Transition States during Conformational Folding/Unfolding of Proteins

A previously introduced kinetic-rate constant (k/k(0)) method, where k and k(0) are the folding (unfolding) rate constants in the mutant and the wild-type forms, respectively, of a protein, has been applied to obtain qualitative information about structure in the transition state ensemble (TSE) of bovine pancreatic ribonuclease A (RNase A), which contains four native disulfide bonds. The method compares the folding (unfolding) kinetics of RNase A, with and without a covalent crosslink and tests whether the crosslinked residues are associated in the folding (unfolding) transition state (TS) of the noncrosslinked version. To confirm that the fifth disulfide bond has not introduced a significant structural perturbation, we solved the crystal structure of the V43C-R85C mutant to 1.6 A resolution. Our findings suggest that residues Val43 and Arg85 are not associated, and that residues Ala4 and Val118 may form nonnative contacts, in the folding (unfolding) TSE of RNase A.

A study of Barstar folding events using boundary value simulations

From extensive biophysical studies of protein folding, two competing mechanisms emerged: hydrophobic collapse and the framework model. Our protein of choice is Barstar—a barnase inhibitor. The approximation algorithm we used to study Barstar folding trajectories is called SDEL—stochastic difference equation in length. Using the native structure as the final boundary value and a collection of unfolded structures as the varying initial boundary value, SDEL calculates an ensemble of least action pathways between these boundaries. The results are atomically detailed folding pathways, with as many intermediate structures as you request in the input. We generated 12 pathways, starting from a structurally wide selection of unfolded conformations. Using the protein's radius of gyration as our primary reaction coordinate, we tracked H-bonds, dihedral angles, native and non-native contacts, and energy along the folding pathways. This paper will follow our findings, with special emphasis on pinpointing hydrophobic collapse as a more appropriate mechanism for Barstar. Comparison with pathway predictions for Barstar using experimental techniques will also be discussed.

NMR characterizations of an amyloidogenic conformational ensemble of the PI3K SH3 domain

Amyloid formation is associated with structural changes of native polypeptides to monomeric intermediate states and their self-assembly into insoluble aggregates. Characterizations of the amyloidogenic intermediate state are, therefore, of great importance in understanding the early stage of amyloidogenesis. Here, we present NMR investigations of the structural and dynamic properties of the acid-unfolded amyloidogenic intermediate state of the phosphatidylinositol 3-kinase (PI3K) SH3 domain--a model peptide. The monomeric amyloidogenic state of the SH3 domain studied at pH 2.0 (35 degrees C) was shown to be substantially disordered with no secondary structural preferences. (15)N NMR relaxation experiments indicated that the unfolded polypeptide is highly flexible on a subnanosecond timescale when observed under the amyloidogenic condition (pH 2.0, 35 degrees C). However, more restricted motions were detected in residues located primarily in the beta-strands as well as in a loop in the native fold. In addition, nonnative long-range interactions were observed between the residues with the reduced flexibility by paramagnetic relaxation enhancement (PRE) experiments. These indicate that the acid-unfolded state of the SH3 domain adopts a partly folded conformation through nonnative long-range contacts between the dynamically restricted residues at the amyloid-forming condition.

Reversible Thermal Denaturation of a 60-kDa Genetically Engineered β-Sheet Polypeptide

A de novo 687-amino-acid residue polypeptide with a regular 32-amino-acid repeat sequence, (GA) 3GY(GA) 3GE(GA) 3GH(GA) 3GK, forms large β-sheet assemblages that exhibit remarkable folding properties and, as well, form fibrillar structures. This construct is an excellent tool to explore the details of β-sheet formation yielding intimate folding information that is otherwise difficult to obtain and may inform folding studies of naturally occurring materials. The polypeptide assumes a fully folded antiparallel β-sheet/turn structure at room temperature, and yet is completely and reversibly denatured at 125°C, adopting a predominant polyproline II conformation. Deep ultraviolet Raman spectroscopy indicated that melting/refolding occurred without any spectroscopically distinct intermediates, yet the relaxation kinetics depend on the initial polypeptide state, as would be indicative of a non-two-state process. Thermal denaturation and refolding on cooling appeared to be monoexponential with characteristic times of ∼1 and ∼60 min, respectively, indicating no detectable formation of hairpin-type nuclei in the millisecond timescale that could be attributed to nonlocal “nonnative” interactions. The polypeptide folding dynamics agree with a general property of β-sheet proteins, i.e., initial collapse precedes secondary structure formation. The observed folding is much faster than expected for a protein of this size and could be attributed to a less frustrated free-energy landscape funnel for folding. The polypeptide sequence suggests an important balance between the absence of strong nonnative contacts (salt bridges or hydrophobic collapse) and limited repulsion of charged side chains.

Identification of a Collapsed Intermediate with Non-native Long-range Interactions on the Folding Pathway of a Pair of Fyn SH3 Domain Mutants by NMR Relaxation Dispersion Spectroscopy

Recent 15N and 13C spin-relaxation dispersion studies of fast-folding mutants of the Fyn SH3 domain have established that folding proceeds through a low-populated on-pathway intermediate (I) where the central β-sheet is at least partially formed, but without interactions between the NH 2- and COOH-terminal β-strands that exist in the folded state (F). Initial studies focused on mutants where Gly48 is replaced; in an effort to establish whether this intermediate is a general feature of Fyn SH3 folding a series of 15N relaxation experiments monitoring the folding of Fyn SH3 mutants N53P/V55L and A39V/N53P/V55L are reported here. For these mutants as well, folding proceeds through an on-pathway intermediate with similar features to those observed for G48M and G48V Fyn SH3 domains. However, the 15N chemical shifts extracted for the intermediate indicate pronounced non-native contacts between the NH 2 and COOH-terminal regions not observed previously. The kinetic parameters extracted for the folding of A39V/N53P/V55L Fyn SH3 from the three-state folding model F ↔ I ↔ U are in good agreement with folding and unfolding rates extrapolated to zero denaturant obtained from stopped-flow experiments analyzed in terms of a simplified two-state folding reaction. The folding of the triple mutant was studied over a wide range of temperatures, establishing that there is no difference in heat capacities between F and I states. This confirms a compact folding intermediate structure, which is supported by the 15N chemical shifts of the I state extracted from the dispersion data. The temperature-dependent relaxation data simplifies data analysis because at low temperatures (<25 °C) the unfolded state (U) is negligibly populated relative to I and F. A comparison between parameters extracted at low temperatures where the F ↔ I exchange model is appropriate with those from the more complex, three-state model at higher temperatures has been used to validate the protocol for analysis of three-site exchange relaxation data.

Analyses of simulations of three-dimensional lattice proteins in comparison with a simplified statistical mechanical model of protein folding

Folding and unfolding simulations of three-dimensional lattice proteins were analyzed using a simplified statistical mechanical model in which their amino acid sequences and native conformations were incorporated explicitly. Using this statistical mechanical model, under the assumption that only interactions between amino acid residues within a local structure in a native state are considered, the partition function of the system can be calculated for a given native conformation without any adjustable parameter. The simulations were carried out for two different native conformations, for each of which two foldable amino acid sequences were considered. The native and non-native contacts between amino acid residues occurring in the simulations were examined in detail and compared with the results derived from the theoretical model. The equilibrium thermodynamic quantities (free energy, enthalpy, entropy, and the probability of each amino acid residue being in the native state) at various temperatures obtained from the simulations and the theoretical model were also examined in order to characterize the folding processes that depend on the native conformations and the amino acid sequences. Finally, the free energy landscapes were discussed based on these analyses.

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

The villin headpiece subdomain (HP36) is the smallest naturally occurring protein that folds cooperatively. The protein folds on a microsecond time scale. Its small size and very rapid folding have made it a popular target for biophysical studies of protein folding. Temperature-dependent one-dimensional (1D) NMR studies of the full-length protein together with CD and 1D NMR studies of the 21-residue peptide fragment (HP21) derived from HP36 have shown that there is significant structure in the unfolded state of HP36 and have demonstrated that HP21 is a good model of these interactions. Here, we characterized the model peptide HP21 in detail by two-dimensional NMR. Strongly upfield shifted C(alpha) protons, the magnitude of the 3J(NH,alpha) coupling constants, and the pattern of backbone-backbone and backbone-side chain NOEs indicate that the ensemble of structures populated by HP21 contains alpha-helical structure and native as well as non-native hydrophobic contacts. The hydrogen-bonded secondary structure inferred from the NOEs is, however, not sufficient to confer significant protection against amide H-D exchange. These studies indicate that there is significant secondary structure and hydrophobic clustering in the unfolded state of HP36. The implications for the folding of HP36 are discussed.

Engineering Diverse Changes in β-Turn Propensities in the N-Terminal β-Hairpin of Ubiquitin Reveals Significant Effects on Stability and Kinetics but a Robust Folding Transition State

Using the N-terminal 17-residue beta-hairpin of ubiquitin as a "host" for mutational studies, we have investigated the influence of the beta-turn sequence on protein stability and folding kinetics by replacing the native G-bulged turn (TLTGK) with more flexible analogues (TG3K and TG5K) and a series of four-residue type I' beta-turn sequences, commonly found in beta-hairpins. Although a statistical analysis of type I' turns demonstrates residue preferences at specific sites, the frequency of occurrence appears to only broadly correlate with experimentally determined protein stabilities. The subsequent engineering of context-dependent non-native tertiary contacts involving turn residues is shown to produce large changes in stability. Relatively few point mutations have been described that probe secondary structure formation in ubiquitin in a manner that is independent of tertiary contacts. To this end, we have used the more rigorous rate-equilibrium free energy relationship (Leffler analysis), rather than the two-point phi value analysis, to show for a family of engineered beta-turn mutants that stability (range of approximately 20 kJ/mol) and folding kinetics (190-fold variation in refolding rate) are linearly correlated (alpha(f) = 0.74 +/- 0.08). The data are consistent with a transition state that is robust with regard to a wide range of statistically favored and disfavored beta-turn mutations and implicate a loosely assembled beta-hairpin as a key template in transition state stabilization with the beta-turn playing a central role.