Protein folding kinetics beyond the phi value: using multiple amino acid substitutions to investigate the structure of the SH3 domain folding transition state

Mega-scale experimental analysis of protein folding stability in biology and design

Advances in DNA sequencing and machine learning are providing insights into protein sequences and structures on an enormous scale1. However, the energetics driving folding are invisible in these structures and remain largely unknown2. The hidden thermodynamics of folding can drive disease3,4, shape protein evolution5-7 and guide protein engineering8-10, and new approaches are needed to reveal these thermodynamics for every sequence and structure. Here we present cDNA display proteolysis, a method for measuring thermodynamic folding stability for up to 900,000 protein domains in a one-week experiment. From 1.8 million measurements in total, we curated a set of around 776,000 high-quality folding stabilities covering all single amino acid variants and selected double mutants of 331 natural and 148 de novo designed protein domains 40-72 amino acids in length. Using this extensive dataset, we quantified (1) environmental factors influencing amino acid fitness, (2) thermodynamic couplings (including unexpected interactions) between protein sites, and (3) the global divergence between evolutionary amino acid usage and protein folding stability. We also examined how our approach could identify stability determinants in designed proteins and evaluate design methods. The cDNA display proteolysis method is fast, accurate and uniquely scalable, and promises to reveal the quantitative rules for how amino acid sequences encode folding stability.

Mutational Analysis of Protein Folding Transition States: Phi Values

The analysis of protein folding reactions by monitoring the kinetic effects of specifically designed single-point mutations, the so-termed phi-value analysis, has been a favorite technique to experimentally probe the mechanisms of protein folding. The idea behind phi-value analysis is that the effects that mutations have on the folding and unfolding rate constants report on the energetic/structural features of the folding transition state ensemble (TSE), which is the highest point in the free energy surface connecting the native and unfolded states, and thus the rate limiting step that ultimately defines the folding mechanism. For single-domain, two-state folding proteins, the general procedure to perform the phi-value analysis of protein folding is relatively simple to implement in the lab. Once the mutations have been produced and purified, the researcher needs to follow a few specific guidelines to perform the experiments and to analyze the data so produced. In this chapter, a step-by-step description of how to measure and interpret the effects induced by site-directed mutations on the folding and unfolding rate constants of a protein of interest is provided. Some possible solutions to the most typical problems that arise when performing phi-value analysis in the lab are also provided.

Effects of Topology and Sequence in Protein Folding Linked via Conformational Fluctuations

Experiments have compared the folding of proteins with different amino acid sequences but the same basic structure, or fold. Results indicate that folding is robust to sequence variations for proteins with some nonlocal folds, such as all-β, whereas the folding of more local, all-α proteins typically exhibits a stronger sequence dependence. Here, we use a coarse-grained model to systematically study how variations in sequence perturb the folding energy landscapes of three model sequences with 3α, 4β + α, and β-barrel folds, respectively. These three proteins exhibit folding features in line with experiments, including expected rank order in the cooperativity of the folding transition and stability-dependent shifts in the location of the free-energy barrier to folding. Using a generalized-ensemble simulation approach, we determine the thermodynamics of around 2000 sequence variants representing all possible hydrophobic or polar single- and double-point mutations. From an analysis of the subset of stability-neutral mutations, we find that folding is perturbed in a topology-dependent manner, with the β-barrel protein being the most robust. Our analysis shows, in particular, that the magnitude of mutational perturbations of the transition state is controlled in part by the size or "width" of the underlying conformational ensemble. This result suggests that the mutational robustness of the folding of the β-barrel protein is underpinned by its conformationally restricted transition state ensemble, revealing a link between sequence and topological effects in protein folding.

Protein folding while chaperone bound is dependent on weak interactions

It is generally assumed that protein clients fold following their release from chaperones instead of folding while remaining chaperone-bound, in part because binding is assumed to constrain the mobility of bound clients. Previously, we made the surprising observation that the ATP-independent chaperone Spy allows its client protein Im7 to fold into the native state while continuously bound to the chaperone. Spy apparently permits sufficient client mobility to allow folding to occur while chaperone bound. Here, we show that strengthening the interaction between Spy and a recently discovered client SH3 strongly inhibits the ability of the client to fold while chaperone bound. The more tightly Spy binds to its client, the more it slows the folding rate of the bound client. Efficient chaperone-mediated folding while bound appears to represent an evolutionary balance between interactions of sufficient strength to mediate folding and interactions that are too tight, which tend to inhibit folding.

Spy is an ATP independent chaperone that allows folding of its client protein Im7 while continuously bound to Spy. Here the authors employ kinetics measurements to study the folding of another Spy client protein SH3 and find that Spy’s ability to allow a client to fold while bound is inversely related to how strongly it interacts with that client.

High-Resolution Mapping of the Folding Transition State of a WW Domain

Fast-folding WW domains are among the best-characterized systems for comparing experiments and simulations of protein folding. Recent microsecond-resolution experiments and long duration (totaling milliseconds) single-trajectory modeling have shown that even mechanistic changes in folding kinetics due to mutation can now be analyzed. Thus, a comprehensive set of experimental data would be helpful to benchmark the predictions made by simulations. Here, we use T-jump relaxation in conjunction with protein engineering and report mutational Φ-values (Φ(M)) as indicators for folding transition-state structure of 65 side chain, 7 backbone hydrogen bond, and 6 deletion and /or insertion mutants within loop 1 of the 34-residue hPin1 WW domain. Forty-five cross-validated consensus mutants could be identified that provide structural constraints for transition-state structure within all substructures of the WW domain fold (hydrophobic core, loop 1, loop 2, β-sheet). We probe the robustness of the two hydrophobic clusters in the folding transition state, discuss how local backbone disorder in the native-state can lead to non-classical Φ(M)-values (Φ(M) > 1) in the rate-determining loop 1 substructure, and conclusively identify mutations and positions along the sequence that perturb the folding mechanism from loop 1-limited toward loop 2-limited folding.

Assessing the effect of dynamics on the closed-loop protein-folding hypothesis

