|
|
||||||||

)8 barrel motif of the alpha subunit of tryptophan synthaseDepartment of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, Massachusetts 01605, USA
(RECEIVED December 4, 2006; FINAL REVISION April 4, 2007; ACCEPTED April 4, 2007)
| Abstract |
|---|
|
|
|---|

)8, TIM barrel protein, the alpha subunit of tryptophan synthase (
TS), was probed by mutational analysis. The F19D46 and I97D124 hydrogen bonds link the N terminus of a
-strand with the C terminus of the succeeding antiparallel
-helix, and the A103D130 hydrogen bond links the N terminus of an
-helix with the C terminus of the succeeding antiparallel
-strand, forming clamps for the respective 
or 
hairpins. The individual replacement of these aspartic acid side chains with alanine leads to what appear to be closely related partially folded structures with significantly reduced far-UV CD ellipticity and thermodynamic stability. Comparisons with the effects of eliminating another main-chainside-chain hydrogen bond, G26S33, and two electrostatic side-chainside-chain hydrogen bonds, D38H92 and D112H146, all in the same N-terminal folding unit of
TS, demonstrated a unique role for the clamp interactions in stabilizing the native barrel conformation. Because neither the asparagine nor glutamic acid variant at position 46 can completely reproduce the spectroscopic, thermodynamic, or kinetic folding properties of aspartic acid, both size and charge are crucial to its unique role in the clamp hydrogen bond. Kinetic studies suggest that the three clamp hydrogen bonds act in concert to stabilize the transition state leading to the fully folded TIM barrel motif.
Keywords: circular dichroism; protein folding; site-directed mutagenesis; thermodynamic and kinetic mechanisms; 
and 
hairpins
| Introduction |
|---|
|
|
|---|
12 kcal mol1 (Horovitz et al. 1990; Serrano et al. 1992; Myers and Pace 1996; Ibarra-Molero et al. 2004). The small contribution to stability for main-chainside-chain interactions is underscored by the relatively rapid exchange of their amide hydrogens with solvent through local fluctuations within the native ensemble of conformers (Englander and Kallenbach 1983; Krishna et al. 2004). Thus, the results of a native-state hydrogen exchange (HX) experiment on the alpha subunit of tryptophan synthase (
TS) are quite striking.
TS is a 29-kDa TIM barrel protein in which the eight parallel
-strands alternate in sequence with amphipathic antiparallel helices that dock on the hydrophobic barrel. H-bonds between
1 and
8 and van der Waals interactions between nonpolar side chains in
1,
1,
8, and
8'/
8 close the barrel (Fig. 1A). Relevant to the present issue, the 1H-15N heteronuclear single quantum coherence (HSQC) NMR spectrum of uniformly 15N-labeled protein (Vadrevu et al. 2003) provided a vehicle for demonstrating that the main-chain amide hydrogens of two residues, F19 and I97, are very resistant to exchange with solvent deuterium (R. Vadrevu, Y. Wu, and C. R. Matthews, in prep.). Inspection of the crystal structure of
TS (Hyde et al. 1988; Nishio et al. 2005) shows that both are involved in long-range main-chainside-chain H-bonds, F19D46 and I97D124 (Fig. 1B). The resistance of the amide hydrogen of F19 to exchange for several days and that of I97 for several months in 2H2O at pH 7.8 and 25°C implies that their amide hydrogens become exposed to solvent in rare high energy states (Bai et al. 1995), not through local fluctuations in the native basin. By implication, the F19D46 and I97D124 H-bonds might play unusually important roles in the structure and stability of both native and partially folded states in
TS.
|
TS that is distinct from the relatively minor perturbations for the other variants. | Results |
|---|
|
|
|---|
TS (Fig. 1A) shows that the F19D46, I97D124, and A103D130 H-bonds each bracket a sequential pair of elements of secondary structure. Each of these three long-range H-bond interactions involves a main-chain amide hydrogen donor and a side-chain oxygen acceptor in aspartic acid. The F19D46 H-bond links the N terminus of
1 with the C terminus of
1, the I97D124 H-bond links the N terminus of
3 with the C terminus of
3, and the A103D130 H-bond links the N terminus of
3 with the C terminus of
4 (Fig. 1B). All three interactions span precisely 27 residues in sequence, and, by analogy, serve as structural clamps for 
hairpins, F19D46 and I97D124, or an 
hairpin, A103D130, in
TS.
The uniqueness of these structural clamp interactions can be assessed by comparisons with other noncovalent interactions involving short-range or surface H-bonds. The G26S33 H-bond braces a tight turn between
1 and
1 that contains the only cis prolyl peptide bond in
TS between D27 and P28. The D38H92 and D112H146 salt bridges on the surface of the barrel link the middle of
1 to the loop between
2 and
2 and the N terminus of
3 with the loop between
4 and
5, respectively. (In the higher resolution structure of
TS from Salmonella typhimurium [Hyde et al. 1988], D112 appears to form a second salt bridge with R145.)
Alanine replacements at D46, D124, or D130 significantly perturb the secondary and tertiary structure of
TS
The effects of replacing D46, D124, or D130 with alanine on the secondary and tertiary structure of the
TS variants were determined by collecting far- and near-UV CD spectra (Fig. 2). The far-UV CD spectra show a maximum mean residue ellipticity (MRE) at 195 nm and minima at 208 nm and 222 nm, demonstrating that all of the variants retain components of the
-helices and
-strands found in the wild-type (WT) protein. However, the significant decrease in the MRE by
40% for the variants (Fig. 2AC) at all three wavelengths indicate that each individual replacement causes a substantial loss of secondary structure.
|
TS, indicate altered chiral environments for one or more of the seven tyrosines and 12 phenylalanines (
TS does not contain tryptophan). The nearly featureless spectrum of D130A
TS between 270 nm and 290 nm suggests that the tyrosine side chains are mobile. Small peaks near 260 nm imply that one or more of the phenylalanines experience chiral packing environments.
These data show that the elimination of any one of the F19D46, I97D124, or A103D130 H-bonds and/or the loss of the negative charge when aspartic acid is replaced by alanine have a substantial impact on the secondary and tertiary structure of
TS. The very similar response in the far-UV CD spectrum for all three variants (Fig. 2D) and two of the three variants in the near-UV CD spectrum (Fig. 2H) suggest that the loss of either the F19D46 or the I97D124 noncovalent interactions lead to similar partially folded structures. These folds retain a substantial fraction of the secondary structure and maintain tertiary packing around some of the aromatic side chains. The contrasting near-UV CD spectrum for D130A
TS shows that at least some aspect of the tertiary structure of this variant differs from those of the other two clamp variants.
Equilibrium folding properties of the clamp-deletion variants
Urea denaturation experiments were performed to test the contribution of the three main-chain amideaspartic acid side-chain H-bonds to the thermodynamic properties of
TS. The dependence of the MRE at 222 nm on the urea concentration for each variant is compared to WT
TS in Figure 3. The equilibrium unfolding curve of the WT protein displays two overlapping sigmoidal transitions that have previously been assigned to interconversions between the native, N, and intermediate, I1, states and between the I1 and the unfolded-like, I2, and unfolded, U, states (Bilsel et al. 1999). The far-UV CD spectra of I2 and U are identical, permitting the application of a three-state model,
, to fit the equilibrium data. The I1 state for WT
TS reaches a maximal population of
60% at 3 M urea (Gualfetti et al. 1999), and its presence is evident in the change in slope of the transition curve at that urea concentration (Fig. 3).
|
TS. This fitting procedure was justified by the appearance of a limited native baseline below 0.5 M urea and a small shoulder near 2 M urea. As shown in Table 1, the free energy difference between the N and I1 states is greatly reduced, as is the urea-dependence of this free energy difference, that is, the m value, from the values for WT
TS. The
TS and WT protein are 1.17 ± 0.39 and 7.19 ± 0.58 kcal mol1, respectively; the respective m values are 1.19 ± 0.21 and 2.85 ± 0.24 kcal mol1 M1. This three-state fit of the equilibrium unfolding data for D130A was confirmed by a complementary kinetic measurement of the stability of the native state (see below).
|

