Role of a solvent-exposed aromatic cluster in the folding of Escherichia coli CspA

Mutational Analysis of Protein Folding Transition States: Phi Values

The analysis of protein folding reactions by monitoring the kinetic effects of specifically designed single-point mutations, the so-termed phi-value analysis, has been a favorite technique to experimentally probe the mechanisms of protein folding. The idea behind phi-value analysis is that the effects that mutations have on the folding and unfolding rate constants report on the energetic/structural features of the folding transition state ensemble (TSE), which is the highest point in the free energy surface connecting the native and unfolded states, and thus the rate limiting step that ultimately defines the folding mechanism. For single-domain, two-state folding proteins, the general procedure to perform the phi-value analysis of protein folding is relatively simple to implement in the lab. Once the mutations have been produced and purified, the researcher needs to follow a few specific guidelines to perform the experiments and to analyze the data so produced. In this chapter, a step-by-step description of how to measure and interpret the effects induced by site-directed mutations on the folding and unfolding rate constants of a protein of interest is provided. Some possible solutions to the most typical problems that arise when performing phi-value analysis in the lab are also provided.

Y-Box Binding Proteins in mRNP Assembly, Translation, and Stability Control

Y-box binding proteins (YB proteins) are DNA/RNA-binding proteins belonging to a large family of proteins with the cold shock domain. Functionally, these proteins are known to be the most diverse, although the literature hardly offers any molecular mechanisms governing their activities in the cell, tissue, or the whole organism. This review describes the involvement of YB proteins in RNA-dependent processes, such as mRNA packaging into mRNPs, mRNA translation, and mRNA stabilization. In addition, recent data on the structural peculiarities of YB proteins underlying their interactions with nucleic acids are discussed.

A conserved and buried edge-to-face aromatic interaction in small ubiquitin-like modifier (SUMO) has a role in SUMO stability and function

Aromatic amino acids buried at a protein's core are often involved in mutual paired interactions. Ab initio energy calculations have highlighted that the conformational orientations and the effects of substitutions are important for stable aromatic interactions among aromatic rings, but studies in the context of a protein's fold and function are elusive. Small ubiquitin-like modifier (SUMO) is a common post-translational modifier that affects diverse cellular processes. Here, we report that a highly conserved aromatic triad of three amino acids, Phe36-Tyr51-Phe64, is a unique SUMO signature that is absent in other ubiquitin-like homologous folds. We found that a specific edge-to-face conformation between the Tyr51-Phe64 pair of interacting aromatics is vital to the fold and stability of SUMO. Moreover, the noncovalent binding of SUMO-interacting motif (SIM) at the SUMO surface was critically dependent on the paired aromatic interactions buried at the core. NMR structural studies revealed that perturbation of the Tyr51-Phe64 conformation disrupts several long-range tertiary contacts in SUMO, leading to a heterogeneous and dynamic protein with attenuated SUMOylation both in vitro and in cells. A subtle perturbation of the edge-to-face conformation by a Tyr to Phe substitution significantly decreased stability, SUMO/SIM affinity, and the rate of SUMOylation. Our results highlight that absolute co-conservation of specific aromatic pairs inside the SUMO protein core has a role in its stability and function.

Spatial ranges of driving forces are a key determinant of protein folding cooperativity and rate diversity

The physical basis of two-state-like folding transitions and the tremendous diversity in folding rates is elucidated by directly simulating the folding kinetics of 52 representative proteins. Relative to the results from a common modeling approach, the diversity of the simulated folding rates can be increased from ~10(2.1) to the experimental ~10(6.0) by a modest decrease in the spatial range of the attractive potential. The required theoretical range is consistent with desolvation physics and is notably much more permissive than that needed for two-state-like homopolymer collapse.

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.

Is there an en route folding intermediate for Cold shock proteins?

Cold shock proteins (Csps) play an important role in cold shock response of a diverse number of organisms ranging from bacteria to humans. Numerous studies of the Csp from various species showed that a two-state folding mechanism is conserved and the transition state (TS) appears to be very compact. However, the atomic details of the folding mechanism of Csp remain unclear. This study presents the folding mechanism of Csp in atomic detail using an all-atom Go model-based simulations. Our simulations predict that there may exist an en route intermediate, in which β strands 1-2-3 are well ordered and the contacts between β1 and β4 are almost developed. Such an intermediate might be too unstable to be detected in the previous fluorescence energy transfer experiments. The transition state ensemble has been determined from the P(fold) analysis and the TS appears even more compact than the intermediate state.

