Protein Science
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Published online before print January 4, 2005, 10.1110/ps.041094205
Protein Science (2005), 14:351-359. Published by Cold Spring Harbor Laboratory Press. Copyright © 2005 The Protein Society
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
ps.041094205v1
14/2/351    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Kim, J.
Right arrow Articles by Blaber, M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Kim, J.
Right arrow Articles by Blaber, M.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Reddit   Add to Technorati  
What's this?

Sequence swapping does not result in conformation swapping for the {beta}4/{beta}5 and {beta}8/{beta}9 {beta}-hairpin turns in human acidic fibroblast growth factor

Jaewon Kim, Jihun Lee, Stephen R. Brych, Timothy M. Logan and Michael Blaber

Kasha Institute of Molecular Biophysics and Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4380, USA

Reprint requests to: Michael Blaber, Kasha Institute of Molecular Biophysics, Florida State University, Tallahassee, FL 32306-4380, USA; e-mail: blaber{at}sb.fsu.edu; fax: (850) 644-7244.

(RECEIVED September 1, 2004; FINAL REVISION September 25, 2004; ACCEPTED September 27, 2004)


    Abstract
 TOP
 Abstract
 Introduction
 Results
 Discussion
 Materials and methods
 References
 
The {beta}-turn is the most common type of nonrepetitive structure in globular proteins, comprising ~25% of all residues; however, a detailed understanding of effects of specific residues upon {beta}-turn stability and conformation is lacking. Human acidic fibroblast growth factor (FGF-1) is a member of the {beta}-trefoil superfold and contains a total of five {beta}-hairpin structures (antiparallel {beta}-sheets connected by a reverse turn). {beta}-Turns related by the characteristic threefold structural symmetry of this superfold exhibit different primary structures, and in some cases, different secondary structures. As such, they represent a useful system with which to study the role that turn sequences play in determining structure, stability, and folding of the protein. Two turns related by the threefold structural symmetry, the {beta}4/{beta}5 and {beta}8/{beta}9 turns, were subjected to both sequence-swapping and poly-glycine substitution mutations, and the effects upon stability, folding, and structure were investigated. In the wild-type protein these turns are of identical length, but exhibit different conformations. These conformations were observed to be retained during sequence-swapping and glycine substitution mutagenesis. The results indicate that the {beta}-turn structure at these positions is not determined by the turn sequence. Structural analysis suggests that residues flanking the turn are a primary structural determinant of the conformation within the turn.

Keywords: {beta}-hairpin; fibroblast growth factor; folding kinetics; stability; turn conformation

Abbreviations: FGF-1, human acidic fibroblast growth factor • NOE, nuclear Overhauser effect • NOESY, nuclear Overhauser spectroscopy • HSQC, heteronuclear single quantum correlation • TROSY, transverse-relaxation optimized spectroscopy • TOCSY, total correlation spectroscopy • ADA, N-(2-acetamido)iminodiacetic acid • DTT, dithiothreitol • GuHCl, guanidinium hydrochloride • RT, room temperature (25°C) • FGFR, fibroblast growth factor receptor • the single-letter amino acid code is utilized in the description of FGF-1 mutations • r.m.s., root-mean-square

Article published online ahead of print. Article and publication date are at http://www.proteinscience.org/cgi/doi/10.1110/ps.041094205.


    Introduction
 TOP
 Abstract
 Introduction
 Results
 Discussion
 Materials and methods
 References
 
Different types of secondary structures are recognized in proteins: helix, sheet, turns, and loops. Among these, turns and loops have received the least attention due to wide conformational variability and lack of simple model systems. Several successfully designed {beta}-hairpins have been reported, including monomeric {beta}-hairpins (Stanger and Gellman 1998; de Alba et al. 1999a) and turns associated with a three-stranded antiparallel {beta}-sheet (Kortemme et al. 1998; Schenck and Gellman 1998; de Alba et al. 1999b). With these model systems, a close relationship was observed between turn stability and turn sequence. Furthermore, the effects of different amino acids, at specific locations within particular types of {beta}-turns, upon stability and folding are being elucidated (Gibbs et al. 2002; Blandl et al. 2003; Kim et al. 2003; Rotondi and Gierasch 2003). However, more than 50% of turns have at least one residue in common with related turn types (Hutchinson and Thornton 1994); thus, the larger structural context within which turn sequences exist appears to be key to understanding their effects upon structure, thermodynamics, and folding kinetics.

Human acidic fibroblast growth factor (FGF-1) has a total of 12 {beta}-strands, forming a series of five {beta}-hairpin structures characteristic of the {beta}-trefoil superfold (Fig. 1Go; Murzin et al. 1992; Orengo et al. 1994). Among these five turns, {beta}4/{beta}5 (i.e., connecting {beta}-strands 4 and 5) and {beta}8/{beta}9 are related by the threefold tertiary structure symmetry inherent in the {beta}-trefoil architecture. Turn {beta}4/{beta}5 has been argued to contribute to different binding specificities of FGF homologous factors (Olsen et al. 2003). This turn in FGF-1 has two overlapping conformations: a type IV turn (an approximate 180° turn in the main chain direction, but residues i+1 and i+2 within the turn do not exhibit characteristic {phi}, {psi} angles) spanning residue positions 48–51 (i.e., residues i through i+3 of the turn) and a type II turn (i+1 {phi}~–60° {psi} ~+120°, i+2 {phi}~+90° {psi} ~0°) spanning residue positions 50–53 (Hutchinson and Thornton 1996). However, based on the crystal structure (Blaber et al. 1996) residue positions 49–52 within turn {beta}4/{beta}5 correspond to residues 90–93 within turn {beta}8/{beta}9 (Fig. 2Go). In this regard, there are no insertions or deletions when comparing the main chain regions that comprise turns {beta}4/{beta}5 and {beta}8/{beta}9, yet these turns exhibit significantly different turn conformations. Residue positions 90–93 in turn {beta}8/{beta}9 are defined as a type I turn (type I {beta}-turn plus G1-type {beta}-bulge; i+1 {phi}~ –60° {psi} ~–30°, i+2 {phi}~ –90° {psi}~0°). The importance of this particular turn upon the protein stability and function has been discussed in several reports. An extensive contact between turn {beta}8/{beta}9 and fibroblast growth factor receptor (FGFR) has been described (Pellegrini et al. 2000). Alternate conformations for the {beta}8/{beta}9 turn have been reported, and have been proposed as a possible basis of various functionalities (Hecht et al. 2001; Fernandez-Tornero et al. 2003); however, a type I turn conformation for this region appears to be adopted in the absence of ligand interactions or crystal packing influences (Kim et al. 2002).



