Folding versus aggregation: polypeptide conformations on competing pathways

Protein aggregation has now become recognised as an important and generic aspect of protein energy landscapes. Since the discovery that numerous human diseases are caused by protein aggregation, the biophysical characterisation of misfolded states and their aggregation mechanisms has received increased attention. Utilising experimental techniques and computational approaches established for the analysis of protein folding reactions has ensured rapid advances in the study of pathways leading to amyloid fibrils and amyloid-related aggregates. Here we describe recent experimental and theoretical advances in the elucidation of the conformational properties of dynamic, heterogeneous and/or insoluble protein ensembles populated on complex, multidimensional protein energy landscapes. We discuss current understanding of aggregation mechanisms in this context and describe how the synergy between biochemical, biophysical and cell-biological experiments are beginning to provide detailed insights into the partitioning of non-native species between protein folding and aggregation pathways.     

Improved Gō-like models demonstrate the robustness of protein folding mechanisms towards non-native interactions

The use of simple theoretical models has provided a considerable contribution to our present understanding of the means by which proteins adopt their native fold from the plethora of available unfolded states. A common assumption in building computationally tractable models has been the neglect of stabilizing non-native interactions in the class of models described as "Gō-like." The focus of this study is the characterization of the folding of a number of proteins via a Gō-like model, which aims to map a maximal amount of information reflecting the protein sequence onto a "minimalist" skeleton. This model is shown to contain sufficient information to reproduce the folding transition states of a number of proteins, including topologically analogous proteins that fold via different transition states. Remarkably, these models also demonstrate consistency with the general features of folding transition states thought to be stabilized by non-native interactions. This suggests that native interactions are the primary determinant of most protein folding transition states, and that non-native interactions lead only to local structural perturbations. A prediction is also included for an asymmetrical folding transition state of bacteriophage lambda protein W, which has yet to be subjected to experimental characterization.     

Molten globule versus variety of intermediates: influence of anions on pH-denatured apomyoglobin.

The molten globule state was shown to be the third thermodynamic state of protein molecules in addition to their native and unfolded states. On the other hand, it was reported that optical and hydrodynamic properties of pH-denatured apomyoglobin depend on the nature of anions added to the protein solution. This observation was used to conclude that there are many 'partly folded' intermediates between the native and unfolded states rather than one distinct molten globule state. However, little is known on the structures of pH-denatured apomyoglobin in the presence of different anions. Two tyrosine residues in horse apomyoglobin have been successively modified by the reaction with tetranitromethane. This approach was employed to measure the distances between tryptophans and modified tyrosines in different states of apomyoglobin by the method of direct energy transfer. Experimental data show that the distance between the middle of the A-helix and the beginning of the G-helix and/or the end of the H-helix in 'anion-induced' states are very close to those in the native holo- and apomyoglobins. This suggests that the AGH helical complex, being the most structured part of apomyoglobin in the molten globule state, exists also in pH-denatured apomyoglobin in the presence of different anions. Consequently, all non-native forms of apomyoglobin studied so far share the common important feature of its native structure.     

Structural examination of phi-value analysis in protein folding

Protein folding intermediates and transition states are commonly characterized using a protein engineering procedure (Phi-value analysis) based on several assumptions, including (1) intermediates and transition states have native-like conformations and (2) single mutations from larger hydrophobic residues to smaller ones do not perturb their structures. Although Phi-value analysis has been widely used, these assumptions have not been tested to date because of the lack of high-resolution structures of intermediates and transition states. We recently have determined the structure of a folding intermediate for a four-helix bundle protein (Rd-apocytochrome b(562)) using NMR. The intermediate has the N-terminal helix unfolded. The other three helices fold in a native-like topology with extensive non-native hydrophobic interactions. Here, we have determined the Phi values for 14 hydrophobic core residues, including those with significant non-native interactions. All of the Phi values are in the normal range from 0 to 1, indicating that these non-native interactions cannot be identified by the common Phi-value analysis, and therefore, the first assumption is not valid for this intermediate. We also determined the structure of a mutant (F65A) of the intermediate and found that the structure of the intermediate is not perturbed by the mutation, supporting the second assumption. Together, these results suggest that Phi-value analysis may be valid for characterizing the energetics of the interactions between the mutated residue and others, but not for determining the detailed structures of intermediates and transition states because non-native interactions may exist and may not be identifiable by the common Phi-value analysis.     

