Correspondence between anomalous m- and DeltaCp-values in protein folding

Proteins folding according to a classical two-state system characteristically show V-shaped chevron plots. We have previously interpreted the symmetrically curved chevron plot of the protein U1A as denaturant-dependent movements in the position of the transition state ensemble (TSE). S6, a structural analog of U1A, shows a classical V-shaped chevron plot indicative of straightforward two-state kinetics, but the mutant LA30 has a curved unfolding limb, which is most consistent with TSE mobility. The kinetic m-values (derivatives of the rate constants with respect to denaturant concentration) in themselves depend on denaturant concentration. To obtain complementary information about putative mobile TSEs, we have carried out a thermodynamic analysis of the three proteins, based on data for refolding and unfolding over the range 10 degrees C to 70 degrees C. The data at all temperatures can be fitted to two-state model systems. Importantly, for all three proteins the activation heat capacities are, within error, identical to the heat capacities measured in independent experiments under equilibrium conditions. Although the equilibrium heat capacities are essentially invariant with regard to denaturant concentration, the activation heat capacities, similar to the structurally equivalent kinetic m-values, show marked denaturant dependence. Furthermore, the values of beta++ at different denaturant concentrations measured by m-values and by heat capacity values are very similar. These observations are consistent with significant transition state movements within the framework of two-state folding. The basis for TSE movement appears to be enthalpic rather than entropic, suggesting that the binding energy of denaturant-protein interactions is a major determinant of the response of energy landscape contours to changing environments.     

An osmolyte mitigates the destabilizing effect of protein crowding

Most theories predict that macromolecular crowding stabilizes globular proteins, but recent studies show that weak attractive interactions can result in crowding-induced destabilization. Osmolytes are ubiquitous in biology and help protect cells against stress. Given that dehydration stress adds to the crowded nature of the cytoplasm, we speculated that cells might use osmolytes to overcome the destabilization caused by the increased weak interactions that accompany desiccation. We used NMR-detected amide proton exchange experiments to measure the stability of the test protein chymotrypsin inhibitor 2 under physiologically relevant crowded conditions in the presence and absence of the osmolyte glycine betaine. The osmolyte overcame the destabilizing effect of the cytosol. This result provides a physiologically relevant explanation for the accumulation of osmolytes by dehydration-stressed cells.     

Structural thermodynamics of protein preferential solvation: osmolyte solvation of proteins, aminoacids, and peptides

Protein stability and solubility depend strongly on the presence of osmolytes, because of the protein preference to be solvated by either water or osmolyte. It has traditionally been assumed that only this relative preference can be measured, and that the individual solvation contributions of water and osmolyte are inaccessible. However, it is possible to determine hydration and osmolyte solvation (osmolation) separately using Kirkwood-Buff theory, and this fact has recently been utilized by several researchers. Here, we provide a thermodynamic assessment of how each surface group on proteins contributes to the overall hydration and osmolation. Our analysis is based on transfer free energy measurements with model-compounds that were previously demonstrated to allow for a very successful prediction of osmolyte-dependent protein stability. When combined with Kirkwood-Buff theory, the Transfer Model provides a space-resolved solvation pattern of the peptide unit, amino acids, and the folding/unfolding equilibrium of proteins in the presence of osmolytes. We find that the major solvation effects on protein side-chains originate from the osmolytes, and that the hydration mostly depends on the size of the side-chain. The peptide backbone unit displays a much more variable hydration in the different osmolyte solutions. Interestingly, the presence of sucrose leads to simultaneous accumulation of both the sugar and water in the vicinity of peptide groups, resulting from a saccharide accumulation that is less than the accumulation of water, a net preferential exclusion. Only the denaturing osmolyte, urea, obeys the classical solvent exchange mechanism in which the preferential interaction with the peptide unit excludes water.     

Folding recombinant spider‐silk in H2O: Effect of osmolytes on the solution conformation of a 15‐repeat spider‐silk mimetic

