Multi-state unfolding of the alpha subunit of tryptophan synthase, a TIM barrel protein: insights into the secondary structure of the stable equilibrium intermediates by hydrogen exchange mass spectrometry

The urea-induced unfolding of the alpha subunit of tryptophan synthase (alphaTS) from Escherichia coli, an eight-stranded (beta/alpha)(8) TIM barrel protein, has been shown to involve two stable equilibrium intermediates, I1 and I2, well populated at approximately 3 M and 5 M urea, respectively. The characterization of the I1 intermediate by circular dichroism (CD) spectroscopy has shown that I1 retains a significant fraction of the native ellipticity; the far-UV CD signal for the I2 species closely resembles that of the fully unfolded form. To obtain detailed insight into the disruption of secondary structure in the urea-induced unfolding process, a hydrogen exchange-mass spectrometry study was performed on alphaTS. The full-length protein was destabilized in increasing concentration of urea, the amide hydrogen atoms were pulse-labeled with deuterium, the labeled samples were quenched in acid and the products were analyzed by electrospray ionization mass spectrometry. Consistent with the CD results, the I1 intermediate protects up to approximately 129 amide hydrogen atoms against exchange while the I2 intermediate offers no protection. Electrospray ionization mass spectrometry analysis of the peptic fragments derived from alphaTS labeled at 3 M urea indicates that most of the region between residues 12-130, which constitutes the first four beta strands and three alpha helices, (beta/alpha)(1-3)beta(4), is structured. The (beta/alpha)(1-3)beta(4) module appears to represent the minimum sub-core of stability of the I1 intermediate. A 4+2+2 folding model is proposed as a likely alternative to the earlier 6+2 folding mechanism for alphaTS.     

An obligatory intermediate controls the folding of the alpha-subunit of tryptophan synthase, a TIM barrel protein

The proposed kinetic folding mechanism of the alpha-subunit of tryptophan synthase (alphaTS), a TIM barrel protein, displays multiple unfolded and intermediate forms which fold through four parallel pathways to reach the native state. To obtain insight into the secondary structure that stabilizes a set of late, highly populated kinetic intermediates, the refolding of urea-denatured alphaTS from Escherichia coli was monitored by pulse-quench hydrogen exchange mass spectrometry. Following dilution from 8 M urea, the protein was pulse-labeled with deuterium, quenched with acid and mass analyzed by electrospray ionization mass spectrometry (ESI-MS). Hydrogen bonds that form prior to the pulse of deuterium offer protection against exchange and, therefore, retain protons at the relevant amide bonds. Consistent with the proposed refolding model, an intermediate builds up rapidly and decays slowly over the first 100 seconds of folding. ESI-MS analysis of the peptic fragments derived from alphaTS mass-labeled and quenched after two seconds of refolding indicates that the pattern of protection of the backbone amide hydrogens in this transient intermediate is very similar to that observed previously for the equilibrium intermediate of alphaTS highly populated at 3 M urea. The protection observed in a contiguous set of beta-strands and alpha-helices in the N terminus implies a significant role for this sub-domain in directing the folding of this TIM barrel protein.     

A cis-prolyl peptide bond isomerization dominates the folding of the alpha subunit of Trp synthase,a TIM barrel protein

The cis/trans isomerization of prolyl peptide bonds has been suggested to dominate the folding of the alpha subunit of tryptophan synthase from Escherichia coli (alphaTS). To test the role of the unique cis isomer between Asp27 and Pro28, the folding properties of P28A, P28G and G(3)P28G, a three-glycine insertion mutant between Asp27 and Gly28, were investigated using urea as a denaturant. Circular dichroism analysis demonstrated that none of the mutations perturb the secondary structure significantly, although the aromatic side-chain packing is altered for P28A and P28G. All three mutant proteins inherited the three-state thermodynamic behavior observed in wild-type alphaTS, ensuring that the fundamental features of the energy surface are intact. Kinetic studies showed that neither alanine nor glycine substitutions at Pro28 results in the elimination of any slow-refolding phases. By contrast, the G(3)P28G mutant eliminates the fastest of the slow-refolding phases and one of the two unfolding phases. Double-jump experiments on G(3)P28G confirm the assignment of the missing refolding phase to the isomerization of the Asp27-Pro28 peptide bond. These results imply that the local stability conveyed by the tight, overlapping turns containing the cis peptide bond is sufficient to favor the cis isomer for several non-prolyl residues. The free energy required to drive the isomerization reaction is provided by the formation of the stable intermediate, demonstrating that the acquisition of structure and stability is required to induce subsequent rate-limiting steps in the folding of alphaTS.     

