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

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

The intestinal fatty acid binding protein: the role of turns in fast and slow folding processes

The intestinal fatty acid binding protein is one of a family of proteins that are composed of two beta-sheets surrounding a large interior cavity into which the ligand binds. Glycine residues occur in many of the turns between adjacent antiparallel beta-strands. In previous work, the effect of replacing these glycine residues with valine has been examined with stopped flow instrumentation using intrinsic tryptophan fluorescence spectroscopy [Kim and Frieden (1998) Protein Sci. 7, 1821-1828]. To resolve the burst phase missing in the stopped flow measurements, these valine mutants have been reexamined with sub-millisecond continuous flow instrumentation. Some of the glycine residues have also been replaced with proline, and the folding reactions of these proline mutants have been compared with those of their valine counterparts. In all cases, the stability of the protein is decreased, but some turns appear to be more critical for final structure stabilization than others. Surprisingly, the rate constants observed for all the mutants measured by sub-millisecond continuous flow methods are quite similar (1400-3000 s(-1)), and in all the mutants, there is a shift in the fluorescence emission maximum from that of the unfolded protein to lower wavelengths, suggesting some collapse of the unfolded state within 200 micros. In contrast to the rate constants observed for the initial folding events measured by the sub-millisecond continuous flow method, the rate constants for the slower phase observed in the stopped flow instrument vary widely for the different mutants. The latter step appears to be related to side chain stabilization rather than secondary structure formation. It is also shown that the ligand binds tightly only to the native protein and not to any intermediate forms.     

Hierarchical folding of intestinal fatty acid binding protein

Intestinal fatty acid binding protein (IFABP) is a member of the lipid binding protein family, members of which have a clam shell type of motif formed by two five-stranded beta-sheets. Understanding the folding mechanism of these proteins has been hindered by the presence of an unresolved burst phase. By initiating the reaction with a sub-millisecond mixer and following its progression by Trp fluorescence, we discovered three distinct phases in the folding reaction of the W6Y mutant of IFABP from which we postulate the following sequence of events. The first phase (k(1) > 10 000 s(-1)) involves collapse of the polypeptide chain around a hydrophobic core. During the second phase (k(2) approximately 1500 s(-1)), beta-strands B-G, mostly located on the top half of the clam shell structure, propagate from this hydrophobic core. It is followed by the final phase (k(3) approximately 5 s(-1)) involving the formation of the last three beta-strands on the bottom half of the clam shell and the establishment of the native hydrogen bonding network throughout the protein molecule.     

Consequences of single-site mutations in the intestinal fatty acid binding protein

The intestinal fatty acid binding protein (IFABP) is a small (15 kDa) protein consisting mostly of 10 antiparallel beta-strands (A-J) and a small helical region that serves as a portal for the ligand. Two beta-sheet structures (strands A-E and F-J) surround a cavity into which the ligand binds. In this work, we investigated how changes in the side chains of specific residues are propagated through the structure. To determine what these changes were and how they relate to changes in stability, (15)N chemical shift perturbations were measured and compared to those of the wild-type protein. Seven mutations, five of which change either valine or leucine to glycine, have been examined. All these mutants were less stable than wild-type IFABP, suggesting some structural changes. For five of the mutants, the data suggest that destabilization of a small region of the protein propagates throughout the structure, resulting in an overall decrease in stability. In two (Leu38Gly and Leu89Gly), the loss of cooperativity in the equilibrium denaturation curves suggests that the destabilization of one region may not be transmitted to other regions in a cooperative manner. It is shown that the effect of mutating hydrophobic residues is much greater than that observed upon mutation of a solvent-exposed polar residue.     

Partially unfolded forms and non-two-state folding of a beta-sandwich: FHA domain from Arabidopsis receptor kinase-associated protein phosphatase