The closed-loop (loop-n-lock) hypothesis of protein folding suggests that loops of about 25 residues, closed through interactions between the loop ends (locks), play an important role in protein structure. Coarse-grain elastic network simulations, and examination of loop lengths in a diverse set of proteins, each supports a bias towards loops of close to 25 residues in length between residues of high stability. Previous studies have established a correlation between total contact distance (TCD), a metric of sequence distances between contacting residues (cf. contact order), and the log-folding rate of a protein. In a set of 43 proteins, we identify an improved correlation (r² = 0.76), when the metric is restricted to residues contacting the locks, compared to the equivalent result when all residues are considered (r² = 0.65). This provides qualified support for the hypothesis, albeit with an increased emphasis upon the importance of a much larger set of residues surrounding the locks. Evidence of a similar-sized protein core/extended nucleus (with significant overlap) was obtained from TCD calculations in which residues were successively eliminated according to their hydrophobicity and connectivity, and from molecular dynamics simulations. Our results suggest that while folding is determined by a subset of residues that can be predicted by application of the closed-loop hypothesis, the original hypothesis is too simplistic; efficient protein folding is dependent on a considerably larger subset of residues than those involved in lock formation.

Testing the diffusing boundary model for the helix-coil transition in peptides

The dynamics of peptide α-helices have been studied extensively for many years, and the kinetic mechanism of the helix-coil dynamics has been discussed controversially. Recent experimental results have suggested that equilibrium helix-coil dynamics are governed by movement of the helix/coil boundary along the peptide chain, which leads to slower unfolding kinetics in the helix center compared with the helix ends and position-independent helix formation kinetics. We tested this diffusion of boundary model in helical peptides of different lengths by triplet-triplet energy transfer measurements and compared the data with simulations based on a kinetic linear Ising model. The results show that boundary diffusion in helical peptides can be described by a classical, Einstein-type, 1D diffusion process with a diffusion coefficient of 2.7⋅10(7) (amino acids)(2)/s or 6.1⋅10(-9) cm(2)/s. In helices with a length longer than about 40 aa, helix unfolding by coil nucleation in a helical region occurs frequently in addition to boundary diffusion. Boundary diffusion is slowed down by helix-stabilizing capping motifs at the helix ends in agreement with predictions from the kinetic linear Ising model. We further tested local and nonlocal effects of amino acid replacements on helix-coil dynamics. Single amino acid replacements locally affect folding and unfolding dynamics with a ϕf-value of 0.35, which shows that interactions leading to different helix propensities for different amino acids are already partially present in the transition state for helix formation. Nonlocal effects of amino acid replacements only influence helix unfolding (ϕf = 0) in agreement with a diffusing boundary mechanism.

A "Link-Psi" strategy using crosslinking indicates that the folding transition state of ubiquitin is not very malleable

Using a combined crosslinking-ψ analysis strategy, we examine whether the structural content of the transition state of ubiquitin can be altered. A synthetic dichloroacetone crosslink is first introduced across two β strands. Whether the structural content in the transition state ensemble has shifted towards the region containing the crosslink is probed by remeasuring the ψ value at another region (ψ identifies the degree to which an inserted bi-Histidine metal ion binding site is formed in the transition state). For sites around the periphery of the obligate transition state nucleus, we find that the resulting changes in ψ values are near or at our detection limit, thereby indicating that the structural content of the transition state has not measurably changed upon crosslinking. This work demonstrates the utility of the simultaneous application of crosslinking and ψ-analysis for examining potential transition state heterogeneity in globular proteins.

A residue in helical conformation in the native state adopts a β-strand conformation in the folding transition state despite its high and canonical Φ-value

Delineating structures of the transition states in protein folding reactions has provided great insight into the mechanisms by which proteins fold. The most common method for obtaining this information is Φ-value analysis, which is carried out by measuring the changes in the folding and unfolding rates caused by single amino acid substitutions at various positions within a given protein. Canonical Φ-values range between 0 and 1, and residues displaying high values within this range are interpreted to be important in stabilizing the transition state structure, and to elicit this stabilization through native-like interactions. Although very successful in defining the general features of transition state structures, Φ-value analysis can be confounded when non-native interactions stabilize this state. In addition, direct information on backbone conformation within the transition state is not provided. In the work described here, we have investigated structure formation at a conserved β-bulge (with helical conformation) in the Fyn SH3 domain by characterizing the effects of substituting all natural amino acids at one position within this structural motif. By comparing the effects on folding rates of these substitutions with database-derived local structure propensity values, we have determined that this position adopts a non-native backbone conformation in the folding transition state. This result is surprising because this position displays a high and canonical Φ-value of 0.7. This work emphasizes the potential role of non-native conformations in folding pathways and demonstrates that even positions displaying high and canonical Φ-values may, nevertheless, adopt a non-native conformation in the transition state.

Differential ordering of the protein backbone and side chains during protein folding revealed by site-specific recombinant infrared probes

The time scale for ordering of the polypeptide backbone relative to the side chains is a critical issue in protein folding. The interplay between ordering of the backbone and ordering of the side chains is particularly important for the formation of β-sheet structures, as the polypeptide chain searches for the native stabilizing cross-strand interactions. We have studied these issues in the N-terminal domain of protein L9 (NTL9), a model protein with mixed α/β structure. We have developed a general approach for introducing site-specific IR probes for the side chains (azide) and backbone ((13)C═(18)O) using recombinant protein expression. Temperature-jump time-resolved IR spectroscopy combined with site-specific labeling enables independent measurement of the respective backbone and side-chain dynamics with single residue resolution. We have found that side-chain ordering in a key region of the β-sheet structure occurs on a slower time scale than ordering of the backbone during the folding of NTL9, likely as a result of the transient formation of non-native side-chain interactions.

Dry molten globule intermediates and the mechanism of protein unfolding

