Repeat-protein folding: new insights into origins of cooperativity, stability, and topology

Predicting repeat protein folding kinetics from an experimentally determined folding energy landscape

The Notch ankyrin domain is a repeat protein whose folding has been characterized through equilibrium and kinetic measurements. In previous work, equilibrium folding free energies of truncated constructs were used to generate an experimentally determined folding energy landscape (Mello and Barrick, Proc Natl Acad Sci USA 2004;101:14102-14107). Here, this folding energy landscape is used to parameterize a kinetic model in which local transition probabilities between partly folded states are based on energy values from the landscape. The landscape-based model correctly predicts highly diverse experimentally determined folding kinetics of the Notch ankyrin domain and sequence variants. These predictions include monophasic folding and biphasic unfolding, curvature in the unfolding limb of the chevron plot, population of a transient unfolding intermediate, relative folding rates of 19 variants spanning three orders of magnitude, and a change in the folding pathway that results from C-terminal stabilization. These findings indicate that the folding pathway(s) of the Notch ankyrin domain are thermodynamically selected: the primary determinants of kinetic behavior can be simply deduced from the local stability of individual repeats.

Comparisons of the conformational stability of cyclin-dependent kinase (CDK) 4-interacting ankyrin repeat (AR) proteins

Ankyrin repeat (AR) proteins are one of the most abundant repeat protein classes in nature, and they are involved in numerous physiological processes through mediating protein/protein interactions. The repetitive and modular architecture of these AR proteins may lead to biochemical and biophysical properties distinct from those of globular proteins. It has been demonstrated that like most globular proteins, AR proteins exhibit a two-state, cooperative transition in chemical- and heat-induced unfolding. However, the biophysical characteristics underlying such cooperative unfolding remain to be further investigated. In the present study, we evaluated the conformational stability of a group of cyclin-dependent kinase (CDK) 4-interacting AR proteins, P16, P18, IkappaBalpha, gankyrin, and their truncated mutants under different conditions, including the presence of denaturants, temperature, and pH. Our results showed that the first four N-terminal ARs are required to form a potent and stable CDK4 modulator. Moreover, in spite of their similarities in skeleton structure, CDK4 binding, and cooperative unfolding, P16, P18, IkappaBalpha, and gankyrin exhibited considerably different biophysical properties with regard to the conformational stability, and these differences mainly resulted from the discrepancies in the primary sequence of the relatively conserved AR motifs. Our results also demonstrated that these sequence discrepancies are able to influence the function of AR proteins to a certain extent. Overall, our results provide important insights into understanding the biophysical properties of AR proteins.

A theoretical model for the mechanical unfolding of repeat proteins

We consider the mechanical stretching of a polypeptide chain formed by multiple interacting repeats. The folding thermodynamics and the interactions among the repeats are described by the Ising model. Unfolded repeats act as soft entropic springs, whereas folded repeats respond to a force as stiffer springs. We show that the resulting force-extension curve may exhibit a pronounced force maximum corresponding to the unfolding of the first repeat. This event is followed by the unfolding of the remaining repeats, which takes place at a lower force. As the protein extension is increased, the force-extension curve of a sufficiently long repeat protein displays a plateau, where the force remains nearly constant and the protein unfolds sequentially so that the number of unfolded repeats is proportional to the extension. Such a sequential mechanical unfolding mechanism is displayed even by the repeat proteins whose thermal denaturation is highly cooperative, provided that they are long enough. By contrast, the unfolding of short repeat progressions can be cooperative.

Invivo folding efficiencies for mutants of the P22 tailspike beta-helix protein correlate with predicted stability changes

Parallel beta-helices are among the simplest repetitive structural elements in proteins. The folding behavior of beta-helix proteins has been studied intensively, also to gain insight on the formation of amyloid fibrils, which share the parallel beta-helix as a central structural motif. An important system for investigating beta-helix folding is the tailspike protein from the Salmonella bacteriophage P22. The central domain of this protein is a right-handed parallel beta-helix with 13 windings. Extensive mutational analyses of the P22 tailspike protein have revealed two main phenotypes: temperature-sensitive-folding (tsf) mutations that reduce the folding efficiency at elevated temperatures, and global suppressor (su) mutations that increase the tailspike folding efficiency. A central question is whether these phenotypes can be understood from changes in the protein stability induced by the mutations. Experimental determination of the protein stability is complicated by the nearly irreversible trimerization of the folded tailspike protein. Here, we present calculations of stability changes with the program FoldX, focusing on a recently published extensive data set of 145 singe-residue alanine mutants. We find that the calculated stability changes are correlated with the experimentally measured invivo folding efficiencies. In addition, we determine the free-energy landscape of the P22 tailspike protein in a nucleation-propagation model to explore the folding mechanism of this protein, and obtain a processive folding route on which the protein nucleates in the N-terminal region of the helix.