FHA domains adopt a beta-sandwich fold with 11 strands. The first evidence of partially unfolded forms of a beta-sandwich is derived from native-state hydrogen exchange (NHX) of the forkhead-associated (FHA) domain from kinase-associated protein phosphatase from Arabidopsis. The folding kinetics of this FHA domain indicate that EX2 behavior prevails at pH 6.3. In the chevron plot, rollover in the folding arm and bends in the unfolding arm suggest folding intermediates. NHX of this FHA domain suggests a core of six most stable beta-strands and two loops, characterized by rare global unfolding events. Flanking this stable core are beta-strands and recognition loops with less stability, termed subglobal motifs. These suggest partially unfolded forms (near-native intermediates) with two levels of stability. The spatial separation of the subglobal motifs on the flanks suggests possible parallelism in their folding as additional beta-strands align with the stable core of six strands. Intermediates may contribute to differences in stabilities and m-values suggested by NHX or kinetics relative to chemical denaturation. Residual structure in the unfolded regime is suggested by superprotection of beta-strand 6 and by GdmCl-dependence of adjustments in amide NMR spectra and residual optical signal. The global folding stability depends strongly on pH, with at least 3 kcal/mol more stability at pH 7.3 than at pH 6.3. This FHA domain is hypothesized to fold progressively with initial hydrophobic collapse of its stable six-stranded core followed by addition of less stable flanking beta-strands and ordering of recognition loops.     

TEM-1 beta-lactamase folds in a nonhierarchical manner with transient non-native interactions involving the C-terminal region

The conformational stability and kinetics of refolding and unfolding of the W290F mutant of TEM-1 beta-lactamase have been determined as a function of guanidinium chloride concentration. The activity and spectroscopic properties of the mutant enzyme did not differ significantly from those of the wild type, indicating that the mutation has only a very limited effect on the structure of the protein. The stability of the folded protein is reduced, however, by 5-10 kJ mol-1 relative to that of the molten globule intermediate (H), but the values of the folding rate constants are unchanged, suggesting that Trp-290 becomes organized in its nativelike environment only after the rate-limiting step; i.e., the C-terminal region of the enzyme folds very late. In contrast to the significant increase in fluorescence intensity seen in the dead time (3-4 ms) of refolding of the wild-type protein, no corresponding burst phase was observed with the mutant enzyme, enabling the burst phase to be attributed specifically to the C-terminal Trp-290. This residue is suggested to be buried in a nonpolar environment from which it has to escape during subsequent folding steps. With both proteins, fast early collapse leads to a folding intermediate in which the C-terminal region of the polypeptide chain is trapped in a non-native structure, consistent with a nonhierarchical folding process.     

Engineering a compact non-native state of intestinal fatty acid-binding protein

The last three C-terminal residues (129-131) of intestinal fatty acid-binding protein (IFABP) participate in four main-chain hydrogen bonds and two electrostatic interactions to sequentially distant backbone and side-chain atoms. To assess if these interactions are involved in the final adjustment of the tertiary structure during folding, we engineered an IFABP variant truncated at residue 128. An additional mutation, Trp-6-->Phe, was introduced to simplify the conformational analysis by optical methods. Although the changes were limited to a small region of the protein surface, they resulted in an IFABP with altered secondary and tertiary structure. Truncated IFABP retains some cooperativity, is monomeric, highly compact, and has the molecular dimensions and shape of the native protein. Our results indicated that residues 129-131 are part of a crucial conformational determinant in which several long-range interactions, essential for the acquisition of the native state, are established. This work suggests that carefully controlled truncation can populate equilibrium non-native states under physiological conditions. These non-native states hold a great promise as experimental models for protein folding.     

Structural and sequence features of two residue turns in beta‐hairpins

Beta‐turns in beta‐hairpins have been implicated as important sites in protein folding. In particular, two residue β‐turns, the most abundant connecting elements in beta‐hairpins, have been a major target for engineering protein stability and folding. In this study, we attempted to investigate and update the structural and sequence properties of two residue turns in beta‐hairpins with a large data set. For this, 3977 beta‐turns were extracted from 2394 nonhomologous protein chains and analyzed. First, the distribution, dihedral angles and twists of two residue turn types were determined, and compared with previous data. The trend of turn type occurrence and most structural features of the turn types were similar to previous results, but for the first time Type II turns in beta‐hairpins were identified. Second, sequence motifs for the turn types were devised based on amino acid positional potentials of two‐residue turns, and their distributions were examined. From this study, we could identify code‐like sequence motifs for the two residue beta‐turn types. Finally, structural and sequence properties of beta‐strands in the beta‐hairpins were analyzed, which revealed that the beta‐strands showed no specific sequence and structural patterns for turn types. The analytical results in this study are expected to be a reference in the engineering or design of beta‐hairpin turn structures and sequences. Proteins 2014; 82:1721–1733. © 2014 Wiley Periodicals, Inc.     

