Protein topology affects the appearance of intermediates during the folding of proteins with a flavodoxin-like fold

The topology of a native protein influences the rate with which it is formed, but does topology affect the appearance of folding intermediates and their specific role in kinetic folding as well? This question is addressed by comparing the folding data recently obtained on apoflavodoxin from Azotobacter vinelandii with those available on all three other alpha-beta parallel proteins the kinetic folding mechanism of which has been studied, i.e. Anabaena apoflavodoxin, Fusarium solani pisi cutinase and CheY. Two kinetic folding intermediates, one on-pathway and the other off-pathway, seem to be present during the folding of proteins with an alpha-beta parallel, also called flavodoxin-like, topology. The on-pathway intermediate lies on a direct route from the unfolded to the native state of the protein involved. The off-pathway intermediate needs to unfold to allow the production of native protein. Available simulation data of the folding of CheY show the involvement of two intermediates with characteristics that resemble those of the two intermediates experimentally observed. Apparently, protein topology governs the appearance and kinetic roles of protein folding intermediates during the folding of proteins that have a flavodoxin-like fold.     

Topological Frustration in βα-Repeat Proteins: Sequence Diversity Modulates the Conserved Folding Mechanisms of α/β/α Sandwich Proteins

The thermodynamic hypothesis of Anfinsen postulates that structures and stabilities of globular proteins are determined by their amino acid sequences. Chain topology, however, is known to influence the folding reaction, in that motifs with a preponderance of local interactions typically fold more rapidly than those with a larger fraction of nonlocal interactions. Together, the topology and sequence can modulate the energy landscape and influence the rate at which the protein folds to the native conformation. To explore the relationship of sequence and topology in the folding of βα-repeat proteins, which are dominated by local interactions, we performed a combined experimental and simulation analysis on two members of the flavodoxin-like, α/β/α sandwich fold. Spo0F and the N-terminal receiver domain of NtrC (NT-NtrC) have similar topologies but low sequence identity, enabling a test of the effects of sequence on folding. Experimental results demonstrated that both response-regulator proteins fold via parallel channels through highly structured submillisecond intermediates before accessing their cis prolyl peptide bond-containing native conformations. Global analysis of the experimental results preferentially places these intermediates off the productive folding pathway. Sequence-sensitive Gō-model simulations conclude that frustration in the folding in Spo0F, corresponding to the appearance of the off-pathway intermediate, reflects competition for intra-subdomain van der Waals contacts between its N- and C-terminal subdomains. The extent of transient, premature structure appears to correlate with the number of isoleucine, leucine, and valine (ILV) side chains that form a large sequence-local cluster involving the central β-sheet and helices α2, α3, and α4. The failure to detect the off-pathway species in the simulations of NT-NtrC may reflect the reduced number of ILV side chains in its corresponding hydrophobic cluster. The location of the hydrophobic clusters in the structure may also be related to the differing functional properties of these response regulators. Comparison with the results of previous experimental and simulation analyses on the homologous CheY argues that prematurely folded unproductive intermediates are a common property of the βα-repeat motif.     

The energy landscape of modular repeat proteins: topology determines folding mechanism in the ankyrin family

Proteins consisting of repeating amino acid motifs are abundant in all kingdoms of life, especially in higher eukaryotes. Repeat-containing proteins self-organize into elongated non-globular structures. Do the same general underlying principles that dictate the folding of globular domains apply also to these extended topologies? Using a simplified structure-based model capturing a perfectly funneled energy landscape, we surveyed the predicted mechanism of folding for ankyrin repeat containing proteins. The ankyrin family is one of the most extensively studied classes of non-globular folds. The model based only on native contacts reproduces most of the experimental observations on the folding of these proteins, including a folding mechanism that is reminiscent of a nucleation propagation growth. The confluence of simulation and experimental results suggests that the folding of non-globular proteins is accurately described by a funneled energy landscape, in which topology plays a determinant role in the folding mechanism.     

Kinetic traps in the folding of beta alpha-repeat proteins: CheY initially misfolds before accessing the native conformation