Specific structure appears at the N terminus in the sub-millisecond folding intermediate of the alpha subunit of tryptophan synthase, a TIM barrel protein

Competing views of the products of sub-millisecond folding reactions observed in many globular proteins have been ascribed either to the formation of discrete, partially folded states or to the random collapse of the unfolded chain under native-favoring conditions. To test the validity of these alternative interpretations for the stopped-flow burst-phase reaction in the (betaalpha)8, TIM barrel motif, a series of alanine replacements were made at five different leucine or isoleucine residues in the alpha subunit of tryptophan synthase (alphaTS) from Escherichia coli. This protein has been proposed to fold, in the sub-millisecond time range, to an off-pathway intermediate with significant stability and approximately 50% of the far-UV circular dichroism (CD) signal of the native conformation. Individual alanine replacements at any of three isoleucine or leucine residues in either alpha1, beta2 or beta3 completely eliminate the off-pathway species. These variants, within 5 ms, access an intermediate whose properties closely resemble those of an on-pathway equilibrium intermediate that is highly populated at moderate urea concentrations in wild-type alphaTS. By contrast, alanine replacements for leucine residues in either beta4 or beta6 destabilize but preserve the off-pathway, burst-phase species. When considered with complementary thermodynamic and kinetic data, this mutational analysis demonstrates that the sub-millisecond appearance of CD signal for alphaTS reflects the acquisition of secondary structure in a distinct thermodynamic state, not the random collapse of an unfolded chain. The contrasting results for replacements in the contiguous alpha1/beta2/beta3 domain and the C-terminal beta4 and beta6 strands imply a heterogeneous structure for the burst-phase species. The alpha1/beta2/beta3 domain appears to be tightly packed, and the C terminus appears to behave as a molten-globule-like structure whose folding is tightly coupled to that of the alpha1/beta2/beta3 domain.     

Probing minimal independent folding units in dihydrofolate reductase by molecular dissection.

Molecular dissection was employed to identify minimal independent folding units in dihydrofolate reductase (DHFR) from Escherichia coli. Eight overlapping fragments of DHFR, spanning the entire sequence and ranging in size from 36 to 123 amino acids, were constructed by chemical cleavage. These fragments were designed to examine the effect of tethering multiple elements of secondary structure on folding and to test if the secondary structural domains represent autonomous folding units. CD and fluorescence spectroscopy demonstrated that six fragments containing up to a total of seven alpha-helices or beta-strands and, in three cases, the adenine binding domain (residues 37-86), are largely disordered. A stoichiometric mixture of the two fragments comprising the large discontinuous domain, 1-36 and 87-159, also showed no evidence for folding beyond that observed for the isolated fragments. A fragment containing residues 1-107 appears to have secondary and tertiary structure; however, spontaneous self-association made it impossible to determine if this structure solely reflects the behavior of the monomeric form. In contrast, a monomeric fragment spanning residues 37-159 possesses significant secondary and tertiary structure. The urea-induced unfolding of fragment 37-159 in the presence of 0.5 M ammonium sulfate was found to be a well-defined, two-state process. The observation that fragment 37-159 can adopt a stable native fold with unique, aromatic side-chain packing is quite striking because residues 1-36 form an integral part of the structural core of the full-length protein.     