The folding of a recombinant spider silk protein‐polymer in the presence of the tri‐methylamine osmolytes TMANO and Betaine in 80% H₂O is reported. Circular dichroism measurements (CD) reveal an increase in α‐helical secondary structure with increasing osmolyte concentrations, as determined by an increase in ellipticity at 222 nm. Consistent with this observation, the signal for random coil sampling, observed at 205 nm, is greatly reduced with increasing trimethylamine. Fluorescence spectra of a single tyrosine positioned within the conserved 33‐amino acid repeat primary sequence (of the spider‐silk mimetic) complements the conformational changes observed by CD. Importantly, there is a correlation between the number of Alkyl‐groups (CH₃‐) on the amine of the osmolyte and enhanced helicity of the 15‐repeat silk‐mimetic for the osmolytes tested, ie TMANO, Betaine, Sarcosine and Glycine. These preliminary results are applicable to storing and processing recombinant silk sequences in H₂O, an important mile‐stone for widespread use of recombinant silk polymers.     

Analysis of the stability of multimeric proteins by effective DeltaG and effective m-values

Analyzing the stability of a multimeric protein is challenging because of the intrinsic difficulty in handling the mathematical model for the folded multimer-unfolded monomer equilibrium. To circumvent this problem, we introduce the concept of effective stability, DeltaGeff (= -RTlnKeff), where Keff is the equilibrium constant expressed in monomer units. Analysis of the denaturant effect on DeltaGeff gives new insight into the stability of multimeric proteins. When a multimeric protein is mostly folded, the dependence of effective stability on denaturant concentration (effective m-value) is simply the m-value of its monomeric unit. However, when the protein is mostly unfolded, its stability depends on denaturant concentration with the m-value of its multimeric form. We also find that the effective m-value at the Cm is a good approximation of the apparent m-value determined by fitting the equilibrium unfolding data from multimeric proteins with a two-state monomer model. Moreover, when the m-value of a monomeric unit is estimated from its size, the effective stability of a multimeric protein can be determined simply from Cm and this estimated m-value. These simple and intuitive approaches will allow a facile analysis of the stability of multimeric proteins. These analyses are also applicable for high-throughput analysis of protein stability on a proteomic scale.     

BPPred: a Web-based computational tool for predicting biophysical parameters of proteins

We exploit the availability of recent experimental data on a variety of proteins to develop a Web-based prediction algorithm (BPPred) to calculate several biophysical parameters commonly used to describe the folding process. These parameters include the equilibrium m-values, the length of proteins, and the changes upon unfolding in the solvent-accessible surface area, in the heat capacity, and in the radius of gyration. We also show that the knowledge of any one of these quantities allows an estimate of the others to be obtained, and describe the confidence limits with which these estimations can be made. Furthermore, we discuss how the kinetic m-values, or the Beta Tanford values, may provide an estimate of the solvent-accessible surface area and the radius of gyration of the transition state for protein folding. Taken together, these results suggest that BPPred should represent a valuable tool for interpreting experimental measurements, as well as the results of molecular dynamics simulations.     

Cooperative folding in a multi-domain protein

Most protein domains are found in multi-domain proteins, yet most studies of protein folding have concentrated on small, single-domain proteins or on isolated domains from larger proteins. Spectrin domains are small (106 amino acid residues), independently folding domains consisting of three long alpha-helices. They are found in multi-domain proteins with a number of spectrin domains in tandem array. Structural studies have shown that in these arrays the last helix of one domain forms a continuous helix with the first helix of the following domain. It has been demonstrated that a number of spectrin domains are stabilised by their neighbours. Here we investigate the molecular basis for cooperativity between adjacent spectrin domains 16 and 17 from chicken brain alpha-spectrin (R16 and R17). We show that whereas the proteins unfold as a single cooperative unit at 25 degrees C, cooperativity is lost at higher temperatures and in the presence of stabilising salts. Mutations in the linker region also cause the cooperativity to be lost. However, the cooperativity does not rely on specific interactions in the linker region alone. Most mutations in the R17 domain cause a decrease in cooperativity, whereas proteins with mutations in the R16 domain still fold cooperatively. We propose a mechanism for this behaviour.     