Specific non-native hydrophobic interactions in a hidden folding intermediate: implications for protein folding

Structures of intermediates and transition states in protein folding are usually characterized by amide hydrogen exchange and protein engineering methods and interpreted on the basis of the assumption that they have native-like conformations. We were able to stabilize and determine the high-resolution structure of a partially unfolded intermediate that exists after the rate-limiting step of a four-helix bundle protein, Rd-apocyt b(562), by multidimensional NMR methods. The intermediate has partial native-like secondary structure and backbone topology, consistent with our earlier native state hydrogen exchange results. However, non-native hydrophobic interactions exist throughout the structure. These and other results in the literature suggest that non-native hydrophobic interactions may occur generally in partially folded states. This can alter the interpretation of mutational protein engineering results in terms of native-like side chain interactions. In addition, since the intermediate exists after the rate-limiting step and Rd-apocyt b(562) folds very rapidly (k(f) approximately 10(4) s(-1)), these results suggest that non-native hydrophobic interactions, in the absence of topological misfolding, are repaired too rapidly to slow folding and cause the accumulation of folding intermediates. More generally, these results illustrate an approach for determining the high-resolution structure of folding intermediates.     

A native disulfide stabilizes non-native helical structures in partially folded states of equine β-lactoglobulin

Equine β-lactoglobulin (ELG) assumes non-native helices during refolding and in partially folded states. Previously, circular dichroism (CD) combined with site-directed mutagenesis identified helical regions in the acid- and cold-denatured states of ELG. It is also known that a fragment of ELG, CHIBL (residues 88-142), has a structure similar to that of the cold-denatured state. For the study reported herein, the structure of a shorter fragment, CHIBLΔF (residues 97-142), was investigated by CD and nuclear magnetic resonance spectroscopy. The secondary chemical shifts clearly showed that non-native α-helices are present in two different regions, residues 98-107 and 114-135, and are connected by a native disulfide bond. The CD spectra of two peptides that correspond to the helical regions are characterized by weak helical signatures, and the sum of their CD spectra is nearly the same as the spectrum of disulfide-reduced CHIBLΔF. Therefore, the non-native helices are stabilized by the disulfide, and non-native helix formation may occur only during the refolding of the disulfide-intact protein. Supporting this conclusion is the observation that tear lipocalin, a homologue of ELG that lacks the disulfide, does not form non-native helices during folding.     

Computational methods for generating models of denatured and partially folded proteins

Partially folded and denatured proteins can give important insights into protein folding, misfolding, and aggregation. Such non-native states of proteins are however very difficult to characterise in detail as they are dynamic, heterogeneous systems comprising of ensembles of interconverting conformers. This article describes methods that produce models for non-native proteins in atomic detail. A variety of molecular dynamics based protocols are discussed together with some recent procedures that include restraints from experimental data. These models provide an important framework for interpreting experimental data from studies of non-native states using nuclear magnetic resonance spectroscopy, fluorescence, circular dichroism, and small angle scattering techniques.     

Electrostatic interactions in the denatured state and in the transition state for protein folding: effects of denatured state interactions on the analysis of transition state structure

The development of electrostatic interactions during the folding of the N-terminal domain of the ribosomal protein L9 (NTL9) is investigated by pH-dependent rate equilibrium free energy relationships. We show that Asp8, among six acidic residues, is involved in non-native, electrostatic interactions with K12 in the transition state for folding as well as in the denatured state. The perturbed native state pK(a) of D8 (pK(a) = 3.0) appears to be maintained through non-native interactions in both the transition state and the denatured state. Mutational effects on the stability of the transition state for protein (un)folding are often analyzed in respect to change in ground states. Thus, the interpretation of transition state analysis critically depends on an understanding of mutational effects on both the native and denatured state. Increasing evidence for structurally biased denatured states under physiological conditions raises concerns about possible denatured state effects on folding studies. We show that the structural interpretation of transition state analysis can be altered dramatically by denatured state effects.     

Structural and dynamic characteristics of a partially folded state of ubiquitin revealed by hydrogen exchange mass spectrometry

Structural and dynamic properties of a partially folded conformation (A-state) of ubiquitin are studied using amide hydrogen exchange in solution (HDX) and mass spectrometric detection. A clear distinction between the native state of the protein and the A-state can be made when HDX is carried out in a semicorrelated regime. Convoluted exchange patterns are interpreted with the aid of HDX simulations in a three-state system (highly structured, partially unstructured, and fully unstructured states). The data clearly indicate a highly dynamic character of the non-native state. Furthermore, combination of HDX and protein ion fragmentation in the gas phase [by means of collision-induced dissociation (CAD)] is used to evaluate the conformational stability of various protein segments specifically in the molten globular state. Chain flexibility appears to be distributed very unevenly in this non-native conformation. The highest degree of structural disorder is displayed by the C-terminal segment (Gly(53)-Gly(76)), which was previously suggested to form a transient alpha-helix. The least dynamic segment of ubiquitin in the A-state is Thr(9)-Glu(18) (which was previously suggested to form a stable nativelike beta-strand), with the adjacent segments exhibiting somewhat diminished conformational stability. The study also demonstrates the power of mass spectrometry as a tool in providing conformer-specific information about the structure and dynamics of both native and non-native protein states coexisting in solution under equilibrium.     

Engineering a compact non-native state of intestinal fatty acid-binding protein

The last three C-terminal residues (129-131) of intestinal fatty acid-binding protein (IFABP) participate in four main-chain hydrogen bonds and two electrostatic interactions to sequentially distant backbone and side-chain atoms. To assess if these interactions are involved in the final adjustment of the tertiary structure during folding, we engineered an IFABP variant truncated at residue 128. An additional mutation, Trp-6-->Phe, was introduced to simplify the conformational analysis by optical methods. Although the changes were limited to a small region of the protein surface, they resulted in an IFABP with altered secondary and tertiary structure. Truncated IFABP retains some cooperativity, is monomeric, highly compact, and has the molecular dimensions and shape of the native protein. Our results indicated that residues 129-131 are part of a crucial conformational determinant in which several long-range interactions, essential for the acquisition of the native state, are established. This work suggests that carefully controlled truncation can populate equilibrium non-native states under physiological conditions. These non-native states hold a great promise as experimental models for protein folding.     

Effects of a helix substitution on the folding mechanism of bovine alpha-lactalbumin

The structure, stability, and unfolding-refolding kinetics of a chimeric protein, in which the amino acid sequence of the flexible loop region (residues 105-110) comes from equine lysozyme and the remainder of the sequence comes from bovine alpha-lactalbumin were studied by circular dichroism spectroscopy and stopped-flow measurements, and the results were compared with those of bovine alpha-lactalbumin. The substitution of the flexible loop in bovine alpha-lactalbumin with the helix D of equine lysozyme destabilizes the molten globule state, although the native state is significantly stabilized by substitution of the flexible loop region. The kinetic refolding and unfolding experiments showed that the chimeric protein refolds significantly faster and unfolds substantially slower than bovine alpha-lactalbumin. To characterize the transition state between the molten globule and the native states, we investigated the guanidine hydrochloride concentration dependence of the rate constants of refolding and unfolding. Despite the significant differences in the stabilities of both the molten globule and native states between the chimeric protein and bovine alpha-lactalbumin, the free energy level of the transition state is not affected by the amino acid substitution in the flexible loop region. Our results suggest that the destabilization in the molten globule state of the chimeric protein is caused by the disruption of the non-native interaction in the flexible loop region and that the disruption of the non-native interaction reduces the free energy barrier of refolding. We conclude that the non-native interaction in the molten globule state may act as a kinetic trap for the folding of alpha-lactalbumin.     

Protein folding assisted by chaperones

Molecular chaperones are one of the most important cell defense mechanisms against protein aggregation and misfolding. These specialized proteins bind non-native states of other proteins and assist them in reaching a correctly folded and functional conformation. Chaperones also participate in protein translocation by membranes, in the stabilization of unstable protein conformers and regulatory factors, in the delivery of substrates for proteolysis and in the recovery of proteins from aggregates.     

A collapsed state functions to self-chaperone RNA folding into a native ribonucleoprotein complex

Most large RNAs achieve their active, native structures only as complexes with one or more cofactor proteins. By varying the Mg(2+) concentration, the catalytic core of the bI5 group I intron RNA can be manipulated into one of three states, expanded, collapsed or native, or into balanced equilibria between these states. Under near-physiological conditions, the bI5 RNA folds rapidly to a collapsed but non-native state. Hydroxyl radical footprinting demonstrates that assembly with the CBP2 protein cofactor chases the RNA from the collapsed state to the native state. In contrast, CBP2 also binds to the RNA in the expanded state to form many non-native interactions. This structural picture is reinforced by functional splicing experiments showing that RNA in an expanded state forms a non-productive, kinetically trapped complex with CBP2. Thus, rapid folding to the collapsed state functions to self-chaperone bI5 RNA folding by preventing premature interaction with its protein cofactor. This productive, self-chaperoning role for RNA collapsed states may be especially important to avert misassembly of large multi-component RNA-protein machines in the cell.     

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.     

Mechanism of chaperone function in small heat shock proteins: dissociation of the HSP27 oligomer is required for recognition and binding of destabilized T4 lysozyme

Mammalian small heat shock proteins (sHSP) form polydisperse and dynamic oligomers that undergo equilibrium subunit exchange. Current models of their chaperone activity hypothesize that recognition and binding of protein non-native states involve changes in the oligomeric state. The equivalent thermodynamic representation is a set of three coupled equilibria that includes the sHSP oligomeric equilibrium, the substrate folding equilibrium, and the equilibrium binding between the sHSP and the substrate non-native states. To test this hypothesis and define the binding-competent oligomeric state of human Hsp27, we have perturbed the two former equilibria and quantitatively determined the consequences on binding. The substrate is a set of T4 lysozyme (T4L) mutants that bind under conditions that favor the folded state over the unfolded state by 10(2)-10(4)-fold. The concentration-dependent oligomer equilibrium of Hsp27 was perturbed by mutations that alter the relative stability of two major oligomeric states including phosphorylation-mimicking mutations that result in the dissociation to a small multimer over a wide range of concentrations. Correlation of binding isotherms with size exclusion chromatography analysis of the Hsp27 oligomer equilibrium demonstrates that the multimer is the binding-competent state. Binding occurs through two modes, each characterized by different affinity and number of binding sites, and results in T4L.Hsp27 complexes of different hydrodynamic properties. Mutants of the Hsp27 phosphorylation mimic that reverse the reduction in oligomer size also reduce the extent of T4L binding. Taken together, these results suggest a central role for the oligomeric equilibrium in regulating the chaperone activity of sHSP. The mutants identify sequence features important for modulating this equilibrium.     

Chaperonin-mediated protein folding: fate of substrate polypeptide

Chaperonins are megadalton ring assemblies that mediate essential ATP-dependent assistance of protein folding to the native state in a variety of cellular compartments, including the mitochondrial matrix, the eukaryotic cytosol, and the bacterial cytoplasm. Structural studies of the bacterial chaperonin, GroEL, both alone and in complex with its co-chaperonin, GroES, have resolved the states of chaperonin that bind and fold non-native polypeptides. Functional studies have resolved the action of ATP binding and hydrolysis in driving the GroEL-GroES machine through its folding-active and binding-active states, respectively. Yet the exact fate of substrate polypeptide during these steps is only poorly understood. For example, while binding involves multivalent interactions between hydrophobic side-chains facing the central cavity of GroEL and exposed hydrophobic surfaces of the non-native protein, the structure of any polypeptide substrate while bound to GroEL remains unknown. It is also unclear whether binding to an open GroEL ring is accompanied by structural changes in the non-native substrate, in particular whether there is an unfolding action. As a polypeptide-bound ring becomes associated with GroES, do the large rigid-body movements of the GroEL apical domains serve as another source of a potential unfolding action? Regarding the encapsulated folding-active state, how does the central cavity itself influence the folding trajectory of a substrate? Finally, how do GroEL and GroES serve, as recently recognized, to assist the folding of substrates too large to be encapsulated inside the machine? Here, such questions are addressed with the findings available to date, and means of further resolving the states of chaperonin-associated polypeptide are discussed.     

Unfolding studies of tissue transglutaminase