The βα-repeat class of proteins, represented by the (βα)₈ barrel and the α/β/α sandwich, are among the most common structural platforms in biology. Previous studies on the folding mechanisms of these motifs have revealed or suggested that the initial event involves the submillisecond formation of a kinetically trapped species that must at least partially unfold before productive folding to the respective native conformation can occur. To test the generality of these observations, CheY, a bacterial response regulator, was subjected to an extensive analysis of its folding reactions. Although earlier studies had proposed the formation of an off-pathway intermediate, the data available were not sufficient to rule out an alternative on-pathway mechanism. A global analysis of single- and double-jump kinetic data, combined with equilibrium unfolding data, was used to show that CheY folds and unfolds through two parallel channels defined by the state of isomerization of a prolyl peptide bond in the active site. Each channel involves a stable, highly structured folding intermediate whose kinetic properties are better described as the properties of an off-pathway species. Both intermediates subsequently flow through the unfolded state ensemble and adopt the native cis-prolyl isomer prior to forming the native state. Initial collapse to off-pathway folding intermediates is a common feature of the folding mechanisms of βα-repeat proteins, perhaps reflecting the favored partitioning to locally determined substructures that cannot directly access the native conformation. Productive folding requires the dissipation of these prematurely folded substructures as a prelude to forming the larger-scale transition state that leads to the native conformation. Results from Gō-modeling studies in the accompanying paper elaborate on the topological frustration in the folding free-energy landscape of CheY.     

The underlying mechanism for the diversity of disulfide folding pathways

The disulfide folding pathway of bovine pancreatic trypsin inhibitor (BPTI) is characterized by the predominance of folding intermediates with native-like structures. Our laboratory has recently analyzed the folding pathway(s) of four 3-disulfide-containing proteins, including hirudin, potato carboxypeptidase inhibitor, epidermal growth factor, and tick anticoagulant peptide. Their folding mechanism(s) differ from that of BPTI by 1) a higher degree of heterogeneity of 1- and 2-disulfide intermediates and 2) the presence of 3-disulfide scrambled isomers as folding intermediates. To search for the underlying causes of these diversities, we conducted kinetic analyses of the reductive unfolding of these five proteins. The experiment of reductive unfolding was designed to evaluate the relative stability and interdependence of disulfide bonds in the native protein. It is demonstrated here that among these five proteins, there exists a striking correlation between the mechanism(s) of reductive unfolding and that of oxidative folding. Those proteins with their native disulfide bonds reduced in a collective and simultaneous manner exhibit both a high degree of heterogeneity of folding intermediates and the accumulation of scrambled isomers along the folding pathway. A sequential reduction of the native disulfide bonds is associated with the presence of predominant intermediates with native- like structures.     

Coevolution of function and the folding landscape: correlation with density of native contacts

The relationship between the folding landscape and function of evolved proteins is explored by comparison of the folding mechanisms for members of the flavodoxin fold. CheY, Spo0F, and NtrC have unrelated functions and low sequence homology but share an identical topology. Recent coarse-grained simulations show that their folding landscapes are uniquely tuned to properly suit their respective biological functions. Enhanced packing in Spo0F and its limited conformational dynamics compared to CheY or NtrC lead to frustration in its folding landscape. Simulation as well as experimental results correlate with the local density of native contacts for these and a sample of other proteins. In particular, protein regions of low contact density are observed to become structured late in folding; concomitantly, these dynamic regions are often involved in binding or conformational rearrangements of functional importance. These observations help to explain the widespread success of Gomacr;-like coarse-grained models in reproducing protein dynamics.     

Protein aggregation kinetics in an Escherichia coli strain overexpressing a Salmonella typhimurium CheY mutant gene.

The tendency of recombinant protein in bacteria to partition into soluble and insoluble forms is attributed, in general, to a kinetic competition between protein folding and aggregation. However, little experimental work has actually been performed in vivo on the kinetics and mechanisms of protein folding and aggregation. Results are presented here from radiolabeling experiments which monitored the kinetics of recombinant protein aggregation in actively growing cultures. The strain used was an Escherichia coli strain overexpressing a Salmonella typhimurium CheY mutant gene. The rate of CheY aggregation was found to be time dependent in that the tendency of CheY to aggregate was greater for newly translated molecules, i.e., those translated within the previous several minutes, than for molecules translated less recently. CheY protein molecules that were translated less recently continued to aggregate for several hours but at a lower rate. The movement of soluble CheY to the insoluble form was enhanced at elevated growth temperatures and inhibited by the presence of chloramphenicol. The latter observation suggests that ongoing translation facilitates the movement of soluble CheY to the insoluble form. The implications of these results for the mechanism of protein aggregation in vivo, i.e., inclusion body formation, are discussed.     

Molecular basis of co-operativity in protein folding. III. Structural identification of cooperative folding units and folding intermediates.