Analysis of the pH-dependent folding and stability of histidine point mutants allows characterization of the denatured state and transition state for protein folding

pH-Dependent studies of the folding kinetics and stability of a set of His to Gln point mutants were used to characterize the denatured state and transition state ensembles for the C-terminal domain of the ribosomal protein L9 (CTL9). CTL9 contains three histidine residues, two of which, H106 and H134, are buried in the native state, while the third, H144, is more exposed. Comparison of the pH-dependent stability calculated using the Tanford-Wyman linkage relationship to the measured values demonstrates that the apparent pK(a) values of the three histidine residues are not significantly perturbed in the denatured state ensemble. Kinetic measurements show that mutation of H134 has a larger effect on the folding process than does mutation of H106 and H144. The Phi-value for H134 is significantly larger than the Phi-values for the other histidine residues, which are near zero at both pH 5.45 and pH 8.0. The Phi-value for H134 is higher, 0.55, at pH 8.0 than at pH 5.45, 0.39. At pH 5.45, H134 is protonated in the unfolded state but deprotonated in the native state, while at pH 8.0 it is deprotonated in both. There is an excellent linear relationship between stability (logK) and folding rates (logk(f)) over the range of pH 5-9 for all mutants. From these plots, the ratio of DeltaQ( not equal)/DeltaQ can be calculated for each mutant. DeltaQ( not equal) is the difference in the number of protons bound to the transition state and to the unfolded state, while DeltaQ represents the difference between folded and denatured state. The linear plots indicate that the relative position of the transition state ensemble as judged by DeltaQ( not equal)/DeltaQ is independent of pH. The linkage analysis is consistent with the Phi-value analysis, showing that H134 is the most critical contributor to the development of pH-dependent interactions, including desolvation effects in the transition state ensemble.     

FRET-based in vivo screening for protein folding and increased protein stability

Fluorescence resonance energy transfer (FRET) was used to establish a novel in vivo screening system that allows rapid detection of protein folding and protein variants with increased thermodynamic stability in the cytoplasm of Escherichia coli. The system is based on the simultaneous fusion of the green fluorescent protein (GFP) to the C terminus of a protein X of interest, and of blue-fluorescent protein (BFP) to the N terminus of protein X. Efficient FRET from BFP to GFP in the ternary fusion protein is observed in vivo only when protein X is folded and brings BFP and GFP into close proximity, while FRET is lost when BFP and GFP are far apart due to unfolding or intracellular degradation of protein X. The screening system was validated by identification of antibody V(L) intradomains with increased thermodynamic stabilities from expression libraries after random mutagenesis, bacterial cell sorting, and colony screening.     

Functional differences in yeast protein disulfide isomerases

PDI1 is the essential gene encoding protein disulfide isomerase in yeast. The Saccharomyces cerevisiae genome, however, contains four other nonessential genes with homology to PDI1: MPD1, MPD2, EUG1, and EPS1. We have investigated the effects of simultaneous deletions of these genes. In several cases, we found that the ability of the PDI1 homologues to restore viability to a pdi1-deleted strain when overexpressed was dependent on the presence of low endogenous levels of one or more of the other homologues. This shows that the homologues are not functionally interchangeable. In fact, Mpd1p was the only homologue capable of carrying out all the essential functions of Pdi1p. Furthermore, the presence of endogenous homologues with a CXXC motif in the thioredoxin-like domain is required for suppression of a pdi1 deletion by EUG1 (which contains two CXXS active site motifs). This underlines the essentiality of protein disulfide isomerase-catalyzed oxidation. Most mutant combinations show defects in carboxypeptidase Y folding as well as in glycan modification. There are, however, no significant effects on ER-associated protein degradation in the various protein disulfide isomerase-deleted strains.     

On the precision of experimentally determined protein folding rates and phi-values