What have we learned from the studies of two-state folders, and what are the unanswered questions about two-state protein folding?

Small proteins with globular structures often fold by simple all-or-none mechanisms, both in an equilibrium and a kinetic sense, despite the very large number of partly folded conformations available. This type of 'two-state' folding will be discussed in terms of experimental tests, underlying molecular mechanisms, and limits to two-state behavior. Factors that appear to be important for two-state folding include topology (sequence distance of contacts in the native structure), molecular cooperativity and local energy distribution. Because their local stability distributions and cooperativities can be dissected and analyzed separately from topological features, recent studies of the folding of symmetric proteins will be discussed as a means to better understand the origins of two-state folding.

Biological regulation via ankyrin repeat folding

By mimicking the phosphorylation of p19(INK4d), a tumor suppressor containing five ankyrin repeats, the native state could be destabilized to such an extent that only a partially folded state is populated at physiological temperature. This partly folded state, which mimics an on-pathway folding intermediate lacking structure in ankyrin repeats 1 and 2, is more rapidly ubiquitinated than the parent construct. Thus, phosphorylation of p19(INK4d) is likely to regulate cell-cycle progression through both biochemical (proteasomal) and biophysical (folding and binding to cyclin-dependent kinases) mechanisms.

Analysis of repeat-protein folding using nearest-neighbor statistical mechanical models

The linear "Ising" model, which has been around for nearly a century, treats the behavior of linear arrays of repetitive, interacting subunits. Linear "repeat-proteins" have only been described in the last decade or so, and their folding energies have only been characterized very recently. Owing to their repetitive structures, linear repeat-proteins are particularly well suited for analysis by the nearest-neighbor Ising formalism. After briefly describing the historical origins and applications of the Ising model to biopolymers, and introducing repeat protein structure, this chapter will focus on the application of the linear Ising model to repeat proteins. When applied to homopolymers, the model can be represented and applied in a fairly simplified form. When applied to heteropolymers, where differences in energies among individual subunits (i.e. repeats) must be included, some (but not all) of this simplicity is lost. Derivations of the linear Ising model for both homopolymer and heteropolymer repeat-proteins will be presented. With the increased complexity required for analysis of heteropolymeric repeat proteins, the ability to resolve different energy terms from experimental data can be compromised. Thus, a simple matrix approach will be developed to help inform on the degree to which different thermodynamic parameters can be extracted from a particular set of unfolding curves. Finally, we will describe the application of these models to analyze repeat-protein folding equilibria, focusing on simplified repeat proteins based on "consensus" sequence information.

Characterization of the unfolded state of repeat proteins

The unfolded state ensemble of proteins has been described as a structurally featureless state. While this approach is supported by the fact that many unfolded proteins follow the scaling law behavior of a random coil, there is evidence that the unfolded states of various proteins are stabilized by native or non-native interactions. Recently, the existence of extensive non-native structure was reported for a repeat protein, which resulted in a scaling law exponent that is significantly smaller than that of a random polymer [Cortajarena et al., J. Mol. Biol. 382(1), 203-212 (2008)]. It was concluded that the high compactness of this protein stems from a significant fraction of interacting PP(II) helical segments in the unfolded state. In this study, we aim at providing possible molecular understanding of this anomalous compactness of the unfolded state and to investigate its origin. Using a hierarchy of computational models, we ask whether in general the unfolded state of a repeat protein is likely to be intrinsically more compact than the unfolded state of globular proteins, or whether this phenomenon depends mostly on the occurrence of a specific sequence that promotes PP(II) conformations. Our results suggest that the formation of the PP(II) conformation is indeed essential, yet the recurring sequence of repeat proteins promotes the interactions between these PP(II) segments and the formation of non-native interactions in the unfolded state.

Studying the folding of multidomain proteins

There have been relatively few detailed comprehensive studies of the folding of protein domains (or modules) in the context of their natural covalently linked neighbors. This is despite the fact that a significant proportion of the proteome consists of multidomain proteins. In this review we highlight some key experimental investigations of the folding of multidomain proteins to draw attention to the difficulties that can arise in analyzing such systems. The evidence suggests that interdomain interactions can significantly affect stability, folding, and unfolding rates. However, preliminary studies suggest that folding pathways are unaffected-to this extent domains can be truly considered to be independent folding units. Nonetheless, it is clear that interactions between domains cannot be ignored, in particular when considering the effects of mutations.