NMR analysis of a kinetically trapped intermediate of a disulfide-deficient mutant of the starch-binding domain of glucoamylase

A thermally unfolded disulfide-deficient mutant of the starch-binding domain of glucoamylase refolds into a kinetically trapped metastable intermediate when subjected to a rapid lowering of temperature. We attempted to characterise this intermediate using multidimensional NMR spectroscopy. The (1)H-(15)N heteronuclear single quantum coherence spectrum after a rapid temperature decrease (the spectrum of the intermediate) showed good chemical shift dispersion but was significantly different from that of the native state, suggesting that the intermediate adopts a nonnative but well-structured conformation. Large chemical shift changes for the backbone amide protons between the native and the intermediate states were observed for residues in the β-sheet consisting of strands 2, 3, 5, 6, and 7 as well as in the C-terminal region. These residues were found to be in close proximity to aromatic residues, suggesting that the chemical shift changes are mainly due to ring current shifts caused by the aromatic residues. The two-dimensional nuclear Overhauser enhancement (NOE) spectroscopy experiments showed that the intermediate contained substantial, native-like NOE connectivities, although there were fewer cross peaks in the spectrum of the intermediate compared with that of the native state. It was also shown that there were native-like interresidue NOEs for residues buried in the protein, whereas many of the NOE cross peaks were lost for the residues involved in a surface-exposed aromatic cluster. These results suggest that, in the intermediate, the aromatic cluster at the surface is structurally less organised, whereas the interior of the protein has relatively rigid, native-like side-chain packing.

The effect of context on the folding of β-hairpins

Small β-hairpin peptides have been widely used as models for the folding of β-sheets. But how applicable is the folding of such models to β-structure in larger proteins with conventional hydrophobic cores? Here we present multiple unfolding simulations of three such proteins that contain the WW domain double hairpin β-sheet motif: cold shock protein A (CspA), cold shock protein B (CspB) and glucose permease IIA domain. We compare the behavior of the free motif in solution and in the context of proteins of different size and architecture. Both Csp proteins lost contacts between the double-hairpin motif and the protein core as the first step of unfolding and proceeded to unfold with loss of the third β-strand, similar to the isolated WW domain. The glucose permease IIA domain is a larger protein and the contacts between the motif and the core were not lost as quickly. Instead the unfolding pathway of glucose permease IIA followed a different pathway with β1 pulling away from the sheet first. Interestingly, when the double hairpin motif was excised from the glucose permease IIA domain and simulated in isolation in water it unfolded by the same pathway as the WW domain, indicating that it is tertiary interactions with the protein that alter the motif's unfolding not a sequence dependent effect on its intrinsic unfolding behavior. With respect to the unfolding of the hairpins, there was no consistent order to the loss of hydrogen bonds between the β-strands in the hairpins in any of the systems. Our results show that while the folding behavior of the isolated WW domain is generally consistent with the double hairpin motif's behavior in the cold shock proteins, it is not the case for the glucose permease IIA domain. So, one must be cautious in extrapolating findings from model systems to larger more complicated proteins where tertiary interactions can overwhelm intrinsic behavior.

Contributions of aromatic pairs to the folding and stability of long-lived human γD-crystallin

Human γD-crystallin (HγD-Crys) is a highly stable protein that remains folded in the eye lens for the majority of an individual's lifetime. HγD-Crys exhibits two homologous crystallin domains, each containing two Greek key motifs with eight β-strands. Six aromatic pairs (four Tyr/Tyr, one Tyr/Phe and one Phe/Phe) are present in the β-hairpin sequences of the Greek keys. Ultraviolet damage to the aromatic residues in lens crystallins may contribute to the genesis of cataract. Mutant proteins with these aromatic residues substituted with alanines were constructed and expressed in E. coli. All mutant proteins except F115A and F117A had lower thermal stability than the WT protein. In equilibrium experiments in guanidine hydrochloride (GuHCl), all mutant proteins had lower thermodynamic stability than the WT protein. N-terminal domain (N-td) substitutions shifted the N-td transition to lower GuHCl concentration, but the C-terminal domain (C-td) transition remained unaffected. C-td substitutions led to a more cooperative unfolding/refolding process, with both the N-td and C-td transitions shifted to lower GuHCl concentration. The aromatic pairs conserved for each Greek key motif (Greek key pairs) had larger contributions to both thermal stability and thermodynamic stability than the other pairs. Aromatic-aromatic interaction was estimated as 1.5-2.0 kcal/mol. In kinetic experiments, N-td substitutions accelerated the early phase of unfolding, while C-td substitutions accelerated the late phase, suggesting independent domain unfolding. Only substitutions of the second Greek key pair of each crystallin domain slowed refolding. The second Greek keys may provide nucleation sites during the folding of the double-Greek-key crystallin domains.

