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.     

Early formation of a beta hairpin during folding of staphylococcal nuclease H124L as detected by pulsed hydrogen exchange

Pulsed hydrogen exchange methods were used to follow the formation of structure during the refolding of acid-denatured staphylococcal nuclease containing a stabilizing Leu substitution at position 124 (H124L SNase). The protection of more than 60 backbone amide protons in uniformly (15)N-labeled H124L SNase was monitored as a function of refolding time by heteronuclear two-dimensional NMR spectroscopy. As found in previous studies of staphylococcal nuclease, partial protection was observed for a subset of amide protons even at the earliest folding time point (10 msec). Protection indicative of marginally stable hydrogen-bonded structure in an early folding intermediate was observed at over 30 amide positions located primarily in the beta-barrel and to a lesser degree in the alpha-helical domain of H124L SNase. To further characterize the folding intermediate, protection factors for individual amide sites were measured by varying the pH of the labeling pulse at a fixed refolding time of 16 msec. Protection factors >5.0 were observed only for amide positions in a beta-hairpin formed by strands 2 and 3 of the beta-barrel domain and a single site near the C-terminus. The results indicate that formation of stable hydrogen-bonded structure in a core region of the beta-sheet is among the earliest structural events in the folding of SNase and may serve as a nucleation site for further structure formation.     

Hydrogen exchange in BPTI variants that do not share a common disulfide bond.

Bovine pancreatic trypsin inhibitor (BPTI) is stabilized by 3 disulfide bonds, between cysteines 30-51, 5-55, and 14-38. To better understand the influence of disulfide bonds on local protein structure and dynamics, we have measured amide proton exchange rates in 2 folded variants of BPTI, [5-55]Ala and [30-51; 14-38]V5A55, which share no common disulfide bonds. These proteins resemble disulfide-bonded intermediates that accumulate in the BPTI folding pathway. Essentially the same amide hydrogens are protected from exchange in both of the BPTI variants studied here as in native BPTI, demonstrating that the variants adopt fully folded, native-like structures in solution. However, the most highly protected amide protons in each variant differ, and are contained within the sequences of previously studied peptide models of related BPTI folding intermediates containing either the 5-55 or the 30-51 disulfide bond.     

Folding subdomains of thioredoxin characterized by native-state hydrogen exchange

Native-state hydrogen exchange (HX) studies, used in conjunction with NMR spectroscopy, have been carried out on Escherichia coli thioredoxin (Trx) for characterizing two folding subdomains of the protein. The backbone amide protons of only the slowest-exchanging 24 amino acid residues, of a total of 108 amino acid residues, could be followed at pH 7. The free energy of the opening event that results in an amide hydrogen exchanging with solvent (DeltaG(op)) was determined at each of the 24 amide hydrogen sites. The values of DeltaG(op) for the amide hydrogens belonging to residues in the helices alpha(1), alpha(2), and alpha(4) are consistent with them exchanging with the solvent only when the fully unfolded state is sampled transiently under native conditions. The denaturant-dependences of the values of DeltaG(op) provide very little evidence that the protein samples partially unfolded forms, lower in energy than the unfolded state. The amide hydrogens belonging to the residues in the beta strands, which form the core of the protein, appear to have higher values of DeltaG(op) than amide hydrogens belonging to residues in the helices, suggesting that they might be more stable to exchange. This apparently higher stability to HX of the beta strands might be either because they exchange out their amide hydrogens in a high energy intermediate preceding the globally unfolded state, or, more likely, because they form residual structure in the globally unfolded state. In either case, the central beta strands-beta(3,) beta(2), and beta(4)-would appear to form a cooperatively folding subunit of the protein. The native-state HX methodology has made it possible to characterize the free energy landscape that Trx can sample under equilibrium native conditions.     

Structural and kinetic characterization of early folding events in beta-lactoglobulin