p25alpha is flexible but natively folded and binds tubulin with oligomeric stoichiometry

p25alpha is a 219-residue protein which stimulates aberrant tubulin polymerization and is implicated in a variety of other functions. The protein has unusual secondary structure involving significant amounts of random coil, and binding to microtubules is accompanied by a large structural change, suggesting a high degree of plasticity. p25alpha has been proposed to be natively unfolded, so that folding is coupled to interaction with its physiological partners. Here we show that recombinant human p25alpha is folded under physiological conditions, since it has a well structured and solvent-sequestered aromatic environment and considerable chemical shift dispersion of amide and aliphatic protons. With increasing urea concentrations, p25alpha undergoes clear spectral changes suggesting significant loss of structure. p25alpha unfolds cooperatively in urea according to a simple two-state transition with a stability in water of approximately 5 kcal/mol. The protein behaves as a monomer and refolds with a transient on-pathway folding intermediate. However, high sensitivity to proteolytic attack and abnormal gel filtration migration behavior suggests a relatively extended structure, possibly organized in distinct domains. A deletion mutant of p25alpha lacking residues 3-43 also unfolds cooperatively and with similar stability, suggesting that the N-terminal region is largely unstructured. Both proteins undergo significant loss of structure when bound to monomeric tubulin. The stoichiometry of binding is estimated to be 3-4 molecules of tubulin per p25alpha and is not significantly affected by the deletion of residues 3-43. In conclusion, we dismiss the proposal that p25alpha is natively unfolded, although the protein is relatively flexible. This flexibility may be linked to its tubulin-binding properties.     

Folding mechanism of indole-3-glycerol phosphate synthase from Sulfolobus solfataricus: a test of the conservation of folding mechanisms hypothesis in (beta(alpha))(8) barrels

As a test of the hypothesis that folding mechanisms are better conserved than sequences in TIM barrels, the equilibrium and kinetic folding mechanisms of indole-3-glycerol phosphate synthase (sIGPS) from the thermoacidophilic archaebacterium Sulfolobus solfataricus were compared to the well-characterized models of the alpha subunit of tryptophan synthase (alphaTS) from Escherichia coli. A multifaceted approach combining urea denaturation and far-UV circular dichroism, tyrosine fluorescence total intensity, and tyrosine fluorescence anisotropy was employed. Despite a sequence identity of only 13%, a stable intermediate (I) in sIGPS was found to be similar to a stable intermediate in alphaTS in terms of its thermodynamic properties and secondary structure. Kinetic experiments revealed that the fastest detectable folding event for sIGPS involves a burst-phase (<5ms) reaction that leads directly to the stable intermediate. The slower of two subsequent phases reflects the formation/disruption of an off-pathway dimeric form of I. The faster phase reflects the conversion of I to the native state and is limited by folding under marginally stable conditions and by isomerization or rearrangement under strongly folding conditions. By contrast, alphaTS is thought to fold via an off-pathway burst-phase intermediate whose unfolding controls access to a set of four on-pathway intermediates that comprise the stable equilibrium intermediate. At least three proline isomerization reactions are known to limit their interconversions and lead to a parallel channel mechanism. The simple sequential mechanism deduced for sIGPS reflects the dominance of the on-pathway burst-phase intermediate and the absence of prolyl residues that partition the stable intermediate into kinetically distinguishable species. Comparison of the results for sIGPS and alphaTS demonstrates that the thermodynamic properties and the final steps of the folding reaction are better conserved than the early events. The initial events in folding appear to be more sensitive to the sequence differences between the two TIM barrel proteins.     

Molecular dissection of the folding mechanism of the alpha subunit of tryptophan synthase: an amino-terminal autonomous folding unit controls several rate-limiting steps in the folding of a single domain protein.