View larger version (38K):
[in this window]
[in a new window]
 
Figure 1. (Upper panel) Relaxed stereo-ribbon diagram of FGF-1 (view looking down the threefold axis of tertiary structure symmetry) and showing the locations of the {beta}-strands comprising the {beta}4/{beta}5 (gray) and {beta}8/{beta}9 (dark gray) {beta}-hairpin turns. The locations of the N and C termini are also indicated. (Lower panel) Primary structure of FGF-1 (single-letter code) aligned according to the threefold tertiary structure symmetry and indicating the locations of the {beta}4, {beta}5, {beta}8, and {beta}9 {beta}-strands.

 


View larger version (54K):
[in this window]
[in a new window]
 
Figure 2. Relaxed stereo diagram of an overlay of the {beta}4/{beta}5 and {beta}8/{beta}9 turn regions. For purposes of clarity, only the main chain atoms are shown. The hydrogen bond interactions between the adjacent {beta}-strands are also indicated.

 
Turns {beta}4/{beta}5 and {beta}8/{beta}9 have different primary structures within the structurally heterogeneous region (S50, V51, G52 in turn {beta}4/{beta}5, and E91, N92, H93 in turn {beta}8/{beta}9) (Fig. 1Go), and an important question is whether the different turn conformations are a consequence of the heterogeneous primary structures within these regions. Furthermore, do the different turn conformations exert differential contributions toward the stability and folding of the protein? In an attempt to answer these questions, we constructed a series of mutations designed to exchange the primary structures of both turns. Furthermore, we also constructed a series of glycine mutations within these turn regions to evaluate the nature of strain inherent within the turns. Folding kinetics, stability studies, and X-ray crystallographic and NMR analyses were used to quantify the effects of these mutations. The results indicate that the residues within the turns are not the primary determinant of the turn structure, but do have a significant influence on protein stability and folding. The primary determinants of the turn structure reside within the local environment, and in particular, are postulated to involve residue positions that flank the turns.


    Results
 TOP
 Abstract
 Introduction
 Results
 Discussion
 Materials and methods
 References
 
Substitution of the {beta}8/{beta}9 turn primary structure into the {beta}4/{beta}5 turn
Mutations were constructed that stepwise substituted the {beta}8/{beta}9 turn primary structure into the {beta}4/{beta}5 turn, including V51N and G52H point mutants, S50E/V51N double mutant, and S50E/V51N/G52H triple mutant. The V51N point mutant and the S50E/V51N double mutant had no effect on protein stability (Table 1Go) or folding kinetics (Table 2Go). However, the G52H point mutation destabilized the protein by 5.5 kJ/mol, essentially due to a fourfold reduction in the folding rate (Table 2Go). The S50E/V51N/G52H triple mutant destabilized the protein by 3.6 kJ/mol, again, primarily due to a decrease (approximately twofold) in the rate of folding.


View this table:
[in this window]
[in a new window]
 
Table 1. Thermodynamic parameters for WT FGF-1 and mutant proteins determined from isothermal equilibrium denaturation
 

View this table:
[in this window]
[in a new window]
 
Table 2. Folding and unfolding kinetic parameters for WT FGF-1 and mutant proteins
 
Substitution of the {beta}4/{beta}5 turn primary structure into the {beta}8/{beta}9 turn
Mutations were constructed that stepwise substituted the {beta}4/{beta}5 turn primary structure into the {beta}8/{beta}9 turn, including a previously described H93G point mutant (Kim et al. 2002), an E91S/N92V double mutant, and an E91S/N92V/H93G triple mutant. The E91S/N92V double mutant destabilized the protein by 5.7 kJ/mol (Table 1Go). This destabilizing effect is due to both a fivefold reduction in the rate of folding and a threefold increase in the rate of unfolding (Table 2Go). The E91S/N92V/H93G triple mutant destabilized the protein by 2.8 kJ/mol, and is due primarily to a twofold increase in the rate of unfolding with only a slight decrease in the rate of folding (Table 2Go). Relative to the H93G mutant, the triple mutant E91S/N92V/H93G was destabilized by a substantial 9.6 kJ/mol, due to both an 18-fold reduction in the rate of folding and an approximately fourfold increase in the rate of unfolding (Table 2Go).