We have defined the structural and dynamic properties of an early folding intermediate of beta-lactoglobulin known to contain non-native alpha-helical structure. The folding of beta-lactoglobulin was monitored over the 100 micros--10 s time range using ultrarapid mixing techniques in conjunction with fluorescence detection and hydrogen exchange labeling probed by heteronuclear NMR. An initial increase in Trp fluorescence with a time constant of 140 micros is attributed to formation of a partially helical compact state. Within 2 ms of refolding, well protected amide protons indicative of stable hydrogen bonded structure were found only in a domain comprising beta-strands F, G and H, and the main alpha-helix, which was thus identified as the folding core of beta-lactoglobulin. At the same time, weak protection (up to approximately 10-fold) of amide protons in a segment spanning residues 12--21 is consistent with formation of marginally stable non-native alpha-helices near the N-terminus. Our results indicate that efficient folding, despite some local non-native structural preferences, is insured by the rapid formation of a native-like alpha/beta core domain.     

Microscopic stability of cold shock protein A examined by NMR native state hydrogen exchange as a function of urea and trimethylamine N-oxide

Native state hydrogen exchange of cold shock protein A (CspA) has been characterized as a function of the denaturant urea and of the stabilizing agent trimethylamine N-oxide (TMAO). The structure of CspA has five strands of beta-sheet. Strands beta1-beta4 have strongly protected amide protons that, based on experiments as a function of urea, exchange through a simple all-or-none global unfolding mechanism. By contrast, the protection of amide protons from strand beta5 is too weak to measure in water. Strand beta5 is hydrogen bonded to strands beta3 and beta4, both of which afford strong protection from solvent exchange. Gaussian network model (GNM) simulations, which assume that the degree of protection depends on tertiary contact density in the native structure, accurately predict the strong protection observed in strands beta1-beta4 but fail to account for the weak protection in strand beta5. The most conspicuous feature of strand beta5 is its low sequence hydrophobicity. In the presence of TMAO, there is an increase in the protection of strands beta1-beta4, and protection extends to amide protons in more hydrophilic segments of the protein, including strand beta5 and the loops connecting the beta-strands. TMAO stabilizes proteins by raising the free energy of the denatured state, due to highly unfavorable interactions between TMAO and the exposed peptide backbone. As such, the stabilizing effects of TMAO are expected to be relatively independent of sequence hydrophobicity. The present results suggest that the magnitude of solvent exchange protection depends more on solvent accessibility in the ensemble of exchange susceptible conformations than on the strength of hydrogen-bonding interactions in the native structure.     

The apomyoglobin folding pathway revisited: structural heterogeneity in the kinetic burst phase intermediate

Extensive analysis of accurate quench-flow hydrogen exchange results indicates that the burst phase kinetic intermediate in the folding of apomyoglobin (apoMb) from urea is structurally heterogeneous. The structural variability is associated with the partial folding of the E helix during the burst phase (<6.4ms) of the folding process. Analysis of the effects of exchange-out of amide proton labels during the labeling pulse ( approximately pH 10) of the quench-flow process indicates that three of the amide protons in the E helix are in fact largely protected in the burst phase of folding, while the remainder of the E helix has a substantial complement of amide protons that show biphasic kinetics, i.e. are protected partly during the burst phase and partly during the slow phase of folding. The locations of these amide protons can be used to map the sites of structural heterogeneity in the kinetic molten globule. These sites include, besides the E helix, the ends of the A and B helices and part of the C helix. Our results give significant support to the hypothesis that the kinetic molten globule intermediate of apoMb is native-like.     

Conformational stability of ribonuclease T1 determined by hydrogen-deuterium exchange.