Phi-values, a relatively direct probe of transition-state structure, are an important benchmark in both experimental and theoretical studies of protein folding. Recently, however, significant controversy has emerged regarding the reliability with which phi-values can be determined experimentally: Because phi is a ratio of differences between experimental observables it is extremely sensitive to errors in those observations when the differences are small. Here we address this issue directly by performing blind, replicate measurements in three laboratories. By monitoring within- and between-laboratory variability, we have determined the precision with which folding rates and phi-values are measured using generally accepted laboratory practices and under conditions typical of our laboratories. We find that, unless the change in free energy associated with the probing mutation is quite large, the precision of phi-values is relatively poor when determined using rates extrapolated to the absence of denaturant. In contrast, when we employ rates estimated at nonzero denaturant concentrations or assume that the slopes of the chevron arms (mf and mu) are invariant upon mutation, the precision of our estimates of phi is significantly improved. Nevertheless, the reproducibility we thus obtain still compares poorly with the confidence intervals typically reported in the literature. This discrepancy appears to arise due to differences in how precision is calculated, the dependence of precision on the number of data points employed in defining a chevron, and interlaboratory sources of variability that may have been largely ignored in the prior literature.     

Stabilization of yeast hexokinase A by polyol osmolytes: correlation with the physicochemical properties of aqueous solutions

Osmolytes of the polyol series are known to accumulate in biological systems under stress and stabilize the structures of a wide variety of proteins. While increased surface tension of aqueous solutions has been considered an important factor in protein stabilization effect, glycerol is an exception, lowering the surface tension of water. To clarify this anomalous effect, the effect of a series of polyols on the thermal stability of a highly thermolabile two domain protein yeast hexokinase A has been investigated by differential scanning calorimetry and by monitoring loss in the biological activity of the enzyme as a function of time. A larger increase in the T(m) of domain 1 compared with that of domain 2, varying linearly with the number of hydroxyl groups in polyols, has been observed, sorbitol being the best stabilizer against both thermal as well as urea denaturation. Polyols help retain the activity of the enzyme considerably and a good correlation of the increase in T(m) (DeltaT(m)) and the retention of activity with the increase in the surface tension of polyol solutions, except glycerol, which breaks this trend, has been observed. However, the DeltaT(m) values show a linear correlation with apparent molal heat capacity and volume of aqueous polyol solutions including glycerol. These results suggest that while bulk solution properties contribute significantly to protein stabilization, interfacial properties are not always a good indicator of the stabilizing effect. A subtle balance of various weak binding and exclusion effects of the osmolytes mediated by water further regulates the stabilizing effect. Understanding these aspects is critical in the rational design of stable protein formulations.     

Probing the folding intermediate of Rd-apocyt b562 by protein engineering and infrared T-jump

Small proteins often fold in an apparent two-state manner with the absence of detectable early-folding intermediates. Recently, using native-state hydrogen exchange, intermediates that exist after the rate-limiting transition state have been identified for several proteins. However, little is known about the folding kinetics from these post-transition intermediates to their corresponding native states. Herein, we have used protein engineering and a laser-induced temperature-jump (T-jump) technique to investigate this issue and have applied it to Rd-apocyt b(562) , a four-helix bundle protein. Previously, it has been shown that Rd-apocyt b(562) folds via an on-pathway hidden intermediate, which has only the N-terminal helix unfolded. In the present study, a double mutation (V16G/I17A) in the N-terminal helix of Rd-apocyt b(562) was made to further increase the relative population of this intermediate state at high temperature by selectively destabilizing the native state. In the circular dichroism thermal melting experiment, this mutant showed apparent two-state folding behavior. However, in the T-jump experiment, two kinetic phases were observed. Therefore, these results are in agreement with the idea that a folding intermediate is populated on the folding pathway of Rd-apocyt b(562) . Moreover, it was found that the exponential growth rate of the native state from this intermediate state is roughly (25 microsec)(-1) at 65 degrees C.     

The folding landscape of Streptomyces griseus protease B reveals the energetic costs and benefits associated with evolving kinetic stability