clamp-deletion variants, D46A and D124A
TS, display broad rather featureless transition curves (Fig. 3) that could not be uniquely fit to a three-state model. In fact, the appearance of the unfolding profiles made it unclear whether the N, I1, and I2/U thermodynamic states were still being sampled in these two variants. As will be shown below, kinetic folding studies can be used to demonstrate the presence of stable folded states in the absence of denaturant and to estimate the stabilities of D46A and D124A
TS.
Kinetic folding properties of the clamp-deletion variants
To obtain insights into the role of the clamp interactions on the kinetic and thermodynamic properties of
TS, a comprehensive set of unfolding and refolding jumps at various final urea concentrations was performed on the clamp-deletion variants. A distinctive feature of the native state of WT
TS is its rate-limiting role in the kinetic unfolding reaction. WT
TS unfolds through two slow phases whose time constants decrease exponentially as the urea concentration is increased (Fig. 4A). It has previously been shown that these two exponential phases correspond to the independent unfolding of a major and a minor native conformer (Bilsel et al. 1999). These conformers reflect the native cis, 93%, and nonnative trans, 7%, isomers for the D27P28 peptide bond (Wu and Matthews 2002a). The subsequent unfolding reactions of the corresponding on-pathway I1 intermediates are far faster, rendering the pair of
reactions as the rate-limiting steps in unfolding for
TS. If the
TS variants, in the absence of denaturant, occupy a stable thermodynamic state related to the native state of the WT protein, they might be expected to display a vestige of the rate-limiting unfolding reactions.
|
TS all occupy distinct thermodynamic states that retain the major cis and minor trans isomers at D27P28 and the rate-limiting barriers to unfolding experienced by WT
TS.
The refolding of WT
TS detected by stopped-flow CD (SF-CD) and stopped-flow fluorescence (SF-FL) proceeds through a submillisecond burst-phase reaction followed by a hundreds of milliseconds reaction that accelerates as the urea concentration is increased (Fig. 4A). This "backtracking" phase reflects the unfolding of a misfolded off-pathway species and controls access to the subsequent on-pathway species, I1. cistrans isomerization at three key prolines, P28, P217, and P261, ultimately limits folding to the major and minor native species evident in the unfolding reaction (Bilsel et al. 1999; Wu and Matthews 2002a,b). These rate-limiting isomerization reactions are the source of the three slow refolding phases in the 10300 sec time range for the WT protein under strongly folding conditions (Fig. 4A).
By contrast, all three variants are missing two of the three slow urea-independent refolding phases observed in the WT protein (Fig. 4AD). Only two phases were detected under strongly folding conditions for D124A (Fig. 4C) and D130A (Fig. 4D)
TS. The relaxation time of the faster phase is weakly dependent on the final urea concentration and coincident with the backtracking phase in the WT protein. The relaxation time of the slower refolding phase was nearly independent of the urea concentration and in the time range of phases previously assigned to cis/trans isomerization reactions among a set of four I1 intermediates in parallel folding channels (Wu and Matthews 2002a). In addition to these two refolding phases, D46A also displays a refolding phase in the seconds time range that connects smoothly with an additional unfolding reaction. This additional phase has also been detected for the L48A variant (Wu et al. 2006), and is thought to reflect the interconversion of the I1 and I2 states. Finally, the measurement of the amplitude of the submillisecond CD signal at 222 nm as a function of the final urea concentration confirmed the presence of a burst-phase intermediate for all three variants (Fig. 4FH).
Taken together, the kinetic refolding data demonstrate that all of the clamp-deletion variants retain the submillisecond off-pathway misfolded intermediate and one of the proline isomerization reactions between two I1-like intermediates. Because the pair of unfolding phases for all three variants so closely agree in relaxation time and amplitude with those previously assigned to isomerization at P28, it is reasonable to suppose that the remaining isomerization reaction in refolding for all three clamp-deletion variants reflects the isomerization at P28. By implication, the P217 and P261 cis/trans isomerization reactions no longer limit access to the native conformation in the three clamp-deletion variants. The C terminus of the clamp-deletion variants, including residues P217 and P261, does not appear to be well folded or at least not thermodynamically coupled to the folded N terminus in their respective native states.
Thermodynamic parameters extracted from the kinetic experiments
The kinetic unfolding experiment that verified the presence of stable thermodynamic states for the clamp variants can be modified to determine the stability of their folded states. Because the amplitudes of the unfolding reactions are directly proportional to the populations of the folded states from which these reactions initiate, the stability of the folded states can be measured by monitoring the amplitudes of the unfolding reactions (to constant urea concentration) as a function of the initial urea concentration. The amplitudes of the slow unfolding phases are expected to decrease in a sigmoidal fashion as the initial urea concentration is increased into the zone where the native states become depopulated.
This procedure was validated on WT
TS (Fig. 5A), where a fit of the sigmoidal decrease in ellipticity at increasing initial urea concentrations to a two-state model yielded a stability of 7.06 ± 0.59 kcal mol1. This value is in excellent agreement with that for the
transition from the standard equilibrium titration experiment on the WT protein, 7.19 ± 0.58 kcal mol1 (Fig. 3; Table 1). The urea dependence of the stability, the m value, is also in excellent agreement, 2.90 ± 0.24 vs. 2.85 ± 0.24 kcal mol1 M1 (Table 1). Conformation of the kinetic measurement of stability was provided by the D130a
TS variant. The traditional equilibrium titration data (Fig. 3) yielded a stability of 1.17 ± 0.39 kcal mol1 and an m value 1.19 ± 0.21 kcal mol1 M1 (Table 1). The kinetic measurement yielded a stability of 1.05 ± 0.17 kcal mol1 and an m value 1.25 ± 0.16 kcal mol1 M1 (Table 1).
|
TS clamp-deletion variants were then determined in a similar fashion. Figure 5B displays a representative example of the amplitude of the major unfolding phase for D46A
TS as a function of the initial urea concentration. The stability of the folded state for D46A
TS, 1.98 ± 0.45 kcal mol1, is 5.21 kcal mol1 less than that of the N state for the WT protein; the m value is also significantly reduced, 0.78 ± 0.17 vs. 2.85 ± 0.24 kcal mol1 M1 (Table 1). The stability of the folded state for D124A
TS is decreased to 2.53 ± 0.40 kcal mol1. The m value for the
TS variant is also significantly decreased (Table 1).
Although the absence of a prominent shoulder in the MRE titration data for the D46A (Fig. 4F) and D124A (Fig. 4G) variants leaves the existence of an I1-like state ambiguous, the folded states for both proteins, as measured by the kinetic amplitude experiment, disappear; for example, at 34 M urea for D46A
TS (Fig. 5B), before the unfolded state is completely populated; for example, at 56 M urea for D46A
TS (Fig. 3). This discrepancy in the ellipticities implies the presence of a stable intermediate between 3 and 6 M urea. The stabilities of the intermediates for D46A and D124A
TS can be estimated by fitting the unfolding data for these clamp-deletion variants (Fig. 3) to a three-state
model, for which the parameters for the
reaction are fixed to those obtained from the kinetic unfolding experiments. The thermodynamic parameters obtained for the
transition are shown in Table 1.
The stabilities of the I1 states for all three clamp-deletion variants (relative to their I2/U states) and the m values for the
transition are within error of the values for the I1 state relative to the I2/U state for WT
TS. Because the removal of the clamp interaction by replacing aspartic acid with alanine has no significant effect on the stabilities of the intermediate states, it appears that each of the clamp interactions is broken in the intermediate of the WT protein.
Uniqueness of the clamp interactions
All three clamp interactions are located in the N-terminal half of
TS, which is known to serve as a core of stability in both the N and I1 states (Zitzewitz and Matthews 1999; Rojsajjakul et al. 2004). To determine whether the surprisingly large perturbations of structure and stability are related to the participation of these clamp H-bonds in this stable core, the S33A, D38A, and D112A variants of
TS were constructed. S33 forms an H-bond with the main-chain amide nitrogen of G26, D38 forms an electrostatic interaction with H92, and D112 forms an electrostatic interaction with H146. The heavy atom distances for these three interactions are 3.00, 2.75, and 2.61 Å, respectively, comparable to 2.92, 2.61, and 2.85 Å for the F19D46, I97D124, and A103D130 clamp H-bonds (Fig. 1B). The solvent accessible surface areas for the G26S33, D38H92, and D112H146 interactions are 16, 89, and 90 Å2, also comparable to those for the F19D46, I97D124, and A103D130 clamp H-bonds, 28, 19, and 80 Å2, respectively.
The D38A and D112A variants display a similar reduction in the far-UV CD signal as the clamp-deletion variants (Fig. 6A,B). However, both retain the distinctive near-UV CD spectrum for the WT protein (Fig. 6D,E), which suggests that they remain a similar tertiary structure as WT protein. The reductions in the far-UV CD signal, therefore, likely reflect the local unfolding of the surface helices and loops linked by these salt bridges. The far-UV (Fig. 6C) and near-UV CD (Fig. 6F) spectra of the S33A variant are both very similar to those for WT
TS. The unfolding transition curves for all three variants reveal native baselines and the biphasic behavior characteristic of a three-state process (Fig. 7A). Fits of these data to a three-state model show that the perturbations on stability for D38A, D112A, and S33A
TS are less than those of their clamp-deletion counterparts for the
transition and negligible for the
transition (Table 1). Thus, the effects on the structure and/or stability of the native state accompanying the loss of the 
or 
clamps in D46A, D124A, and D130A
TS are quantitatively different than those of other main-chainside-chain or side-chainside-chain electrostatic interactions in the N-terminal folding unit of
TS.
|
|
TS variant experiences a 20% reduction in the far-UV CD signal and retains the three-state equilibrium unfolding profile (Fig. 7B). The free energy change for the
TS (Table 1). The absence of a significant perturbation in the free energy change for the
TS. Kinetic analysis showed that D46N
TS retains the two unfolding phases observed for the WT protein; however, two of the three proline isomerization-limited refolding reactions are absent (data not shown). Thus, the isosteric replacement of aspartic acid with asparagine at position 46 in the N-terminal domain appears to preclude the proper folding and/or docking of the C-terminal region containing P217 and P261.
The ellipticity at 222 nm for the D46E
TS variant is smaller than that for any other variant examined, including that for D46A
TS (Fig. 7B). D46E
TS retains the pair of slow unfolding reactions characteristic of the native conformation of
TS (data not shown), and the kinetic measurement of the free energy change for the
transition is also greatly destabilized, 2.15 ± 0.46 kcal mol1 (Table 1). Similar to D46A and D46N, the free energy difference for the
reaction is within error of that for WT
TS. Neither asparagine nor glutamic acid can completely recapitulate the spectroscopic, thermodynamicm, and kinetic folding properties of aspartic acid at position 46.
| Discussion |
|---|
|
|
|---|
TS
or 
clamp interactions in
TS, F19D46, I97D124, or A103D130, results in surprisingly large decreases in the stability of the native state relative to the equilibrium intermediate, I1. Equally surprising was the appearance of what appear to be closely related partially folded states with substantially diminished secondary structure and altered tertiary structure for all three clamp-deletion variants. This tentative conclusion is based on (1) virtually identical far-UV CD spectra for all three clamp-deletion variants, and identical near-UV CD spectra for the D46A and D124A
TS variants; (2) similar reductions in the stability of the N state, 4.666.02 kcal mol1; and (3) a pair of unfolding reactions with identical relaxation times and relative amplitudes for all three proteins. The contrasting near-UV CD spectrum for D130A
TS (Fig. 2G) may reflect the peculiar and localized behavior of Y102, adjacent to the peptide bond containing the amide hydrogen of A103 with which D130 forms an H-bond. The observation of a common urea-independent slow refolding reaction, in the absence of the prolyl isomerization reactions for P217 and P261, implies that the C-terminal region of
TS cannot properly fold and/or dock on the N-terminal domain in the native states of all three clamp-deletion variants.
Other variants in which either aspartic acid, which serves as the anionic partner in surface electrostatic interactions, D38H92 and D112H146, or serine, which serves as an H-bond acceptor for a nearby main-chain amide nitrogen donor, G26S33, is replaced with alanine show smaller decreases in stability and lesser effects on the secondary and/or tertiary structure. The comparable range of H-bond lengths and solvent-accessible surface areas with the clamp-deletion variants suggest that other factors, possibly the concerted disruption of all three clamp H-bonds (see below), are responsible for the dramatically different effects. The observation that neither asparagine nor glutamic acid can completely recover the structural, thermodynamic, or kinetic folding properties of WT
TS implies, at least for this position, that both the formal charge and the precise positioning of that charge are crucial for full acquisition of the TIM barrel structure and stability.
The role of the clamp interactions in the folding mechanism of
TS
The retention of the burst-phase ellipticity and the hundreds of milliseconds unfolding reaction under refolding conditions, that is, the backtracking of the off-pathway burst-phase intermediate to the productive folding pathway, for all three clamp-deletion variants (Fig. 4FH) demonstrate that none of the clamps are essential for the misfolding reaction. The observation of a single urea-independent subsequent refolding reaction attributed to the trans
cis isomerization of the peptide main chain at D27P28 shows that the prolyl isomerizations at P217 and P261 are no longer rate limiting in folding for any of the clamp-deletion variants. By contrast, the retention of the cis isomer for P28 in the native state of all three variants implies that the tight turn between
1 and
1, responsible for stabilizing this higher energy form, is preserved in the absence of the clamps.
These results can be explained by examination of the properties of a stable C-terminal truncated version of
TS, 1188
TS (Zitzewitz and Matthews 1999). 1188
TS retains the burst-phase reaction, the early backtracking reaction, a single slow folding reaction, and two slow unfolding reactions. All but the unmeasured burst-phase reaction have relaxation times similar to full-length
TS. Although the secondary and tertiary structures of the clamp-deletion variants are somewhat disrupted compared to the WT 1188 fragment (with the clamp H-bonds), the similarities of the kinetic folding properties of 1188
TS (Zitzewitz and Matthews 1999) and the clamp-deletion variants suggest that the C-terminal region of each of the clamp-deletion variants is not well folded under native-favoring conditions. Apparently, the loss of any of the clamp interactions in the N terminus alters the structure of the N terminus so as to preclude the folding and docking of the C terminus. An earlier fragmentation analysis of
TS (Higgins et al. 1979) showed that the C-terminal region, residues 189268, is poorly folded in isolation. The spontaneous recovery of structure and function for
TS when these two fragments are complemented showed that the N-terminal region serves as the template for the folding and docking of the C-terminal region. If the template structure is altered, as appears to be the case for all three clamp-deletion variants, the C-terminal region cannot fold and dock on the N-terminal region.
It is also interesting to note that the two slow unfolding relaxation times exhibited by all three variants are within error of those for WT
TS at high denaturant concentration. Given the substantial loss in the stability of their respective native states, the similar destabilization of the unfolding transition states for all three variants means that all three clamp interactions in the WT
TS must be maintained in both unfolding transition states. Because the
reaction is reversible, these clamps must therefore also be very native-like in the transition states for the refolding reactions. The absence of the clamps in the I1 state, as implied by the lack of perturbation in the nature of
when aspartic acid is replaced by alanine (Table 1), implies that the clamps first appear in the final, rate-limiting transition state for folding, after the U
I2 and I2
I1 reactions are complete. By acting in concert to stabilize the transition state of WT
TS, the clamps direct the folding to the fully formed TIM barrel structure rather than to the partially folded states observed for the clamp variants. The stabilization of the N terminus of the native state for WT
TS by the clamps presumably enhances the folding of the C terminus by providing a properly structured template on which the C terminus can fold.
The closely related partially folded structures implied by the nearly identical far-UV CD spectra (Fig. 2D) and unfolding rates (Fig. 4BD) for the clamp-deletion variants suggest that the loss of a single clamp results in the rupture of the two remaining clamps. Thus, the unusually large decrease in stability observed for each clamp-deletion variant (Table 1) does not reflect the strength of a single side-chainmain-chain H-bond. This conjecture is supported by the substantial perturbations in the far- and near-UV CD spectra that imply the loss of numerous noncovalent interactions in these variants.
Structural implication of strong protection against HX for F19 NH and I97 NH
The absence of the F19D46 and I97D124 clamp interactions in the stable I1 folding intermediate for
TS (Table 1) begs the question as to the source of the unusually strong protection of these amide hydrogens against HX with solvent. If I1 were that source, one would have expected exchange at both of these positions in a few hours at 25°C and pH 7.8 via an EX1 mechanism (Englander and Kallenbach 1983; Bilsel et al. 1999). The implication of protection sufficient to retard exchange for days and months is either that these amide hydrogens have an alternative H-bond acceptor partner and/or are buried in nonnative nonpolar environments that inhibit exchange with solvent in the I1 state. Exchange must occur through other, higher energy states on the folding free energy surface, for example, the I2 state that lacks discernable secondary structure. Either implied explanation means that the N termini of at least two
-strands,
1 and
3, are involved in a very stable but nonnative structure(s) that differentiate the I1 and the N state. It is interesting to speculate that both aspartic acid H-bond donors are replaced by carbonyl oxygens in adjacent strands, that is,
2 and
4, whose lengths are extended to enable these additional stabilizing interactions. Whatever the source of these nonnative stabilizing interactions, they must be disrupted to enable formation if the clamp H-bonds are required to access the transition state leading to the fully folded TIM barrel.
Clamp interactions define the boundaries of 
and 
hairpins
The clamp interactions probed in this analysis are analogous to helix-capping interactions previously identified by Rose (Presta and Rose 1988). H-bond acceptor side chains near the N termini of helices can form H-bonds with unsatisfied main-chain amide hydrogen donors within but near the N termini of the helices. These H-bonds are typically sheltered from solvent by a pair of nonpolar side chains that presumably serve to increase the strength of the H-bond. Inspection of the local environments of each of the 
and 
clamp interactions in
TS shows hydrophobic clusters in those regions: A103 has van der Waals contacts with V128 and V131, I97 has contacts with L85 and I88, and F19 has contacts with V259 and M262. As for the helixcapping interactions, screening of the solvent by hydrophobic clusters would be expected to strengthen the clamp H-bond.
In addition to enhancing the stability of the native conformation, helixcapping interactions are thought to limit the length of the helices and, thereby, promote the formation of loops and turns that link the helices to neighboring elements of secondary structure. The clamp interactions may serve a similar role for 
and 
hairpins by delineating and reinforcing the boundaries of its constituent
strand and
helix. The specificity conferred by the formation of a clamp interaction would ensure the proper register of the strand and helix forming the hairpin. The identification of these long-range clamp interactions in the amino acid sequence of TIM barrel proteins of unknown structure could, therefore, assist in the prediction of their structures from the sequence. Nonpolar side chains at the N termini of predicted
strands or
helices that precede polar H-bond acceptor side chains at the C termini of subsequent predicted helices or strands would be candidates for the clamp H-bonds.
Broader involvement of clamp interactions
Similar 
and 
hairpin clamps exist in other TIM barrel proteins, including indole-3-glycerol phosphate synthase (IGPS) from Sulfolubus solfataricus (I107D128, T129E155), the iolI protein from Bacillus subtilis (L106N72), the cyclase moiety of imidazole-3-glycerol phosphate synthase (HisF) from Thermotoga maritime (F77D98), N-(5'-phosphoribosyl)anthranilate isomerase (PRAI) from Escherichia coli (I355D374), and the founding member of the motif, triosephosphate isomerase (TIM) from chicken (I206D227). Expansion of the search to other types of (
)n structures, for example, the
/
/
sandwich flavodoxin-fold, the Rossmann-fold, and the leucine-rich repeat-fold revealed the presence of one or more clamps in members of these motifs as well (S. Kathuria, unpubl.). Further experimentation is required to determine if the enhancement of stability gained by linking consecutive elements of secondary structure in
TS through these types of long-range noncovalent interactions is a general property of the TIM barrel and other (
)n classes of proteins.
| Materials and Methods |
|---|
|
|
|---|
Protein expression and purification
TS variant proteins were expressed in CB149 cells and purified as described previously (Wu and Matthews 2002b). Purity was demonstrated by the appearance of a single band by Coomassie blue-stained PAGE.
Circular dichroism
CD spectroscopy was employed to monitor the secondary structure and the tertiary structure near aromatic side chains. Spectra were obtained on a Jasco Model J-810 spectropolarimeter equipped with a thermoelectric cell holder. Far-UV CD data were collected from 280 nm to 185 nm at a scan rate of 50 nm/min and at 1-nm intervals using a 0.1-cm pathlength cell, a bandwidth of 2.5 nm, and an averaging time of 8 sec. Three replicate spectra were collected and averaged. The protein concentration was
5 µM. Near-UV CD data were collected from 320 nm to 250 nm at 5 nm/min using a 0.5-cm path length cell, and the protein concentration was 50150 µM. The temperature was maintained at 25°C with a computer-controlled Peltier system.
Equilibrium experiments
The stability of the
TS variants was measured by urea denaturation as described previously (Wu and Matthews 2002b) in a buffer containing 10 mM potassium phosphate, pH 7.8, 0.2 mM K2EDTA, and 1 mM
ME. A Hamilton 540B automatic titrator was used to prepare the samples to enhance the precision of the measurements. The samples were incubated overnight at 25°C to ensure equilibration.
Kinetic experiments
CD manual-mixing kinetic experiments were performed on a Jasco Model J-810 spectropolarimeter equipped with a thermoelectric cell holder using a 1-cm pathlength cell, a bandwidth of 2.5 nm, and an averaging time of 1 sec. The dead time of the experiments was 3 sec and the instrument response time was about 5 sec. The change in ellipticity as a function of time was monitored at 222 nm. Stopped-flow CD experiments were performed with an Aviv model 202SF stopped-flow CD spectrometer. The cell path length was 1.0 mm, and the dead time for mixing was 5 msec. Stopped-flow fluorescence experiments were performed on an Applied Photophysics SX-18MV stopped-flow fluorometer. The excitation wavelength was 280 nm, and the emission above 300 nm was detected through a cutoff filter.
Kinetic unfolding experiments to determine the stability of the folded states of the clamp-deletion variants were performed by jumping from different initial urea concentration (05.6 M) to a final concentration 6 M urea. Protein samples were first equilibrated in the initial urea concentration overnight and jumped to corresponding urea buffer solutions by a 1:10 dilution. The final protein concentration was around 5 µM.
Data analysis
Equilibrium data were fit to a three-state model (Bilsel et al. 1999), and the kinetic data were fit to a sum of exponentials and a constant, as described previously (Bilsel et al. 1999). All thermodynamic and kinetic folding data were fit using Savuka version 6.2, an in-house, nonlinear, least-squares program.
| Footnotes |
|---|
Abbreviations: H-bond, hydrogen bond;
TS,
subunit of tryptophan synthase from Escherichia coli; HSQC, heteronuclear single quantum coherence; CD, circular dichroism; MRE, mean residue ellipticity; SF, stopped flow; FL, fluorescence; HX, hydrogen exchange; I1, equilibrium intermediate of
TS populated at
3 M urea; I2, equilibrium intermediate of
TS populated at
5 M urea; NMR, nuclear magnetic resonance; N, native state; U, unfolded state; WT, wild type.
Article and publication are at http://www.proteinscience.org/cgi/doi/10.1110/ps.062704507.
| Acknowledgments |
|---|
|
|
|---|
| References |
|---|
|
|
|---|
Bai, Y., Sosnick, T.R., Mayne, L., and Englander, S.W. 1995. Protein folding intermediates: Native-state hydrogen exchange. Science 269: 192197.
Bilsel, O., Zitzewitz, J.A., Bowers, K.E., and Matthews, C.R. 1999. Folding mechanism of the alpha-subunit of tryptophan synthase, an alpha/beta barrel protein: Global analysis highlights the interconversion of multiple native, intermediate, and unfolded forms through parallel channels. Biochemistry 38: 10181029.[CrossRef][Medline]
Delano, W.L. 2002. The PyMOL molecular graphics system. DeLano Scientific, Palo Alto, CA.
Englander, S.W. and Kallenbach, N.R. 1983. Hydrogen exchange and structural dynamics of proteins and nucleic acids. Q. Rev. Biophys. 16: 521655.[Medline]
Eriksson, A.E., Baase, W.A., Zhang, X.J., Heinz, D.W., Blaber, M., Baldwin, E.P., and Matthews, B.W. 1992. Response of a protein structure to cavity-creating mutations and its relation to the hydrophobic effect. Science 255: 178183.
Gualfetti, P.J., Bilsel, O., and Matthews, C.R. 1999. The progressive development of structure and stability during the equilibrium folding of the alpha subunit of tryptophan synthase from Escherichia coli . Protein Sci. 8: 16231635.[Abstract]
Higgins, W., Fairwell, T., and Miles, E.W. 1979. An active proteolytic derivative of the
subunit of tryptophan synthase. Identification of the site of cleavage and characterization of the fragments. Biochemistry 18: 48274835.[CrossRef][Medline]
Horovitz, A., Serrano, L., Avron, B., Bycroft, M., and Fersht, A.R. 1990. Strength and co-operativity of contributions of surface salt bridges to protein stability. J. Mol. Biol. 216: 10311044.[Medline]
Hyde, C.C., Ahmed, S.A., Padlan, E.A., Miles, E.W., and Davies, D.R. 1988. Three-dimensional structure of the tryptophan synthase alpha 2 beta 2 multienzyme complex from Salmonella typhimurium . J. Biol. Chem. 263: 1785717871.
Ibarra-Molero, B., Zitzewitz, J.A., and Matthews, C.R. 2004. Salt-bridges can stabilize but do not accelerate the folding of the homodimeric coiled-coil peptide GCN4-p1. J. Mol. Biol. 336: 989996.[CrossRef][Medline]
Jackson, S.E., Moracci, M., elMasry, N., Johnson, C.M., and Fersht, A.R. 1993. Effect of cavity-creating mutations in the hydrophobic core of chymotrypsin inhibitor 2. Biochemistry 32: 1125911269.[CrossRef][Medline]
Krishna, M.M., Hoang, L., Lin, Y., and Englander, S.W. 2004. Hydrogen exchange methods to study protein folding. Methods 34: 5164.[CrossRef][Medline]
McDonald, I.K. and Thornton, J.M. 1994. Satisfying hydrogen bonding potential in proteins. J. Mol. Biol. 238: 777793.[CrossRef][Medline]
Myers, J.K. and Pace, C.N. 1996. Hydrogen bonding stabilizes globular proteins. Biophys. J. 71: 20332039.[Medline]
Nishio, K., Morimoto, Y., Ishizuka, M., Ogasahara, K., Tsukihara, T., and Yutani, K. 2005. Conformational changes in the alpha-subunit coupled to binding of the beta 2-subunit of tryptophan synthase from Escherichia coli: Crystal structure of the tryptophan synthase alpha-subunit alone. Biochemistry 44: 11841192.[CrossRef][Medline]
Presta, L.G. and Rose, G.D. 1988. Helix signals in proteins. Science 240: 16321641.
Rojsajjakul, T., Wintrode, P., Vadrevu, R., Matthews, C.Robert, and Smith, D.L. 2004. 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. J. Mol. Biol. 341: 241253.[CrossRef][Medline]
Serrano, L., Kellis Jr, J.T., Cann, P., Matouschek, A., and Fersht, A.R. 1992. The folding of an enzyme. II. Substructure of barnase and the contribution of different interactions to protein stability. J. Mol. Biol. 224: 783804.[CrossRef][Medline]
Shortle, D., Stites, W.E., and Meeker, A.K. 1990. Contributions of the large hydrophobic amino acids to the stability of staphylococcal nuclease. Biochemistry 29: 80338041.[CrossRef][Medline]
Vadrevu, R., Falzone, C.J., and Matthews, C.R. 2003. Partial NMR assignments and secondary structure mapping of the isolated alpha subunit of Escherichia coli tryptophan synthase, a 29-kD TIM barrel protein. Protein Sci. 12: 185191.
Vaughan, C.K., Harryson, P., Buckle, A.M., and Fersht, A.R. 2002. A structural double-mutant cycle: Estimating the strength of a buried salt bridge in barnase. Acta Crystallogr. D Biol. Crystallogr. 58: 591600.[CrossRef][Medline]
Wu, Y. and Matthews, C.R. 2002a. A cis-prolyl peptide bond isomerization dominates the folding of the alpha subunit of Trp synthase, a TIM barrel protein. J. Mol. Biol. 322: 713.[CrossRef][Medline]
Wu, Y. and Matthews, C.R. 2002b. Parallel channels and rate-limiting steps in complex protein folding reactions: Prolyl isomerization and the alpha subunit of Trp synthase, a TIM barrel protein. J. Mol. Biol. 323: 309325.[CrossRef][Medline]
Wu, Y., Vadrevu, R., Kathuria, S., Yang, X., and Matthews, C.R. 2006. A tightly packed hydrophobic cluster directs the formation of an off-pathway sub-millisecond folding intermediate in the alpha subunit of tryptophan synthase, a TIM barrel protein. J. Mol. Biol. 366: 16241638.[CrossRef][Medline]
Zitzewitz, J.A. and Matthews, C.R. 1999. 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. Biochemistry 38: 1020510214.[CrossRef][Medline]
![]()
CiteULike
Connotea
Del.icio.us
Digg
Reddit
Technorati What's this?
| ||||||