The hierarchical partition function formalism for protein folding developed earlier has been extended through the use of three-dimensional polar and apolar contact plots. For each amino acid residue in the protein, these plots indicate the apolar and polar surfaces that are buried from the solvent, the identity of all amino acid residues that contribute to this shielding, and the magnitude of their contributions. These contact plots are then used to examine the distribution of the free energy of stabilization throughout the protein molecule. Analysis of these data allows identification of co-operative folding units and their hierarchical levels, and the identification of partially folded intermediates with a significant probability of being populated. The overall folding/unfolding thermodynamics of 12 globular proteins, for which crystallographic and experimental thermodynamics are available, is predicted within error. An energetic classification of partially folded intermediates is presented and the results compared to those cases for which structural and thermodynamic experimental information is available. Four different types of partially folded states and their structural energies are considered. (1) Local intermediates, in which only a local region of the protein loses secondary and tertiary interactions, while the rest of the protein remains intact. (2) Global intermediates, corresponding to the standard molten globule definition, in which significant secondary structure is maintained but native-like tertiary structure contacts are disrupted. (3) Extended intermediates characterized by the existence of secondary structure elements (e.g. alpha-helices) exposed to solvent. (4) Folding intermediates in proteins with two structural domains. The structure and energetics of folding intermediates of apo-myoglobin, alpha-lactalbumin, phosphoglycerate kinase and arabinose-binding protein are considered in detail.     

Evidence for sequential barriers and obligatory intermediates in apparent two-state protein folding

Many small proteins fold fast and without detectable intermediates. This is frequently taken as evidence against the importance of partially folded states, which often transiently accumulate during folding of larger proteins. To get insight into the properties of free energy barriers in protein folding we analyzed experimental data from 23 proteins that were reported to show non-linear activation free-energy relationships. These non-linearities are generally interpreted in terms of broad transition barrier regions with a large number of energetically similar states. Our results argue against the presence of a single broad barrier region. They rather indicate that the non-linearities are caused by sequential folding pathways with consecutive distinct barriers and a few obligatory high-energy intermediates. In contrast to a broad barrier model the sequential model gives a consistent picture of the folding barriers for different variants of the same protein and when folding of a single protein is analyzed under different solvent conditions. The sequential model is also able to explain changes from linear to non-linear free energy relationships and from apparent two-state folding to folding through populated intermediates upon single point mutations or changes in the experimental conditions. These results suggest that the apparent discrepancy between two-state and multi-state folding originates in the relative stability of the intermediates, which argues for the importance of partially folded states in protein folding.     

Protein Structure Classification Based on Conserved Hydrophobic Residues

Protein folding is frequently guided by local residue interactions that form clusters in the protein core. The interactions between residue clusters serve as potential nucleation sites in the folding process. Evidence postulates that the residue interactions are governed by the hydrophobic propensities that the residues possess. An array of hydrophobicity scales has been developed to determine the hydrophobic propensities of residues under different environmental conditions. In this work, we propose a graph-theory-based data mining framework to extract and isolate protein structural features that sustain invariance in evolutionary-related proteins, through the integrated analysis of five well-known hydrophobicity scales over the 3D structure of proteins. We hypothesize that proteins of the same homology contain conserved hydrophobic residues and exhibit analogous residue interaction patterns in the folded state. The results obtained demonstrate that discriminatory residue interaction patterns shared among proteins of the same family can be employed for both the structural and the functional annotation of proteins. We obtained on the average 90 percent accuracy in protein classification with a significantly small feature vector compared to previous results in the area. This work presents an elaborate study, as well as validation evidence, to illustrate the efficacy of the method and the correctness of results reported.     

Lattice simulations of cotranslational folding of single domain proteins

We use lattice protein models and Monte Carlo simulations to study cotranslational folding of small single domain proteins. We show that the assembly of native structure begins during late extrusion stages, but final formation of native state occurs during de novo folding, when all residues are extruded. There are three main results in our study. First, for the sequences displaying two-state refolding mechanism de novo cotranslational folding pathway differs from that sampled in in vitro refolding. The change in folding pathways is due to partial assembly of native interactions during extrusion that results in different starting conditions for in vitro refolding and for de novo cotranslational folding. For small single domain proteins cotranslational folding is slower than in vitro refolding, but is generally fast enough to be completed before the release from a ribosome. Second, we found that until final stages of biosynthesis cotranslational folding is essentially equilibrium. This observation is explained by low stability of structured states for partially extruded chains. Finally, our data suggest that the proteins, which refold in vitro slowly via intermediates, complete their de novo folding after the release from a ribosome. Comparison of our lattice cotranslational simulations with recent experimental and computational studies is discussed.     