New experimental results show that either gain or loss of close packing can be observed as a discrete step in protein folding or unfolding reactions. This finding poses a significant challenge to the conventional two-state model of protein folding. Results of interest involve dry molten globule (DMG) intermediates, an expanded form of the protein that lacks appreciable solvent. When an unfolding protein expands to the DMG state, side chains unlock and gain conformational entropy, while liquid-like van der Waals interactions persist. Four unrelated proteins are now known to form DMGs as the first step of unfolding, suggesting that such an intermediate may well be commonplace in both folding and unfolding. Data from the literature show that peptide amide protons are protected in the DMG, indicating that backbone structure is intact despite loss of side-chain close packing. Other complementary evidence shows that secondary structure formation provides a major source of compaction during folding. In our model, the major free-energy barrier separating unfolded from native states usually occurs during the transition between the unfolded state and the DMG. The absence of close packing at this barrier provides an explanation for why phi-values, derived from a Brønsted-Leffler plot, depend primarily on structure at the mutational site and not on specific side-chain interactions. The conventional two-state folding model breaks down when there are DMG intermediates, a realization that has major implications for future experimental work on the mechanism of protein folding.

Energetics and mechanisms of folding and flipping the myristoyl switch in the {beta}-trefoil protein, hisactophilin

Myristoylation, the covalent linkage of a saturated, C(14) fatty acyl chain to the N-terminal glycine in a protein, plays a vital role in reversible membrane binding and signaling by the modified proteins. Currently, little is known about the effects of myristoylation on protein folding and stability, or about the energetics and molecular mechanisms of switching involving states with sequestered versus accessible myristoyl group. Our analysis of these effects in hisactophilin, a histidine-rich protein that binds cell membranes and actin in a pH-dependent manner, shows that myristoylation significantly increases hisactophilin stability, while also markedly increasing global protein folding and unfolding rates. The switching between sequestered and accessible states is pH dependent, with an apparent pK(switch) of 6.95, and an apparent free energy change of 2.0 kcal·mol(-1). The myristoyl switch is linked to the reversible uptake of ∼1.5 protons, likely by histidine residues. This pH dependence of switching appears to be the physical basis of the sensitive, pH-dependent regulation of membrane binding observed in vivo. We conclude that an increase in protein stability upon modification and burial of the attached group is likely to occur in numerous proteins modified with fatty acyl or other hydrophobic groups, and that the biophysical effects of such modification are likely to play an important role in their functional switches. In addition, the increased global dynamics caused by myristoylation of hisactophilin reveals a general mechanism whereby hydrophobic moieties can make nonnative interactions or relieve strain in transition states, thereby increasing the rates of interconversion between different states.

Investigation of an anomalously accelerating substitution in the folding of a prototypical two-state protein

The folding rates of two-state single-domain proteins are generally resistant to small-scale changes in amino acid sequence. For example, having surveyed here over 700 single-residue substitutions in 24 well-characterized two-state proteins, we find that the majority (55%) of these substitutions affect folding rates by less than a factor of 2, and that only 9% affect folding rates by more than a factor of 8. Among those substitutions that significantly affect folding rates, we find that accelerating substitutions are an order of magnitude less common than those that decelerate the process. One of the most extreme outliers in this data set, an arginine-to-phenylalanine substitution at position 48 (R48F) of chymotrypsin inhibitor 2 (CI2), accelerates the protein's folding rate by a factor of 36 relative to that of the wild-type protein and is the most accelerating substitution reported to date in a two-state protein. In order to better understand the origins of this anomalous behavior, we have characterized the kinetics of multiple additional substitutions at this position. We find that substitutions at position 48 in CI2 fall into two distinct classes. The first, comprising residues that ablate the charge of the wild-type arginine but retain the hydrophobicity of its alkane chain, accelerate folding by at least 10-fold. The second class, comprising all other residues, produces folding rates within a factor of two of the wild-type rate. A significant positive correlation between hydrophobicity and folding rate across all of the residues we have characterized at this position suggests that the hydrophobic methylene units of the wild-type arginine play a significant role in stabilizing the folding transition state. Likewise, studies of the pH dependence of the histidine substitution indicate a strong correlation between folding rate and charge state. Thus, mutations that ablate the arginine's positive charge while retaining the hydrophobic contacts of its methylene units tend to dramatically accelerate folding. Previous studies have suggested that arginine 48 plays an important functional role in CI2, which may explain why it is highly conserved despite the anomalously large deceleration it produces in the folding of this protein.

What lessons can be learned from studying the folding of homologous proteins?

The studies of the folding of structurally related proteins have proved to be a very important tool for investigating protein folding. Here we review some of the insights that have been gained from such studies. Our highlighted studies show just how such an investigation should be designed and emphasise the importance of the synergy between experiment and theory. We also stress the importance of choosing the right system carefully, exploiting the excellent structural and sequence databases at our disposal.

Mutational investigation of protein folding transition states by Phi-value analysis and beyond: lessons from SH3 domain folding

Understanding how proteins adopt their unique native structures requires a complete structural characterization of the rate-limiting transition state(s) along the folding pathway. By definition, transition states are not significantly populated and are only accessible via folding kinetics studies. In this respect, interpreting the kinetic effects of amino acid substitutions (especially to Ala) via Phi-value analysis is the most common method to probe the structure of these transient, yet important states. A critical review of the key assumptions required for rigorous interpretation of Phi values reveals that a multiple substitution strategy in which a position of interest is mutated to a variety of amino acids, and not exclusively to Ala, provides the best means to characterize folding transition states. This approach has proven useful in revealing non-native interactions and (or) conformations in folding transition states. Moreover, by simultaneously examining the folding kinetics of multiple substitutions made at a single surface-exposed position using the Brønsted analysis the backbone conformation in a folding transition state can be investigated. For folding equilibria with exchange rates on the order of milliseconds, the kinetic parameters for Phi-value analysis can be obtained from NMR relaxation dispersion experiments, under fully native conditions, along with a wealth of high-resolution structural information about the states in exchange (native, denatured, and intermediate states that populate the pathway). This additional structural information, which is not readily obtained through stopped-flow based methods, can significantly facilitate the interpretation of Phi values because it often reports on the validity of the assumptions required for a rigorous interpretation of Phi values.

The osmolyte trimethylamine-N-oxide stabilizes the Fyn SH3 domain without altering the structure of its folding transition state