Activation of tissue transglutaminase by calcium involves a conformational change which allows exposition of the active site to the substrate via movements of domains 3 and 4 that lead to an increase of the inter-domain distance. The inhibitor GTP counteracts these changes. Here we investigate the possible existence of non-native conformational states still compatible with the enzyme activity produced by chemical and thermal perturbations. The results indicate that chemical denaturation is reversible at low guanidine concentrations but irreversible at high concentrations of guanidine. Indeed, at low guanidine concentrations tissue TG-ase exists in a non-native state which is still affected by the ligands as in the native form. In contrast, thermal unfolding is always irreversible, with aggregation and protein self-crosslinkage in the presence of calcium. DSC thermograms of the native protein in the absence of ligands consist of two partly overlapped transitions, which weaken in the presence of calcium and merge together and strengthen in the presence of GTP. Overall, the present work shows, for the first time, the reversible denaturation of a TG-ase isoenzyme and suggests the possibility that also in in vivo, the enzyme may acquire non-native conformations relevant to its patho-physiological functions.     

Using scores derived from statistical coupling analysis to distinguish correct and incorrect folds in de-novo protein structure prediction

Distinguishing native from non-native folds remains a challenging problem for protein structure prediction. We describe a method, SCA-distance scoring, based on results from statistical coupling analysis which discriminates between native and non-native folds produced by a de novo protein structure prediction method for four out of five test proteins. The method is particularly good at discriminating non-native folds which are close in RMSD to the true fold but contain a change in an internal structural element. SCA-distance scoring is a useful addition to the tools available for distinguishing native from non-native folds in protein structure prediction.     

NOE data demonstrating a compact unfolded state for an SH3 domain under non-denaturing conditions.

The N-terminal SH3 domain of drk (drkN SH3) is unstable, existing in equilibrium between a folded state (Fexch) and an unfolded state (Uexch) under non-denaturing buffer conditions. Using a15N/2H-labeled sample, long range amide NOEs can be observed in the Uexchstate as a result of reduced relaxation, in some cases correlating protons over 40 residues apart. These long range NOEs disappear upon addition of 2 M guanidinium chloride, demonstrating that there are substantial differences between the Uexchand the guanidine denatured states. Calculations using the long range NOEs of the Uexchstate yield highly compact structures having non-native turns and a non-native buried tryptophan residue. These structures agree with experimental stopped-flow fluorescence data and analytical ultracentrifugation results. Since protein stability depends on the structural and dynamic properties of both the folded and unfolded states, this study provides insights into the stability of the drkN SH3 domain. These results provide the first strong NOE-based evidence for compact unfolded states of proteins and suggest that some unfolded states under physiological conditions have specific interactions leading to compact structures.     

Denaturation-renaturation of the fibrin-stabilizing factor XIII a-chain isolated from human placenta. Properties of the native and reconstituted protein

The denaturation-renaturation transition between the native and unfolded states of the dimeric blood coagulation factor XIIIa has been examined by far-UV circular dichroism, fluorescence spectroscopy, activity measurements, sedimentation equilibrium analysis, and size exclusion high performance liquid chromatography. Guanidine hydrochloride and urea-dependent denaturation in the absence and in the presence of 5mM dithioerythritol or glutathione (5mM GSH) exhibit biphasic transitions. The first stage represents a sharp transition characterized by a change in secondary structure without subunit dissociation. This step is accompanied by the irreversible loss of biological activity. The second transition reflects the dissociation and complete unfolding of the protein to a random coil. After loss of biological activity no reactivation can be accomplished under any of the following conditions: (i) denaturation and renaturation under reducing or non-reducing conditions, (ii) variation of the protein concentration and temperature, (iii) addition of specific ligands (Ca2+, substrate), (iv) presence of stabilizing and/or destabilizing agents. Attempts to renature the protein under standard conditions (0.1 M Tris/HCl pH 7.5-9.0, 5mM DTE, 5mM EDTA) lead to refolding intermediates which exhibit a strong tendency to aggregate. A soluble product of reconstitution can be obtained by refolding at low protein concentration, low temperature, and in the presence of small amounts of destabilizing agents such as arginine or urea in the renaturation buffer at pH 7.5 to 9. The spectroscopic and hydrodynamic characterization of the partially reconstituted (non-native inactive) protein shows that partially reconstituted factor XIIIa exhibits the fluorescence properties and the dimeric structure of the native protein.(ABSTRACT TRUNCATED AT 250 WORDS)