The hydrogen-deuterium exchange kinetics of 37 backbone amide residues in RNase T1 have been monitored at 25, 40, 45, and 50 degrees C at pD 5.6 and at 40 and 45 degrees C at pD 6.6. The hydrogen exchange rate constants of the hydrogen-bonded residues varied over eight orders of magnitude at 25 degrees C with 13 residues showing exchange rates consistent with exchange occurring as a result of global unfolding. These residues are located in strands 2-4 of the central beta-pleated sheet. The residues located in the alpha-helix and the remaining strands of the beta-sheet exhibited exchange behaviors consistent with exchange occurring due to local structural fluctuations. For several residues at 25 degrees C, the global free energy change calculated from the hydrogen exchange data was over 2 kcal/mol greater than the free energy of unfolding determined from urea denaturation experiments. The number of residues showing this unexpected behavior was found to increase with temperature. This apparent inconsistency can be explained quantitatively if the cis-trans isomerization of the two cis prolines, Pro-39 and Pro-55, is taken into account. The cis-trans isomerization equilibrium calculated from kinetic data indicates the free energy of the unfolded state will be 2.6 kcal/mol higher at 25 degrees C when the two prolines are cis rather than trans (Mayr LM, Odefey CO, Schutkowski M, Schmid FX. 1996. Kinetic analysis of the unfolding and refolding of ribonuclease T1 by a stopped-flow double-mixing technique. Biochemistry 35: 5550-5561). The hydrogen exchange results are consistent with the most slowly exchanging hydrogens exchanging from a globally higher free energy unfolded state in which Pro-55 and Pro-39 are still predominantly in the cis conformation. When the conformational stabilities determined by hydrogen exchange are corrected for the proline isomerization equilibrium, the results are in excellent agreement with those from an analysis of urea denaturation curves.     

Distributions of beta sheets in proteins with application to structure prediction

We recently developed the Rosetta algorithm for ab initio protein structure prediction, which generates protein structures from fragment libraries using simulated annealing. The scoring function in this algorithm favors the assembly of strands into sheets. However, it does not discriminate between different sheet motifs. After generating many structures using Rosetta, we found that the folding algorithm predominantly generates very local structures. We surveyed the distribution of beta-sheet motifs with two edge strands (open sheets) in a large set of non-homologous proteins. We investigated how much of that distribution can be accounted for by rules previously published in the literature, and developed a filter and a scoring method that enables us to improve protein structure prediction for beta-sheet proteins. Proteins 2002;48:85-97.     

Three-dimensional structure of Escherichia coli asparagine synthetase B: a short journey from substrate to product

Asparagine synthetase B catalyzes the assembly of asparagine from aspartate, Mg(2+)ATP, and glutamine. Here, we describe the three-dimensional structure of the enzyme from Escherichia colidetermined and refined to 2.0 A resolution. Protein employed for this study was that of a site-directed mutant protein, Cys1Ala. Large crystals were grown in the presence of both glutamine and AMP. Each subunit of the dimeric protein folds into two distinct domains. The N-terminal region contains two layers of antiparallel beta-sheet with each layer containing six strands. Wedged between these layers of sheet is the active site responsible for the hydrolysis of glutamine. Key side chains employed for positioning the glutamine substrate within the binding pocket include Arg 49, Asn 74, Glu 76, and Asp 98. The C-terminal domain, responsible for the binding of both Mg(2+)ATP and aspartate, is dominated by a five-stranded parallel beta-sheet flanked on either side by alpha-helices. The AMP moiety is anchored to the protein via hydrogen bonds with O(gamma) of Ser 346 and the backbone carbonyl and amide groups of Val 272, Leu 232, and Gly 347. As observed for other amidotransferases, the two active sites are connected by a tunnel lined primarily with backbone atoms and hydrophobic and nonpolar amino acid residues. Strikingly, the three-dimensional architecture of the N-terminal domain of asparagine synthetase B is similar to that observed for glutamine phosphoribosylpyrophosphate amidotransferase while the molecular motif of the C-domain is reminiscent to that observed for GMP synthetase.     

Apparent local stability of the secondary structure of Azotobacter vinelandii holoflavodoxin II as probed by hydrogen exchange: implications for redox potential regulation and flavodoxin folding.