Like most extracellular bacterial proteases, Streptomyces griseus protease B (SGPB) and alpha-lytic protease (alphaLP) are synthesized with covalently attached pro regions necessary for their folding. In this article, we characterize the folding free energy landscape of SGPB and compare it to the folding landscapes of alphaLP and trypsin, a mammalian homolog that folds independently of its zymogen peptide. In contrast to the thermodynamically stable native state of trypsin, SGPB and alphaLP fold to native states that are thermodynamically marginally stable or unstable, respectively. Instead, their apparent stability arises kinetically, from unfolding free energy barriers that are both large and highly cooperative. The unique unfolding transitions of SGPB and alphaLP extend their functional lifetimes under highly degradatory conditions beyond that seen for trypsin; however, the penalty for evolving kinetic stability is remarkably large in that each factor of 2.4-8 in protease resistance is accompanied by a cost of ~10(5) in the spontaneous folding rate and ~5-9 kcal/mole in thermodynamic stability. These penalties have been overcome by the coevolution of increasingly effective pro regions to facilitate folding. Despite these costs, kinetic stability appears to be a potent mechanism for developing native-state properties that maximize protease longevity.     

Folding and association of an extremely stable dimeric protein from Sulfolobus islandicus

ORF56 is a plasmid-encoded protein from Sulfolobus islandicus, which probably controls the copy number of the pRN1 plasmid by binding to its own promotor. The protein showed an extremely high stability in denaturant, heat, and pH-induced unfolding transitions, which can be well described by a two-state reaction between native dimers and unfolded monomers. The homodimeric character of native ORF56 was confirmed by analytical ultracentrifugation. Far-UV circular dichroism and fluorescence spectroscopy gave superimposable denaturant-induced unfolding transitions and the midpoints of both heat as well as denaturant-induced unfolding depend on the protein concentration supporting the two-state model. This model was confirmed by GdmSCN-induced unfolding monitored by heteronuclear 2D NMR spectroscopy. Chemical denaturation was accomplished by GdmCl and GdmSCN, revealing a Gibbs free energy of stabilization of -85.1 kJ/mol at 25 degrees C. Thermal unfolding was possible only above 1 M GdmCl, which shifted the melting temperature (t(m)) below the boiling point of water. Linear extrapolation of t(m) to 0 M GdmCl yielded a t(m) of 107.5 degrees C (5 microM monomer concentration). Additionally, ORF56 remains natively structured over a remarkable pH range from pH 2 to pH 12. Folding kinetics were followed by far-UV CD and fluorescence after either stopped-flow or manual mixing. All kinetic traces showed only a single phase and the two probes revealed coincident folding rates (k(f), k(u)), indicating the absence of intermediates. Apparent first-order refolding rates depend linearly on the protein concentration, whereas the unfolding rates do not. Both lnk(f) and lnk(u) depend linearly on the GdmCl concentration. Together, folding and association of homodimeric ORF56 are concurrent events. In the absence of denaturant ORF56 refolds fast (7.0 x 10(7)M(-1)s(-1)) and unfolds extremely slowly (5.7 year(-1)). Therefore, high stability is coupled to a slow unfolding rate, which is often observed for proteins of extremophilic organisms.     

An exclusively nuclear RNA-binding protein affects asymmetric localization of ASH1 mRNA and Ash1p in yeast

The localization of ASH1 mRNA to the distal tip of budding yeast cells is essential for the proper regulation of mating type switching in Saccharomyces cerevisiae. A localization element that is predominantly in the 3'-untranslated region (UTR) can direct this mRNA to the bud. Using this element in the three-hybrid in vivo RNA-binding assay, we identified a protein, Loc1p, that binds in vitro directly to the wild-type ASH1 3'-UTR RNA, but not to a mutant RNA incapable of localizing to the bud nor to several other mRNAs. LOC1 codes for a novel protein that recognizes double-stranded RNA structures and is required for efficient localization of ASH1 mRNA. Accordingly, Ash1p gets symmetrically distributed between daughter and mother cells in a loc1 strain. Surprisingly, Loc1p was found to be strictly nuclear, unlike other known RNA-binding proteins involved in mRNA localization which shuttle between the nucleus and the cytoplasm. We propose that efficient cytoplasmic ASH1 mRNA localization requires a previous interaction with specific nuclear factors.     