Folding of elongated proteins: conventional or anomalous?

The complexity of the mechanisms by which proteins fold has been shown by many studies to be governed by their native-state topologies. This was manifested in the ability of the native topology-based model to capture folding mechanisms and the success of folding rate predictions based on various topological measures, such as the contact order. However, while the finer details of topological complexity have been thoroughly examined and related to folding kinetics, simpler characteristics of the protein, such as its overall shape, have been largely disregarded. In this study, we investigated the folding of proteins with an unusual elongated geometry that differs substantially from the common globular structure. To study the effect of the elongation degree on the folding kinetics, we used repeat proteins, which become more elongated as they include more repeating units. Some of these have apparently anomalous experimental folding kinetics, with rates that are often less than expected on the basis of rates for globular proteins possessing similar topological complexity. Using experimental folding rates and a larger set of rates obtained from simulations, we have shown that as the protein becomes increasingly elongated, its folding kinetics becomes slower and deviates more from the rate expected on the basis of topology measures fitted for globular proteins. The observed slow kinetics is a result of a more complex pathway in which stable intermediates composed of several consecutive repeats can appear. We thus propose a novel measure, an elongation-sensitive contact order, that takes into account both the extent of elongation and the topological complexity of the protein. This new measure resolves the apparent discrimination between the folding of globular and elongated repeat proteins. Our study extends the current capabilities of folding-rate predictions by unifying the kinetics of repeat and globular proteins.

Thermodynamics, kinetics, and salt dependence of folding of YopM, a large leucine-rich repeat protein

Small globular proteins have many contacts between residues that are distant in primary sequence. These contacts create a complex network between sequence-distant segments of secondary structure, which may be expected to promote the cooperative folding of globular proteins. Although repeat proteins, which are composed of tandem modular units, lack sequence-distant contacts, several of considerable length have been shown to undergo cooperative two-state folding. To explore the limits of cooperativity in repeat proteins, we have studied the unfolding of YopM, a leucine-rich repeat (LRR) protein of over 400 residues. Despite its large size and modular architecture (15 repeats), YopM equilibrium unfolding is highly cooperative, and shows a very strong dependence on the concentration of urea. In contrast, kinetic studies of YopM folding indicate a mechanism that includes one or more transient intermediates. The urea dependence of the folding and unfolding rates suggests a relatively small transition state ensemble. As with the urea dependence, we have found an extreme dependence of the free energy of unfolding on the concentration of salt. This salt dependence likely results from general screening of a large number of unfavorable columbic interactions in the folded state, rather than from specific cation binding.

Navigating the RNA folding landscape

Analysis of individual RNA folding reactions reveals that, as in proteins, cooperative interactions selectively drive RNA toward its biologically active, native conformation. This new work establishes a platform for future investigations of the physical principles underlying the assembly of large RNA enzymes.

The leucine-rich repeat domain of Internalin B folds along a polarized N-terminal pathway

The leucine-rich repeat domain of Internalin B is composed of seven tandem leucine-rich repeats, which each contain a short beta strand connected to a 3(10) helix by a short turn, and an N-terminal alpha-helical capping motif. To determine whether folding proceeds along a single, discrete pathway or multiple, parallel pathways, and to map the structure of the transition state ensemble, we examined the effects of destabilizing substitutions of conserved residues in each repeat. We find that, despite the structural redundancy among the repeats, folding proceeds through an N-terminal transition state ensemble in which the extent of structure formation is biased toward repeats one and two and includes both local and interrepeat interactions. Our results suggest that the N-terminal capping motif serves to polarize the folding pathway by acting as a fast-growing nucleus onto which consecutive repeats fold in the transition state ensemble, and highlight the importance of sequence-specific interactions in pathway selection.

Folding landscapes of ankyrin repeat proteins: experiments meet theory

Nearly 6% of eukaryotic protein sequences contain ankyrin repeat (AR) domains, which consist of several repeats and often function in binding. AR proteins show highly cooperative folding despite a lack of long-range contacts. Both theory and experiment converge to explain that formation of the interface between elements is more favorable than formation of any individual repeat unit. IkappaBalpha and Notch both undergo partial folding upon binding perhaps influencing the binding free energy. The simple architecture, combined with identification of consensus residues that are important for stability, has enabled systematic perturbation of the energy landscape by single point mutations that affect stability or by addition of consensus repeats. The folding energy landscapes appear highly plastic, with small perturbations re-routing folding pathways.

The energy landscapes of repeat-containing proteins: topology, cooperativity, and the folding funnels of one-dimensional architectures

Repeat-proteins are made up of near repetitions of 20- to 40-amino acid stretches. These polypeptides usually fold up into non-globular, elongated architectures that are stabilized by the interactions within each repeat and those between adjacent repeats, but that lack contacts between residues distant in sequence. The inherent symmetries both in primary sequence and three-dimensional structure are reflected in a folding landscape that may be analyzed as a quasi-one-dimensional problem. We present a general description of repeat-protein energy landscapes based on a formal Ising-like treatment of the elementary interaction energetics in and between foldons, whose collective ensemble are treated as spin variables. The overall folding properties of a complete "domain" (the stability and cooperativity of the repeating array) can be derived from this microscopic description. The one-dimensional nature of the model implies there are simple relations for the experimental observables: folding free-energy (DeltaG(water)) and the cooperativity of denaturation (m-value), which do not ordinarily apply for globular proteins. We show how the parameters for the "coarse-grained" description in terms of foldon spin variables can be extracted from more detailed folding simulations on perfectly funneled landscapes. To illustrate the ideas, we present a case-study of a family of tetratricopeptide (TPR) repeat proteins and quantitatively relate the results to the experimentally observed folding transitions. Based on the dramatic effect that single point mutations exert on the experimentally observed folding behavior, we speculate that natural repeat proteins are "poised" at particular ratios of inter- and intra-element interaction energetics that allow them to readily undergo structural transitions in physiologically relevant conditions, which may be intrinsically related to their biological functions.

Direct measurement of tertiary contact cooperativity in RNA folding

All structured biological macromolecules must overcome the thermodynamic folding problem to populate a unique functional state among a vast ensemble of unfolded and alternate conformations. The exploration of cooperativity in protein folding has helped reveal and distinguish the underlying mechanistic solutions to this folding problem. Analogous dissections of RNA tertiary stability remain elusive, however, despite the central biological importance of folded RNA molecules and the potential to reveal fundamental properties of structured macromolecules via comparisons of protein and RNA folding. We report a direct quantitative measure of tertiary contact cooperativity in a folded RNA. We precisely measured the stability of an independently folding P4-P6 domain from the Tetrahymena thermophila group I intron by single molecule fluorescence resonance energy transfer (smFRET). Using wild-type and mutant RNAs, we found that cooperativity between the two tertiary contacts enhances P4-P6 stability by 3.2 +/- 0.2 kcal/mol.

Rerouting the folding pathway of the Notch ankyrin domain by reshaping the energy landscape

The modular nature of repeat proteins has made them a successful target for protein design. Ankyrin repeat, TPR, and leucine rich repeat domains that have been designed solely on consensus information have been shown to have higher thermostability than their biological counterparts. We have previously shown that we can reshape the energy landscape of a repeat protein by adding multiple C-terminal consensus ankyrin repeats to the five N-terminal repeats of the Notch ankyrin domain. Here we explore how the folding mechanism responds to reshaping of the energy landscape. We have used analogous substitutions of a conserved alanine with glycine in each repeat to determine the distribution of structure in the transition state ensembles of constructs containing one (Nank1-5C1) and two consensus (Nank1-5C2) ankyrin repeats. Whereas folding of the wild-type Notch ankyrin domain is slowed by substitutions in its central repeats, (1) folding of Nank1-5C1 and Nank1-5C2 is slowed by substitutions in the C-terminal repeats. Thus, the addition of C-terminal stabilizing repeats shifts the transition state ensemble toward the C-terminal repeats, rerouting the folding pathway of the ankyrin repeat domain. These findings indicate that, for the Notch ankyrin domain, folding pathways are selected based on local energetics.

Folding thermodynamics and kinetics of the leucine-rich repeat domain of the virulence factor Internalin B

Although the folding of alpha-helical repeat proteins has been well characterized, much less is known about the folding of repeat proteins containing beta-sheets. Here we investigate the folding thermodynamics and kinetics of the leucine-rich repeat (LRR) domain of Internalin B (InlB), an extracellular virulence factor from the bacterium Lysteria monocytogenes. This domain contains seven tandem leucine-rich repeats, of which each contribute a single beta-strand that forms a continuous beta-sheet with neighboring repeats, and an N-terminal alpha-helical capping motif. Despite its modular structure, InlB folds in an equilibrium two-state manner, as reflected by the identical thermodynamic parameters obtained by monitoring its sigmoidal urea-induced unfolding transition by different spectroscopic probes. Although equilibrium two-state folding is common in alpha-helical repeat proteins, to date, InlB is the only beta-sheet-containing repeat protein for which this behavior is observed. Surprisingly, unlike other repeat proteins exhibiting equilibrium two-state folding, InlB also folds by a simple two-state kinetic mechanism lacking intermediates, aside from the effects of prolyl isomerization on the denatured state. However, like other repeat proteins, InlB also folds significantly more slowly than expected from contact order. When plotted against urea, the rate constants for the fast refolding and single unfolding phases constitute a linear chevron that, when fitted with a kinetic two-state model, yields thermodynamic parameters matching those observed for equilibrium folding. Based on these kinetic parameters, the transition state is estimated to comprise 40% of the total surface area buried upon folding, indicating that a large fraction of the native contacts are formed in the rate-limiting step to folding.

Pre-folding IkappaBalpha alters control of NF-kappaB signaling

Transcription complex components frequently show coupled folding and binding but the functional significance of this mode of molecular recognition is unclear. IkappaBalpha binds to and inhibits the transcriptional activity of NF-kappaB via its ankyrin repeat (AR) domain. The beta-hairpins in ARs 5-6 in IkappaBalpha are weakly-folded in the free protein, and their folding is coupled to NF-kappaB binding. Here, we show that introduction of two stabilizing mutations in IkappaBalpha AR 6 causes ARs 5-6 to fold cooperatively to a conformation similar to that in NF-kappaB-bound IkappaBalpha. Free IkappaBalpha is degraded by a proteasome-dependent but ubiquitin-independent mechanism, and this process is slower for the pre-folded mutants both in vitro and in cells. Interestingly, the pre-folded mutants bind NF-kappaB more weakly, as shown by both surface plasmon resonance and isothermal titration calorimetry in vitro and immunoprecipitation experiments from cells. One consequence of the weaker binding is that resting cells containing these mutants show incomplete inhibition of NF-kappaB activation; they have significant amounts of nuclear NF-kappaB. Additionally, the weaker binding combined with the slower rate of degradation of the free protein results in reduced levels of nuclear NF-kappaB upon stimulation. These data demonstrate clearly that the coupled folding and binding of IkappaBalpha is critical for its precise control of NF-kappaB transcriptional activity.

Mapping the energy landscape of repeat proteins using NMR-detected hydrogen exchange

Repeat proteins contain tandem arrays of a simple structural motif. In contrast to globular proteins, repeat proteins are stabilized only by interactions between residues that are relatively close together in the sequence, with no "long-range" interactions. Our work focuses on the tetratricopeptide repeat (TPR), a 34 amino acid helix-turn-helix motif found in tandem arrays in many natural proteins. Earlier, we reported the design and characterization of a series of consensus TPR (CTPR) proteins, which are built as arrays of multiple tandem copies of a 34 amino acid consensus sequence. Here, we present the results of extensive hydrogen exchange (HX) studies of the folding-unfolding behavior of two CTPR proteins (CTPR2 and CTPR3). We used HX to detect and characterize partially folded species that are populated at low frequency in the nominally folded state. We show that for both proteins the equilibrium folding-unfolding transition is non-two-state, but sequential, with the outermost helices showing a significantly higher probability than inner helices of being unfolded. We show that the experimentally observed unfolding behavior is consistent with the predictions of a simple Ising model, in which individual helices are treated as "spin-equivalents". The results that we present have general implications for our understanding of the thermodynamic properties of repeat proteins.

Structural insights into an equilibrium folding intermediate of an archaeal ankyrin repeat protein

Repeat proteins are widespread in nature, with many of them functioning as binding molecules in protein-protein recognition. Their simple structural architecture is used in biotechnology for generating proteins with high affinities to target proteins. Recent folding studies of ankyrin repeat (AR) proteins revealed a new mechanism of protein folding. The formation of an intermediate state is rate limiting in the folding reaction, suggesting a scaffold function of this transient state for intrinsically less stable ARs. To investigate a possible common mechanism of AR folding, we studied the structure and folding of a new thermophilic AR protein (tANK) identified in the archaeon Thermoplasma volcanium. The x-ray structure of the evolutionary much older tANK revealed high homology to the human CDK inhibitor p19(INK4d), whose sequence was used for homology search. As for p19(INK4d), equilibrium and kinetic folding analyses classify tANK to the family of sequential three-state folding proteins, with an unusual fast equilibrium between native and intermediate state. Under equilibrium conditions, the intermediate can be populated to >90%, allowing characterization on a residue-by-residue level using NMR spectroscopy. These data clearly show that the three C-terminal ARs are natively folded in the intermediate state, whereas native cross-peaks for the rest of the molecule are missing. Therefore, the formation of a stable folding unit consisting of three ARs is the necessary rate-limiting step before AR 1 and 2 can assemble to form the native state.