Intermediates in the folding equilibrium of repeat proteins from the TPR family

In recent decades, advances in computational methods and experimental biophysical techniques have improved our understanding of protein folding. Although some of these advances have been remarkable, the structural variability of globular proteins usually encountered makes it difficult to extract general features of their folding processes. To overcome this difficulty, experimental and computational studies of the folding of repeat (or modular) proteins are of interest. Because their native structures can be described as linear arrays of the same, repeated, supersecondary structure unit, it is possible to seek  a possibly independent behavior of the different modules without taking into account the intrinsic stability associated with different secondary structure motifs. In this work we have used a Monte Carlo-based simulation to study the folding equilibrium of four repeat proteins belonging to the tetratricopeptide repeat family. Our studies provide new insights into their energy profiles, enabling investigation about the existence of intermediate states and their relative stabilities. We have also performed structural analyses to describe the structure of these intermediates, going through the vast number of conformations obtained from the simulations. In this way, we have tried to identify the regions of each protein in which the modular structure yields a different behavior and, more specifically, regions of the proteins that can stay folded when the rest of the chain has been thermally denatured.     

Secondary structure length as a determinant of folding rate of proteins with two- and three-state kinetics

We present a simple method for determining the folding rates of two- and three-state proteins from the number of residues in their secondary structures (secondary structure length). The method is based on the hypothesis that two- and three-state foldings share a common pattern. Three-state proteins first condense into metastable intermediates, subsequent forming of alpha-helices, turns, and beta-sheets at slow rate-limiting step. The folding rate of such proteins anticorrelate with the length of these beta-secondary structures. It is also assumed that in two-state folding, rapidly folded alpha-helices and turns may facilitate formation of fleeting unobservable intermediates and thus show two-state behavior. There is an inverse relationship between the folding rate and the length of beta-sheets and loops. Our study achieves 94.0 and 88.1% correlations with folding rates determined experimentally for 21 three- and 38 two-state proteins, respectively, suggesting that protein-folding rates are determined by the secondary structure length. The kinetic kinds are selected on the basis of a competitive formation of hydrophobic collapse and alpha-structure in early intermediates.     

Subdomain competition, cooperativity, and topological frustration in the folding of CheY

The folding of multidomain proteins often proceeds in a hierarchical fashion with individual domains folding independent of one another. A large single-domain protein, however, can consist of multiple modules whose folding may be autonomous or interdependent in ways that are unclear. We used coarse-grained simulations to explore the folding landscape of the two-subdomain bacterial response regulator CheY. Thermodynamic and kinetic characterization shows the landscape to be highly analogous to the four-state landscape reported for another two-subdomain protein, T4 lysozyme. An on-pathway intermediate structured in the more stable nucleating subdomain was observed, as were transient states frustrated in off-pathway contacts prematurely structured in the weaker subdomain. Local unfolding, or backtracking, was observed in the frustrated state before the native conformation could be reached. Nonproductive frustration was attributable to competition for van der Waals contacts between the two subdomains. In an accompanying article, stopped-flow kinetic measurements support an off-pathway burst-phase intermediate, seemingly consistent with our prediction of early frustration in the folding landscape of CheY. Comparison of the folding mechanisms for CheY, T4 lysozyme, and interleukin-1β leads us to postulate that subdomain competition is a general feature of large single-domain proteins with multiple folding modules.     

1H NMR and CD evidence of the folding of the isolated ribonuclease 50-61 fragment

In our search for potential folding intermediates we have prepared and characterized the fragment of RNase A corresponding to residues 50-61. Proton chemical shift variations with temperature, addition of stabilizing (TFE) or denaturing agents (urea) provide a strong experimental basis for concluding that in aqueous solution this RNase fragment forms an alpha-helix structure similar to that in the intact RNase A crystal. This conclusion lends strong support to the idea that elements of secondary structure (mainly alpha-helices) can be formed in the absence of tertiary interactions and act as nucleation centers in the protein folding process.     

Topological and energetic factors: what determines the structural details of the transition state ensemble and "en-route" intermediates for protein folding? An investigation for small globular proteins