The alpha subunit of tryptophan synthase (alphaTS) from Escherichia coli is a 268-residue 8-stranded beta/alpha barrel protein. Two autonomous folding units, comprising the first six strands (residues 1-188) and the last two strands (residues 189-268), have been previously identified in this single structural domain protein by tryptic digestion [Higgins, W., Fairwell, T., and Miles, E. W. (1979) Biochemistry 18, 4827-4835]. The larger, amino-terminal fragment, alphaTS(1-188), was overexpressed and independently purified, and its equilibrium and kinetic folding properties were studied by absorbance, fluorescence, and near- and far-UV circular dichroism spectroscopies. The native state of the fragment unfolds cooperatively in an apparent two-state transition with a stability of 3.98 +/- 0.19 kcal mol(-1) in the absence of denaturant and a corresponding m value of 1.07 +/- 0.05 kcal mol(-1) M(-1). Similar to the full-length protein, the unfolding of the fragment shows two kinetic phases which arise from the presence of two discrete native state populations. Additionally, the fragment exhibits a significant burst phase in unfolding, indicating that a fraction of the folded state ensemble under native conditions has properties similar to those of the equilibrium intermediate populated at 3 M urea in full-length alphaTS. Refolding of alphaTS(1-188) is also complex, exhibiting two detectable kinetic phases and a burst phase that is complete within 5 ms. The two slowest isomerization phases observed in the refolding of the full-length protein are absent in the fragment, suggesting that these phases reflect contributions from the carboxy-terminal segment. The folding mechanism of alphaTS(1-188) appears to be a simplified version of the mechanism for the full-length protein [Bilsel, O., Zitzewitz, J. A., Bowers, K.E, and Matthews, C. R.(1999) Biochemistry 38, 1018-1029]. Four parallel channels in the full-length protein are reduced to a pair of channels that most likely reflect a cis/trans proline isomerization reaction in the amino-terminal fragment. The off- and on-pathway intermediates that exist for both full-length alphaTS and alphaTS(1-188) may reflect the preponderance of local interactions in the beta/alpha barrel motif.     

Proline replacements and the simplification of the complex,parallel channel folding mechanism for the alpha subunit of Trp synthase,a TIM barrel protein

The kinetic folding mechanism for the alpha subunit of tryptophan synthase (alphaTS) from Escherichia coli involves four parallel channels whose inter-conversions are controlled by three cis/trans prolyl isomerization reactions (tau(1), tau(2) and tau(3)). A previous mutational analysis of all 19 proline positions, including the unique cis Asp27-Pro28 peptide bond, revealed that the G(3)P28G, P78A or P96A mutations selectively eliminated the fast, tau(1) (ten seconds), folding phase, while the P217M and P261A mutations eliminated the medium, tau(2) (40 seconds) and the slow, tau(3) ( approximately 300 seconds) folding phases, respectively. To further elucidate the role of these proline residues and to simplify the folding mechanism, a series of double and triple mutants were constructed at these critical positions, and comprehensive kinetic and thermodynamic experiments were performed. Although it was not possible to construct a stable system that was free of proline isomerization constraints, a double mutant variant, G(3)P28G/P217M, in which the refolding of more than 90% of the unfolded protein is not limited by proline isomerization reactions was identified. Further, long-range interactions between several of these residues appear to be a crucial part of the cooperative network of structure that stabilizes the TIM barrel motif for alphaTS.     

Parallel channels and rate-limiting steps in complex protein folding reactions: prolyl isomerization and the alpha subunit of Trp synthase,a TIM barrel protein