As a first step to determine the folding pathway of a protein with an alpha/beta doubly wound topology, the 1H, 13C, and 15N backbone chemical shifts of Azotobacter vinelandii holoflavodoxin II (179 residues) have been determined using multidimensional NMR spectroscopy. Its secondary structure is shown to contain a five-stranded parallel beta-sheet (beta2-beta1-beta3-beta4-beta5) and five alpha-helices. Exchange rates for the individual amide protons of holoflavodoxin were determined using the hydrogen exchange method. The amide protons of 65 residues distributed throughout the structure of holoflavodoxin exchange slowly at pH* 6.2 [kex < 10(-5) s(-1)] and can be used as probes in future folding studies. Measured exchange rates relate to apparent local free energies for transient opening. We propose that the amide protons in the core of holoflavodoxin only exchange by global unfolding of the apo state of the protein. The results obtained are discussed with respect to their implications for flavodoxin folding and for modulation of the flavin redox potential by the apoprotein. We do not find any evidence that A. vinelandii holoflavodoxin II is divided into two subdomains based on its amide proton exchange rates, as opposed to what is found for the structurally but not sequentially homologous alpha/beta doubly wound protein Che Y.     

Temperature-dependent folding pathways of Pin1 WW domain: an all-atom molecular dynamics simulation of a Gō model

We study the folding thermodynamics and kinetics of the Pin1 WW domain, a three-stranded beta-sheet protein, by using all-atom (except nonpolar hydrogens) discontinuous molecular dynamics simulations at various temperatures with a Gō model. The protein exhibits a two-state folding kinetics near the folding transition temperature. A good agreement between our simulations and the experimental measurements by the Gruebele group has been found, and the simulation sheds new insights into the structure of transition state, which is hard to be straightforwardly captured in experiments. The simulation also reveals that the folding pathways at approximately the transition temperature and at low temperatures are much different, and an intermediate state at a low temperature is predicted. The transition state of this small beta-protein at its folding transition temperature has a well-established hairpin 1 made of beta1 and beta2 strands while its low-temperature kinetic intermediate has a formed hairpin 2 composed of beta2 and beta3 strands. Theoretical results are compared with other simulation results as well as available experimental data. This study confirms that specific side-chain packing in an all-atom Gō model can yield a reasonable prediction of specific folding kinetics for a given protein. Different folding behaviors at different temperatures are interpreted in terms of the interplay of entropy and enthalpy in folding process.     

Situations of gamma-turns in proteins. Their relation to alpha-helices, beta-sheets and ligand binding sites

Gamma-turns occur as one of two possible enantiomers with regard to their main-chain structure, called classic and inverse. Of these, inverse ones are more common. Unlike other hydrogen bonds, those in inverse gamma-turns include a large proportion that are weak. If such hydrogen bonds are included, these turns may be said to be abundant in proteins. A significant number of inverse gamma-turns, usually weak ones, exist as consecutive turns in a structural feature, now called the 2.2(7)-helix, proposed for polypeptides as long ago as 1943 by Huggins. Most of these features occur within strands of beta-sheet. The less-weak inverse gamma-turns fall into several structural subgroups. They are frequently situated directly at either end of alpha-helices or of strands of beta-sheet, or adjacent to certain loop motifs. In general, they are well conserved during evolution and some are found at key positions in proteins. One occurs in the first hypervariable loop in the heavy chain of immunoglobulins.     

Detection and characterization of an early folding intermediate of T4 lysozyme using pulsed hydrogen exchange and two-dimensional NMR.

