Folding of elongated proteins: conventional or anomalous?

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.

A disorder-induced domino-like destabilization mechanism governs the folding and functional dynamics of the repeat protein IκBα

The stability of the repeat protein IκBα, a transcriptional inhibitor in mammalian cells, is critical in the functioning of the NF-κB signaling module implicated in an array of cellular processes, including cell growth, disease, immunity and apoptosis. Structurally, IκBα is complex, with both ordered and disordered regions, thus posing a challenge to the available computational protocols to model its conformational behavior. Here, we introduce a simple procedure to model disorder in systems that undergo binding-induced folding that involves modulation of the contact map guided by equilibrium experimental observables in combination with an Ising-like Wako-Saitô-Muñoz-Eaton model. This one-step procedure alone is able to reproduce a variety of experimental observables, including ensemble thermodynamics (scanning calorimetry, pre-transitions, m-values) and kinetics (roll-over in chevron plot, intermediates and their identity), and is consistent with hydrogen-deuterium exchange measurements. We further capture the intricate distance-dynamics between the domains as measured by single-molecule FRET by combining the model predictions with simple polymer physics arguments. Our results reveal a unique mechanism at work in IκBα folding, wherein disorder in one domain initiates a domino-like effect partially destabilizing neighboring domains, thus highlighting the effect of symmetry-breaking at the level of primary sequences. The offshoot is a multi-state and a dynamic conformational landscape that is populated by increasingly partially folded ensembles upon destabilization. Our results provide, in a straightforward fashion, a rationale to the promiscuous binding and short intracellular half-life of IκBα evolutionarily engineered into it through repeats with variable stabilities and expand the functional repertoire of disordered regions in proteins.

Complex energy landscape of a giant repeat protein

Here, we reveal a remarkable complexity in the unfolding of giant HEAT-repeat protein PR65/A, a molecular adaptor for the heterotrimeric PP2A phosphatases. The repeat array ruptures at multiple sites, leading to intermediate states with noncontiguous folded subdomains. There is a dominant sequence of unfolding, which reflects a nonuniform stability distribution across the repeat array and can be rationalized by theoretical models accounting for heterogeneous contact density in the folded structure. Unfolding of certain intermediates is, however, competitive, leading to parallel unfolding pathways. The low-stability, central repeats sample unfolded conformations under physiological conditions, suggesting how folding directs function: certain regions present rigid motifs for molecular recognition, whereas others have the flexibility with which to broaden the search area, as in the fly-casting mechanism. Partial unfolding of PR65/A also impacts catalysis by altering the proximity of bound catalytic subunit and substrate. Thus, the repeat array orchestrates the assembly and activity of PP2A.

Stabilization of a protein conferred by an increase in folded state entropy

Entropic stabilization of native protein structures typically relies on strategies that serve to decrease the entropy of the unfolded state. Here we report, using a combination of experimental and computational approaches, on enhanced thermodynamic stability conferred by an increase in the configurational entropy of the folded state. The enhanced stability is observed upon modifications of a loop region in the enzyme acylphosphatase and is achieved despite significant enthalpy losses. The modifications that lead to increased stability, as well as those that result in destabilization, however, strongly compromise enzymatic activity, rationalizing the preservation of the native loop structure even though it does not provide the protein with maximal stability or kinetic foldability.

Coarse-grained models for protein folding and aggregation

Coarse-grained models for protein folding and aggregation are used to explore large dimension scales and timescales that are inaccessible to all-atom models in explicit aqueous solution. Combined with enhanced configuration search methods, these simplified models with various levels of granularity offer the possibility to determine equilibrium structures, compare folding kinetics and thermodynamics with experiments for single proteins and understand the dynamic assembly of amyloid proteins leading to neurodegenerative diseases. I shall describe recent progress in developing such models, and discuss their potentials and limitations in probing the folding and misfolding of proteins with computer simulations.

