An Overlapping Region between the Two Terminal Folding Units of the Outer Surface Protein A (OspA) Controls Its Folding Behavior

Although many naturally occurring proteins consist of multiple domains, most studies on protein folding to date deal with single-domain proteins or isolated domains of multi-domain proteins. Studies of multi-domain protein folding are required for further advancing our understanding of protein folding mechanisms. Borrelia outer surface protein A (OspA) is a β-rich two-domain protein, in which two globular domains are connected by a rigid and stable single-layer β-sheet. Thus, OspA is particularly suited as a model system for studying the interplays of domains in protein folding. Here, we studied the equilibria and kinetics of the urea-induced folding–unfolding reactions of OspA probed with tryptophan fluorescence and ultraviolet circular dichroism. Global analysis of the experimental data revealed compelling lines of evidence for accumulation of an on-pathway intermediate during kinetic refolding and for the identity between the kinetic intermediate and a previously described equilibrium unfolding intermediate. The results suggest that the intermediate has the fully native structure in the N-terminal domain and the single layer β-sheet, with the C-terminal domain still unfolded. The observation of the productive on-pathway folding intermediate clearly indicates substantial interactions between the two domains mediated by the single-layer β-sheet. We propose that a rigid and stable intervening region between two domains creates an overlap between two folding units and can energetically couple their folding reactions.     

The protein folding 'speed limit'

How fast can a protein possibly fold? This question has stimulated experimentalists to seek fast folding proteins and to engineer them to fold even faster. Proteins folding at or near the speed limit are prime candidates for all-atom molecular dynamics simulations. They may also have no free energy barrier, allowing the direct observation of intermediate structures on the pathways from the unfolded to the folded state. Both experimental and theoretical approaches predict a speed limit of approximately N/100micros for a generic N-residue single-domain protein, with alpha proteins folding faster than beta or alphabeta. The predicted limits suggest that most known ultrafast folding proteins can be engineered to fold more than ten times faster.     

Dynamics, energetics, and structure in protein folding

For many decades, protein folding experimentalists have worked with no information about the time scales of relevant protein folding motions and without methods for estimating the height of folding barriers. Protein folding experiments have been interpreted using chemical models in which the folding process is characterized as a series of equilibria between two or more distinct states that interconvert with activated kinetics. Accordingly, the information to be extracted from experiments was circumscribed to apparent equilibrium constants and relative folding rates. Recent developments are changing this situation dramatically. The combination of fast-folding experiments with the development of analytical methods more closely connected to physical theory reveals that folding barriers in native conditions range from minimally high (approximately 14RT for the very slow folder AcP) to nonexistent. While slow-folding (i.e., > or = 1 ms) single-domain proteins are expected to fold in a two-state fashion, microsecond-folding proteins should exhibit complex behavior arising from crossing marginal or negligible folding barriers. This realization opens a realm of exciting opportunities for experimentalists. The free energy surface of a protein with a marginal (or no) barrier can be mapped using equilibrium experiments, which could resolve energetic factors from structural factors in folding. Kinetic experiments on these proteins provide the unique opportunity to measure folding dynamics directly. Furthermore, the complex distributions of time-dependent folding behaviors expected for these proteins might be accessible to single-molecule measurements. Here, we discuss some of these recent developments in protein folding, emphasizing aspects that can serve as a guide for experimentalists interested in exploiting this new avenue of research.     

Why do proteins divide into domains? Insights from lattice model simulations

It is known that larger globular proteins are built from domains, relatively independent structural units. A domain size seems to be limited, and a single domain consists of from few tens to a couple of hundred amino acids. Based on Monte Carlo simulations of a reduced protein model restricted to the face centered simple cubic lattice, with a minimal set of short-range and long-range interactions, we have shown that some model sequences upon the folding transition spontaneously divide into separate domains. The observed domain sizes closely correspond to the sizes of real protein domains. Short chains with a proper sequence pattern of the hydrophobic and polar residues undergo a two-state folding transition to the structurally ordered globular state, while similar longer sequences follow a multistate transition. Homopolymeric (uniformly hydrophobic) chains and random heteropolymers undergo a continuous collapse transition into a single globule, and the globular state is much less ordered. Thus, the factors responsible for the multidomain structure of proteins are sufficiently long polypeptide chain and characteristic, protein-like, sequence patterns. These findings provide some hints for the analysis of real sequences aimed at prediction of the domain structure of large proteins.     

Class-specific correlations between protein folding rate, structure-derived, and sequence-derived descriptors