Two-dimensional 1H-15N NMR techniques combined with pulsed hydrogen-deuterium exchange have been used to characterize the folding pathway of T4 lysozyme. In the unfolded state, there is little differential protection of the various amides from hydrogen exchange. In the native folded structure, 84 amides of the 164 residues are sufficiently spectrally resolved and protected from solvent exchange to serve as probes of the folding pathway. These probes are located in both the N-terminal and C-terminal domains of the native folded structure of the protein. The studies described here show that at least one intermediate is formed early during refolding at low denaturant concentrations. This intermediate (or intermediates) forms very rapidly (within the 10-ms temporal resolution of our mixing device) under the conditions used and is completed at least 10 times faster than the overall folding event. The intermediate(s) protect(s) from exchange a subset of amides in the N-terminal and C-terminal regions of the protein. In the final folded states these protected regions correspond to two alpha-helices and a beta-sheet region. These amides are protected from exchange by factors between 20 and 200 as compared to the fully unfolded protein. Protection of this magnitude is consistent with the formation of somewhat exposed secondary structure in these regions and could represent a "molten globule"-like or a "framework"-like structure for the intermediate(s) in which specific parts of the sequence form isolated secondary structures that are not stabilized by extensive tertiary interactions.     

Hydrogen bonding in globular proteins.

A global census of the hydrogen bonds in 42 X-ray-elucidated proteins was taken and the following demographic trends identified: (1) Most hydrogen bonds are local, i.e. between partners that are close in sequence, the primary exception being hydrogen-bonded ion pairs. (2) Most hydrogen bonds are between backbone atoms in the protein, an average of 68%. (3) All proteins studied have extensive hydrogen-bonded secondary structure, an average of 82%. (4) Almost all backbone hydrogen bonds are within single elements of secondary structure. An approximate rule of thirds applies: slightly more than one-third (37%) form i----i--3 hydrogen bonds, almost one-third (32%) form i----i--4 hydrogen bonds, and slightly less than one-third (26%) reside in paired strands of beta-sheet. The remaining 5% are not wholly within an individual helix, turn or sheet. (5) Side-chain to backbone hydrogen bonds are clustered at helix-capping positions. (6) An extensive network of hydrogen bonds is present in helices. (7) To a close approximation, the total number of hydrogen bonds is a simple function of a protein's helix and sheet content. (8) A unique quantity, termed the reduced number of hydrogen bonds, is defined as the maximum number of hydrogen bonds possible when every donor:acceptor pair is constrained to be 1:1. This quantity scales linearly with chain length, with 0.71 reduced hydrogen bond per residue. Implications of these results for pathways of protein folding are discussed.     

Interplay between hydrophobic cluster and loop propensity in beta-hairpin formation

Autonomously folding beta-hairpins have recently emerged as powerful tools for elucidating the origins of antiparallel beta-sheet folding preferences. Analysis of such model systems has suggested four potential sources of beta-sheet stability: (1) the conformational propensity of the loop segment that connects adjacent strands; (2) favorable contacts between side-chains on adjacent strands; (3) interstrand hydrogen bonds; and (4) the intrinsic beta-sheet propensities of the strand residues. We describe the design and analysis of a series of isomeric 20 residue peptides in which factors (1)-(4) are identical. Differences in beta-hairpin formation within this series demonstrate that these four factors, individually, are not sufficient to explain beta-sheet stability. In agreement with the prediction of a simple statistical mechanical model for beta-hairpin formation, our results show that the separation between the loop segment and an interstrand cluster of hydrophobic side-chains strongly influences beta-hairpin size and stability, with a smaller separation leading to greater stability.     

Protein folding kinetics by combined use of rapid mixing techniques and NMR observation of individual amide protons