Modulation of folding kinetics of repeat proteins: interplay between intra- and interdomain interactions

Repeat proteins have unique elongated structures that, unlike globular proteins, are quite modular. Despite their simple one-dimensional structure, repeat proteins exhibit intricate folding behavior with a complexity similar to that of globular proteins. Therefore, repeat proteins allow one to quantify fundamental aspects of the biophysics of protein folding. One important feature of repeat proteins is the interfaces between the repeating units. In particular, the distribution of stabilities within and between the repeats was previously suggested to affect their folding characteristics. In this study, we explore how the interface affects folding kinetics and cooperativity by investigating two families of repeat proteins, namely, the Ankyrin and tetratricopeptide repeat proteins, which differ in the number of interfacial contacts that are formed between their units as well as in their folding behavior. By using simple topology-based models, we show that modulating the energetic strength of the interface relative to that of the repeat itself can drastically change the protein stability, folding rate, and cooperativity. By further dissecting the interfacial contacts into several subsets, we isolated the effects of each of these groups on folding kinetics. Our study highlights the importance of interface connectivity in determining the folding behavior.

Modulation of the multistate folding of designed TPR proteins through intrinsic and extrinsic factors

Tetratricopeptide repeats (TPRs) are a class of all alpha-helical repeat proteins that are comprised of 34-aa helix-turn-helix motifs. These stack together to form nonglobular structures that are stabilized by short-range interactions from residues close in primary sequence. Unlike globular proteins, they have few, if any, long-range nonlocal stabilizing interactions. Several studies on designed TPR proteins have shown that this modular structure is reflected in their folding, that is, modular multistate folding is observed as opposed to two-state folding. Here we show that TPR multistate folding can be suppressed to approximate two-state folding through modulation of intrinsic stability or extrinsic environmental variables. This modulation was investigated by comparing the thermodynamic unfolding under differing buffer regimes of two distinct series of consensus-designed TPR proteins, which possess different intrinsic stabilities. A total of nine proteins of differing sizes and differing consensus TPR motifs were each thermally and chemically denatured and their unfolding monitored using differential scanning calorimetry (DSC) and CD/fluorescence, respectively. Analyses of both the DSC and chemical denaturation data show that reducing the total stability of each protein and repeat units leads to observable two-state unfolding. These data highlight the intimate link between global and intrinsic repeat stability that governs whether folding proceeds by an observably two-state mechanism, or whether partial unfolding yields stable intermediate structures which retain sufficient stability to be populated at equilibrium.

From artificial antibodies to nanosprings: the biophysical properties of repeat proteins

In this chapter we review recent studies of repeat proteins, a class of proteins consisting of tandem arrays of small structural motifs that stack approximately linearly to produce elongated structures. We discuss the observation that, despite lacking the long-range tertiary interactions that are thought to be the hallmark of globular protein stability, repeat proteins can be as stable and as co-orperatively folded as their globular counterparts. The symmetry inherent in the structures of repeat arrays, however, means there can be many partly folded species (whether it be intermediates or transition states) that have similar stabilities. Consequently they do have distinct folding properties compared with globular proteins and these are manifest in their behaviour both at equilibrium and under kinetic conditions. Thus, when studying repeat proteins one appears to be probing a moving target: a relatively small perturbation, by mutation for example, can result in a shift to a different intermediate or transition state. The growing literature on these proteins illustrates how their modular architecture can be adapted to a remarkable array of biological and physical roles, both in vivo and in vitro. Further, their simple architecture makes them uniquely amenable to redesign-of their stability, folding and function-promising exciting possibilities for future research.

Accurate prediction of protein folding rates from sequence and sequence-derived residue flexibility and solvent accessibility