Trimethylamine-N-oxide (TMAO) is a naturally occurring osmolyte that stabilizes proteins against denaturation. Although the impact of TMAO on the folding thermodynamics of many proteins has been well characterized, far fewer studies have investigated its effects on protein folding kinetics. In particular, no previous studies have used Phi-value analysis to determine whether TMAO may alter the structure of the folding transition state. Here we have measured the effects on folding kinetics of 16 different amino acid substitutions distributed across the structure of the Fyn SH3 domain both in the presence and absence of TMAO. The folding and unfolding rates in TMAO, on average, improved to equivalent degrees, with a twofold increase in the protein folding rate accompanied by a twofold decrease in the unfolding rate. Importantly, TMAO caused little alteration to the Phi-values of the mutants tested, implying that this compound minimally perturbs the folding transition state structure. Furthermore, the solvent accessibility of the transition state was not altered as reflected in an absence of a TMAO-induced change in the denaturant beta(T) (D) factors. Through TMAO-induced folding studies, a beta(T) (TMAO) factor of 0.5 was calculated for this compound, suggesting that the protein backbone, which is the target of action of TMAO, is 50% exposed in the transition state as compared to the native state. This finding is consistent with the equivalent effects of TMAO on the folding and unfolding rates. Through thermodynamic analysis of mutants, we also discovered that the stabilizing effect of TMAO is lessened with increasing temperature.

Quantifying the structural requirements of the folding transition state of protein A and other systems

The B-domain of protein A is a small three-helix bundle that has been the subject of considerable experimental and theoretical investigation. Nevertheless, a unified view of the structure of the transition-state ensemble (TSE) is still lacking. To characterize the TSE of this surprisingly challenging protein, we apply a combination of psi analysis (which probes the role of specific side-chain to side-chain contacts) and kinetic H/D amide isotope effects (which measures hydrogen-bond content), building upon previous studies using mutational phi analysis (which probes the energetic influence of side-chain substitutions). The second helix is folded in the TSE, while helix formation appears just at the carboxy and amino termini of the first and third helices, respectively. The experimental data suggest a homogenous yet plastic TS with a native-like topology. This study generalizes our earlier conclusion, based on two larger alpha/beta proteins, that the TSEs of most small proteins achieve approximately 70% of their native state's relative contact order. This high percentage limits the degree of possible TS heterogeneity and requires a reevaluation of the structural content of the TSE of other proteins, especially when they are characterized as small or polarized.

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.

Side chain burial and hydrophobic core packing in protein folding transition states

A critical step in the folding pathway of globular proteins is the formation of a tightly packed hydrophobic core. Several mutational studies have addressed the question of whether tight packing interactions are present during the rate-limiting step of folding. In some of these investigations, substituted side chains have been assumed to form native-like interactions in the transition state when the folding rates of mutant proteins correlate with their native-state stabilities. Alternatively, it has been argued that side chains participate in nonspecific hydrophobic collapse when the folding rates of mutant proteins correlate with side-chain hydrophobicity. In a reanalysis of published data, we have found that folding rates often correlate similarly well, or poorly, with both native-state stability and side-chain hydrophobicity, and it is therefore not possible to select an appropriate transition state model based on these one-parameter correlations. We show that this ambiguity can be resolved using a two-parameter model in which side chain burial and the formation of all other native-like interactions can occur asynchronously. Notably, the model agrees well with experimental data, even for positions where the one-parameter correlations are poor. We find that many side chains experience a previously unrecognized type of transition state environment in which specific, native-like interactions are formed, but hydrophobic burial dominates. Implications of these results to the design and analysis of protein folding studies are discussed.

Computational design of the Fyn SH3 domain with increased stability through optimization of surface charge charge interactions

Computational design of surface charge-charge interactions has been demonstrated to be an effective way to increase both the thermostability and the stability of proteins. To test the robustness of this approach for proteins with predominantly beta-sheet secondary structure, the chicken isoform of the Fyn SH3 domain was used as a model system. Computational analysis of the optimal distribution of surface charges showed that the increase in favorable energy per substitution begins to level off at five substitutions; hence, the designed Fyn sequence contained four charge reversals at existing charged positions and one introduction of a new charge. Three additional variants were also constructed to explore stepwise contributions of these substitutions to Fyn stability. The thermodynamic stabilities of the variants were experimentally characterized using differential scanning calorimetry and far-UV circular dichroism spectroscopy and are in very good agreement with theoretical predictions from the model. The designed sequence was found to have increased the melting temperature, DeltaT (m) = 12.3 +/- 0.2 degrees C, and stability, DeltaDeltaG(25 degrees C) = 7.1 +/- 2.2 kJ/mol, relative to the wild-type protein. The experimental data suggest that a significant increase in stability can be achieved through a very small number of amino acid substitutions. Consistent with a number of recent studies, the presented results clearly argue for a seminal role of surface charge-charge interactions in determining protein stability and suggest that the optimization of surface interactions can be an attractive strategy to complement algorithms optimizing interactions in the protein core to further enhance protein stability.

Plasticity within the obligatory folding nucleus of an immunoglobulin-like domain

A number of β-sandwich immunoglobulin-like domains have been shown to fold using a set of structurally equivalent residues that form a folding nucleus deep within the core of the protein. Formation of this nucleus is sufficient to establish the complex Greek key topology of the native state. These nucleating residues are highly conserved within the immunoglobulin superfamily, but are less well conserved in the fibronectin type III (fnIII) superfamily, where the requirement is simply to have four interacting hydrophobic residues. However, there are rare examples where this nucleation pattern is absent. In this study, we have investigated the folding of a novel member of the fnIII superfamily whose nucleus appears to lack one of the four buried hydrophobic residues. We show that the folding mechanism is unaltered, but the folding nucleus has moved within the hydrophobic core.

Protein folding kinetics provides a context-independent assessment of beta-strand propensity in the Fyn SH3 domain