Glycine substitutions in the {beta}4/{beta}5 and {beta}8/{beta}9 turn regions
Three different combination glycine mutations (E49G, S50G/V51G, and E49G/S50G/V51G) were created within the {beta}4/{beta}5 turn. All three mutants were marginally more stable than the wild-type protein, with {Delta}{Delta}G values ranging from –0.4 to –1.3 kJ/mol (Table 1Go). The unfolding rate for all three mutants was essentially unaffected, and the observed modest increases in stability are due primarily to increases in the rate of folding (Table 2Go). In reference to the G52H background, the introduction of glycine residues at positions 50 and 51 resulted in a slight increase in stability of –2.1 kJ/mol, due primarily to a twofold increase in the rate of folding. E91G, N92G, and H93G point mutations as well as an E91G/N92G double mutant in FGF-1 were characterized as part of a prior study of the {beta}8/{beta}9 loop region (Kim et al. 2002), and a combination triple glycine mutation, E91G/N92G/H93G, is reported here. The introduction of three glycine substitutions in the {beta}8/{beta}9 turn resulted in an increase in stability of –2.6 kJ/mol in comparison to the wild-type protein. As with the glycine mutations in the {beta}4/{beta}5 turn, this increase in stability is due almost exclusively to an increase in the rate of folding, while the rate of unfolding is largely unaffected.

X-ray structure determinations
The V51N point mutation, S50G/V51G double mutant, and S50E/V51N double mutant yielded diffraction quality crystals, each in space group C2221, and produced high-resolution data sets with good completion (Table 3Go). Each structure refined to similar values of Rfree (21.4%–22.7%) and Rcryst (18.5%–19.7%) and with appropriate stereochemistry. The 2Fo – Fc difference electron density was unambiguous at the mutation site(s), and the mutant structures could be accurately modeled in each case.


View this table:
[in this window]
[in a new window]
 
Table 3. Crystallographic data collection and refinement statistics
 
NMR analysis
The backbone HN, N, CA, CB, and CO resonances of wild-type FGF were assigned (with the exception of residue E81) from correlations obtained in triple resonance NMR spectra, and are in excellent agreement with previous reports (Pineda-Lucena et al. 1994; Ogura et al. 1999). The wild-type assignments were used to guide the backbone assignments for the E91S/N92V/H93G triple mutant, and assignments were obtained for all residues except E81, Y74, and Y94. Figure 3Go compares the backbone HN-HN NOE pattern obtained in a 3D 15N-separated NOESY-HSQC spectrum for the residues found in the {beta}4/{beta}5 and {beta}8/{beta}9 turns within the E91S/N92V/H93G triple mutant. The NOE cross-peaks between residues 52 and 53, combined with the lack of cross-peaks between residues 51 and 52, are consistent with a type II turn for residue positions 50–53 (as identified from the X-ray structure [Blaber et al. 1996]). Despite the sequence identity within the {beta}8/{beta}9 turn region for the E91S/N92V/H93G triple mutant, the NOE cross-peaks between residues 91 and 92 and between 92 and 93 confirm a type I turn for residue positions 90–93 (identical to the wild-type {beta}8/{beta}9 turn).



View larger version (36K):
[in this window]
[in a new window]
 
Figure 3. A composite plot of several slices taken from the 3D 15N-separated NOESY-HSQC spectrum of the E91S/N92V/H93G triple mutant. The left panel shows the HN-HN NOEs observed for residues V51-G52-G53 of turn {beta}4/{beta}5; the right panel shows the HN-HN NOEs observed for residues E91S-N92V-H93G in turn {beta}8/{beta}9. The lines indicate the location of sequential NOEs observed for each residue. The NOE cross-peaks between residues 52 and 53, combined with the lack of cross-peaks between residues 51 and 52, are consistent with a type II turn for residue positions 50–53. Despite the sequence identity within the {beta}8/{beta}9 turn region for the E91S/N92V/H93G triple mutant, the NOE cross-peaks between residues 91 and 92, and between 92 and 93, confirm a type I turn for residue positions 90–93.

 

    Discussion
 TOP
 Abstract
 Introduction
 Results
 Discussion
 Materials and methods
 References
 
While the residues comprising the {beta}4/{beta}5 and {beta}8/{beta}9 turn regions (i.e., positions 50–52 and 91–93, respectively) differ in conformation, the adjacent {beta}-sheet residues (43–49, 53–59 and 84–90, 94–100) exhibit a high degree of structural identity, and the main chain atoms overlay with an r.m.s. deviation of 0.6 Å (Fig. 2Go). The substitution of the {beta}4/{beta}5 turn primary structure into the {beta}8/{beta}9 turn, and the converse substitution of the {beta}8/{beta}9 turn primary structure into the {beta}4/{beta}5 turn, were observed to destabilize the protein in both cases. The degree of destabilization was of similar magnitude (2.8 and 3.6 kJ/mol, respectively), although it was achieved by differential effects on folding and unfolding rates (i.e., the mutation of the {beta}4/{beta}5 turn destabilized primarily through a decrease in the folding rate, whereas the mutation of the {beta}8/{beta}9 turn destabilized primarily through an increase in the unfolding rate). These results indicate that the primary structure within these turns, in both cases, is optimized to interact with the local structural environment.