Small single-domain proteins that fold by simple two-state kinetics have been shown to exhibit a wide variation in their folding rates. It has been proposed that folding mechanisms in these proteins are largely determined by the native-state topology, and a significant correlation between folding rate and measures of the average topological complexity, such as relative contact order (RCO), has been reported. We perform a statistical analysis of folding rate and RCO in all three major structural classes (alpha, beta, and alpha/beta) of small two-state proteins and of RCO in groups of analogous and homologous small single-domain proteins with the same topology. We also study correlation between folding rate and the average physicochemical properties of amino acid sequences in two-state proteins. Our results indicate that 1) helical proteins have statistically distinguishable, class-specific folding rates; 2) RCO accounts for essentially all the variation of folding rate in helical proteins, but for only a part of the variation in beta-sheet-containing proteins; and 3) only a small fraction of the protein topologies studied show a topology-specific RCO. We also report a highly significant correlation between the folding rate and average intrinsic structural propensities of protein sequences. These results suggest that intrinsic structural propensities may be an important determinant of the rate of folding in small two-state proteins.     

Maximal stabilities of reversible two-state proteins

The hydrophobic effect is the major force driving protein folding. Around room temperature, small organic solutes and hydrophobic amino acids have low solubilities in water and the hydrophobic effect is the strongest. These facts suggest that globular proteins should be maximally stable around room temperature. While this fundamental paradigm has been expected, it has not actually been shown to hold. Toward this goal, we have collected and analyzed experimental thermodynamic data for 31 proteins that show reversible two-state folding <--> unfolding transitions at or near neutral pH. Twenty-six of these are unique, and 20 of the 26 are maximally stable around room temperature irrespective of their structural properties, the melting temperature, or the living temperatures of their source organisms. Their average temperature of maximal stability is 293 +/- 8 K (20 +/- 8 degrees C). These proteins differ in size, fold, and number of domains, hydrophobic folding units, and oligomeric states. They derive from the cold-loving psychrophiles, from mesophiles, and from thermophiles. Analysis of the single-domain proteins present in this set shows that the variations in their thermodynamic parameters are correlated in a way which may explain the adaptation of the proteins to the living temperatures of the organisms from which they derive. The average energetic contribution of the individual amino acids toward protein stability decreases with an increase in protein size, suggesting that there may be an upper limit for protein maximal thermodynamic stability. For the remaining proteins, deviation of the maximal stability temperatures from room temperature may be due to greater uncertainties in their heat capacity change (DeltaC(p)) values, a weaker hydrophobic effect, and/or a stronger electrostatic contribution.     

Cunning simplicity of a hierarchical folding

A hierarchic scheme of protein folding does not solve the Levinthal paradox since it cannot provide a simultaneous explanation for major features observed for protein folding: (i) folding within non-astronomical time, (ii) independence of the native structure on large variations in the folding rates of given protein under different conditions, and (iii) co-existence, in a visible quantity, of only the native and the unfolded molecules during folding of moderate size (single-domain) proteins. On the contrary, a nucleation mechanism can account for all these major features simultaneously and resolves the Levinthal paradox.     

The folding pathway of T4 lysozyme: the high-resolution structure and folding of a hidden intermediate

Folding intermediates have been detected and characterized for many proteins. However, their structures at atomic resolution have only been determined for two small single domain proteins: Rd-apocytochrome b(562) and engrailed homeo domain. T4 lysozyme has two easily distinguishable but energetically coupled domains: the N and C-terminal domains. An early native-state hydrogen exchange experiment identified an intermediate with the C-terminal domain folded and the N-terminal domain unfolded. We have used a native-state hydrogen exchange-directed protein engineering approach to populate this intermediate and demonstrated that it is on the folding pathway and exists after the rate-limiting step. Here, we determined its high-resolution structure and the backbone dynamics by multi-dimensional NMR methods. We also characterized the folding behavior of the intermediate using stopped-flow fluorescence, protein engineering, and native-state hydrogen exchange. Unlike the folding intermediates of the two single-domain proteins, which have many non-native side-chain interactions, the structure of the hidden folding intermediate of T4 lysozyme is largely native-like. It folds like many small single domain proteins. These results have implications for understanding the folding mechanism and evolution of multi-domain proteins.     

Contact order revisited: influence of protein size on the folding rate

Guided by the recent success of empirical model predicting the folding rates of small two-state folding proteins from the relative contact order (CO) of their native structures, by a theoretical model of protein folding that predicts that logarithm of the folding rate decreases with the protein chain length L as L(2/3), and by the finding that the folding rates of multistate folding proteins strongly correlate with their sizes and have very bad correlation with CO, we reexamined the dependence of folding rate on CO and L in attempt to find a structural parameter that determines folding rates for the totality of proteins. We show that the Abs_CO = CO x L, is able to predict rather accurately folding rates for both two-state and multistate folding proteins, as well as short peptides, and that this Abs_CO scales with the protein chain length as L(0.70 +/- 0.07) for the totality of studied single-domain proteins and peptides.     