Structural database-derived propensities for amino acids to adopt particular local protein structures, such as alpha-helix and beta-strand, have long been recognized and effectively exploited for the prediction of protein secondary structure. However, the experimental verification of database-derived propensities using mutagenesis studies has been problematic, especially for beta-strand propensities, because local structural preferences are often confounded by non-local interactions arising from formation of the native tertiary structure. Thus, the overall thermodynamic stability of a protein is not always altered in a predictable manner by changes in local structural propensity at a single position. In this study, we have undertaken an investigation of the relationship between beta-strand propensity and protein folding kinetics. By characterizing the effects of a wide variety of amino acid substitutions at two different beta-strand positions in an SH3 domain, we have found that the observed changes in protein folding rates are very well correlated to beta-strand propensities for almost all of the substitutions examined. In contrast, there is little correlation between propensities and unfolding rates. These data indicate that beta-strand conformation is well formed in the structured portion of the SH3 domain transition state, and that local structure propensity strongly influences the stability of the transition state. Since the transition state is known to be packed more loosely than the native state and likely lacks many of the non-local stabilizing interactions seen in the native state, we suggest that folding kinetics studies may generally provide an effective means for the experimental validation of database-derived local structural propensities.

The folding pathway of an FF domain: characterization of an on-pathway intermediate state under folding conditions by (15)N, (13)C(alpha) and (13)C-methyl relaxation dispersion and (1)H/(2)H-exchange NMR spectroscopy

The FF domain from the human protein HYPA/FBP11 folds via a low-energy on-pathway intermediate (I). Elucidation of the structure of such folding intermediates and denatured states under conditions that favour folding are difficult tasks. Here, we investigated the millisecond time-scale equilibrium folding transition of the 71-residue four-helix bundle wild-type protein by (15)N, (13)C(alpha) and methyl(13)C Carr-Purcell-Meiboom-Gill (CPMG) NMR relaxation dispersion experiments and by (1)H/(2)H-exchange measurements. The relaxation data for the wild-type protein fitted a simple two-site exchange process between the folded state (F) and I. Destabilization of F in mutants A17G and Q19G allowed the detection of the unfolded state U by (15)N CPMG relaxation dispersion. The dispersion data for these mutants fitted a three-site exchange scheme, U<-->I<-->F, with I populated higher than U. The kinetics and thermodynamics of the folding reaction were obtained via temperature and urea-dependent relaxation dispersion experiments, along with structural information on I from backbone (15)N, (13)C(alpha) and side-chain methyl (13)C chemical shifts, with further information from protection factors for the backbone amide groups from (1)H/(2)H-exchange. Notably, helices H1-H3 are at least partially formed in I, while helix H4 is largely disordered. Chemical shift differences for the methyl (13)C nuclei suggest a paucity of stable, native-like hydrophobic interactions in I. These data are consistent with Phi-analysis of the rate-limiting transition state between I and F. The combination of relaxation dispersion and Phi data can elucidate whole experimental folding pathways.

The folding and evolution of multidomain proteins

Analyses of genomes show that more than 70% of eukaryotic proteins are composed of multiple domains. However, most studies of protein folding focus on individual domains and do not consider how interactions between domains might affect folding. Here, we address this by analysing the three-dimensional structures of multidomain proteins that have been characterized experimentally and observe that where the interface is small and loosely packed, or unstructured, the folding of the domains is independent. Furthermore, recent studies indicate that multidomain proteins have evolved mechanisms to minimize the problems of interdomain misfolding.

Small Proteins Fold Through Transition States With Native-like Topologies

The folding pathway of common-type acyl phosphatase (ctAcP) is characterized using psi-analysis, which identifies specific chain-chain contacts using bi-histidine (biHis) metal-ion binding sites. In the transition state ensemble (TSE), the majority of the protein is structured with a near-native topology, only lacking one beta-strand and an alpha-helix. psi-Values are zero or unity for all sites except one at the amino terminus of helix H2. This fractional psi-value remains unchanged when three metal ions of differing coordination geometries are used, indicating this end of the helix experiences microscopic heterogeneity through fraying in the TSE. Ubiquitin, the other globular protein characterized using psi-analysis, also exhibits a single consensus TSE structure. Hence, the TSE of both proteins have converged to a single configuration, albeit one that contains some fraying at the periphery. Models of the TSE of both proteins are created using all-atom Langevin dynamics simulations using distance constraints derived from the experimental psi-values. For both proteins, the relative contact order of the TS models is approximately 80% of the native value. This shared value viewed in the context of the known correlation between contact order and folding rates, suggests that other proteins will have a similarly high fraction of the native contact order. This constraint greatly limits the range of possible configurations at the rate-limiting step.

Characterization of the hydrodynamic properties of the folding transition state of an SH3 domain by magnetization transfer NMR spectroscopy

Protein folding kinetic data have been obtained for the marginally stable N-terminal Src homology 3 domain of the Drosophila protein drk (drkN SH3) in an investigation of the hydrodynamic properties of its folding transition state. Due to the presence of NMR resonances of both folded and unfolded states at equilibrium, kinetic data can be derived from NMR magnetization transfer techniques under equilibrium conditions. Kinetic analysis as a function of urea (less than approximately 1 M) and glycerol enables determination of alpha values, measures of the energetic sensitivity of the transition state to the perturbation relative to the end states of the protein folding reaction (the folded and unfolded states). Both end states have previously been studied experimentally by NMR spectroscopic and other biophysical methods in great detail and under nondenaturing conditions. Combining these results with the kinetic folding data obtained here, we can characterize the folding transition state without requiring empirical models for the unfolded state structure. We are thus able to give a reliable measure of the solvent-accessible surface area of the transition state of the drkN SH3 domain (4730 +/- 360 A(2)) based on urea titration data. Glycerol titration data give similar results and additionally demonstrate that folding of this SH3 domain is dependent on solvent viscosity, which is indicative of at least partial hydration of the transition state. Because SH3 domains appear to fold by a common folding mechanism, the data presented here provide valuable insight into the transition states of the drkN and other SH3 domains.

Phi-analysis of the folding of the B domain of protein A using multiple optical probes

We examined the co-operativity of ultra-fast folding of a protein and whether the Phi-value analysis of its transition state depended on the location of the optical probe. We incorporated in turn a tryptophan residue into each of the three helices of the B domain of Protein A. Each Trp mutant of the three-helix bundle protein was used as a pseudo-wild-type parent for Phi-analysis in which the intrinsic Trp fluorescence probed the formation of each helix during the transition state. Apart from local effects in the immediate vicinity of the probe, the three separate sets of Phi-values were in excellent agreement, demonstrating the overall co-operativity of folding and the robustness of the Phi-analysis. The transition state of folding of Protein A contains the second helix being well formed with many stabilizing tertiary hydrophobic interactions. In contrast, the first and the third helices are more poorly structured in the transition state. The mechanism of folding thus involves the concurrent formation of secondary and tertiary interactions, and is towards the nucleation-condensation extreme in the nucleation-condensation-framework continuum of mechanism, with helix 2 being the nucleus. We provide an error analysis of Phi-values derived purely from the kinetics of two-state chevron plots.

Do proteins with similar folds have similar transition state structures? A diffuse transition state of the 169 residue apoflavodoxin

Apoflavodoxin from Anabaena PCC 7119 is a 169 residue globular protein of known structure and energetics. Here, we present a comprehensive Phi-value analysis to characterize the structure of its transition state. A total of 34 non-disruptive mutations are made throughout the structure and a range of Phi-values from zero to one are observed. In addition, a small set of eight aliphatic small-to-large mutations have been introduced in the hydrophobic core of the protein and they have been analyzed to investigate the feasibility of stabilizing the unfolding transition state by creating new non-native interactions. We find that the transition state of apoflavodoxin (so far the largest protein subjected to Phi-analysis) is diffuse and that it can be stabilized by unspecific hydrophobic interactions that can speed up the folding reaction. The data gathered on the apoflavodoxin transition state are compared with results from experimental studies in other proteins to revisit the relationship between the native state topology and transition state structure.

On the precision of experimentally determined protein folding rates and phi-values

Phi-values, a relatively direct probe of transition-state structure, are an important benchmark in both experimental and theoretical studies of protein folding. Recently, however, significant controversy has emerged regarding the reliability with which phi-values can be determined experimentally: Because phi is a ratio of differences between experimental observables it is extremely sensitive to errors in those observations when the differences are small. Here we address this issue directly by performing blind, replicate measurements in three laboratories. By monitoring within- and between-laboratory variability, we have determined the precision with which folding rates and phi-values are measured using generally accepted laboratory practices and under conditions typical of our laboratories. We find that, unless the change in free energy associated with the probing mutation is quite large, the precision of phi-values is relatively poor when determined using rates extrapolated to the absence of denaturant. In contrast, when we employ rates estimated at nonzero denaturant concentrations or assume that the slopes of the chevron arms (mf and mu) are invariant upon mutation, the precision of our estimates of phi is significantly improved. Nevertheless, the reproducibility we thus obtain still compares poorly with the confidence intervals typically reported in the literature. This discrepancy appears to arise due to differences in how precision is calculated, the dependence of precision on the number of data points employed in defining a chevron, and interlaboratory sources of variability that may have been largely ignored in the prior literature.

Cofactor effects on the protein folding reaction: acceleration of alpha-lactalbumin refolding by metal ions

About 30% of proteins require cofactors for their proper folding. The effects of cofactors on the folding reaction have been investigated with alpha-lactalbumin as a model protein and metal ions as cofactors. Metal ions accelerate the refolding of alpha-lactalbumin by lessening the energy barrier between the molten globule state and the transition state, mainly by decreasing the difference of entropy between the two states. These effects are linked to metal ion binding to the protein in the native state. Hence, relationships between the metal affinities for the intermediate states and those for the native state are observed. Some residual specificity for the calcium ion is still observed in the molten globule state, this specificity getting closer in the transition state to that of the native state. The comparison between kinetic and steady-state data in association with the Phi value method indicates the binding of the metal ions on the unfolded state of alpha-lactalbumin. Altogether, these results provide insight into cofactor effects on protein folding. They also suggest new possibilities to investigate the presence of residual native structures in the unfolded state of protein and the effects of such structures on the protein folding reaction and on protein stability.

Improved pairwise alignments of proteins in the Twilight Zone using local structure predictions

MOTIVATION

In recent years, advances have been made in the ability of computational methods to discriminate between homologous and non-homologous proteins in the 'twilight zone' of sequence similarity, where the percent sequence identity is a poor indicator of homology. To make these predictions more valuable to the protein modeler, they must be accompanied by accurate alignments. Pairwise sequence alignments are inferences of orthologous relationships between sequence positions. Evolutionary distance is traditionally modeled using global amino acid substitution matrices. But real differences in the likelihood of substitutions may exist for different structural contexts within proteins, since structural context contributes to the selective pressure.

RESULTS

HMMSUM (HMMSTR-based substitution matrices) is a new model for structural context-based amino acid substitution probabilities consisting of a set of 281 matrices, each for a different sequence-structure context. HMMSUM does not require the structure of the protein to be known. Instead, predictions of local structure are made using HMMSTR, a hidden Markov model for local structure. Alignments using the HMMSUM matrices compare favorably to alignments carried out using the BLOSUM matrices or structure-based substitution matrices SDM and HSDM when validated against remote homolog alignments from BAliBASE. HMMSUM has been implemented using local Dynamic Programming and with the Bayesian Adaptive alignment method.

Fine structure analysis of a protein folding transition state; distinguishing between hydrophobic stabilization and specific packing

Developing a detailed understanding of the structure and energetics of protein folding transition states is a key step in describing the folding process. The phi-value analysis approach allows the energetic contribution of side-chains to be mapped out by comparing wild-type with individual mutants where conservative changes are introduced. Studies where multiple substitutions are made at individual sites are much rarer but are potentially very useful for understanding the contribution of each element of a side-chain to transition state formation, and for distinguishing the relative importance of specific packing versus hydrophobic interactions. We have made a series of conservative mutations at multiple buried sites in the N-terminal domain of L9 in order to assess the relative importance of specific side-chain packing versus less specific hydrophobic stabilization of the transition state. A total of 28 variants were prepared using both naturally occurring and non-naturally occurring amino acids at six sites. Analysis of the mutants by NMR and CD showed no perturbation of the structure. There is no correlation between changes in hydrophobicity and changes in stability. In contrast, there is excellent linear correlation between the hydrophobicity of a side-chain and the log of the folding rate, ln(k(f)). The correlation between ln(k(f)) and the change in hydrophobicity holds even for substitutions that change the shape and/or size of a side-chain significantly. For most sites, the correlation with the logarithm of the unfolding rate, ln(k(u)), is much worse. Mutants with more hydrophobic amino acid substitutions fold faster, and those with less hydrophobic amino acid substitutions fold slower. The results show that hydrophobic interactions amongst core residues are an important driving force for forming the transition state, and are more important than specific tight packing interactions. Finally, a number of substitutions lead to negative phi-values and the origin of these effects are described.