Recent experimental results suggest that the native fold, or topology, plays a primary role in determining the structure of the transition state ensemble, at least for small, fast-folding proteins. To investigate the extent of the topological control of the folding process, we studied the folding of simplified models of five small globular proteins constructed using a Go-like potential to retain the information about the native structures but drastically reduce the energetic frustration and energetic heterogeneity among residue-residue native interactions. By comparing the structure of the transition state ensemble (experimentally determined by Phi-values) and of the intermediates with those obtained using our models, we show that these energetically unfrustrated models can reproduce the global experimentally known features of the transition state ensembles and "en-route" intermediates, at least for the analyzed proteins. This result clearly indicates that, as long as the protein sequence is sufficiently minimally frustrated, topology plays a central role in determining the folding mechanism.     

Are the same or different amino acid residues responsible for correct and incorrect protein folding?

It has been shown for 20 proteins that amino acid residues included into the protein folding nucleus, determined experimentally, are often involved in the theoretically determined amyloidogenic fragments. For 18 proteins, Φ-values indicative of the extent of residue involvement into the folding nucleus are on average higher for amino acid residues within amyloidogenic regions. Amyloidogenic fragments were predicted for 20 proteins by two methods chosen from four on the basis of comparison of prediction of amyloidogenic regions known from experimental data. Since theoretical folding nuclei are detected by the protein three-dimensional structure and amyloidogenic regions by the protein chain primary structure, the detected regularity makes possible predictions of folding nucleation sites on the basis of amino acid sequence.     

Conservation of folding pathways in evolutionarily distant globin sequences

To test the hypothesis that the folding pathways of evolutionarily related proteins with similar three-dimensional structures but widely different sequences should be similar, the folding pathway of apoleghemoglobin has been characterized using stopped-flow circular dichroism, heteronuclear NMR pulse labeling techniques and mass spectrometry. The pathway of folding was found to differ significantly from that of a protein of the same family, apomyoglobin, although both proteins appear to fold through helical burst phase intermediates. For leghemoglobin, the burst phase intermediate exhibits stable helical structure in the G and H helices, together with a small region in the center of the E helix. The A and B helices are not stabilized until later stages of the folding process. The structure of the burst phase folding intermediate thus differs from that of apomyoglobin, in which stable helical structure is formed in the A, B, G and H helix regions.     

Structural resemblance between the families of bacterial signal-transduction proteins and of G proteins revealed by graph theoretical techniques

The first application of a novel technique for the identification of common folding motifs in proteins is presented. Using techniques derived from graph theory, developed in order to compare secondary structure motifs in proteins, we have established that there is a striking resemblance in the tertiary fold of the Salmonella typhimurium Che Y chemotaxis protein and that of the GDP-binding domain of Escherichia coli elongation factor Tu (EF Tu). These two protein structures are representatives of two major macromolecular classes: CheY is a signal-transduction protein with sequence homologies to a wide range of bacterial proteins involved in regulation of chemotaxis, membrane synthesis and sporulation; whilst EF Tu is one of a family of guanosine-nucleotide-binding proteins which include the ras oncogene proteins and signal-transducing G proteins. The similarity we have found extends far beyond the previously recognized resemblances of each protein's fold to that of a generic nucleotide-binding domain. The lack of significant sequence homology between the two classes of proteins may mean that the common fold of the two proteins constitutes a particularly stable folding motif. However, an alternative possibility is that the strong three-dimensional structural resemblance may be indicative of a remote shared common ancestry between the bacterial signal-transduction proteins and the GDP-binding proteins.     

A detailed unfolding pathway of a (beta/alpha)8-barrel protein as studied by molecular dynamics simulations

The (beta/alpha)(8)-barrel is the most common protein fold. Similar structural properties for folding intermediates of (beta/alpha)(8)-barrel proteins involved in tryptophan biosynthesis have been reported in a number of experimental studies; these intermediates have the last two beta-strands and three alpha-helices partially unfolded, with other regions of the polypeptide chain native-like in conformation. To investigate the detailed folding/unfolding pathways of these (beta/alpha)(8)-barrel proteins, temperature-induced unfolding simulations of N-(5'-phosphoribosyl)anthranilate isomerase from Escherichia coli were carried out using a special-purpose parallel computer system. Unfolding simulations at five different temperatures showed a sequential unfolding pathway comprised of several events. Early events in unfolding involved disruption of the last two strands and three helices, producing an intermediate ensemble similar to those detected in experimental studies. Then, denaturation of the first two betaalpha units and separation of the sixth strand from the fifth took place independently. The remaining central betaalphabetaalphabeta module persisted the longest during all simulations, suggesting an important role for this module as the incipient folding scaffold. Our simulations also predicted the presence of a nucleation site, onto which several hydrophobic residues condensed forming the foundation for the central betaalphabetaalphabeta module.