Cavity formation before stable hydrogen bonding in the folding of a beta-clam protein

Polyethylene glycol perturbs the unfolding of CRABP I: A correlation between experimental and theoretical approach

The importance of macromolecules paves the way towards a detailed molecular level investigation as all most all cellular processes occurring at the interior of cells in the form of proteins, enzymes, and other biological molecules are significantly affected because of their crowding. Thus, exploring the role of crowding environment on the denaturation and renaturation kinetics of protein molecules is of great importance. Here, CRABP I (cellular retinoic acid binding protein I) is employed as a model protein along with different molecular weights of Polyethylene glycol (PEG) as molecular crowders. The experimental evaluations are done by accessing the protein secondary structure analysis using circular dichroism (CD) spectroscopy and unfolding kinetics using intrinsic fluorescence of CRABP I at 37 °C to mimic the in vivo crowding environment. The unfolding kinetics results indicated that both PEG 2000 and PEG 4000 act as stabilizers by retarding the unfolding kinetic rates. Both kinetic and stability outcomes presented the importance of crowding environment on stability and kinetics of CRABP I. The molecular dynamics (MD) studies revealed that thirteen PEG 2000 molecules assembled during the 500 ns simulation, which increases the stability and percentage of β-sheet. The experimental findings are well supported by the molecular dynamics simulation results.

Kinetic versus thermodynamic control of mutational effects on protein homeostasis: A perspective from computational modeling and experiment

The effect of mutations in individual proteins on protein homeostasis, or "proteostasis," can in principle depend on the mutations' effects on the thermodynamics or kinetics of folding, or both. Here, we explore this issue using a computational model of in vivo protein folding that we call FoldEcoSlim. Our model predicts that kinetic versus thermodynamic control of mutational effects on proteostasis hinges on the relationship between how fast a protein's folding reaction reaches equilibrium and a critical time scale that characterizes the lifetime of a protein in its environment: for rapidly dividing bacteria, this time scale is that of cell division; for proteins that are produced in heterologous expression systems, this time scale is the amount of time before the protein is harvested; for proteins that are synthesized in and then exported from the eukaryotic endoplasmic reticulum, this time scale is that of protein secretion, and so forth. This prediction was validated experimentally by examining the expression yields of the wild type and several destabilized mutants of a model protein, the mouse ortholog of cellular retinoic acid-binding protein 1.

Local and non-local topological information in the denatured state ensemble of a β-barrel protein

The folding of predominantly β-sheet proteins is complicated by the presence of a large number of non-local interactions in their native states, which increase the ruggedness of their folding energy landscapes. However, forming non-local contacts early in folding or even in the unfolded state can smooth the energy landscape and facilitate productive folding. We report that several sequence regions of a β-barrel protein, cellular retinoic acid-binding protein 1 (CRABP1), populate native-like secondary structure to a significant extent in the denatured state in 8 M urea. In addition, we provide evidence for both local and non-local interactions in the denatured state of CRABP1. NMR chemical shift perturbations (CSPs) under denaturing conditions upon substitution of single residues by mutation support the presence of several non-local interactions in topologically key sites, arguing that the denatured state is conformationally restricted and contains topological information for the native fold. Among the most striking non-local interactions are those between the N- and C-terminal regions, which are involved in closure of the native β-barrel. In addition, CSPs support the presence of two features in the denatured state: a major hydrophobic cluster involving residues from various parts of the sequence and a native-like interaction similar to one identified in previous studies as forming early in folding (Budyak et al., Structure 21, 476 [2013]). Taken together, our data support a model in which transient structures involving nonlocal interactions prime early folding interactions in CRABP1, determine its barrel topology, and may protect this predominantly β-sheet protein against aggregation.

Energy landscapes of functional proteins are inherently risky

Evolutionary pressure for protein function leads to unavoidable sampling of conformational states that are at risk of misfolding and aggregation. The resulting tension between functional requirements and the risk of misfolding and/or aggregation in the evolution of proteins is becoming more and more apparent. One outcome of this tension is sensitivity to mutation, in which only subtle changes in sequence that may be functionally advantageous can tip the delicate balance toward protein aggregation. Similarly, increasing the concentration of aggregation-prone species by reducing the ability to control protein levels or compromising protein folding capacity engenders increased risk of aggregation and disease. In this Perspective, we describe examples that epitomize the tension between protein functional energy landscapes and aggregation risk. Each case illustrates how the energy landscapes for the at-risk proteins are sculpted to enable them to perform their functions and how the risks of aggregation are minimized under cellular conditions using a variety of compensatory mechanisms.

Effect of charged amino acid side chain length on lateral cross-strand interactions between carboxylate-containing residues and lysine analogues in a β-hairpin

β-Sheets are one of the fundamental three-dimensional building blocks for protein structures. Oppositely charged amino acids are frequently observed directly across one another in antiparallel sheet structures, suggesting the importance of cross-strand ion pairing interactions. Despite the apparent electrostatic nature of ion pairing interactions, the charged amino acids Asp, Glu, Arg, Lys have different numbers of hydrophobic methylenes linking the charged functionality to the backbone. Accordingly, the effect of charged amino acid side chain length on cross-strand ion pairing interactions at lateral non-hydrogen bonded positions was investigated in a β-hairpin motif. The negatively charged residues with a carboxylate (Asp, Glu, Aad in increasing length) were incorporated at position 4, and the positively charged residues with an ammonium (Dap, Dab, Orn, Lys in increasing length) were incorporated at position 9. The fraction folded population and folding free energy were derived from the chemical shift deviation data. Double mutant cycle analysis was used to determine the interaction energy for the potential lateral ion pairs. Only the Asp/Glu-Dap interactions with shorter side chains and the Aad-Orn/Lys interactions with longer side chains exhibited stabilizing energetics, mostly relying on electrostatics and hydrophobics, respectively. This suggested the need for length matching of the interacting residues to stabilize the β-hairpin motif. A survey of a nonredundant protein structure database revealed that the statistical sheet pair propensity followed the trend Asp-Lys < Glu-Lys, also implying the need for length matching of the oppositely charged residues.