Nucleation and the transition state of the SH3 domain

We present a verified computational model of the SH3 domain transition state (TS) ensemble. This model was built for three separate SH3 domains using experimental phi-values as structural constraints in all-atom protein folding simulations. While averaging over all conformations incorrectly considers non-TS conformations as transition states, quantifying structures as pre-TS, TS, and post-TS by measurement of their transmission coefficient ("probability to fold", or p(fold)) allows for rigorous conclusions regarding the structure of the folding nucleus and a full mechanistic analysis of the folding process. Through analysis of the TS, we observe a highly polarized nucleus in which many residues are solvent-exposed. Mechanistic analysis suggests the hydrophobic core forms largely after an early nucleation step. SH3 presents an ideal system for studying the nucleation-condensation mechanism and highlights the synergistic relationship between experiment and simulation in the study of protein folding.

Probing the interactions between the folding elements early in the folding of Escherichia coli dihydrofolate reductase by systematic sequence perturbation analysis

One of the necessary conditions for a protein to be foldable is the presence of a complete set of folding elements (FEs) that are short contiguous peptide segments distributed over an amino acid sequence. Previous studies indicated the FE assembly model of protein folding, in which the FEs interact with each other and coalesce to form an intermediate(s) early in the folding reaction. This suggests that a clue to the understanding of the determinants of protein foldability can be found by investigating how the FEs interact with each other early in the folding and thereby elucidating roles of the FEs in protein folding. To reveal the formation process of FE-FE interactions, we studied the early folding events of Escherichia coli dihydrofolate reductase (DHFR) utilizing systematic sequence perturbation analysis. Here, systematic single amino acid substitutions were introduced inside of the FEs (W30X in FE2, V40X in FE3, N59X in FE4, and I155X in FE10; X refers to various amino acid residues), and their kinetic refolding reactions were measured by stopped-flow circular dichroism and fluorescence. We show that the interactions around Trp30 and Ile155 are formed in the burst phase intermediate, while those around Val40 and Asn59 are formed in the transition state of the subsequent folding phase (tau5-phase) and in much later processes, respectively. These and previous results suggest that FE2 and FE10, and also FE1 and FE7, involved in the loop subdomain of DHFR, interact with each other within a millisecond time range, while the stable FE3-FE4 interactions are formed in the later processes. This may highlight the important roles of the FEs mainly inside of the loop subdomain in formation of the burst phase intermediate having a hydrophobic cluster and native-like overall topology and in acquisition of the foldability of DHFR.

Comparison of the transition states for folding of two Ig-like proteins from different superfamilies

In the "fold approach" proteins with a similar fold but different sequences are compared in order to investigate the relationship between native state structure and folding behaviour. Here we compare the properties of the transition states for folding of TI I27, the 27th immunoglobulin domain from human cardiac titin, and that of TNfn3, the third fibronectin type III domain from human tenascin. Experimental phi-values were used as restraints in molecular dynamics simulations to determine the structures that make up the transition state ensembles (TSEs) for folding of the two proteins. The restrained simulations that we present allow a detailed structural comparison of the two TSEs to be made. Further calculations show explicitly that for both proteins the formation of the interactions involving the residues in the folding nucleus is sufficient for the establishment of the topology of the Ig-like fold. We found that, although the folding nuclei of the two proteins are similar, the packing of the folding nucleus of TI I27 is much tighter than that of TNfn3, reflecting the higher experimental phi-values and beta(T) (Tanford Beta) of TI I27. These results suggest that the folding nucleus can be significantly deformed to accommodate extensive sequence variation while conserving the same folding mechanism.

Differences in the folding transition state of ubiquitin indicated by phi and psi analyses

We compare the folding transition state (TS) of ubiquitin previously identified by using psi analysis to that determined by using analysis. Both methods attempt to identify interactions and their relative populations at the rate-limiting step for folding. The TS ensemble derived from psi analysis has an extensive native-like chain topology, with a four-stranded beta-sheet network and a portion of the major helix. According to analysis, however, the TS is much smaller and more polarized, with only a local helix/hairpin motif. We find that structured regions can have values far from unity, the canonical value for such sites, because of structural relaxation of the TS. Consequently, these sites may be incorrectly interpreted as contributing little to the structure of the TS. These results stress the need for caution when interpreting and drawing conclusions from analysis alone and highlight the need for more specific tools for examining the structure and energetics of the TS ensemble.

Context-dependent effects of proline residues on the stability and folding pathway of ubiquitin

Substitution of trans-proline at three positions in ubiquitin (residues 19, 37 and 38) produces significant context-dependent effects on protein stability (both stabilizing and destabilizing) that reflect changes to a combination of parameters including backbone flexibility, hydrophobic interactions, solvent accessibility to polar groups and intrinsic backbone conformational preferences. Kinetic analysis of the wild-type yeast protein reveals a predominant fast-folding phase which conforms to an apparent two-state folding model. Temperature-dependent studies of the refolding rate reveal thermodynamic details of the nature of the transition state for folding consistent with hydrophobic collapse providing the overall driving force. Brønsted analysis of the refolding and unfolding rates of a family of mutants with a variety of side chain substitutions for P37 and P38 reveals that the two prolines, which are located in a surface loop adjacent to the C terminus of the main alpha-helix (residues 24-33), are not significantly structured in the transition state for folding and appear to be consolidated into the native structure only late in the folding process. We draw a similar conclusion regarding position 19 in the loop connecting the N-terminal beta-hairpin to the main alpha-helix. The proline residues of ubiquitin are passive spectators in the folding process, but influence protein stability in a variety of ways.

Formation of the folding nucleus of an SH3 domain investigated by loosely coupled molecular dynamics simulations