A part of ice nucleation protein exhibits the ice-binding ability

We generated a recombinant 96-residue polypeptide corresponding to a sequence Tyr176-Gly273 of ice nucleation protein from Pseudomonas syringae (denoted INP96). INP96 exhibited an ability to shape an ice crystal, whose morphology is highly similar to the hexagonal-bipyramid generally identified for antifreeze protein. INP96 also showed a non-linear, concentration-dependent retardation of ice growth. Additionally, circular dichroism and NMR measurements suggested a local structural construction in INP96, which undergoes irreversible thermal denaturation. These data imply that a part of INP constructs a unique structure so as to interact with the ice crystal surfaces.     

The osmolyte trimethylamine‐N‐oxide stabilizes the Fyn SH3 domain without altering the structure of its folding transition state

Trimethylamine‐N‐oxide (TMAO) is a naturally occurring osmolyte that stabilizes proteins against denaturation. Although the impact of TMAO on the folding thermodynamics of many proteins has been well characterized, far fewer studies have investigated its effects on protein folding kinetics. In particular, no previous studies have used Φ‐value analysis to determine whether TMAO may alter the structure of the folding transition state. Here we have measured the effects on folding kinetics of 16 different amino acid substitutions distributed across the structure of the Fyn SH3 domain both in the presence and absence of TMAO. The folding and unfolding rates in TMAO, on average, improved to equivalent degrees, with a twofold increase in the protein folding rate accompanied by a twofold decrease in the unfolding rate. Importantly, TMAO caused little alteration to the Φ‐values of the mutants tested, implying that this compound minimally perturbs the folding transition state structure. Furthermore, the solvent accessibility of the transition state was not altered as reflected in an absence of a TMAO‐induced change in the denaturant β factors. Through TMAO‐induced folding studies, a β factor of 0.5 was calculated for this compound, suggesting that the protein backbone, which is the target of action of TMAO, is 50% exposed in the transition state as compared to the native state. This finding is consistent with the equivalent effects of TMAO on the folding and unfolding rates. Through thermodynamic analysis of mutants, we also discovered that the stabilizing effect of TMAO is lessened with increasing temperature.     

Expression of the sucrose binding protein from soybean: renaturation and stability of the recombinant protein

The sucrose binding protein (SBP) belongs to the cupin family of proteins and is structurally related to vicilin-like storage proteins. In this investigation, a SBP isoform (GmSBP2/S64) was expressed in E. coli and large amounts of the protein accumulated in the insoluble fraction as inclusion bodies. The renatured protein was studied by circular dichroism (CD), intrinsic fluorescence, and binding of the hydrophobic probes ANS and Bis-ANS. The estimated content of secondary structure of the renatured protein was consistent with that obtained by theoretical modeling with a large predominance of beta-strand structure (42%) over the alpha-helix (9.9%). The fluorescence emission maximum of 303 nm for SBP2 indicated that the fluorescent tryptophan was completely buried within a highly hydrophobic environment. We also measured the equilibrium dissociation constant (K(d)) of sucrose binding by fluorescence titration using the refolded protein. The low sucrose binding affinity (K(d)=2.79+/-0.22 mM) of the renatured protein was similar to that of the native protein purified from soybean seeds. Collectively, these results indicate that the folded structure of the renatured protein was similar to the native SBP protein. As a member of the bicupin family of proteins, which includes highly stable seed storage proteins, SBP2 was fairly stable at high temperatures. Likewise, it remained folded to a similar extent in the presence or absence of 7.6M urea or 6.7 M GdmHCl. The high stability of the renatured protein may be a reminiscent property of SBP from its evolutionary relatedness to the seed storage proteins.