Early folding events protect aggregation-prone regions of a β-rich protein

Protein folding and aggregation inevitably compete with one another. This competition is even keener for proteins with frustrated landscapes, such as those rich in β structure. It is interesting that, despite their rugged energy landscapes and high β sheet content, intracellular lipid-binding proteins (iLBPs) appear to successfully avoid aggregation, as they are not implicated in aggregation diseases. In this study, we used a canonical iLBP, cellular retinoic acid-binding protein 1 (CRABP1), to understand better how folding is favored over aggregation. Analysis of folding kinetics of point mutants reveals that the folding pathway of CRABP1 involves early barrel closure. This folding mechanism protects sequences in CRABP1 that comprise cores of aggregates as identified by nuclear magnetic resonance. The amino acid conservation pattern in other iLBPs suggests that early barrel closure may be a general strategy for successful folding and minimization of aggregation. We suggest that folding mechanisms in general may incorporate steps that disfavor aggregation.

The Role of Aromatic-Aromatic Interactions in Strand-Strand Stabilization of β-Sheets

Aromatic-aromatic interactions have long been believed to play key roles in protein structure, folding, and binding functions. However, we still lack full understanding of the contributions of aromatic-aromatic interactions to protein stability and the timing of their formation during folding. Here, using an aromatic ladder in the β-barrel protein, cellular retinoic acid-binding protein 1 (CRABP1), as a case study, we find that aromatic π stacking plays a greater role in the Phe65-Phe71 cross-strand pair, while in another pair, Phe50-Phe65, hydrophobic interactions are dominant. The Phe65-Phe71 pair spans β-strands 4 and 5 in the β-barrel, which lack interstrand hydrogen bonding, and we speculate that it compensates energetically for the absence of strand-strand backbone interactions. Using perturbation analysis, we find that both aromatic-aromatic pairs form after the transition state for folding of CRABP1, thus playing a role in the final stabilization of the β-sheet rather than in its nucleation as had been earlier proposed. The aromatic interaction between strands 4 and 5 in CRABP1 is highly conserved in the intracellular lipid-binding protein (iLBP) family, and several lines of evidence combine to support a model wherein it acts to maintain barrel structure while allowing the dynamic opening that is necessary for ligand entry. Lastly, we carried out a bioinformatics analysis and found 51 examples of aromatic-aromatic interactions across non-hydrogen-bonded β-strands outside the iLBPs, arguing for the generality of the role played by this structural motif.

Early turn formation and chain collapse drive fast folding of the major cold shock protein CspA of Escherichia coli

The folding mechanism of the β-sheet protein CspA, the major cold shock protein of Escherichia coli, was previously reported to be a concerted, two-state process. We have reexamined the folding of CspA using multiple spectroscopic probes of the equilibrium transition and laser-induced temperature jump (T-jump) to achieve better time resolution of the kinetics. Equilibrium temperature-dependent Fourier transform infrared (1634 cm(-1)) and tryptophan fluorescence measurements reveal probe-dependent thermal transitions with midpoints (T(m)) of 66 ± 1 and 61 ± 1 °C, respectively. Singular-value decomposition analysis with global fitting of the temperature-dependent infrared (IR) difference spectra reveals two spectral components with distinct melting transitions with different midpoints. T-jump relaxation measurements of CspA probed by IR and fluorescence spectroscopy show probe-dependent multiexponential kinetics characteristic of non-two-state folding. The frequency-dependent IR transients all show biphasic relaxation with average time constants of 50 ± 7 and 225 ± 25 μs at a T(f) of 77 °C and almost equal amplitudes. Similar biphasic kinetics are observed using Trp fluorescence of the wild-type protein and the Y42W and T68W mutants, with comparable lifetimes. All of these observations support a model for the folding of CspA through a compact intermediate state. The transient IR and fluorescence spectra are consistent with a diffuse intermediate having β-turns and substantial β-sheet structure. The loop β3-β4 structure is likely not folded in the intermediate state, allowing substantial solvent penetration into the barrel structure.

A career pathway in protein folding: from model peptides to postreductionist protein science

This review discusses the inherent challenge of linking "reductionist" approaches to decipher the information encoded in protein sequences with burgeoning efforts to explore protein folding in native environments-"postreductionist" approaches. Because the invitation to write this article came as a result of my selection to receive the 2010 Dorothy Hodgkin Award of the Protein Society, I use examples from my own work to illustrate the evolution from the reductionist to the postreductionist perspective. I am incredibly honored to receive the Hodgkin Award, but I want to emphasize that it is the combined effort, creativity, and talent of many students, postdoctoral fellows, and collaborators over several years that has led to any accomplishments on which this selection is based. Moreover, I do not claim to have unique insight into the topics discussed here; but this writing opportunity allows me to illustrate some threads in the evolution of protein folding research with my own experiences and to point out to those embarking on careers how the twists and turns in anyone's scientific path are influenced and enriched by the scientific context of our research. The path my own career has taken thus far has been shaped by the timing of discoveries in the field of protein science; together with our contemporaries, we become part of a knowledge evolution. In my own case, this has been an epoch of great discovery in protein folding and I feel very fortunate to have participated in it.

Macromolecular crowding remodels the energy landscape of a protein by favoring a more compact unfolded state

The interior of cells is highly crowded with macromolecules, which impacts all physiological processes. To explore how macromolecular crowding may influence cellular protein folding, we interrogated the folding landscape of a model beta-rich protein, cellular retinoic acid-binding protein I (CRABP I), in the presence of an inert crowding agent (Ficoll 70). Urea titrations revealed a crowding-induced change in the water-accessible polar amide surface of its denatured state, based on an observed ca. 15% decrease in the change in unfolding free energy with respect to urea concentration (the m-value), and the effect of crowding on the equilibrium stability of CRABP I was less than our experimental error (i.e., < or = 1.2 kcal/mol). Consequently, we directly probed the effect of crowding on the denatured state of CRABP I by measuring side-chain accessibility using iodide quenching of tryptophan fluorescence and chemical modification of cysteines. We observed that the urea-denatured state is more compact under crowded conditions, and the observed extent of reduction of the m-value by crowding agent is fully consistent with the extent of reduction of the accessibility of the Trp and Cys probes, suggesting a random and nonspecific compaction of the unfolded state. The thermodynamic consequences of crowding-induced compaction are discussed. In addition, over a wide range of Ficoll concentration, crowding significantly retarded the unfolding kinetics of CRABP I without influencing the urea dependence of the unfolding rate, arguing for no appreciable change in the nature of the transition state. Our results demonstrate how macromolecular crowding may influence protein folding by effects on both the unfolded state ensemble and unfolding kinetics.

Slow formation of aggregation-resistant beta-sheet folding intermediates

Protein folding has been studied extensively for decades, yet our ability to predict how proteins reach their native state from a mechanistic perspective is still rudimentary at best, limiting our understanding of folding-related processes in vivo and our ability to manipulate proteins in vitro. Here, we investigate the in vitro refolding mechanism of a large beta-helix protein, pertactin, which has an extended, elongated shape. At 55 kDa, this single domain, all-beta-sheet protein allows detailed analysis of the formation of beta-sheet structure in larger proteins. Using a combination of fluorescence and far-UV circular dichroism spectroscopy, we show that the pertactin beta-helix refolds remarkably slowly, with multiexponential kinetics. Surprisingly, despite the slow refolding rates, large size, and beta-sheet-rich topology, pertactin refolding is reversible and not complicated by off-pathway aggregation. The slow pertactin refolding rate is not limited by proline isomerization, and 30% of secondary structure formation occurs within the rate-limiting step. Furthermore, site-specific labeling experiments indicate that the beta-helix refolds in a multistep but concerted process involving the entire protein, rather than via initial formation of the stable core substructure observed in equilibrium titrations. Hence pertactin provides a valuable system for studying the refolding properties of larger, beta-sheet-rich proteins, and raises intriguing questions regarding the prevention of aggregation during the prolonged population of partially folded, beta-sheet-rich refolding intermediates. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

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.

Transporter-to-trap conversion: a disulfide bond formation in cellular retinoic acid binding protein I mutant triggered by retinoic acid binding irreversibly locks the ligand inside the protein

Transport proteins must bind their ligands reversibly to enable release at the point of delivery, while irreversible binding is usually associated with the extreme cases of ligand sequestration. Protein conformational dynamics is an important modulator of binding kinetics, as increased flexibility in the regions adjacent to the binding site may facilitate both association and dissociation processes. Ligand entry to, and exit from, the internal binding site of the cellular retinoic acid binding protein I (CRABP I) occurs via a flexible portal region, which functions as a dynamic aperture. We designed and expressed a CRABP I mutant (A35C/T57C), in which a small-scale conformational switch caused by the ligand binding event triggers formation of a disulfide bond in the portal region, thereby arresting structural fluctuations and effectively locking the ligand inside the binding cavity. At the same time, no formation of the disulfide bond is observed in the apo form of the mutant, and most characteristics of the mutant, including protein stability, are very similar to those of the wild-type protein in the absence of retinoic acid. The mutation does not alter the kinetics of retinoic acid binding to the protein, although the disulfide formation makes the binding effectively irreversible, as suggested by the absence of retinoic acid transfer from the holo form of the mutant to lipid vesicles in the absence of a reducing agent. Taken together, these data suggest that the disulfide bond formation in the portal region arrests large-scale structural fluctuations, which are required for retinoic acid release from the protein. The unique properties of the CRABP I mutant described in this work can be used to inspire and guide a design of nanodevices for multiple tasks ranging from sequestering small-molecule toxins in both tissue and circulation to nutrient deprivation of pathogens.

Apo-nitrophorin 4 at atomic resolution

The nitrophorins from Rhodnius prolixus, the kissing bug, are heme-containing proteins used for the transport of nitric oxide to aide the insect in obtaining a blood meal. The Rhodnius nitrophorins display an eight-stranded antiparallel beta-barrel motif, typical of lipocalins, with a histidine-linked heme in the open end of the barrel. Heme is stabilized in the ferric state and highly distorted, displaying a ruffled conformation that may be of importance in the setting of the reduction potential. To help in understanding the means by which the protein matrix, an inherently soft material, is able to distort the heme from its low-energy planar conformation, we have determined the crystal structure of apo-nitrophorin 4-1.1 A resolution. Removal of the heme from nitrophorin 4 has very little effect on its structure: The heme binding cavity remains open and the loops near the cavity entrance respond to lower pH in the same manner as the intact protein. We conclude that the general stability of the lipocalin fold and apparent rigidity of the beta-barrel provide the means for distorting the heme cofactor.

From the test tube to the cell: exploring the folding and aggregation of a beta-clam protein

A crucial challenge in present biomedical research is the elucidation of how fundamental processes like protein folding and aggregation occur in the complex environment of the cell. Many new physico-chemical factors like crowding and confinement must be considered, and immense technical hurdles must be overcome in order to explore these processes in vivo. Understanding protein misfolding and aggregation diseases and developing therapeutic strategies to these diseases demand that we gain mechanistic insight into behaviors and misbehaviors of proteins as they fold in vivo. We have developed a fluorescence approach using FlAsH labeling to study the thermodynamics of folding of a model beta-rich protein, cellular retinoic acid binding protein (CRABP) in Escherichia coli cells. The labeling approach has also enabled us to follow aggregation of a modified version of CRABP and chimeras between CRABP and huntingtin exon 1 with its glutamine repeat tract. In this article, we review our recent results using FlAsH labeling to study in-vivo folding and present new observations that hint at fundamental differences between the thermodynamics and kinetics of protein folding in vivo and in vitro.

Evolutionary coupling of structural and functional sequence information in the intracellular lipid-binding protein family

We have mined the evolutionary record for the large family of intracellular lipid-binding proteins (iLBPs) by calculating the statistical coupling of residue variations in a multiple sequence alignment using methods developed by Ranganathan and coworkers (Lockless and Ranganathan, Science 1999:286;295-299). The 213 sequences analyzed have a wide range of ligand-binding functions as well as highly divergent phylogenetic origins, assuring broad sampling of sequence space. Emerging from this analysis were two major clusters of coupled residues, which when mapped onto the structure of a representative iLBP under study in our laboratory, cellular retinoic-acid binding protein I, are largely contiguous and provide useful points of comparison to available data for the folding of this protein. One cluster comprises a predominantly hydrophobic core away from the ligand-binding site and likely represents key structural information for the iLBP fold. The other cluster includes the portal region where ligand enters its binding site, regions of the ligand-binding cavity, and the region where the 10-stranded beta-barrel characteristic of this family closes (between strands 1' and 10). Linkages between these two clusters suggest that evolutionary pressures on this family constrain structural and functional sequence information in an interdependent fashion. The necessity of the structure to wrap around a hydrophobic ligand confounds the typical sequestration of hydrophobic side chains. Additionally, ligand entry and exit require these structures to have a capacity for specific conformational change during binding and release. We conclude that an essential and structurally apparent separation of local and global sequence information is conserved throughout the iLBP family.

Fluorine-19 NMR studies on the acid state of the intestinal fatty acid binding protein

The intestinal fatty acid binding protein (IFABP) is composed of two beta-sheets with a large hydrophobic cavity into which ligands bind. After eight 4-(19)F-phenylalanines were incorporated into the protein, the acid state of both apo- and holo-IFABP (at pH 2.8 and 2.3) was characterized by means of (1)H NMR diffusion measurements, circular dichroism, and (19)F NMR. Diffusion measurements show a moderately increased hydrodynamic radius while near- and far-UV CD measurements suggest that the acid state has substantial secondary structure as well as persistent tertiary interactions. At pH 2.8, these tertiary interactions have been further characterized by (19)F NMR and show an NOE cross-peak between residues that are located on different beta-strands. Side chain conformational heterogeneity on the millisecond time scale was captured by phase-sensitive (19)F-(19)F NOESY. At pH 2.3, native NMR peaks are mostly gone, but the protein can still bind fatty acid to form the holoprotein. An exchange cross-peak of one phenylalanine in the holoprotein is attributed to increased motional freedom of the fatty acid backbone caused by the slight opening of the binding pocket at pH 2.8. In the acid environment Phe128 and Phe17 show dramatic line broadening and chemical shift changes, reflecting greater degrees of motion around these residues. We propose that there is a separation of specific regions of the protein that gives rise to the larger radius of hydration. Temperature and urea unfolding studies indicate that persistent hydrophobic clusters are nativelike and may account for the ability of ligand to bind and induce nativelike structure, even at pH 2.3.

Extended polyglutamine tracts cause aggregation and structural perturbation of an adjacent beta barrel protein

Formation of fibrillar intranuclear inclusions and related neuropathologies of the CAG-repeat disorders are linked to the expansion of a polyglutamine tract. Despite considerable effort, the etiology of these devastating diseases remains unclear. Although polypeptides with glutamine tracts recapitulate many of the observed characteristics of the gene products with CAG repeats, such as in vitro and in vivo aggregation and toxicity in model organisms, extended polyglutamine segments have also been reported to structurally perturb proteins into which they are inserted. Additionally, the sequence context of a polyglutamine tract has recently been shown to modulate its propensity to aggregate. These findings raise the possibility that indirect influences of the repeat tract on adjacent protein domains are contributory to pathologies. Destabilization of an adjacent domain may lead to loss of function, as well as favoring non-native structures in the neighboring domain causing them to be prone to intermolecular association and consequent aggregation. To explore these phenomena, we have used chimeras of a well studied globular protein and exon 1 of huntingtin. We find that expansion of the polyglutamine segment beyond the pathological threshold (>35 glutamines) results in structural perturbation of the neighboring protein whether the huntingtin exon is N- or C-terminal. Elongation of the polyglutamine region also substantially increases the propensity of the chimera to aggregate, both in vitro and in vivo, and in vitro aggregation kinetics of a chimera with a 53-glutamine repeat follow a nucleation polymerization mechanism with a monomeric nucleus.

Water and urea interactions with the native and unfolded forms of a beta-barrel protein

A fundamental understanding of protein stability and the mechanism of denaturant action must ultimately rest on detailed knowledge about the structure, solvation, and energetics of the denatured state. Here, we use (17)O and (2)H magnetic relaxation dispersion (MRD) to study urea-induced denaturation of intestinal fatty acid-binding protein (I-FABP). MRD is among the few methods that can provide molecular-level information about protein solvation in native as well as denatured states, and it is used here to simultaneously monitor the interactions of urea and water with the unfolding protein. Whereas CD shows an apparently two-state transition, MRD reveals a more complex process involving at least two intermediates. At least one water molecule binds persistently (with residence time >10 nsec) to the protein even in 7.5 M urea, where the large internal binding cavity is disrupted and CD indicates a fully denatured protein. This may be the water molecule buried near the small hydrophobic folding core at the D-E turn in the native protein. The MRD data also provide insights about transient (residence time <1 nsec) interactions of urea and water with the native and denatured protein. In the denatured state, both water and urea rotation is much more retarded than for a fully solvated polypeptide. The MRD results support a picture of the denatured state where solvent penetrates relatively compact clusters of polypeptide segments.

Hierarchical Formation of Disulfide Bonds in the Immunoglobulin Fc Fragment Is Assisted by Protein-disulfide Isomerase

Antibodies provide an excellent system to study the folding and assembly of all beta-sheet proteins and to elucidate the hierarchy of intra/inter chain disulfide bonds formation during the folding process of multimeric and multidomain proteins. Here, the folding process of the Fc fragment of the heavy chain of the antibody MAK33 was investigated. The Fc fragment consists of the C(H)3 and C(H)2 domains of the immunoglobulin heavy chain, both containing a single S-S bond. The folding process was investigated both in the absence and presence of the folding catalyst protein-disulfide isomerase (PDI), monitoring the evolution of intermediates by electrospray mass spectrometry. Moreover, the disulfide bonds present at different times in the folding mixture were identified by mass mapping to determine the hierarchy of disulfide bond formation. The analysis of the uncatalyzed folding showed that the species containing one intramolecular disulfide predominated throughout the entire process, whereas the fully oxidized Fc fragment never accumulated in significant amounts. This result suggests the presence of a kinetic trap during the Fc folding, preventing the one-disulfide-containing species (1S2H) to reach the fully oxidized protein (2S). The assignment of disulfide bonds revealed that 1S2H is a homogeneous species characterized by the presence of a single disulfide bond (Cys-130-Cys-188) belonging to the C(H)3 domain. When the folding experiments were carried out in the presence of PDI, the completely oxidized species accumulated and predominated at later stages of the process. This species contained the two native S-S bonds of the Fc protein. Our results indicate that the two domains of the Fc fragment fold independently, with a precise hierarchy of disulfide formation in which the disulfide bond, especially, of the C(H)2 domain requires catalysis by PDI.

Sequence and structural analysis of cellular retinoic acid-binding proteins reveals a network of conserved hydrophobic interactions

Proteins in the intracellular lipid-binding protein (iLBP) family show remarkably high structural conservation despite their low-sequence identity. A multiple-sequence alignment using 52 sequences of iLBP family members revealed 15 fully conserved positions, with a disproportionately high number of these (n=7) located in the relatively small helical region. The conserved positions displayed high structural conservation based on comparisons of known iLBP crystal structures. It is striking that the beta-sheet domain had few conserved positions, despite its high structural conservation. This observation prompted us to analyze pair-wise interactions within the beta-sheet region to ask whether structural information was encoded in interacting amino acid pairs. We conducted this analysis on the iLBP family member, cellular retinoic acid-binding protein I (CRABP I), whose folding mechanism is under study in our laboratory. Indeed, an analysis based on a simple classification of hydrophobic and polar amino acids revealed a network of conserved interactions in CRABP I that cluster spatially, suggesting a possible nucleation site for folding. Significantly, a small number of residues participated in multiple conserved interactions, suggesting a key role for these sites in the structure and folding of CRABP I. The results presented here correlate well with available experimental evidence on folding of CRABPs and their family members and suggest future experiments. The analysis also shows the usefulness of considering pair-wise conservation based on a simple classification of amino acids, in analyzing sequences and structures to find common core regions among homologues.

Monitoring protein stability and aggregation in vivo by real-time fluorescent labeling

In vivo fluorescent labeling of an expressed protein has enabled the observation of its stability and aggregation directly in bacterial cells. Mammalian cellular retinoic acid-binding protein I (CRABP I) was mutated to incorporate in a surface-exposed omega loop the sequence Cys-Cys-Gly-Pro-Cys-Cys, which binds specifically to a biarsenical fluorescein dye (FlAsH). Unfolding of labeled tetra-Cys CRABP I is accompanied by enhancement of FlAsH fluorescence, which made it possible to determine the free energy of unfolding of this protein by urea titration in cells and to follow in real time the formation of inclusion bodies by a slow-folding, aggregationprone mutant (FlAsH-labeled P39A tetra-Cys CRABP I). Aggregation in vivo displayed a concentration-dependent apparent lag time similar to observations of protein aggregation in purified in vitro model systems.

Native structural propensity in cellular retinoic acid-binding protein I 64-88: the role of locally encoded structure in the folding of a beta-barrel protein

A central question in protein folding is the relative importance of locally encoded structure and cooperative interactions among residues distant in sequence. We have been exploring this question in a predominantly beta-sheet protein, since beta-structure formation clearly relies on both local and global sequence information. We present evidence that a 24-residue peptide corresponding to two linked hairpins of cellular retinoic acid-binding protein I (CRABP I) adopts significant native structure in aqueous solution. Prior work from our laboratory showed that the two turns contained in this fragment (turns III and IV) had the highest tendency of any of the eight turns in this anti-parallel beta-barrel to fold into native turns. In addition, the primary sequence of these two turns is well conserved throughout the structural family to which CRABP I belongs, and residues in the turns and their associated hairpins participate in a network of conserved long-range interactions. We propose that the strong local-sequence biases within the chain segment comprising turns III and IV favor longer-range interactions that are crucial to the folding and native-state stability of CRABP I, and may play a similar role in related intracellular lipid-binding proteins (iLBPs).

Identification of rare partially unfolded states in equilibrium with the native conformation in an all beta-barrel protein

Human acidic fibroblast growth factor 1 (hFGF-1) is an all beta-barrel protein, and the secondary structural elements in the protein include 12 antiparallel beta-strands arranged into a beta-trefoil fold. In the present study, we investigate the stability of hFGF-1 by hydrogen-deuterium exchange as a function of urea concentration. Urea-induced equilibrium unfolding of hFGF-1 monitored by fluorescence and CD spectroscopy suggests that the protein unfolds by a two-state (native to denatured) mechanism. Hydrogen exchange in hFGF-1, under the experimental conditions used, occurs by the EX2 mechanism. In contrast to the equilibrium unfolding events monitored by optical probes, native state hydrogen exchange data show that the beta-trefoil architecture of hFGF-1 does not behave as a single cooperative unit. There are at least two structurally independent units with differing stabilities in hFGF-1. Beta-strands I, II, III, VI, VII, X, XI, and XII fit into the global unfolding isotherm. By contrast, residues in beta-strands IV, V, VIII, and IX exchange by the subfolding isotherm and could be responsible for the occurrence of high-energy partially unfolded state(s) in hFGF-1. There appears to be a broad continuum of stabilities among the four beta-strands (beta-strands IV, V, VIII, and IX) constituting the subglobal folding unit. The slow exchanging residues in hFGF-1 do not represent the folding nucleus of the protein.

Conserved and nonconserved features of the folding pathway of hisactophilin, a beta-trefoil protein

Based on previous studies of interleukin-1beta (IL-1beta) and both acidic and basic fibroblast growth factors (FGFs), it has been suggested that the folding of beta-trefoil proteins is intrinsically slow and may occur via the formation of essential intermediates. Using optical and NMR-detected quenched-flow hydrogen/deuterium exchange methods, we have measured the folding kinetics of hisactophilin, another beta-trefoil protein that has < 10% sequence identity and unrelated function to IL-1beta and FGFs. We find that hisactophilin can fold rapidly and with apparently two-state kinetics, except under the most stabilizing conditions investigated where there is evidence for formation of a folding intermediate. The hisactophilin intermediate has significant structural similarities to the IL-1beta intermediate that has been observed experimentally and predicted theoretically using a simple, topology-based folding model; however, it appears to be different from the folding intermediate observed experimentally for acidic FGF. For hisactophilin and acidic FGF, intermediates are much less prominent during folding than for IL-1beta. Considering the structures of the different beta-trefoil proteins, it appears that differences in nonconserved loops and hydrophobic interactions may play an important role in differential stabilization of the intermediates for these proteins.

Studies of biomolecular conformations and conformational dynamics by mass spectrometry

In the post-genomic era, a wealth of structural information has been amassed for proteins from NMR and crystallography. However, static protein structures alone are not a sufficient description: knowledge of the dynamic nature of proteins is essential to understand their wide range of functions and behavior during the life cycle from synthesis to degradation. Furthermore, few proteins have the ability to act alone in the crowded cellular environment. Assemblies of multiple proteins governed by complex signaling pathways are often required for the tasks of target recognition, binding, transport, and function. Mass spectrometry has emerged over the past several years as a powerful tool to address many of these questions. Recent improvements in "soft" ionization techniques have enabled researchers to study proteins and biomolecular complexes, both directly and indirectly. Likewise, continuous improvements in instrumental design in recent years have resulted in a dramatic expansion of the m/z range and resolution, enabling observation of large multi-protein assemblies whose structures are retained in the gas phase. In this article, we discuss some of the mass spectrometric techniques applied to investigate the nature of the conformations and dynamical properties that govern protein function.

Examination of the folding of E. coli CspA through tryptophan substitutions

Escherichia coli cold shock protein, CspA, folds very rapidly (time constant, tau = 4 msec) by an apparent two-state mechanism. However, recent time-resolved infrared (IR) temperature-jump experiments indicate that the folding trajectory of CspA may be more complicated. The sole tryptophan of wild-type CspA (Trp11), which is used to monitor the folding process by fluorescence spectroscopy, is located in an unusual aromatic cluster on the surface of CspA within the nucleic acid binding site. To gain a more global picture of the folding kinetics of CspA and to determine if there are any previously undetected intermediates, we have introduced a second tryptophan at three different surface locations in the protein. The three mutations did not significantly alter the tertiary structure of CspA, although two of the substitutions were found to be slightly stabilizing. Two-state folding, as detected by stopped-flow fluorescence spectroscopy, is preserved in all three mutants. These results indicate that the fast folding of CspA is driven by a concerted mechanism.

Folding of intracellular retinol and retinoic acid binding proteins

The folding mechanisms of cellular retinol binding protein II (CRBP II), cellular retinoic acid binding protein I (CRABP I), and cellular retinoic acid binding protein II (CRABP II) were examined. These beta-sheet proteins have very similar structures and higher sequence homologies than most proteins in this diverse family. They have similar stabilities and show completely reversible folding at equilibrium with urea as a denaturant. The unfolding kinetics of these proteins were monitored during folding and unfolding by circular dichroism (CD) and fluorescence. During unfolding, CRABP II showed no intermediates, CRABP I had an intermediate with nativelike secondary structure, and CRBP II had an intermediate that lacked secondary structure. The refolding kinetics of these proteins were more similar. Each protein showed a burst-phase change in intensity by both CD and fluorescence, followed by a single observed phase by both CD and fluorescence and one or two additional refolding phases by fluorescence. The fluorescence spectral properties of the intermediate states were similar and suggested a gradual increase in the amount of native tertiary structure present for each step in a sequential path. However, the rates of folding differed by as much as 3 orders of magnitude and were slower than those expected from the contact order and topology of these proteins. As such, proteins with the same final structure may not follow the same route to the native state.

Denaturation and partial renaturation of a tightly tetramerized DsRed protein under mildly acidic conditions

The red fluorescent protein, DsRed, recently cloned from coral Discosoma sp. has one of the longest fluorescence waves and one of the most complex absorbance spectra among the family of fluorescent proteins. In this work we found that with time DsRed fluorescence decreases under mildly acidic conditions (pH 4.0-4.8) in a pH-dependent manner, and this fluorescence inactivation could be partially recovered by subsequent re-alkalization. The DsRed absorbance and circular dichroism spectra under these conditions revealed that the fluorescence changes were caused by denaturation followed by partial renaturation of the protein. Further, analytical ultracentrifugation determined that native DsRed formed a tight tetramer under various native conditions. Quantitative analysis of the data showed that several distinct states of protein exist during the fluorescence inactivation and recovery, and the inactivation of fluorescence can be caused by protonation of a single ionogenic group in each monomer of DsRed tetramer.

Multiple roles of prolyl residues in structure and folding

To explore the ways that proline residues may influence the conformational options of a polypeptide backbone, we have characterized Pro-->Ala mutants of cellular retinoic acid-binding protein I (CRABP I). While all three Xaa-Pro bonds are in the trans conformation in the native protein and the equilibrium stability of each mutant is similar to that of the parent protein, each has distinct effects on folding and unfolding kinetics. The mutation of Pro105 does not alter the kinetics of folding of CRABP I, which indicates that the flexible loop containing this residue is passive in the folding process. By contrast, replacement of Pro85 by Ala abolishes the observable slow phase of folding, revealing that correct configuration of the 84-85 peptide bond is prerequisite to productive folding. Substitution of Pro39 by Ala yields a protein that folds and unfolds more slowly. Removal of the conformational constraint imposed by the proline ring likely raises the transition state barrier by increasing the entropic cost of narrowing the conformational ensemble. Additionally, the Pro-->Ala mutation removes a helix-termination signal that is important for efficient folding to the native state.

The folding of an immunoglobulin-like Greek key protein is defined by a common-core nucleus and regions constrained by topology

TNfn3, the third fibronectin type III domain of human tenascin, is an immunoglobulin-like protein that is a good model for experimental and theoretical analyses of Greek key folding. The third fibronectin type III domain of human tenascin folds and unfolds in a two-state fashion over a range of temperature and pH values, and in the presence of stabilising salts. Here, we present a high resolution protein engineering analysis of the single rate determining transition state. The 48 mutations report on the contribution of side-chains at 32 sites in the core and loop regions. Three areas in the protein exhibit high Phi-values, indicating that they are partially structured in the transition state. First, a common-core ring of four positions in the central strands B, C, E and F, that are in close contact, form a nucleus of tertiary interactions. The two other regions that appear well-formed are the C' region and the E-F loop. The Phi-values gradually decrease away from these regions such that the very ends of the two terminal strands A and G, have Phi-values of zero. We propose a model for the folding of immunoglobulin-like proteins in which the common-core "ring" forms the nucleus for folding, whilst the C' and E-F regions are constrained by topology to pack early. Folding characteristics of a group of structurally related proteins appear to support this model.

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.

Unfolding dynamics of a beta-sheet protein studied by mass spectrometry

The unfolding dynamics of cellular retinoic acid-binding protein I (CRABP I), an 18 kDa predominantly beta-sheet protein, were studied by monitoring the hydrogen-deuterium (H-D) exchange reaction under various solution conditions. A bimodal charge state distribution was observed when a denaturing agent was added to the protein aqueous solution. These two populations exhibit different kinetics of H-D exchange, with the high charge state ions undergoing very rapid isotope exchange, while the low charge state protein ions exchange cooperatively but at much slower rates. Transiently populated intermediate states were detected indirectly using hydrogen exchange measurement in aqueous solution at various pHs. At pH 2.5 and room temperature, three distinct populations of CRABP I ions exist over an extended period of time, each corresponding to a specific degree of backbone amide hydrogen atom protection. Mass spectral data are complementary to hydrogen exchange measurements by NMR, since the former samples a much faster time-scale of dynamic events in solution.

Folding and association of the antibody domain CH3: prolyl isomerization preceeds dimerization

The simplest naturally occurring model system for studying immunoglobulin folding and assembly is the non-covalent homodimer formed by the C-terminal domains (CH3) of the heavy chains of IgG. Here, we describe the structure of recombinant CH3 dimer as determined by X-ray crystallography and an analysis of the folding pathway of this protein. Under conditions where prolyl isomerization does not contribute to the folding kinetics, formation of the beta-sandwich structure is the rate-limiting step. beta-Sheet formation of CH3 is a slow process, even compared to other antibody domains, while the subsequent association of the folded monomers is fast. After long-time denaturation, the majority of the unfolded CH3 molecules reaches the native state in two serial reactions, involving the re-isomerization of the Pro35-peptide bond to the cis configuration. The species with the wrong isomer accumulate as a monomeric intermediate. Importantly, the isomerization to the correct cis configuration is the prerequisite for dimerization of the CH3 domain. In contrast, in the Fab fragment of the same antibody, prolyl isomerization occurs after dimerization demonstrating that within one protein, comprised of highly homologous domains, both the kinetics of beta-sandwich formation and the stage at which prolyl isomerization occurs during the folding process can be completely different.

Context-dependence of amino acid residue pairing in antiparallel beta-sheets

In an effort to understand the driving forces behind antiparallel beta-sheet assembly, we have investigated the mutational tolerance of four pairs of residues in CspA, the major cold shock protein of E. coli. Two buried pairs and two exposed pairs of neighboring amino acids were separately randomized and the corresponding effects on protein stability were assessed using a protein expression screen. The thermal denaturation of a subset of the recovered proteins was measured by circular dichroism spectroscopy in order to determine the range of stabilities sampled by the expressed mutants. As anticipated, buried sites are substantially less tolerant of substitutions than exposed sites with more than half of the exposed residue combinations giving rise to stably folded proteins. The two exposed residue pairs, however, display different degrees of tolerance to substitution and accept different residue pair combinations. Except for the prohibition of proline from interior strand positions, no obvious correlations of mutant stability with any single parameter such as beta-sheet propensity or hydrophobicity can be detected. Mutant combinations recovered in both orientations (e.g. XY and YX) at a given exposed pair site often show markedly different stabilities, indicating that the local environment plays a substantial role in modulating the pairing preferences of residues in beta-sheets.

Biophysics of the green fluorescent protein

It is almost certainly a truism that interpretation of the fluorescence of a protein matrix-embedded chromophore in terms of the physicochemical character of its environment requires that the tertiary structure of the protein be known to high resolution. This reality derives from the complexity of the photophysics of most fluorescent molecules, complexity that reveals the imperfections of available theory. The accuracy of these dicta is highlighted by the biophysical properties of the green fluorescent protein now being so elegantly elucidated from the application of X-ray crystallography, ultrafast optical spectroscopy, and site-specific mutagenesis. Despite the mass of recent data, however, the physicochemical basis of the green fluorescence cannot be regarded as having been fully defined, nor has the role of protein folding in chromophore formation been solved. In addition, GFP consistently yields surprises typified by the recent experiments of Vanden Bout et al. (1997) on the fluorescence of single molecules mutant (T203F and T2034) GFPs, which showed unique reversible photobleaching and were whimsically termed "blinking molecules" by Moerner (1977). Given the apparent malleability of the GFP sequence and the sensitivity of the chromophore's photophysics to a broad spectrum of physicochemical factors, it is inevitable that additional useful and intriguing biophysical properties will emerge from the study of other mutants. Although on the surface it may seem mundane, determination of the amino acid sequence and tertiary structures of the GFPs from other coelenterates is quite likely to provide very useful insights into the biophysical bases of both protein folding and the green fluorescence per se. Finally, a broader set of spectroscopic techniques need to be applied to the study of GFPs, and future fluorescence examination should include measurements of transient absorption and fluorescence emission anisotropy decays.

Probing the folding pathway of a beta-clam protein with single-tryptophan constructs

BACKGROUND

Cellular retinoic acid binding protein I (CRABPI) is a small, predominantly beta-sheet protein with a simple architecture and no disulfides or cofactors. Folding of mutants containing only one of the three native tryptophans has been examined using stopped-flow fluorescence and circular dichroism at multiple wavelengths.

RESULTS

Within 10 ms, the tryptophan fluorescence of all three mutants shows a blue shift, and stopped-flow circular dichroism shows significant secondary structure content. The local environment of Trp7, a completely buried residue located near the intersection of the N and C termini, develops on a 100 ms time scale. Spectral signatures of the other two tryptophan residues (87 and 109) become native-like in a 1 s kinetic phase.

CONCLUSIONS

Formation of the native beta structure of CRABPI is initiated by rapid hydrophobic collapse, during which local segments of chain adopt significant secondary structure. Subsequently, transient yet specific interactions of amino acid residues restrict the arrangement of the chain topology and initiate long-range associations such as the docking of the N and C termini. The development of native tertiary environments, including the specific packing of the beta-sheet sidechains, occurs in a final, highly cooperative step simultaneous with stable interstrand hydrogen bonding.

Folding mechanism of three structurally similar beta-sheet proteins

The folding mechanism of cellular retinoic acid binding protein I (CRABP I), cellular retinol binding protein II (CRBP II), and intestinal fatty acid binding protein (IFABP) were investigated to determine if proteins with similar native structures have similar folding mechanisms. These mostly beta-sheet proteins have very similar structures, despite having as little as 33% sequence similarity. The reversible urea denaturation of these proteins was characterized at equilibrium by circular dichroism and fluorescence. The data were best fit by a two-state model for each of these proteins, suggesting that no significant population of folding intermediates were present at equilibrium. The native states were of similar stability with free energies (linearly extrapolated to 0 M urea, deltaGH2O) of 6.5, 8.3, and 5.5 kcal/mole for CRABP I, CRBP II, and IFABP, respectively. The kinetics of the folding and unfolding processes for these proteins was monitored by stopped-flow CD and fluorescence. Intermediates were observed during both the folding and unfolding of all of these proteins. However, the overall rates of folding and unfolding differed by nearly three orders of magnitude. Further, the spectroscopic properties of the intermediate states were different for each protein, suggesting that different amounts of secondary and/or tertiary structure were associated with each intermediate state for each protein. These data show that the folding path for proteins in the same structural family can be quite different, and provide evidence for different folding landscapes for these sequences.

Turn scanning by site-directed mutagenesis: application to the protein folding problem using the intestinal fatty acid binding protein

We have systematically mutated residues located in turns between beta-strands of the intestinal fatty acid binding protein (IFABP), and a glycine in a half turn, to valine and have examined the stability, refolding rate constants and ligand dissociation constants for each mutant protein. IFABP is an almost all beta-sheet protein exhibiting a topology comprised of two five-stranded sheets surrounding a large cavity into which the fatty acid ligand binds. A glycine residue is located in seven of the eight turns between the antiparallel beta-strands and another in a half turn of a strand connecting the front and back sheets. Mutations in any of the three turns connecting the last four C-terminal strands slow the folding and decrease stability with the mutation between the last two strands slowing folding dramatically. These data suggest that interactions between the last four C-terminal strands are highly cooperative, perhaps triggered by an initial hydrophobic collapse. We suggest that this trigger is collapse of the highly hydrophobic cluster of amino acids in the D and E strands, a region previously shown to also affect the last stage of the folding process (Kim et al., 1997). Changing the glycine in the strand between the front and back sheets also results in a unstable, slow folding protein perhaps disrupting the D-E strand interactions. For most of the other turn mutations there was no apparent correlation between stability and refolding rate constants. In some turns, the interaction between strands, rather than the turn type, appears to be critical for folding while in others, turn formation itself appears to be a rate limiting step. Although there is no simple correlation between turn formation and folding kinetics, we propose that turn scanning by mutagenesis will be a useful tool for issues related to protein folding.

Kinetic studies of beta-sheet protein folding

New studies have shown that folding of beta-sheet proteins can occur with and without intermediates, with fast to slow refolding rates and late to very late transition states. These experiments demonstrate that, despite early speculation to the contrary, beta-sheet protein folding does not appear to be fundamentally different from that of helical and mixed alpha, beta proteins.