A kinetic folding mechanism for the alpha subunit of tryptophan synthase (alphaTS) from Escherichia coli, involving four parallel channels with multiple native, intermediate and unfolded forms, has recently been proposed. The hypothesis that cis/trans isomerization of several Xaa-Pro peptide bonds is the source of the multiple folding channels was tested by measuring the sensitivity of the three rate-limiting phases (tau(1), tau(2), tau(3)) to catalysis by cyclophilin, a peptidyl-prolyl isomerase. Although the absence of catalysis for the tau(1) (fast) phase leaves its assignment ambiguous, our previous mutational analysis demonstrated its connection to the unique cis peptide bond preceding proline 28. The acceleration of the tau(2) (medium) and tau(3) (slow) refolding phases by cyclophilin demonstrated that cis/trans prolyl isomerization is also the source of these phases. A collection of proline mutants, which covered all of the remaining 18 trans proline residues of alphaTS, was constructed to obtain specific assignments for these phases. Almost all of the mutant proteins retained the complex equilibrium and kinetic folding properties of wild-type alphaTS; only the P217A, P217G and P261A mutations caused significant changes in the equilibrium free energy surface. Both the P78A and P96A mutations selectively eliminated the tau(1) folding phase, while the P217M and P261A mutations eliminated the tau(2) and tau(3) folding phases, respectively. The redundant assignment of the tau(1) phase to Pro28, Pro78 and Pro96 may reflect their mutual interactions in non-random structure in the unfolded state. The non-native cis isomers for Pro217 and Pro261 may destabilize an autonomous C-terminal folding unit, thereby giving rise to kinetically distinct unfolded forms. The nature of the preceding amino acid, the solvent exposure, or the participation in specific elements of secondary structure in the native state, in general, are not determinative of the proline residues whose isomerization reactions can limit folding.     

Partial NMR assignments and secondary structure mapping of the isolated alpha subunit of Escherichia coli tryptophan synthase,a 29-kD TIM barrel protein

The alpha subunit of tryptophan synthase (alphaTS) from S. typhimurium belongs to the triosephosphate isomerase (TIM) or the (beta/alpha)(8) barrel fold, one of the most common structures in biology. To test the conservation of the global fold in the isolated Escherichia coli homolog, we have obtained a majority of the backbone assignments for the 29-kD alphaTS by using standard heteronuclear multidimensional NMR methods on uniformly (15)N- and (15)N/(13)C-labeled protein and on protein selectively (15)N-labeled at key hydrophobic residues. The secondary structure mapped by chemical shift index, nuclear Overhauser enhancements (NOEs), and hydrogen-deuterium (H-D) exchange, and several abnormal chemical shifts are consistent with the conservation of the global TIM barrel fold of the isolated E. coli alphaTS. Because most of the amide protons that are slow to exchange with solvent correspond to the beta-sheet residues, the beta-barrel is likely to play an important role in stabilizing the previously detected folding intermediates for E. coli alphaTS. A similar combination of uniform and selective labeling can be extended to other TIM barrel proteins to obtain insight into the role of the motif in stabilizing what appear to be common partially folded forms.     

Anabaena apoflavodoxin hydrogen exchange: on the stable exchange core of the alpha/beta(21345) flavodoxin-like family

An important issue in modern protein biophysics is whether structurally homologous proteins share common stability and/or folding features. Flavodoxin is an archetypal alpha/beta protein organized in three layers: a central beta-sheet (strand order 21345) flanked by helices 1 and 5 on one side and helices 2, 3, and 4 on the opposite side. The backbone internal dynamics of the apoflavodoxin from Anabaena is analyzed here by the hydrogen exchange method. The hydrogen exchange rates indicate that 46 amide protons, distributed throughout the structure of apoflavodoxin, exchange relatively slowly at pH 7.0 (k(ex) < 10(-1) min(-1)). According to their distribution in the structure, protein stability is highest on the beta-sheet, helix 4, and on the layer formed by helices 1 and 5. The exchange kinetics of Anabaena apoflavodoxin was compared with those of the apoflavodoxin from Azotobacter, with which it shares a 48% sequence identity, and with Che Y and cutinase, two other alpha/beta (21345) proteins with no significant sequence homology with flavodoxins. Both similarities and differences are observed in the cores of these proteins. It is of interest that a cluster of a few structurally equivalent residues in the central beta-strands and in helix 5 is common to the cores.     

Topology and sequence in the folding of a TIM barrel protein: global analysis highlights partitioning between transient off-pathway and stable on-pathway folding intermediates in the complex folding mechanism of a (betaalpha)8 barrel of unknown function from B. subtilis