A method to be used for experimental studies of protein folding introduced by Schmid and Baldwin (J. Mol. Biol. 135: 199-215, 1979), which is based on the competition between amide hydrogen exchange and protein refolding, was extended by using rapid mixing techniques and 1H NMR to provide site-resolved kinetic information on the early phases of protein structure acquisition. In this method, a protonated solution of the unfolded protein is rapidly mixed with a deuterated buffer solution at conditions assuring protein refolding in the mixture. This simultaneously initiates the exchange of unprotected amide protons with solvent deuterium and the refolding of protein segments which can protect amide groups from further exchange. After variable reaction times the amide proton exchange is quenched while folding to the native form continues to completion. By using 1H NMR, the extent of exchange at individual amide sites is then measured in the refolded protein. Competition experiments at variable reaction times or variable pH indicate the time at which each amide group is protected in the refolding process. This technique was applied to the basic pancreatic trypsin inhibitor, for which sequence-specific assignments of the amide proton NMR lines had previously been obtained. For eight individual amide protons located in the beta-sheet and the C-terminal alpha-helix of this protein, apparent refolding rates in the range from 15 s-1 to 60 s-1 were observed. These rates are on the time scale of the fast folding phase observed with optical probes.     

FTIR reveals structural differences between native beta-sheet proteins and amyloid fibrils

The presence of beta-sheets in the core of amyloid fibrils raised questions as to whether or not beta-sheet-containing proteins, such as transthyretin, are predisposed to form such fibrils. However, we show here that the molecular structure of amyloid fibrils differs more generally from the beta-sheets in native proteins. This difference is evident from the amide I region of the infrared spectrum and relates to the distribution of the phi/psi dihedral angles within the Ramachandran plot, the average number of strands per sheet, and possibly, the beta-sheet twist. These data imply that amyloid fibril formation from native beta-sheet proteins can involve a substantial structural reorganization.     

Solution structure of the kringle 4 domain from human plasminogen by 1H nuclear magnetic resonance spectroscopy and distance geometry

Kringle 4 is an autonomous structural and folding domain within the proenzyme plasminogen. Homologous domains are found throughout the blood clotting and fibrinolytic proteins. In this paper, we present the almost complete assignment of the 1H nuclear magnetic resonance (n.m.r.) spectrum of the kringle 4 domain of human plasminogen. A detailed structural analysis has been completed. The sequential pattern of nuclear Overhauser enhancements indicated little regular secondary structure but rather a series of turns and loops connecting beta-strands. A small stretch of antiparallel beta-sheet was identified between the residues 61 to 63 and 71 to 73 and the close proximity of other strands was determined from two-dimensional nuclear Overhauser enhancement spectra. Slowly exchanging amide (NH) resonances were found to be associated with residues of the beta-sheet and neighbouring strands that support the hydrophobic core of the domain. A total of 526 interproton distance constraints and two hydrogen bonds were specified as input to the distance geometry program DISGEO. Tertiary structures were produced that were consistent with the n.m.r. data. The structures were compared with that of our earlier model based on n.m.r. studies and with that of prothrombin fragment 1 determined crystallographically.     

Hydrogen exchange study of canine milk lysozyme: stabilization mechanism of the molten globule

The native state (1)H, (15)N resonance assignment of 123 of the 128 nonproline residues of canine milk lysozyme has enabled measurements of the amide hydrogen exchange of over 70 amide hydrogens in the molten globule state. To elucidate the mechanism of protein folding, the molten globule state has been studied as a model of the folding intermediate state. Lysozyme and alpha-lactalbumin are homologous to each other, but their equilibrium unfolding mechanisms differ. Generally, the folding mechanism of lysozyme obeys a two-state model, whereas that of alpha-lactalbumin follows a three-state model. Exceptions to this rule are equine and canine milk lysozymes, which exhibit a partially unfolded state during the equilibrium unfolding; this state resembles the molten globule state of alpha-lactalbumin but with extreme stability. Study of the molten globules of alpha-lactalbumin and equine milk lysozyme showed that the stabilities of their alpha-helices are similar, despite the differences in the thermodynamic stability of their molten globule states. On the other hand, our hydrogen exchange study of the molten globule of canine milk lysozyme showed that the alpha-helices are more stabilized than in alpha-lactalbumin or equine milk lysozyme and that this enhanced stability is caused by the strengthened cooperative interaction between secondary structure elements. Thus, our results underscore the importance of the cooperative interaction in the stability of the molten globule state.