The experimentally well-established folding mechanism of the src-SH3 domain, and in particular the phi-value analysis of its transition state, represents a sort of testing table for computational investigations of protein folding. Here, parallel molecular dynamics simulations of the src-SH3 domain have been performed starting from denatured conformations. By rescuing and restarting only trajectories approaching the folding transition state, an ensemble of conformations was obtained with a completely structured central beta-sheet and a native-like packing of residues Ile-110, Ala-121, and Ile-132. An analysis of the trajectories shows that there are several pathways leading to the formation of the central beta-sheet whereas its two hairpins form in a different but consistent way.

Low-populated folding intermediates of Fyn SH3 characterized by relaxation dispersion NMR

Many biochemical processes proceed through the formation of functionally significant intermediates. Although the identification and characterization of such species can provide vital clues about the mechanisms of the reactions involved, it is challenging to obtain information of this type in cases where the intermediates are transient or present only at low population. One important example of such a situation involves the folding behaviour of small proteins that represents a model for the acquisition of functional structure in biology. Here we use relaxation dispersion nuclear magnetic resonance (NMR) spectroscopy to identify, for two mutational variants of one such protein, the SH3 domain from Fyn tyrosine kinase, a low-population folding intermediate in equilibrium with its unfolded and fully folded states. By performing the NMR experiments at different temperatures, this approach has enabled characterization of the kinetics and energetics of the folding process as well as providing structures of the intermediates. A general strategy emerges for an experimental determination of the energy landscape of a protein by applying this methodology to a series of mutants whose intermediates have differing degrees of native-like structure.

Phi-value analysis and the nature of protein-folding transition states

Phi values are used to map structures of protein-folding transition states from changes in free energies of denaturation (DeltaDeltaG(D-N)) and activation on mutation. A recent reappraisal proposed that Phi values for DeltaDeltaG(D-N) < 1.7 kcal/mol are artifactual. On discarding such derived Phi values from published studies, the authors concluded that there are no high Phi values in diffuse transition states, which are consequently uniformly diffuse with no evidence for nucleation. However, values of DeltaDeltaG(D-N) > 1.7 kcal/mol are often found for large side chains that make dispersed tertiary interactions, especially in hydrophobic cores that are in the process of being formed in the transition state. Conversely, specific local interactions that probe secondary structure tend to have DeltaDeltaG(D-N) approximately 0.5-2 kcal/mol. Discarding Phi values from lower-energy changes discards the crucial information about local interactions and makes transition states appear uniformly diffuse by overemphasizing the dispersed tertiary interactions. The evidence for the 1.7 kcal/mol cutoff was based on mutations that had been deliberately designed to be unsuitable for Phi-value analysis because they are structurally disruptive. We confirm that reliable Phi values can be derived from the recommended mutations in suitable proteins with 0.6 < DeltaDeltaG(D-N) < 1.7 kcal/mol, and there are many reliable high Phi values. Transition states vary from being rather diffuse to being well formed with islands of near-complete secondary structure. We also confirm that the structures of transition-state ensembles can be perturbed by mutations with DeltaDeltaG(D-N) >> 2 kcal/mol and that protein-folding transition states do move on the energy surface on mutation.

Dramatic acceleration of protein folding by stabilization of a nonnative backbone conformation

Through a mutagenic investigation of Gly-48, a highly conserved position in the Src homology 3 domain, we have discovered a series of amino acid substitutions that are highly destabilizing, yet dramatically accelerate protein folding, some up to 10-fold compared with the wild-type rate. The unique folding properties of these mutants allowed for accurate measurement of their folding and unfolding rates in water with no denaturant by using an NMR spin relaxation dispersion technique. A strong correlation was found between beta-sheet propensity and the folding rates of the Gly-48 mutants, even though Gly-48 lies in an unusual non-beta-strand backbone conformation in the native state. This finding indicates that the accelerated folding rates of the Gly-48 mutants are the result of stabilization of a nonnative beta-strand conformation in the transition-state structure at this position, thus providing the first, to our knowledge, experimentally elucidated example of a mechanism by which folding can occur fastest through a nonnative conformation. We also demonstrate that residues that are most stabilizing in the transition-state structure are most destabilizing in the native state, and also cause the greatest reductions in in vitro functional activity. These data indicate that the unusual native conformation of the Gly-48 position is important for function, and that evolutionary selection for function can result in a domain that folds at a rate far below the maximum possible.

Conformational Heterogeneity of an Equilibrium Folding Intermediate Quantified and Mapped by Scanning Mutagenesis

It is challenging to experimentally define an energy landscape for protein folding that comprises multiple partially unfolded states. Experimental results are often ambiguous as to whether a non-native state is conformationally homogeneous. Here, we tested an approach combining systematic mutagenesis and a Brønsted-like analysis to reveal and quantify conformational heterogeneity of folding intermediate states. Using this method, we resolved an otherwise apparently homogeneous equilibrium folding intermediate of Borrelia burgdorferi OspA into two conformationally distinct species and determined their relative populations. Furthermore, we mapped the structural differences between these intermediate species, which are consistent with the non-native species that we previously proposed based on native-state hydrogen exchange studies. When treated as a single state, the intermediate ensemble exhibited fractional Phi-values for mutations and Hammond-type behaviors that are often observed for folding transition states. We found that a change in relative population of the two species within the intermediate ensemble explains these properties well, suggesting that fractional Phi-values and Hammond-type behaviors exhibited by folding intermediates and transition states may arise more often from conformational heterogeneity than from a single partial structure. Our results are consistent with the presence of multiple minima in a rugged energy landscape predicted from theoretical studies. The method described here provides a promising means to probe a complex folding energy landscape.

Transition states for protein folding have native topologies despite high structural variability

We present a structural analysis of the folding transition states of three SH3 domains. Our results reveal that the secondary structure is not yet fully formed at this stage of folding and that the solvent is only partially excluded from the interior of the protein. Comparison of the members of the transition state ensemble with a database of native folds shows that, despite substantial local variability, the transition state structures can all be classified as having the topology characteristic of an SH3 domain. Our results suggest a mechanism for folding in which the formation of a network of interactions among a subset of hydrophobic residues ensures that the native topology is generated. Such a mechanism enables high fidelity in folding while minimizing the need to establish a large number of specific interactions in the conformational search.