Is protein folding rate dependent on number of folding stages? Modeling of protein folding with ferredoxin-like fold

Statistical analysis of protein folding rates has been done for 84 proteins with available experimental data. A surprising result is that the proteins with multi-state kinetics from the size range of 50–100 amino acid residues (a.a.) fold as fast as proteins with two-state kinetics from the same size range. At the same time, the proteins with two-state kinetics from the size range 101–151 a.a. fold faster than those from the size range 50–100 a.a. Moreover, it turns out unexpectedly that usually in the group of structural homologs from the size range 50–100 a.a., proteins with multi-state kinetics fold faster than those with two-state kinetics. The protein folding for six proteins with a ferredoxin-like fold and with a similar size has been modeled using Monte Carlo simulations and dynamic programming. Good correlation between experimental folding rates, some structural parameters, and the number of Monte Carlo steps has been obtained. It is shown that a protein with multi-state kinetics actually folds three times faster than its structural homologs.     

How cooperative are protein folding and unfolding transitions?

A thermodynamically and kinetically simple picture of protein folding envisages only two states, native (N) and unfolded (U), separated by a single activation free energy barrier, and interconverting by cooperative two‐state transitions. The folding/unfolding transitions of many proteins occur, however, in multiple discrete steps associated with the formation of intermediates, which is indicative of reduced cooperativity. Furthermore, much advancement in experimental and computational approaches has demonstrated entirely non‐cooperative (gradual) transitions via a continuum of states and a multitude of small energetic barriers between the N and U states of some proteins. These findings have been instrumental towards providing a structural rationale for cooperative versus noncooperative transitions, based on the coupling between interaction networks in proteins. The cooperativity inherent in a folding/unfolding reaction appears to be context dependent, and can be tuned via experimental conditions which change the stabilities of N and U. The evolution of cooperativity in protein folding transitions is linked closely to the evolution of function as well as the aggregation propensity of the protein. A large activation energy barrier in a fully cooperative transition can provide the kinetic control required to prevent the accumulation of partially unfolded forms, which may promote aggregation. Nevertheless, increasing evidence for barrier‐less “downhill” folding, as well as for continuous “uphill” unfolding transitions, indicate that gradual non‐cooperative processes may be ubiquitous features on the free energy landscape of protein folding.     

Folding pathway of the b1 domain of protein G explored by multiscale modeling

The understanding of the folding mechanisms of single-domain proteins is an essential step in the understanding of protein folding in general. Recently, we developed a mesoscopic CA-CB side-chain protein model, which was successfully applied in protein structure prediction, studies of protein thermodynamics, and modeling of protein complexes. In this research, this model is employed in a detailed characterization of the folding process of a simple globular protein, the B1 domain of IgG-binding protein G (GB1). There is a vast body of experimental facts and theoretical findings for this protein. Performing unbiased, ab initio simulations, we demonstrated that the GB1 folding proceeds via the formation of an extended folding nucleus, followed by slow structure fine-tuning. Remarkably, a subset of native interactions drives the folding from the very beginning. The emerging comprehensive picture of GB1 folding perfectly matches and extends the previous experimental and theoretical studies.     

Direct Observation of Parallel Folding Pathways Revealed Using a Symmetric Repeat Protein System

Although progress has been made to determine the native fold of a polypeptide from its primary structure, the diversity of pathways that connect the unfolded and folded states has not been adequately explored. Theoretical and computational studies predict that proteins fold through parallel pathways on funneled energy landscapes, although experimental detection of pathway diversity has been challenging. Here, we exploit the high translational symmetry and the direct length variation afforded by linear repeat proteins to directly detect folding through parallel pathways. By comparing folding rates of consensus ankyrin repeat proteins (CARPs), we find a clear increase in folding rates with increasing size and repeat number, although the size of the transition states (estimated from denaturant sensitivity) remains unchanged. The increase in folding rate with chain length, as opposed to a decrease expected from typical models for globular proteins, is a clear demonstration of parallel pathways. This conclusion is not dependent on extensive curve-fitting or structural perturbation of protein structure. By globally fitting a simple parallel-Ising pathway model, we have directly measured nucleation and propagation rates in protein folding, and have quantified the fluxes along each path, providing a detailed energy landscape for folding. This finding of parallel pathways differs from results from kinetic studies of repeat-proteins composed of sequence-variable repeats, where modest repeat-to-repeat energy variation coalesces folding into a single, dominant channel. Thus, for globular proteins, which have much higher variation in local structure and topology, parallel pathways are expected to be the exception rather than the rule.     

Cunning Simplicity of a Hierarchical Folding

Abstract

A hierarchic scheme of protein folding does not solve the Levinthal paradox since it cannot provide a simultaneous explanation for major features observed for protein folding: (i) folding within non-astronomical time, (ii) independence of the native structure on large variations in the folding rates of given protein under different conditions, and (iii) co-existence, in a visible quantity, of only the native and the unfolded molecules during folding of moderate size (single-domain) proteins. On the contrary, a nucleation mechanism can account for all these major features simultaneously and resolves the Levinthal paradox.     

Compactness determines protein folding type

We have demonstrated here that protein compactness, which we define as the ratio of the accessible surface area of a protein to that of the ideal sphere of the same volume, is one of the factors determining the mechanism of protein folding. Proteins with multi-state kinetics, on average, are more compact (compactness is 1.49 +/- 0.02 for proteins within the size range of 101-151 amino acid residues) than proteins with two-state kinetics (compactness is 1.59 +/- 0.03 for proteins within the same size range of 101-151 amino acid residues). We have shown that compactness for homologous proteins can explain both the difference in folding rates and the difference in folding mechanisms.     

Influence of medium- and long-range interactions in different folding types of globular proteins

Recognition of protein fold from amino acid sequence is a challenging task. The structure and stability of proteins from different fold are mainly dictated by inter-residue interactions. In our earlier work, we have successfully used the medium- and long-range contacts for predicting the protein folding rates, discriminating globular and membrane proteins and for distinguishing protein structural classes. In this work, we analyze the role of inter-residue interactions in commonly occurring folds of globular proteins in order to understand their folding mechanisms. In the medium-range contacts, the globin fold and four-helical bundle proteins have more contacts than that of DNA-RNA fold although they all belong to all-alpha class. In long-range contacts, only the ribonuclease fold prefers 4-10 range and the other folding types prefer the range 21-30 in alpha/beta class proteins. Further, the preferred residues and residue pairs influenced by these different folds are discussed. The information about the preference of medium- and long-range contacts exhibited by the 20 amino acid residues can be effectively used to predict the folding type of each protein.     

More compact protein globules exhibit slower folding rates

We have demonstrated that, among proteins of the same size, alpha/beta proteins have on the average a greater number of contacts per residue due to their more compact (more "spherical") structure, rather than due to tighter packing. We have examined the relationship between the average number of contacts per residue and folding rates in globular proteins according to general protein structural class (all-alpha, all-beta, alpha/beta, alpha+beta). Our analysis demonstrates that alpha/beta proteins have both the greatest number of contacts and the slowest folding rates in comparison to proteins from the other structural classes. Because alpha/beta proteins are also known to be the oldest proteins, it can be suggested that proteins have evolved to pack more quickly and into looser structures.     

Investigation of folding/unfolding kinetics of apomyoglobin

Apomyoglobin kinetic and equilibrium unfolding and folding processes were studied at pH 6.2, 11 degrees C by stopped-flow tryptophan fluorescence. There are two distinct consecutive processes in apomyoglobin folding process, namely, the protein fast transition between the unfolded (U) and an intermediate (I) states (U <----> I) and slow transition between the intermediate and the native (N) states (I <----> N). Accumulation of the intermediate state was observed in the wide range of urea concentrations. The presence of the intermediate state was shown even beyond the middle transition on the unfolding limb. The dependence of observed folding/unfolding rates on urea concentration (chevron plot) was obtained. The shape of this dependence was compared with that of two-state proteins, folding from the U to N state.     

Combination of the human prolyl isomerase FKBP12 with unrelated chaperone domains leads to chimeric folding enzymes with high activity

Folding enzymes often use distinct domains for the binding of substrate proteins ("chaperone domains") and for the catalysis of slow folding reactions such as disulfide formation or prolyl isomerization. The human prolyl isomerase FKBP12 is a small single-domain protein without a chaperone domain. Its very low folding activity could previously be increased by inserting the chaperone domain from the homolog SlyD (sensitive-to-lysis protein D) of Escherichia coli. We now inserted three unrelated chaperone domains into human FKBP12: the apical domain of the chaperonin GroEL from E. coli, the chaperone domain of protein disulfide isomerase from yeast, or the chaperone domain of SurA from the periplasm of E. coli. All three conveyed FKBP12 with a high affinity for unfolded proteins and increased its folding activity. Substrate binding and release of the chimeric folding enzymes were found to be very fast. This allows rapid substrate transfer from the chaperone domain to the catalytic domain and ensures efficient rebinding of protein chains that were unable to complete folding. The advantage of having separate sites, first for generic protein binding and then for specific catalysis, explains why our construction of the artificial folding enzymes with foreign chaperone domains was successful.