Given that the {beta}4/{beta}5 and {beta}8/{beta}9 turn regions exhibit different conformations in the wild-type protein, an important question is whether the swapping of primary structure within these turns results in a corresponding conformational change. A crystal structure was obtained for the S50E/V51N double mutant that converts the i+1 and i+2 residues within the {beta}4/{beta}5 turn to the equivalent {beta}8/{beta}9 primary structure. These mutations resulted in no significant structural change to the {beta}4/{beta}5 turn (Fig. 4Go). The H93G point mutation structure was reported previously (Kim et al. 2002), and results in no conformational change to the {beta}8/{beta}9 turn. Thus, the same sequence E(i+1), N(i+2), and G(i+3) introduced into the {beta}4/{beta}5 turn adopts the characteristic overlapping type IV and type II turns, but when introduced into the {beta}8/{beta}9 turn adopts the characteristic type I turn. NMR spectra were used to further analyze the conformation of the {beta}8/{beta}9 turn with the E91S/N92V/H93G triple mutant (i.e., introducing the {beta}4/{beta}5 primary structure into the {beta}8/{beta}9 turn). An observed intense NOE between residue 91 (i+1) and 92 (i+2) eliminates the possibility that the {beta}8/{beta}9 turn adopts a type II or type II' conformation. Likewise, an observed intense NOE between 92 (i+2) and 93 (i+3) eliminates the possibility that a type I' is formed in the {beta}8/{beta}9 turn. Thus, this finding is consistent with the interpretation that the {beta}4/{beta}5 sequence introduced into the {beta}8/{beta}9 turn did not result in a conversion to the {beta}4/{beta}5 turn conformation. Although we were unable to obtain structural data for the S50E/V51N/G52H triple mutant, the available data show that swapping the amino acid sequence of the {beta}4/{beta}5 and {beta}8/{beta}9 turns does not concomitantly convert the turn conformations. Therefore, we conclude that the local environment influences both the stability and the conformation of the {beta}4/{beta}5 and {beta}8/{beta}9 turns.



View larger version (41K):
[in this window]
[in a new window]
 
Figure 4. Relaxed stereo diagram illustrating an overlay of the {beta}4/{beta}5 turn regions for WT* (CPK coloring), V51N (red), S50G/V51G (magenta), and S50E/V51N (green) mutants. None of the sequence substitution, or polyglycine substitutions, within this turn affected the structure of the turn.

 
Glycine substitutions within ordered regions of a protein can destabilize the structure due to an increased entropy of unfolding (Matthews et al. 1987). However, such substitutions at turn positions known to fall outside the preferred positions of the Ramachandran plot can result in an increase in stability (due to the reduction in structural strain associated with eclipsed main chain/main chain and main chain/C{beta} groups) (Hutchinson and Thornton 1994; Kwasigroch et al. 1996; Takano et al. 2001; Kim et al. 2003). The S50G/V51G double mutant and E91G/N92G/H93G triple mutant both resulted in a modest stabilization of the protein (–1.3 and –2.6 kJ/mol, respectively). Thus, a poly-glycine sequence in both turns can effectively substitute for the wild-type sequence (position 52 in the {beta}4/{beta}5 turn is glycine in the wild-type protein). However, despite the similar overall stability effects of such a substitution, there appear to be notable differences in these turn regions with regard to structural strain.

A comparison of the wild-type protein with the S50G/V51G double mutant, and similarly, a comparison of the H93G point mutant with the E91G/N92G/H93G triple mutant, evaluates the effect of glycine mutations at the i+1 and i+2 positions within the {beta}4/{beta}5 and {beta}8/{beta}9 turns when the i+3 position is a glycine. In the case of the {beta}4/{beta}5 turn, the conversion of the i+1 and i+2 positions to glycine results in a modest stabilization of –1.3 kJ/mol, and does not result in any structural changes within the turn region (Fig. 4Go). In the case of the {beta}8/{beta}9 turn, the conversion of the i+1 and i+2 positions to glycine results in a significant destabilization of +4.0 kJ/mol (Table 1Go). A comparison of the H93G and G52H point mutations (i.e., comparing the effects of glycine vs. histidine at the i+3 position within the {beta}4/{beta}5 and {beta}8/{beta}9 turn regions) yields relatively consistent results for both turns, namely a destabilization of 5.5 to 6.6 kJ/mol upon the introduction of a histidine for glycine at this position. Thus, the {beta}8/{beta}9 turn i+3 histidine exhibits significant structural strain (being located in the left-handed {alpha}-helical region of the Ramachandran plot) (Kim et al. 2002), whereas the glycine residue at this symmetry-related position in the {beta}4/{beta}5 turn does not. Conversely, the i+1 and i+2 residues in the {beta}8/{beta}9 turn actively contribute to the stability of the structure (which is destabilized when converted to glycine), whereas similar positions in the {beta}4/{beta}5 turn do not. These results identify a fundamental difference between these two turns: the {beta}4/{beta}5 turn appears generally optimized and to have limited strain associated with it, whereas the {beta}8/{beta}9 turn exhibits an inherent structural strain, and can be further optimized to substantially (i.e., ~6.6 kJ/mol) increase protein stability. Thus, while a simple polyglycine substitution within both turn regions can be accommodated, additional stability is achievable with an appropriate selection of side chains (e.g., E91/N92/G93 in the {beta}8/{beta}9 turn region).

What structural features outside of the immediate turns might be responsible for influencing the structure of the turns? An examination of an overlay of the {beta}4/{beta}5 and {beta}8/{beta}9 turn regions provides some possible answers. If the {beta}8/{beta}9 turn adopted the {beta}4/{beta}5 turn structure, it would result in a 2.9 Å close contact between the main chain carbonyl oxygen of position E91 (equivalent to S50) and the side chain L89 C{delta}1 atom (Fig. 5Go). Associated with this contact is the hydrogen-bonding requirement of the main chain carbonyl group of E91, which could not be satisfied by such an interaction; thus, this conformation appears unfavorable. A more egregious close contact would occur if the {beta}4/{beta}5 turn adopted the {beta}8/{beta}9 turn structure. In this case an ~0.5 Å contact would occur between a side chain C{gamma} atom of V51 (equivalent to N92) and the side chain O{varepsilon}1 atom of E53 (which is involved in a salt bridge with adjacent K100) (Fig. 5Go). The above interactions involve residues in the primary structure that are adjacent to the turn residues (i.e., positions 48/89 and 53/94); thus, the conformation of the residues within the turns appears to be determined not by the primary structure of the turn residues, but rather by the residues flanking the turns. In this regard, the main chain atom overlay of the two turns (Fig. 2Go) indicates that the structural deviations initiate with the {psi} angles of residue positions 49 and 90, and extend through the {psi} angles of residues 53 and 94. These different main chain {phi}, {psi} angles may dictate the observed structural difference within the two turns; however, since the residues at these positions are neither proline or glycine, the differences in main chain angles must be due to side chain interactions. Again, the results suggest that such interactions do not involve the side chains within the immediate turn region.