The relative contributions of chain topology and amino acid sequence in directing the folding of a (betaalpha)(8) TIM barrel protein of unknown function encoded by the Bacillus subtilis iolI gene (IOLI) were assessed by reversible urea denaturation and a combination of circular dichroism, fluorescence and time-resolved fluorescence anisotropy spectroscopy. The equilibrium reaction for IOLI involves, in addition to the native and unfolded species, a stable intermediate with significant secondary structure and stability and self-associated forms of both the native and intermediate states. Global kinetic analysis revealed that the unfolded state partitions between an off-pathway refolding intermediate and the on-pathway equilibrium intermediate early in folding. Comparisons with the folding mechanisms of two other TIM barrel proteins, indole-3-glycerol phosphate synthase from the thermophile Sulfolobus solfataricus (sIGPS) and the alpha subunit of Escherichia coli tryptophan synthase (alphaTS), reveal striking similarities that argue for a dominant role of the topology in both early and late events in folding. Sequence-specific effects are apparent in the magnitudes of the relaxation times and relative stabilities, in the presence of additional monomeric folding intermediates for alphaTS and sIGPS and in rate-limiting proline isomerization reactions for alphaTS.     

On the role of some conserved and nonconserved amino acid residues in the transitional state and intermediate of apomyoglobin folding

The contributions of some amino acid residues in the A, B, G, and H helices to the formation of the folding nucleus and folding intermediate of apomyoglobin were estimated. The effects of point substitutions of Ala for hydrophobic amino acid residues on the structural stability of the native (N) protein and its folding intermediate (I), as well as on the folding/unfolding rates for four mutant apomyoglobin forms, were studied. The equilibrium and kinetic studies of the folding/unfolding rates of these mutant proteins in a wide range of urea concentrations demonstrated that their native state was considerably destabilized as compared with the wild-type protein, whereas the stability of the intermediate state changed moderately. It was shown that the amino acid residues in the A, G, and H helices contributed insignificantly to the stabilization of the apomyoglobin folding nucleus in the rate-limiting I ? N transition, taking place after the formation of the intermediate, whereas the residue of the B helix was of great importance in the formation of the folding nucleus in this transition.     

Mapping the structure of folding cores in TIM barrel proteins by hydrogen exchange mass spectrometry: the roles of motif and sequence for the indole-3-glycerol phosphate synthase from Sulfolobus solfataricus

To test the roles of motif and amino acid sequence in the folding mechanisms of TIM barrel proteins, hydrogen-deuterium exchange was used to explore the structure of the stable folding intermediates for the of indole-3-glycerol phosphate synthase from Sulfolobus solfataricus (sIGPS). Previous studies of the urea denaturation of sIGPS revealed the presence of an intermediate that is highly populated at approximately 4.5 M urea and contains approximately 50% of the secondary structure of the native (N) state. Kinetic studies showed that this apparent equilibrium intermediate is actually comprised of two thermodynamically distinct species, I(a) and I(b). To probe the location of the secondary structure in this pair of stable on-pathway intermediates, the equilibrium unfolding process of sIGPS was monitored by hydrogen-deuterium exchange mass spectrometry. The intact protein and pepsin-digested fragments were studied at various concentrations of urea by electrospray and matrix-assisted laser desorption ionization time-of-flight mass spectrometry, respectively. Intact sIGPS strongly protects at least 54 amide protons from hydrogen-deuterium exchange in the intermediate states, demonstrating the presence of stable folded cores. When the protection patterns and the exchange mechanisms for the peptides are considered with the proposed folding mechanism, the results can be interpreted to define the structural boundaries of I(a) and I(b). Comparison of these results with previous hydrogen-deuterium exchange studies on another TIM barrel protein of low sequence identify, alpha-tryptophan synthase (alphaTS), indicates that the thermodynamic states corresponding to the folding intermediates are better conserved than their structures. Although the TIM barrel motif appears to define the basic features of the folding free energy surface, the structures of the partially folded states that appear during the folding reaction depend on the amino acid sequence. Markedly, the good correlation between the hydrogen-deuterium exchange patterns of sIGPS and alphaTS with the locations of hydrophobic clusters defined by isoleucine, leucine, and valine residues suggests that branch aliphatic side-chains play a critical role in defining the structures of the equilibrium intermediates.     