Folding pattern of the alpha-crystallin domain in alphaA-crystallin determined by site-directed spin labeling

The folding pattern of the alpha-crystallin domain, a conserved protein module encoding the molecular determinants of structure and function in the small heat-shock protein superfamily, was determined in the context of the lens protein alphaA-crystallin by systematic application of site-directed spin labeling. The sequence-specific secondary structure was assigned primarily from nitroxide scanning experiments in which the solvent accessibility and mobility of a nitroxide probe were measured as a function of residue number. Seven beta-strands were identified and their orientation relative to the aqueous solvent determined, thus defining the residues lining the hydrophobic core. The pairwise packing of adjacent strands in the primary structure was deduced from patterns of proximities in nitroxide pairs with one member on the exposed surface of each strand. In addition to identifying supersecondary structures, these proximities revealed that the seven strands are arranged in two beta-sheets. The overall packing of the two sheets was determined by application of the general rules of protein structure and from proximities in nitroxide pairs designed to distinguish between known all beta-sheet folds. Our data are consistent with an immunoglobulin-like fold consisting of two aligned beta-sheets. Comparison of this folding pattern to that of the evolutionary distant alpha-crystallin domain in Methanococcus jannaschii heat-shock protein 16.5 reveals a conserved core structure with the differences sequestered at one edge of the beta-sandwich. A beta-strand deletion in alphaA-crystallin disrupts a subunit interface and allows for a different dimerization motif. Putative substrate binding regions appear to include a buried loop and a buried turn, suggesting that the chaperone function involves a disassembly of the oligomer.     

Tertiary structure-dependence of misfolding substitutions in loops of the maltose-binding protein.

We previously identified and characterized amino acid substitutions in a loop connecting helix I to strand B, the alphaI/betaB loop, of the N-domain that are critical for in vivo folding of the maltose-binding protein (MalE31). The tertiary context-dependence of this mutation in MalE folding was assessed by probing the tolerance of an equivalent alphabeta loop of the C-domain to the same amino acid substitutions (MalE219). Moving the loop mutation from the N- to the C-domain eliminated the in vivo misfolding step that led to the formation of inclusion bodies. In vitro, both loop variants exhibited an important decrease of stability, but their intrinsic tendency to aggregate was well correlated with their periplasmic fates in Escherichia coli. Furthermore, the noncoincidence of the unfolding and refolding transition curves and increase of light scattering during the refolding of MalE31 indicate that a competing off-pathway reaction could occurs on the folding pathway of this variant. These results strongly support the notion that the formation of super-secondary structures of the N-domain is a rate-limiting step in the folding pathway of MalE.     

Structural consequences of an amino acid deletion in the B1 domain of protein G

We describe the NMR structure of a deletion mutant of the B1 IgG-binding domain from Group G Streptococcus. The deletion occurs within the last beta-strand of the protein, where it may potentially have a deleterious effect on the stability of the protein if the protein were not able to conformationally adjust to the perturbation. In particular, the deletion changes the registry of the final three residues in the sheet, forcing a polar Thr to be buried in the interior of the protein and exposing a hydrophobic Val to solvent. The deletion could also potentially create a large cavity in the beta-sheet and force the alpha- and gamma-carboxylates of the C-terminal Glu residue into a partially buried region of the sheet. The structure of the mutant illustrates how the conformation of the protein adjusts to the deletion, thereby mitigating some of the potentially deleterious consequences. Although the elements of secondary structure are retained between the mutant and the wt domain, there are multiple small adjustments in the segments connecting secondary structure elements. In particular, a hydrogen bond between the Glu57 carboxylates and two main chain amides is introduced that alters the conformation in the loop connecting the helix to strand 3. In addition, to minimize hydrophobic surface exposure, the turn connecting strands 1 and 2 folds toward the core so that the molecular volume is decreased.     

Structural events during the refolding of an all beta-sheet protein