Investigation of an anomalously accelerating substitution in the folding of a prototypical two-state protein

The folding rates of two-state single-domain proteins are generally resistant to small-scale changes in amino acid sequence. For example, having surveyed here over 700 single-residue substitutions in 24 well-characterized two-state proteins, we find that the majority (55%) of these substitutions affect folding rates by less than a factor of 2, and that only 9% affect folding rates by more than a factor of 8. Among those substitutions that significantly affect folding rates, we find that accelerating substitutions are an order of magnitude less common than those that decelerate the process. One of the most extreme outliers in this data set, an arginine-to-phenylalanine substitution at position 48 (R48F) of chymotrypsin inhibitor 2 (CI2), accelerates the protein's folding rate by a factor of 36 relative to that of the wild-type protein and is the most accelerating substitution reported to date in a two-state protein. In order to better understand the origins of this anomalous behavior, we have characterized the kinetics of multiple additional substitutions at this position. We find that substitutions at position 48 in CI2 fall into two distinct classes. The first, comprising residues that ablate the charge of the wild-type arginine but retain the hydrophobicity of its alkane chain, accelerate folding by at least 10-fold. The second class, comprising all other residues, produces folding rates within a factor of two of the wild-type rate. A significant positive correlation between hydrophobicity and folding rate across all of the residues we have characterized at this position suggests that the hydrophobic methylene units of the wild-type arginine play a significant role in stabilizing the folding transition state. Likewise, studies of the pH dependence of the histidine substitution indicate a strong correlation between folding rate and charge state. Thus, mutations that ablate the arginine's positive charge while retaining the hydrophobic contacts of its methylene units tend to dramatically accelerate folding. Previous studies have suggested that arginine 48 plays an important functional role in CI2, which may explain why it is highly conserved despite the anomalously large deceleration it produces in the folding of this protein.

The critical role of N- and C-terminal contact in protein stability and folding of a family 10 xylanase under extreme conditions

Background

Stabilization strategies adopted by proteins under extreme conditions are very complex and involve various kinds of interactions. Recent studies have shown that a large proportion of proteins have their N- and C-terminal elements in close contact and suggested they play a role in protein folding and stability. However, the biological significance of this contact remains elusive.

Methodology

In the present study, we investigate the role of N- and C-terminal residue interaction using a family 10 xylanase (BSX) with a TIM-barrel structure that shows stability under high temperature, alkali pH, and protease and SDS treatment. Based on crystal structure, an aromatic cluster was identified that involves Phe4, Trp6 and Tyr343 holding the N- and C-terminus together; this is a unique and important feature of this protein that might be crucial for folding and stability under poly-extreme conditions.

Conclusion

A series of mutants was created to disrupt this aromatic cluster formation and study the loss of stability and function under given conditions. While the deletions of Phe4 resulted in loss of stability, removal of Trp6 and Tyr343 affected in vivo folding and activity. Alanine substitution with Phe4, Trp6 and Tyr343 drastically decreased stability under all parameters studied. Importantly, substitution of Phe4 with Trp increased stability in SDS treatment. Mass spectrometry results of limited proteolysis further demonstrated that the Arg344 residue is highly susceptible to trypsin digestion in sensitive mutants such as ΔF4, W6A and Y343A, suggesting again that disruption of the Phe4-Trp6-Tyr343 (F-W-Y) cluster destabilizes the N- and C-terminal interaction. Our results underscore the importance of N- and C-terminal contact through aromatic interactions in protein folding and stability under extreme conditions, and these results may be useful to improve the stability of other proteins under suboptimal conditions.

Role of the Charge–Charge Interactions in Defining Stability and Halophilicity of the CspB Proteins

Charge-charge interactions on the surface of native proteins are important for protein stability and can be computationally redesigned in a rational way to modulate protein stability. Such computational effort led to an engineered protein, CspB-TB that has the same core as the mesophilic cold shock protein CspB-Bs from Bacillus subtilis, but optimized distribution of charge-charge interactions on the surface. The CspB-TB protein shows an increase in the transition temperature by 20 degrees C relative to the unfolding temperature of CspB-Bs. The CspB-TB and CspB-Bs protein pair offers a unique opportunity to further explore the energetics of charge-charge interactions as the substitutions at the same sequence positions are done in largely similar structural but different electrostatic environments. In particular we addressed two questions. What is the contribution of charge-charge interactions in the unfolded state to the protein stability and how amino acid substitutions modulate the effect of increase in ionic strength on protein stability (i.e. protein halophilicity). To this end, we experimentally measured the stabilities of over 100 variants of CspB-TB and CspB-Bs proteins with substitutions at charged residues. We also performed computational modeling of these protein variants. Analysis of the experimental and computational data allowed us to conclude that the charge-charge interactions in the unfolded state of two model proteins CspB-Bs and CspB-TB are not very significant and computational models that are based only on the native state structure can adequately, i.e. qualitatively (stabilizing versus destabilizing) and semi-quantitatively (relative rank order), predict the effects of surface charge neutralization or reversal on protein stability. We also show that the effect of ionic strength on protein stability (protein halophilicity) appears to be mainly due to the screening of the long-range charge-charge interactions.

Non-linear effects of temperature and urea on the thermodynamics and kinetics of folding and unfolding of hisactophilin

Extensive measurements and analysis of thermodynamic stability and kinetics of urea-induced unfolding and folding of hisactophilin are reported for 5-50 degrees C, at pH 6.7. Under these conditions hisactophilin has moderate thermodynamic stability, and equilibrium and kinetic data are well fit by a two-state transition between the native and the denatured states. Equilibrium and kinetic m values decrease with increasing temperature, and decrease with increasing denaturant concentration. The betaF values at different temperatures and urea concentrations are quite constant, however, at about 0.7. This suggests that the transition state for hisactophilin unfolding is native-like and changes little with changing solution conditions, consistent with a narrow free energy profile for the transition state. The activation enthalpy and entropy of unfolding are unusually low for hisactophilin, as is also the case for the corresponding equilibrium parameters. Conventional Arrhenius and Eyring plots for both folding and unfolding are markedly non-linear, but these plots become linear for constant DeltaG/T contours. The Gibbs free energy changes for structural changes in hisactophilin have a non-linear denaturant dependence that is comparable to non-linearities observed for many other proteins. These non-linearities can be fit for many proteins using a variation of the Tanford model, incorporating empirical quadratic denaturant dependencies for Gibbs free energies of transfer of amino acid constituents from water to urea, and changes in fractional solvent accessible surface area of protein constituents based on the known protein structures. Noteworthy exceptions that are not well fit include amyloidogenic proteins and large proteins, which may form intermediates. The model is easily implemented and should be widely applicable to analysis of urea-induced structural transitions in proteins.

A segment of cold shock protein directs the folding of a combinatorial protein

It has been suggested that protein domains evolved by the non-homologous recombination of building blocks of subdomain size. In earlier work we attempted to recapitulate domain evolution in vitro. We took a polypeptide segment comprising three beta-strands in the monomeric, five-stranded beta-barrel cold shock protein (CspA) of Escherichia coli as a building block. This segment corresponds to a complete exon in homologous eukaryotic proteins and includes residues that nucleate folding in CspA. We recombined this segment at random with fragments of natural proteins and succeeded in generating a range of folded chimaeric proteins. We now present the crystal structure of one such combinatorial protein, 1b11, a 103-residue polypeptide that includes segments from CspA and the S1 domain of the 30S ribosomal subunit of E. coli. The structure reveals a segment-swapped, six-stranded beta-barrel of unique architecture that assembles to a tetramer. Surprisingly, the CspA segment retains its structural identity in 1b11, recapitulating its original fold and deforming the structure of the S1 segment as necessary to complete a barrel. Our work provides structural evidence that (i) random shuffling of nonhomologous polypeptide segments can lead to folded proteins and unique architectures, (ii) many structural features of the segments are retained, and (iii) some segments can act as templates around which the rest of the protein folds.

Structures of the agouti signaling protein