Folding propensities of synthetic peptide fragments covering the entire sequence of phage 434 Cro protein.

The phage 434 Cro protein, the N-terminal domain of its repressor (R1-69) and that of phage lambda (lambda6-85) constitute a group of small, monomeric, single-domain folding units consisting of five helices with striking structural similarity. The intrinsic helix stabilities in lambda6-85 have been correlated to its rapid folding behavior, and a residual hydrophobic cluster found in R1-69 in 7 M urea has been proposed as a folding initiation site. To understand the early events in the folding of 434 Cro, and for comparison with R1-69 and lambda6-85, we examined the conformational behavior of five peptides covering the entire 434 Cro sequence in water, 40% (by volume) TFE/water, and 7 M urea solutions using CD and NMR. Each peptide corresponds to a helix and adjacent residues as identified in the native 434 Cro NMR and crystal structures. All are soluble and monomeric in the solution conditions examined except for the peptide corresponding to the 434 Cro helix 4, which has low water solubility. Helix formation is observed for the 434 Cro helix 1 and helix 2 peptides in water, for all the peptides in 40% TFE and for none in 7 M urea. NMR data indicate that the helix limits in the peptides are similar to those in the native protein helices. The number of side-chain NOEs in water and TFE correlates with the helix content, and essentially none are observed in 7 M urea for any peptide, except that for helix 5, where a hydrophobic cluster may be present. The low intrinsic folding propensities of the five helices could account for the observed stability and folding behavior of 434 Cro and is, at least qualitatively, in accord with the results of the recently described diffusion-collision model incorporating intrinsic helix propensities.     

Benchmarking of TASSER in the ab initio limit

A significant number of protein sequences in a given proteome have no obvious evolutionarily related protein in the database of solved protein structures, the PDB. Under these conditions, ab initio or template-free modeling methods are the sole means of predicting protein structure. To assess its expected performance on proteomes, the TASSER structure prediction algorithm is benchmarked in the ab initio limit on a representative set of 1129 nonhomologous sequences ranging from 40 to 200 residues that cover the PDB at 30% sequence identity and which adopt alpha, alpha + beta, and beta secondary structures. For sequences in the 40-100 (100-200) residue range, as assessed by their root mean square deviation from native, RMSD, the best of the top five ranked models of TASSER has a global fold that is significantly close to the native structure for 25% (16%) of the sequences, and with a correct identification of the structure of the protein core for 59% (36%). In the absence of a native structure, the structural similarity among the top five ranked models is a moderately reliable predictor of folding accuracy. If we classify the sequences according to their secondary structure content, then 64% (36%) of alpha, 43% (24%) of alpha + beta, and 20% (12%) of beta sequences in the 40-100 (100-200) residue range have a significant TM-score (TM-score > or = 0.4). TASSER performs best on helical proteins because there are less secondary structural elements to arrange in a helical protein than in a beta protein of equal length, since the average length of a helix is longer than that of a strand. In addition, helical proteins have shorter loops and dangling tails. If we exclude these flexible fragments, then TASSER has similar accuracy for sequences containing the same number of secondary structural elements, irrespective of whether they are helices and/or strands. Thus, it is the effective configurational entropy of the protein that dictates the average likelihood of correctly arranging the secondary structure elements.     

Intestinal fatty acid binding protein: the folding mechanism as determined by NMR studies

The intestinal fatty acid binding protein is composed of two beta-sheets surrounding a large interior cavity. There is a small helical domain associated with the portal for entry of the ligand into the cavity. Denaturation of the protein has been monitored in a residue-specific manner by collecting a series of two-dimensional (1)H-(15)N heteronuclear single-quantum coherence (HSQC) NMR spectra from 0 to 6.5 M urea under equilibrium conditions. In addition, rates for hydrogen-deuterium exchange have been measured as a function of denaturant concentration. Residual, native-like structure persists around hydrophobic clusters at very high urea concentrations. This residual structure (reflecting only about 2-7% persistence of native-like structure) involves the turns between beta-strands and between the two short helices. If this persistence is assumed to reflect transient native-like structure in these regions of the polypeptide chain, these sites may serve as nucleation sites for folding. The data obtained at different urea concentrations are then analyzed on the basis of peak intensities relative to the intensities in the absence of urea reflecting the extent of secondary structure formation. At urea concentrations somewhat below 6.5 M, specific hydrophobic residues in the C-terminal beta-sheet interact and two strands, the D and E strands in the N-terminal beta-sheet, are stabilized. These latter strands surround one of the turns showing residual structure. With decreasing urea concentrations, the remaining strands are stabilized in a specific order. The early strand stabilization appears to trigger the formation of the remainder of the C-terminal beta-sheet. At low urea concentrations, hydrogen bonds are formed. A pathway is proposed on the basis of the data describing the early, intermediate, and late folding steps for this almost all beta-sheet protein. The data also show that there are regions of the protein which appear to act in a concerted manner at intermediate steps in refolding.     

Identification and characterization of key substructures involved in the early folding events of a (beta/alpha)8-barrel protein as studied by experimental and computational methods

A number of studies have examined the structural properties of late folding intermediates of (beta/alpha)8-barrel proteins involved in tryptophan biosynthesis, whereas there is little information available about the early folding events of these proteins. To identify the contiguous polypeptide segments important to the folding of the (beta/alpha)8-barrel protein Escherichia coli N-(5'-phosphoribosyl)anthranilate isomerase, we structurally characterized fragments and circularly permuted forms of the protein. We also simulated thermal unfolding of the protein using molecular dynamics. Our fragmentation experiments demonstrate that the isolated (beta/alpha)(1-4)beta5 fragment is almost as stable as the full-length protein. The far and near-UV CD spectra of this fragment are indicative of native-like secondary and tertiary structures. Structural analysis of the circularly permutated proteins shows that if the protein is cleaved within the two N-terminal betaalpha modules, the amount of secondary structure is unaffected, whereas, when cleaved within the central (beta/alpha)(3-4)beta5 segment, the protein simply cannot fold. An ensemble of the denatured structures produced by thermal unfolding simulations contains a persistent local structure comprised of beta3, beta4 and beta5. The presence of this three-stranded beta-barrel suggests that it may be an important early-stage folding intermediate. Interactions found in (beta/alpha)(3-4)beta5 may be essential for the early events of ePRAI folding if they provide a nucleation site that directs folding.     

Detection and structure determination of an equilibrium unfolding intermediate of Rd-apocytochrome b562: native fold with non-native hydrophobic interactions

The absence of detectable kinetic and equilibrium folding intermediates by optical probes is commonly taken to indicate that protein folding is a two-state process. However, for some small proteins with apparent two-state behavior, unfolding intermediates have been identified in native-state hydrogen exchange or kinetic unfolding experiments monitored by nuclear magnetic resonance. Rd-apocytochrome b(562), a four-helix bundle, is one such protein. Here, we found another unfolding intermediate for Rd-apocytochrome b(562). It is based on a cooperative transition of (15)N chemical shifts of amide protons as a function of urea concentrations before the global unfolding. We have solved the high-resolution structure of the protein at 2.8 M urea, which is after this cooperative transition but before the global unfolding. All four helices remained intact, but a number of hydrophobic core residues repacked. This intermediate provides a possible structural interpretation for the kinetic unfolding intermediates observed using nuclear magnetic resonance methods for several proteins and has important implications for theoretical studies of protein folding.