View larger version (30K):
[in this window]
[in a new window]
 
Figure 5. Relaxed stereo diagram of a close-up of an overlay of the {beta}4/{beta}5 (light gray) and {beta}8/{beta}9 (dark gray) turns showing the local side chain interactions.

 
Despite the importance of {beta}-hairpin structures in proteins, the principles underlying their formation and stability are not well understood. There are two different views regarding the contributions of residues within turn regions upon conformation. One view is that turn residues can determine the turn conformation (de Alba et al. 1997). The other view is that turn regions in a protein may be able to maintain their conformation despite variations in their primary structure (Searle et al. 1995). Several groups have examined the influence of turn sequences on the formation and stability of {beta}-hairpins (Ramirez-Alvarado et al. 1997; Blanco et al. 1998; de Alba et al. 1999b; Griffiths-Jones et al. 1999; Chen et al. 2001). However, these studies used relatively small peptides, and the effects of adjacent regions outside of the turn (as would be present in larger globular proteins) could not be determined. In the present study, involving two symmetry-related turns within FGF-1, the data indicate that neither glycine substitution, nor sequence swapping, affected the local turn conformation. Furthermore, structural analysis suggests a role for flanking residues in determination of the turn conformation. These results support the general hypothesis of Searle et al. (1995) with regard to turn sequences and structures within proteins. Thus, predicting the stability and structure of turn regions within proteins requires knowledge of the larger structural context. However, the results suggest that a simple polyglycine primary structure will suffice for mesophilic applications. Finally, the results show that symmetry-related turn regions within a symmetric superfold can exhibit substantially different properties with regard to strain and contributions to overall folding kinetics.


    Materials and methods
 TOP
 Abstract
 Introduction
 Results
 Discussion
 Materials and methods
 References
 
Mutagenesis and expression
All studies utilized a synthetic gene for the 140-amino-acid form of human FGF-1 (Gimenez-Gallego et al. 1986; Linemeyer et al. 1990; Ortega et al. 1991; Blaber et al. 1996) with the addition of an amino-terminal six-His tag to facilitate purification (Brych et al. 2001). The QuikChange site-directed mutagenesis protocol (Stratagene) was used to introduce all point mutations, and all mutations were confirmed by nucleic acid sequence analysis (Biomolecular Analysis Synthesis and Sequencing Laboratory, Florida State University). Mutant nomenclature utilizes the single-letter amino acid code (e.g., H93G refers to a glycine substitution for a histidine residue at position 93). All expression and purification followed previously published procedures (Brych et al. 2001). Purified protein was exchanged into 20 mM N-(2-acetamido)iminodiacetic acid (ADA), 100 mM NaCl, 2 mM DTT (pH 6.60) ("ADA buffer") for biophysical studies.

Isothermal equilibration denaturation
Protein samples (40–10 µM) in various concentrations of GuHCl/ADA buffer were allowed to equilibrate overnight at room temperature (298 K). This study made use of the fluorescence signal of the single endogenous tryptophan residue at position 107. This residue is ~90% buried in the native structure (Blaber et al. 1996). Complete details of the instrumentation, data collection, and analysis procedure were reported (Blaber et al. 1999). Briefly, the fluorescence signal of FGF-1 is atypical in that W107 exhibits greater quenching in the native state rather than the denatured state. Excitation at 295 nm provides selective excitation of W107 in comparison with the six tyrosine residues that are present in the structure (Blaber et al. 1999; Brych et al. 2001). Triplicate scans were collected and averaged, and buffer traces were collected and subsequently subtracted from the protein scans. All scans were integrated to quantify the total fluorescence as a function of denaturant concentration. The data were analyzed using the general purpose nonlinear least-squares fitting program DataFit (Oakdale Engineering) implementing a six-parameter, two-state model (Eftink 1994):


(1)

where the denaturant concentration is given by [D]; the native state (0 M denaturant) fluorescence intercept and slope are F0N and SN, respectively; the denatured state fluorescence intercept and slope are F0D and SD, respectively; and the free energy of unfolding function intercept and slope are {Delta}G0 and m, respectively. The {Delta}G0 and m values describe the linear function of the free energy of unfolding as a function of denaturant under isothermal equilibrium conditions (where {Delta}G0 is the {Delta}G value extrapolated to 0 M denaturant, and m is reported as –d{Delta}G/d[D]). The midpoint of the transition, i.e., the denaturant concentration where {Delta}G = 0, is defined as Cm. The effect of a given mutation upon the stability of the protein ({Delta} {Delta}G) was calculated using the method of Pace (Pace and Scholtz 1997):


(2)

where a negative value indicates that the mutation is stabilizing in relationship to the wild-type protein.