Expression of the agouti signaling protein (ASIP) during hair growth produces the red/yellow pigment pheomelanin. ASIP, and its neuropeptide homolog the agouti-related protein (AgRP) involved in energy balance, are novel, paracrine signaling molecules that act as inverse agonists at distinct subsets of melanocortin receptors. Ubiquitous ASIP expression in mice gives rise to a pleiotropic phenotype characterized by a uniform yellow coat color, obesity, overgrowth, and metabolic derangements similar to type II diabetes in humans. Here we report the synthesis and NMR structure of ASIP's active, cysteine-rich, C-terminal domain. ASIP adopts the inhibitor cystine knot fold and, along with AgRP, are the only known mammalian proteins in this structure class. Moreover, ASIP populates two distinct conformers resulting from a cis peptide bond at Pro102-Pro103 and a coexistence of cis/trans isomers of Ala104-Pro105. Pharmacologic studies of Pro-->Ala mutants demonstrate that the minor conformation with two cis peptide bonds is responsible for activity at all MCRs. The loop containing the heterogeneous Ala-Pro peptide bond is conserved in mammals, and suggests that ASIP is either trapped by evolution in this unusual configuration or possesses function outside of strict MCR antagonism.

The importance of explicit chain representation in protein folding models: an examination of Ising-like models

A class of models that represents a protein chain as a sequence of "folded" and "unfolded" residues has recently been used to correlate rates and mechanisms of protein folding with the protein native structure. In order to better understand the conditions under which these "Ising-like" models apply, we compare results from this model to those obtained from an off-lattice model which uses the same potential function. We find that Ising-like models by construction impose folding via a highly sequential nucleation-condensation mechanism, which in turn leads to more rugged energy landscapes, fewer "pathways" to the native state, and in the specific case examined here, the cold shock protein A from Escherichia coli, a qualitative difference in the most likely order of events in folding.

Evidence for sequential barriers and obligatory intermediates in apparent two-state protein folding

Many small proteins fold fast and without detectable intermediates. This is frequently taken as evidence against the importance of partially folded states, which often transiently accumulate during folding of larger proteins. To get insight into the properties of free energy barriers in protein folding we analyzed experimental data from 23 proteins that were reported to show non-linear activation free-energy relationships. These non-linearities are generally interpreted in terms of broad transition barrier regions with a large number of energetically similar states. Our results argue against the presence of a single broad barrier region. They rather indicate that the non-linearities are caused by sequential folding pathways with consecutive distinct barriers and a few obligatory high-energy intermediates. In contrast to a broad barrier model the sequential model gives a consistent picture of the folding barriers for different variants of the same protein and when folding of a single protein is analyzed under different solvent conditions. The sequential model is also able to explain changes from linear to non-linear free energy relationships and from apparent two-state folding to folding through populated intermediates upon single point mutations or changes in the experimental conditions. These results suggest that the apparent discrepancy between two-state and multi-state folding originates in the relative stability of the intermediates, which argues for the importance of partially folded states in protein folding.

Unfolding and double-stranded DNA binding of the cold shock protein homologue Cla h 8 from Cladosporium herbarum

The cloning, purification, and biophysical characterization of the first eukaryotic cold shock protein homologue, Cla h 8, expressed as single functional polypeptide is reported here. It was discovered as a minor allergen of the mold Cladosporium herbarum by phage display using a library selectively enriched for IgE-binding proteins. Based on the sequence homology of Cla h 8 with bacterial cold shock proteins (CSPs), a homology-based computer model of the allergen was computed indicating an all-beta structure of Cla h 8. This major structural feature was confirmed by CD spectroscopy. Despite the structural similarities with bacterial CSPs, the DNA-binding and unfolding behavior of Cla h 8 exhibited unique and previously undescribed characteristics. High affinities of Cla h 8 for single-stranded DNA as well as for double-stranded DNA corresponding to the human Y-box were detected. The affinity for double-stranded DNA increased significantly with decreasing temperature, which was paralleled by an increase in the beta sheet content of the protein. Temperature-dependent fluorescence anisotropy and far-UV CD measurements revealed different unfolding transitions at 28 and at 35.7 degrees C, respectively, indicating a multistate transition, which is uncommon for CSPs. The enhanced affinity for DNA at low temperatures together with the low unfolding transition refer to the functional significance of Cla h 8 at reduced temperatures.

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.