Protein folding rates vary by several orders of magnitude and they depend on the topology of the fold and the size and composition of the sequence. Although recent works show that the rates can be predicted from the sequence, allowing for high-throughput annotations, they consider only the sequence and its predicted secondary structure. We propose a novel sequence-based predictor, PFR-AF, which utilizes solvent accessibility and residue flexibility predicted from the sequence, to improve predictions and provide insights into the folding process. The predictor includes three linear regressions for proteins with two-state, multistate, and unknown (mixed-state) folding kinetics. PFR-AF on average outperforms current methods when tested on three datasets. The proposed approach provides high-quality predictions in the absence of similarity between the predicted and the training sequences. The PFR-AF's predictions are characterized by high (between 0.71 and 0.95, depending on the dataset) correlation and the lowest (between 0.75 and 0.9) mean absolute errors with respect to the experimental rates, as measured using out-of-sample tests. Our models reveal that for the two-state chains inclusion of solvent-exposed Ala may accelerate the folding, while increased content of Ile may reduce the folding speed. We also demonstrate that increased flexibility of coils facilitates faster folding and that proteins with larger content of solvent-exposed strands may fold at a slower pace. The increased flexibility of the solvent-exposed residues is shown to elongate folding, which also holds, with a lower correlation, for buried residues. Two case studies are included to support our findings.

Residue-specific analysis of frustration in the folding landscape of repeat beta/alpha protein apoflavodoxin

Flavodoxin adopts the common repeat beta/alpha topology and folds in a complex kinetic reaction with intermediates. To better understand this reaction, we analyzed a set of Desulfovibrio desulfuricans apoflavodoxin variants with point mutations in most secondary structure elements by in vitro and in silico methods. By equilibrium unfolding experiments, we first revealed how different secondary structure elements contribute to overall protein resistance to heat and urea. Next, using stopped-flow mixing coupled with far-UV circular dichroism, we probed how individual residues affect the amount of structure formed in the experimentally detected burst-phase intermediate. Together with in silico folding route analysis of the same point-mutated variants and computation of growth in nucleation size during early folding, computer simulations suggested the presence of two competing folding nuclei at opposite sides of the central beta-strand 3 (i.e., at beta-strands 1 and 4), which cause early topological frustration (i.e., misfolding) in the folding landscape. Particularly, the extent of heterogeneity in folding nuclei growth correlates with the in vitro burst-phase circular dichroism amplitude. In addition, phi-value analysis (in vitro and in silico) of the overall folding barrier to apoflavodoxin's native state revealed that native-like interactions in most of the beta-strands must form in transition state. Our study reveals that an imbalanced competition between the two sides of apoflavodoxin's central beta-sheet directs initial misfolding, while proper alignment on both sides of beta-strand 3 is necessary for productive folding.

Exploring the folding energy landscape of a series of designed consensus tetratricopeptide repeat proteins

Repeat proteins contain short, tandem arrays of simple structural motifs (20-40 aa). These stack together to form nonglobular structures that are stabilized by short-range interactions from residues close in primary sequence. Unlike globular proteins, they have few, if any, long-range nonlocal stabilizing interactions. One ubiquitous repeat is the tetratricopeptide motif (TPR), a 34-aa helix-turn-helix motif. In this article we describe the folding kinetics of a series of 7 designed TPR proteins that are assembled from arraying identical designed consensus repeats (CTPRan). These range from the smallest 2-repeat protein to a large 10-repeat protein (approximately 350 aa). In particular, we describe how the energy landscape changes with the addition of repeat units. The data reveal that although the CTPRa proteins have low local frustration, their highly symmetric, modular native structure is reflected in their multistate kinetics of unfolding and folding. Moreover, although the initial folding of all CTPRan proteins involves a nucleus with similar solvent accessibility, their subsequent folding to the native structure depends directly on repeat number. This corresponds to an increasingly complex landscape that culminates in CTPRa10 populating a misfolded, off-pathway intermediate. These results extend our current understanding of the malleable folding pathways of repeat proteins and highlight the consequences of adding identical repeats to the energy landscape.

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.