Folding kinetic measurements
Denatured protein samples were prepared by overnight dialysis against 20 mM ADA, 100 mM NaCl, 2 mM DTT (pH 6.60) containing either 2.5 M or 3.0 M GuHCl (Blaber et al. 1999). Initial studies using manual mixing indicated that the relaxation times for the folding process were more appropriate for stopped-flow data collection. All folding kinetic data were collected using a Kintek SF2000 stopped-flow system (Kintek). Folding was initiated by a 1:10 dilution of 40 µM denatured protein into 20 mM ADA, 100 mM NaCl (pH 6.60) with denaturant concentrations varying from 0.25 M GuHCl, in increments of 0.05 M, to the midpoint of denaturation as determined by the above described isothermal equilibrium denaturation measurements. The data collection strategy was designed to span approximately five half-lives, or >97% of the expected fluorescence signal change between the fully denatured and native states.

Unfolding kinetic measurements
Native protein samples were dialyzed against 20 mM ADA 100 mM NaCl, 2 mM DTT (pH 6.60) overnight at 298 K. Unfolding was initiated by a 1:10 dilution of 25 µM protein into 20 mM ADA, 100 mM NaCl (pH 6.60) with concentrations of GuHCl varying between 1.5 M and 5 M. The fluorescence signal associated with protein unfolding was quantified using a Varian Eclipse fluorescence spectrophotometer, with a wavelength of 295 nm for excitation and 350 nm for emission, and maintained at 298 K with a temperature-controlled Peltier cell holder (Varian). The unfolding kinetics exhibited relaxation times that were appropriate for manual mixing techniques. The data collection strategy was designed to span approximately 2–3 half-lives, or >80% of the expected fluorescence signal change between the native and denatured states.

Kinetic analysis
Both folding and unfolding kinetic data were collected in triplicate at each denaturant concentration, with typically six runs per sample and all data being averaged. The kinetic rate constants and amplitudes, as a function of denaturant concentration, were determined from the time-dependent change in fluorescence intensity implementing a single exponential model:


(3)

where I(t) is the intensity of fluorescent signal at time t, A is the corresponding amplitude, k is the observed rate constant for the reaction, and C is a constant corresponding to the asymptotic signal limit.

Folding and unfolding rate constant data were fit to a global function describing the contribution of both rate constants to the observed kinetics as a function of denaturant ("Chevron" plot) as described by Fersht (1999):


(4)

where kf0 and ku0 are the folding and unfolding rate constants, respectively, extrapolated to 0 M denaturant; mkf and mku are the slopes of the linear functions relating ln(kf) and ln(ku), respectively, to denaturant concentration; and D is denaturant concentration.

Crystallization of FGF-1
Purified FGF-1 solutions (~13 mg/mL) were filtered through 0.2- µm filters (Pall Life Sciences) immediately prior to crystallization trials. X-ray diffraction quality crystals of S50G/V51G grew within 1 wk at RT in 10- µL hanging drops by vapor diffusion against 1-mL reservoirs of 3.0–3.2 M formate and 0.9 M ammonium sulfate in crystallization buffer. V51N and S50E/V51N crystals were obtained within 1 wk at RT in 10- µL hanging drops by vapor diffusion against 3.4–3.6 M formate.

Data collection, molecular replacement, and refinement
X-ray diffraction data were collected using a Rigaku RU-H2R rotating anode X-ray source (Rigaku) equipped with Osmic Blue, or Purple, confocal mirrors (MarUSA) coupled with either a Rigaku R-axis IIc image plate or a MarCCD 165 detector. The crystals were mounted using Hampton Research nylon-mounted cryo-turns and frozen in a stream of nitrogen gas at 100K. Diffraction data were indexed, integrated, and scaled using the DENZO software package (Otwinowski 1993; Otwinowski and Minor 1997). His-tagged wild-type FGF-1 was used as the search model in molecular replacement in each case. Model refinement was carried out with CNS (Brunger et al. 1998) using the maximum likelihood target function. Model building and visualization utilized the O molecular graphics program (Jones et al. 1991).

NMR studies
NMR samples of FGF-1 in 100 mM phosphate buffer, 200 mM ammonium sulfate, and 2 mM DTT (pH 6.0), were labeled uniformly with 15N or 13C using conditions described by Chi et al. (2002). The protein was concentrated to 1 mM, dried by lyophilization, and then resuspended by dissolving in 10% D2O containing 10 mM phosphate buffer with 90% H2O. The two- and three-dimensional NMR experiments were performed at 30°C on a 600 MHz Varian UnityPlus spectrometer equipped with waveform generators and pulsed field gradient accessories. Individual spin systems of the wild-type protein were identified via TROSY experiments and then confirmed by 3D HNCACB and CBCA(CO)NH. Individual spin systems of mutant proteins were assigned by HSQC and 3D 15N TOCSY-HSQC experiment.


    Acknowledgments
 
We thank Dr. E.G. Hutchinson, University of Reading, for kindly supplying updated turn classification, and Dr. Ewa Bienkiewicz of the Physical Biochemistry Facility at the Institute of Molecular Biophysics, Florida State University for assistance and helpful discussions. This work was supported by grant MCB 0314740 from the NSF.


    References
 TOP
 Abstract
 Introduction
 Results
 Discussion
 Materials and methods
 References
 
Blaber, M., DiSalvo, J., and Thomas, K.A. 1996. X-ray crystal structure of human acidic fibroblast growth factor. Biochemistry 35: 2086–2094.[CrossRef][Medline]

Blaber, S.I., Culajay, J.F., Khurana, A., and Blaber, M. 1999. Reversible thermal denaturation of human FGF-1 induced by low concentrations of guanidine hydrochloride. Biophys. J. 77: 470–477.[Abstract/Free Full Text]