The refolding kinetics of the 140-residue, all beta-sheet, human fibroblast growth factor (hFGF-1) is studied using a variety of biophysical techniques such as stopped-flow fluorescence, stopped-flow circular dichroism, and quenched-flow hydrogen exchange in conjunction with multidimensional NMR spectroscopy. Urea-induced unfolding of hFGF-1 under equilibrium conditions reveals that the protein folds via a two-state (native <--> unfolded) mechanism without the accumulation of stable intermediates. However, measurement of the unfolding and refolding rates in various concentrations of urea shows that the refolding of hFGF-1 proceeds through accumulation of kinetic intermediates. Results of the quenched-flow hydrogen exchange experiments reveal that the hydrogen bonds linking the N- and C-terminal ends are the first to form during the refolding of hFGF-1. The basic beta-trefoil framework is provided by the simultaneous formation of beta-strands I, IV, IX, and X. The other beta-strands comprising the beta-barrel structure of hFGF-1 are formed relatively slowly with time constants ranging from 4 to 13 s.     

NMR structure of an archaeal homologue of ribonuclease P protein Rpp29

A protein component of the Archaeoglobus fulgidus RNase P was expressed in Escherichia coli, purified, and structurally characterized using multidimensional NMR methods. The dominant structural feature of this 11 kDa protein is a sheet of six antiparallel beta-strands, wrapped around a core of conserved hydrophobic amino acids. Amide proton exchange and (15)N relaxation rate data provide evidence that the first 16 residues of the protein, located before the start of the first beta-strand, and the last 24 residues, located past the end of the last beta-strand, are relatively flexible; this contrasts with the relatively rigid and well-defined structure of the beta-sheet. Amino acid sequence comparisons among a diverse set of species indicate that the A. fulgidus protein is homologous to the human RNase P protein Rpp29, yeast RNase P protein Pop4, and a known archaeal RNase P protein from Methanobacter thermoautotrophicus; conserved hydrophobic residues indicate that the homologous protein in each of these species contains a similar beta-sheet structure. Conserved surface residues located in the loop connecting strands beta2 and beta3, the loop connecting strands beta4 and beta5, and in the flexible N- and C-terminal tails are most likely to have specific interactions with the RNA and other proteins of RNase P. The structural model of an RNase P protein component provided by the present work provides an essential step toward eventually understanding the overall architecture of this complex enzyme and the mechanism by which it performs its functions.     

Roles of beta-turns in protein folding: from peptide models to protein engineering

Reverse turns are a major class of protein secondary structure; they represent sites of chain reversal and thus sites where the globular character of a protein is created. It has been speculated for many years that turns may nucleate the formation of structure in protein folding, as their propensity to occur will favor the approximation of their flanking regions and their general tendency to be hydrophilic will favor their disposition at the solvent-accessible surface. Reverse turns are local features, and it is therefore not surprising that their structural properties have been extensively studied using peptide models. In this article, we review research on peptide models of turns to test the hypothesis that the propensities of turns to form in short peptides will relate to the roles of corresponding sequences in protein folding. Turns with significant stability as isolated entities should actively promote the folding of a protein, and by contrast, turn sequences that merely allow the chain to adopt conformations required for chain reversal are predicted to be passive in the folding mechanism. We discuss results of protein engineering studies of the roles of turn residues in folding mechanisms. Factors that correlate with the importance of turns in folding indeed include their intrinsic stability, as well as their topological context and their participation in hydrophobic networks within the protein's structure.     

Dramatic stabilization of an SH3 domain by a single substitution: roles of the folded and unfolded states

The N-terminal SH3 domain of the Drosophila drk protein (drkN SH3) exists in equilibrium between folded and unfolded states under non-denaturing buffer conditions. In order to examine the origins of this instability, we have made mutations in the domain and characterized the thermodynamics and kinetics of folding. Results of substitutions of negatively charged residues to neutral amino acid residues suggest that the large electrostatic potential of the domain does not play a dominant role in the instability of the domain. Sequence alignment of a large number of SH3 domains reveals that the drkN SH3 domain has a threonine (T22) at a position corresponding to an otherwise highly conserved glycine residue in the diverging beta-turn connecting the beta3 and beta4 strands. Mutation of T22 to glycine results in significant stabilization of the drkN SH3 domain by 2.5 kcal/mole. To further characterize the basis for the stabilization of the T22 mutant relative to wild-type, we made additional mutant proteins with substitutions of residue T22. A strong correlation is seen between protein stability or folding rate and propensity for native beta-turn structure at this position. Correlation of folding rates with AGADIR predictions of non-native helical structure in the diverging turn region, along with our previous NMR evidence for non-native structure in this region of the unfolded state of the drkN SH3 domain, suggests that the free energy of the unfolded state also plays a role in stability. This result highlights the importance of both folded and unfolded states for understanding protein stability.     

Measuring cooperativity in the formation of a three-stranded beta sheet (double hairpin)

Recently validated chemical shift measures of hairpin structuring have been applied to a series of turn mutants of the Schenck-Gellman three-strand beta-sheet model with the aim of measuring the entropic advantage associated with aligning an additional strand onto an existing hairpin versus aligning the same two strands in an initial hairpin formation. In a four-state analysis (unfolded, 2 single hairpins, and the double hairpin fold in equilibrium) a cooperativity index can be defined as the factor by which the equilibrium constant for hairpin formation is improved when one strand is prestructured. This cooperativity index is 2.7 +/- 0.7 for hairpin formation about a stable D-Pro-Gly turn locus and increases to 7.6 +/- 1.2 for an Asn-Gly turn locus. The latter corresponds to a cooperativity induced DeltaDeltaG increment of 4.9 kJ/mol for the folding of a hairpin. Although larger than previous experimental measures of folding cooperativity in three-stranded sheets, the magnitude of this effect (which is considerably less than the TDeltaDeltaS expectation for prestructuring three or more beta-strand residue sites) likely reflects the intrinsic preference of these designed sequences for extended conformations. If similar or larger effects apply to protein beta-sheet folding, it is not surprising that particularly favorable hairpin alignments serve as nucleation sites in protein folding pathways.     

Comparison of the human and worm p53 structures suggests a way for enhancing stability

Maintaining the native conformation is essential for the proper function of tumor suppressor protein p53. However, p53 is a low-stability protein that can easily lose its function upon structural perturbations such as those resulting from missense mutations, leading to the development of cancer. Therefore, it is important to develop strategies to design stable p53 which still maintains its normal function. Here, we compare the stabilities of the human and worm p53 core domains using molecular dynamics simulations. We find that the worm p53 is significantly more stable than the human form. Detailed analysis of the structural fluctuations shows that the stability difference lies in the peripheral structural motifs that contrast in their structural features and flexibility. The most dramatic difference in stability originates from loop L1, from the turn between helix H1 and beta-strand S5, and from the turn that connects beta-strands S7 and S8. Structural analysis shows significant differences for these motifs between the two proteins. Loop L1 lacks secondary structure, and the turns between helix H1 and strand S5 and between strands S7 and S8 are much longer in the human form p53. On the basis of these differences, we designed a mutant by shortening the turn between strands S7 and S8 to enhance the stability. Surprisingly, this mutant was very stable when probed by molecular dynamics simulations. In addition, the stabilization was not localized in the turn region. Loop L1 was also significantly stabilized. Our results show that stabilizing peripheral structural motifs can greatly enhance the stability of the p53 core domain and therefore is likely to be a viable alternative in the development of stable p53. In addition, loop- or turn-related mutants with different stabilities may also be used to probe the relationship between function, a particular structural motif, and its flexibility.     

Context-dependent effects of proline residues on the stability and folding pathway of ubiquitin.

Substitution of trans-proline at three positions in ubiquitin (residues 19, 37 and 38) produces significant context-dependent effects on protein stability (both stabilizing and destabilizing) that reflect changes to a combination of parameters including backbone flexibility, hydrophobic interactions, solvent accessibility to polar groups and intrinsic backbone conformational preferences. Kinetic analysis of the wild-type yeast protein reveals a predominant fast-folding phase which conforms to an apparent two-state folding model. Temperature-dependent studies of the refolding rate reveal thermodynamic details of the nature of the transition state for folding consistent with hydrophobic collapse providing the overall driving force. Br?nsted analysis of the refolding and unfolding rates of a family of mutants with a variety of side chain substitutions for P37 and P38 reveals that the two prolines, which are located in a surface loop adjacent to the C terminus of the main alpha-helix (residues 24-33), are not significantly structured in the transition state for folding and appear to be consolidated into the native structure only late in the folding process. We draw a similar conclusion regarding position 19 in the loop connecting the N-terminal beta-hairpin to the main alpha-helix. The proline residues of ubiquitin are passive spectators in the folding process, but influence protein stability in a variety of ways.