Blanco, F., Ramirez-Alvarado, M., and Serrano, L. 1998. Formation and stability of {beta} hairpin structures in polypeptides. Curr. Opin. Struct. Biol. 8: 107–111.[CrossRef][Medline]

Blandl, T., Cochran, A.G., and Skelton, N.J. 2003. Turn stability in {beta}-hairpin peptides: Investigation of peptides containing 3:5 type I G1 bulge turns. Protein Sci. 12: 237–247.[Abstract/Free Full Text]

Brunger, A.T., Adams, P.D., Clore, G.M., DeLano, W.L., Gros, P., Grosse-Kunstleve, R.W., Jiang, J.-S., Kuszewski, J., Nilges, N., Pannu, N.S., et al. 1998. Crystallography and NMR system (CNS): A new software system for macromolecular structure determination. Acta Crystallogr. D Biol. Crystallogr. 54: 905–921.[CrossRef][Medline]

Brych, S.R., Blaber, S.I., Logan, T.M., and Blaber, M. 2001. Structure and stability effects of mutations designed to increase the primary sequence symmetry within the core region of a {beta}-trefoil. Protein Sci. 10: 2587–2599.[Abstract/Free Full Text]

Chen, P.Y., Lin, C.K., Lee, C.T., Jan, H., and Chan, S.I. 2001. Effects of turn residues in directing the formation of the {beta}-sheet and in the stability of the {beta}-sheet. Protein Sci. 10: 1794–1800.[Abstract/Free Full Text]

Chi, Y.H., Kumar, T.K., Chiu, I.M., and Yu, C. 2002. Identification of rare partially unfolded states in equilibrium with the native conformation in an all {beta}-barrel protein. J. Biol. Chem. 277: 34941–34948.[Abstract/Free Full Text]

de Alba, E., Jimenez, M.A., and Rico, M. 1997. Turn residue sequence determines {beta}-hairpin conformation in designed peptides. J. Am. Chem. Soc. 119: 175–183.[CrossRef]

de Alba, E., Rico, M., and Jimenez, M.A. 1999a. The turn sequence directs {beta}-strand alignment in designed {beta}-hairpins. Protein Sci. 8: 2234–2244.[Abstract]

de Alba, E., Santoro, J., Rico, M., and Jimenez, M.A. 1999b. De novo design of a monomeric three-stranded antiparallel {beta}-sheet. Protein Sci. 8: 854–865.[Abstract]

Eftink, M.R. 1994. The use of fluorescence methods to monitor unfolding transitions in proteins. Biophys. J. 66: 482–501.[Medline]

Fernandez-Tornero, C., Lozano, R.M., Redondo-Horcajo, M., Gomez, A.M., Lopez, J.C., Quesada, E., Uriel, C., Cuevas, P., Romero, A., and Gimenez-Gallego, G. 2003. Leads for development of new naphthalenesulfonate derivatives with enhanced antiangiogenic activity—Crystal structure of acidic fibroblast growth factor in complex with 5-amino-2-naphthalenesulfonate. J. Biol. Chem. 278: 21774–21781.[Abstract/Free Full Text]

Fersht, A.R. 1999. Kinetics of protein folding. W.H. Freeman and Co., New York.

Gibbs, A.C., Bjorndahl, T.C., Hodges, R.S., and Wishart, D.S. 2002. Probing the determinants of type II' {beta}-turn formation in peptides and proteins. J. Am. Chem. Soc. 124: 1203–1209.[CrossRef][Medline]

Gimenez-Gallego, G., Conn, G., Hatcher, V.B., and Thomas, K.A. 1986. The complete amino acid sequence of human brain-derived acidic fibroblast growth factor. Biochem. Biophys. Res. Commun. 128: 611–617.

Griffiths-Jones, S.R., Maynard, A.J., and Searle, M.S. 1999. Dissecting the stability of a {beta}-hairpin peptide that folds in water: NMR and molecular dynamics analysis of the {beta}-turn and {beta}-strand contributions to folding. J. Mol. Biol. 292: 1051–1069.[CrossRef][Medline]

Hecht, H.J., Adar, R., Hofmann, B., Bogin, O., Weich, H., and Yayon, A. 2001. Structure of fibroblast growth factor 9 shows a symmetric dimer with unique receptor- and heparin-binding interfaces. Acta Crystallogr. D Biol. Crystallogr. 57: 378–384.[CrossRef][Medline]

Hutchinson, E.G. and Thornton, J.M. 1994. A revised set of potentials for {beta}-turn formation in proteins. Protein Sci. 3: 2207–2216.[Abstract]

———. 1996. PROMOTIF—A program to identify and analyze structural motifs in proteins. Protein Sci. 5: 212–220.[Abstract]

Jones, T.A., Zou, J.Y., Cowan, S.W., and Kjeldgaard, M. 1991. Improved methods for the building of protein models in electron density maps and the location of errors in these models. Acta Crystallogr. A 47: 110–119.

Kim, J., Blaber, S.I., and Blaber, M. 2002. Alternative type I and I' turn conformations in the {beta}8/{beta}9 {beta}-hairpin of human acidic fibroblast growth factor. Protein Sci. 11: 459–466.[Abstract/Free Full Text]

Kim, J., Brych, S.R., Lee, J., Logan, T.M., and Blaber, M. 2003. Identification of a key structural element for protein folding within {beta}-hairpin turns. J. Mol. Biol. 328: 951–961.[CrossRef][Medline]

Kortemme, T., Ramirez-Alvarado, M., and Serrano, L. 1998. Design of a 20-amino acid, three-stranded {beta}-sheet protein. Science 281: 253–256.[Abstract/Free Full Text]

Kwasigroch, J.M., Chomilier, J., and Mornon, J.P. 1996. A global taxonomy of loops in globular proteins. J. Mol. Biol. 259: 855–872.[CrossRef][Medline]

Linemeyer, D.L., Menke, J.G., Kelly, L.J., Disalvo, J., Soderman, D., Schaeffer, M.-T., Ortega, S., Gimenez-Gallego, G., and Thomas, K.A. 1990. Disulfide bonds are neither required, present, nor compatible with full activity of human recombinant acidic fibroblast growth factor. Growth Factors 3: 287–298.[Medline]

Matthews, B., Nicholson, H., and Becktel, W. 1987. Enhanced protein thermostability from site-directed mutations that decrease the entropy of unfolding. Proc. Natl. Acad. Sci. 84: 6663–6667.[Abstract/Free Full Text]

Murzin, A.G., Lesk, A.M., and Chothia, C. 1992. {beta}-Trefoil fold. Patterns of structure and sequence in the kunitz inhibitors interleukins-1{beta}and 1{alpha} and fibroblast growth factors. J. Mol. Biol. 223: 531–543.[CrossRef][Medline]

Ogura, K., Nagata, K., Hatanaka, H., Habuchi, H., Kimata, K., Tate, S., Ravera, M.W., Jaye, M., Schlessinger, J., and Inagaki, F. 1999. Solution structure of human acidic fibroblast growth factor and interaction with heparin-derived hexasaccharide. J. Biomol. NMR 13: 11–24.[CrossRef][Medline]

Olsen, S.K., Garbi, M., Zampieri, N., Eliseenkova, A.V., Ornitz, D.M., Goldfarb, M., and Mohammadi, M. 2003. Fibroblast growth factor (FGF) homologous factors share structural but not functional homology with FGFs. J. Biol. Chem. 278: 34226–34236.[Abstract/Free Full Text]

Orengo, C.A., Jones, D.T., and Thornton, J.M. 1994. Protein superfamilies and domain superfolds. Nature 372: 631–634.[CrossRef][Medline]

Ortega, S., Schaeffer, M.-T., Soderman, D., DiSalvo, J., Linemeyer, D.L., Gimenez-Gallego, G., and Thomas, K.A. 1991. Conversion of cysteine to serine residues alters the activity, stability, and heparin dependence of acidic fibroblast growth factor. J. Biol. Chem. 266: 5842–5846.[Abstract/Free Full Text]

Otwinowski, Z. 1993. Oscillation data reduction program. In Proceeding of the CCP4 Study Weekend: Data Collection and Processing, 29–30 January 1993 (eds. L. Sawyer et al.), pp. 56–62. SERC Daresbury Laboratory, England.

Otwinowski, Z. and Minor, W. 1997. Processing of x-ray diffraction data collected in oscillation mode. Methods Enzymol. 276: 307–326.

Pace, C.N. and Scholtz, J.M. 1997. Measuring the conformational stability of a protein. In Protein structure: A practical approach (ed. T.E. Creighton), pp. 299–321. Oxford University Press, Washington, DC.

Pellegrini, L., Burke, D.F., von Delft, F., Mulloy, B., and Blundell, T.L. 2000. Crystal structure of fibroblast growth factor receptor ectodomain bound to ligand and heparin. Nature 407: 1029–1034.[CrossRef][Medline]

Pineda-Lucena, A., Jimenez, M.A., Nieto, J.L., Santoro, J., Rico, M., and Gimenez-Gallego, G. 1994. 1H-NMR assignment and solution structure of human acidic fibroblast growth factor activated by inositol hexasulfate. J. Mol. Biol. 242: 81–98.[CrossRef][Medline]

Ramirez-Alvarado, M., Serrano, L., and Blanco, F.J. 1997. Conformational analysis of peptides corresponding to all the secondary structure elements of protein L B1 domain: Secondary structure propensities are not conserved in proteins with the same fold. Protein Sci. 6: 162–174.[Abstract]

Rotondi, K.S. and Gierasch, L.M. 2003. Local sequence information in cellular retinoic acid binding protein I: Specific residue roles in {beta}-turns. Biopolymers 71: 638–651.[CrossRef][Medline]

Schenck, H.L. and Gellman, S.H. 1998. Use of a designed triple-stranded antiparallel {beta}-sheet to probe {beta}-sheet cooperativity in aqueous solution. J. Am. Chem. Soc. 120: 4869–4870.[CrossRef]

Searle, M.S., Williams, D.H., and Packman, L.C. 1995. A short linear peptide derived from the N-terminal sequence of ubiquitin folds into a water-stable non-native {beta}-hairpin. Nat. Struct. Biol. 2: 999–1006.[CrossRef][Medline]

Stanger, H.E. and Gellman, S.H. 1998. Rules for antiparallel {beta}-sheet design: D-Pro-Gly is superior to L-Asn-Gly for {beta}-hairpin nucleation. J. Am. Chem. Soc. 120: 4236–4237.[CrossRef]

Takano, K., Yamagata, Y., and Yutani, K. 2001. Role of non-glycine residues in left-handed helical conformation for the conformational stability of human lysozyme. Proteins 44: 233–243.[CrossRef][Medline]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us   Add to Digg Digg   Add to Reddit Reddit   Add to Technorati Technorati    What's this?



This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow All Versions of this Article:
ps.041094205v1
14/2/351    most recent
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Kim, J.
Right arrow Articles by Blaber, M.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation