Role of the chain termini for the folding transition state of the cold shock protein

Effect of imidazolium based ionic liquids on CO-association dynamics and thermodynamic stability of Ferrocytochrome c

Analysis of kinetic and thermodynamic parameters measured for CO-association reaction of Ferrocytochrome c (Ferrocyt c) under variable concentrations of 1-butyl-3-methylimidazolium with varying anion ([Bmim]X) (X = Cl^-, I^-, Br^-, HSO₄^-) at pH 7 revealed that the low concentration of [Bmim]X (≤0.5 M) constrains the CO-association dynamics of Ferrocyt c and typically follows the order: [Bmim]HSO₄ > [Bmim]Cl > [Bmim]Br > [Bmim]I. At relatively higher concentrations (>0.5), the chaotropic action of [Bmim]⁺ dominates which consequently increases the thermal-fluctuations responsible to denature the protein and thus accelerates the speed of CO-association reaction. Analysis of thermal denaturation curves of Ferrocyt c measured at different concentrations of [Bmim]X revealed that the [Bmim]X decreases the thermodynamic stability of protein and typically follows the order: [Bmim]I > [Bmim]Br > [Bmim]Cl > [Bmim]CH₃COO > [Bmim]HSO₄, demonstrating that the effect of [Bmim]X on thermodynamic stability of protein is not in accordance to Hofmeister series effect of anions because instead of increasing the kosmotropic anion carrying [Bmim]X ([Bmim]CH₃COO and [Bmim]HSO₄) also decreases the thermodynamic stability of protein.

Electrostatic Interactions in Protein Structure, Folding, Binding, and Condensation

Charged and polar groups, through forming ion pairs, hydrogen bonds, and other less specific electrostatic interactions, impart important properties to proteins. Modulation of the charges on the amino acids, e.g., by pH and by phosphorylation and dephosphorylation, have significant effects such as protein denaturation and switch-like response of signal transduction networks. This review aims to present a unifying theme among the various effects of protein charges and polar groups. Simple models will be used to illustrate basic ideas about electrostatic interactions in proteins, and these ideas in turn will be used to elucidate the roles of electrostatic interactions in protein structure, folding, binding, condensation, and related biological functions. In particular, we will examine how charged side chains are spatially distributed in various types of proteins and how electrostatic interactions affect thermodynamic and kinetic properties of proteins. Our hope is to capture both important historical developments and recent experimental and theoretical advances in quantifying electrostatic contributions of proteins.

The entropic nature of protein thermal stabilization

We performed thermodynamic analysis of temperature-induced unfolding of mesophilic and thermophilic proteins. It was shown that the variability in protein thermostability associated with pH-dependent unfolding or linked to the substitution of amino acid residues on the protein surface is evidence of the governing role of the entropy factor. Numerical values of conformational components in enthalpy, entropy and free energy which characterize protein unfolding in the "gas phase" were obtained. Based on the calculated absolute values of entropy and free energy, a model of protein unfolding is proposed in which the driving force is the conformational entropy of native protein, as an energy of the heat motion (T·S(NC)) increasing with temperature and acting as an factor devaluating the energy of intramolecular weak bonds in the transition state.

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.

The origin of nonmonotonic complex behavior and the effects of nonnative interactions on the diffusive properties of protein folding

We present a method for calculating the configurational-dependent diffusion coefficient of a globular protein as a function of the global folding process. Using a coarse-grained structure-based model, we determined the diffusion coefficient, in reaction coordinate space, as a function of the fraction of native contacts formed Q for the cold shock protein (TmCSP). We find nonmonotonic behavior for the diffusion coefficient, with high values for the folded and unfolded ensembles and a lower range of values in the transition state ensemble. We also characterized the folding landscape associated with an energetically frustrated variant of the model. We find that a low-level of frustration can actually stabilize the native ensemble and increase the associated diffusion coefficient. These findings can be understood from a mechanistic standpoint, in that the transition state ensemble has a more homogeneous structural content when frustration is present. Additionally, these findings are consistent with earlier calculations based on lattice models of protein folding and more recent single-molecule fluorescence measurements.

What lessons can be learned from studying the folding of homologous proteins?

The studies of the folding of structurally related proteins have proved to be a very important tool for investigating protein folding. Here we review some of the insights that have been gained from such studies. Our highlighted studies show just how such an investigation should be designed and emphasise the importance of the synergy between experiment and theory. We also stress the importance of choosing the right system carefully, exploiting the excellent structural and sequence databases at our disposal.

Mapping protein collapse with single-molecule fluorescence and kinetic synchrotron radiation circular dichroism spectroscopy

We have used the combination of single-molecule Förster resonance energy transfer and kinetic synchrotron radiation circular dichroism experiments to probe the conformational ensemble of the collapsed unfolded state of the small cold shock protein CspTm under near-native conditions. This regime is physiologically most relevant but difficult to access experimentally, because the equilibrium signal in ensemble experiments is dominated by folded molecules. Here, we avoid this problem in two ways. One is the use of single-molecule Förster resonance energy transfer, which allows the separation of folded and unfolded subpopulations at equilibrium and provides information on long-range intramolecular distance distributions. From experiments with donor and acceptor chromophores placed at different positions within the chain, we find that the distance distributions in unfolded CspTm agree surprisingly well with a Gaussian chain not only at high concentrations of denaturant, where the polypeptide chain is expanded, but also at low denaturant concentrations, where the chain is collapsed. The second, complementary approach is synchrotron radiation circular dichroism spectroscopy of collapsed unfolded molecules transiently populated with a microfluidic device that enables rapid mixing. The results indicate a beta-structure content of the collapsed unfolded state of approximately 20% compared with the folded protein. This suggests that collapse can induce secondary structure in an unfolded state without interfering with long-range distance distributions characteristic of a random coil, which were previously found only for highly expanded unfolded proteins.

Key Role of Coulombic Interactions for the Folding Transition State of the Cold Shock Protein

The cold shock protein CspB shows a five-stranded beta-sheet structure, and it folds rapidly via a native-like transition state. A previous Phi value analysis showed that most of the residues with Phi values close to one reside in strand beta1, and two of them, Lys5 and Lys7 are partially exposed charged residues. To elucidate how coulombic interactions of these two residues contribute to the energetic organisation of the folding transition state we performed comparative folding experiments in the presence of an ionic denaturant (guanidinium chloride) and a non-ionic denaturant (urea) and a double-mutant analysis. Lys5 contributes 6.6 kJ mol(-1) to the stability of the transition state, and half of it originates from screenable coulombic interactions. Lys7 contributes 5.3 kJ mol(-1), and 3.4 kJ mol(-1) of it are screened by salt. In the folded protein Lys7 interacts with Asp25, and the screenable coulombic interaction between these two residues is fully formed in the transition state. This suggests that long-range coulombic interactions such as those originating from Lys5 and Lys7 of CspB can be important for organizing and stabilizing native-like structure early in protein folding.

Similarity and difference in the unfolding of thermophilic and mesophilic cold shock proteins studied by molecular dynamics simulations

Molecular dynamics simulations were performed to unfold a homologous pair of thermophilic and mesophilic cold shock proteins at high temperatures. The two proteins differ in just 11 of 66 residues and have very similar structures with a closed five-stranded antiparallel beta-barrel. A long flexible loop connects the N-terminal side of the barrel, formed by three strands (beta1-beta3), with the C-terminal side, formed by two strands (beta4-beta5). The two proteins were found to follow the same unfolding pathway, but with the thermophilic protein showing much slower unfolding. Unfolding started with the melting of C-terminal strands, leading to exposure of the hydrophobic core. Subsequent melting of beta3 and the beta-hairpin formed by the first two strands then resulted in unfolding of the whole protein. The slower unfolding of the thermophilic protein could be attributed to ion pair formation of Arg-3 with Glu-46, Glu-21, and the C-terminal. These ion pairs were also found to be important for the difference in folding stability between the pair of proteins. Thus electrostatic interactions appear to play similar roles in the difference in folding stability and kinetics between the pair of proteins.

Specificity of the initial collapse in the folding of the cold shock protein

The two-state folding reaction of the cold shock protein from Bacillus caldolyticus (Bc-Csp) is preceded by a rapid chain collapse. A fast shortening of intra-protein distances was revealed by Förster resonance energy transfer (FRET) measurements with protein variants that carried individual pairs of donor and acceptor chromophores at various positions along the polypeptide chain. Here we investigated the specificity of this rapid compaction. Energy transfer experiments that probed the stretching of strand beta2 and the close approach between the strands beta1 and beta2 revealed that the beta1-beta2 hairpin is barely formed in the collapsed form, although it is native-like in the folding transition state of Bc-Csp. The time course of the collapse could not be resolved by pressure or temperature jump experiments, indicating that the collapsed and extended forms are not separated by an energy barrier. The co-solute (NH4)2SO4 stabilizes both native Bc-Csp and the collapsed form, which suggests that the large hydrated SO4(2-) ions are excluded from the surface of the collapsed form in a similar fashion as they are excluded from folded Bc-Csp. Ethylene glycol increases the stability of proteins because it is excluded preferentially from the backbone, which is accessible in the unfolded state. The collapsed form of Bc-Csp resembles the unfolded form in its interaction with ethylene glycol, suggesting that in the collapsed form the backbone is still accessible to water and small molecules. Our results thus rule out that the collapsed form is a folding intermediate with native-like chain topology. It is better described as a mixture of compact conformations that belong to the unfolded state ensemble. However, some of its structural elements are reminiscent of the native protein.

Tryptophan rotamers as evidenced by X-ray, fluorescence lifetimes, and molecular dynamics modeling

We investigated the native-state dynamics of the Bacillus caldolyticus cold-shock protein mutant Bc-Csp L66E, using fluorescence and appropriate molecular dynamics methods. Two fluorescence lifetimes were found, the amplitudes of which agree very well with tryptophan rotamer populations, obtained from parallel tempering calculations. Rotamer lifetimes were predicted by transition-state theory from high-temperature simulations. Transition pathways were extracted from the transition rates between individual rotameric states. The molecular dynamics also reveal the loop fluctuations in the native state.

The contribution of the residues from the main hydrophobic core of ribonuclease A to its pressure-folding transition state

The role of hydrophobic interactions established by the residues that belong to the main hydrophobic core of ribonuclease A in its pressure-folding transition state was investigated using the Phi-value method. The folding kinetics was studied using pressure-jump techniques both in the pressurization and depressurization directions. The ratio between the folding activation volume and the reaction volume (beta p-value), which is an index of the compactness or degree of solvation of the transition state, was calculated. All the positions analyzed presented fractional Phi f-values, and the lowest were those corresponding to the most critical positions for the ribonuclease A stability. The structure of the transition state of the hydrophobic core of ribonuclease A, from the point of view of formed interactions, is a relatively, uniformly expanded form of the folded structure with a mean Phi f-value of 0.43. This places it halfway between the folded and unfolded states. On the other hand, for the variants, the average of beta p-values is 0.4, suggesting a transition state that is 40% native-like. Altogether the results suggest that the pressure-folding transition state of ribonuclease A looks like a collapsed globule with some secondary structure and a weakened hydrophobic core. A good correlation was found between the Phi f-values and the Deltabeta p-values. Although the nature of the transition state inferred from pressure-induced folding studies and the results of the protein engineering method have been reported to be consistent for other proteins, to the best of our knowledge this is the first direct comparison using a set of mutants.

Pressure as a tool to study protein-unfolding/refolding processes: The case of ribonuclease A

This paper gives an overview of the application of high-pressure to study the folding/unfolding processes of proteins using Ribonuclease A as a model protein. A particular focus is the study of pressure-equilibrium unfolding and folding kinetics using variants and the information obtained by comparing these with the wild-type enzyme.

Unfolding studies on soybean agglutinin and concanavalin a tetramers: a comparative account

The unfolding pathway of two very similar tetrameric legume lectins soybean agglutinin (SBA) and Concanavalin A (ConA) were determined using GdnCl-induced denaturation. Both proteins displayed a reversible two-state unfolding mechanism. The analysis of isothermal denaturation data provided values for conformational stability of the two proteins. It was found that the DeltaG of unfolding of SBA was much higher than ConA at all the temperatures at which the experiments were done. ConA had a T(g) 18 degrees C less than SBA. The higher conformational stability of SBA in comparison to ConA is largely due to substantial differences in their degrees of subunit interactions. Ionic interactions at the interface of the two proteins especially at the noncanonical interface seem to play a significant role in the observed stability differences between these two proteins. Furthermore, SBA is a glycoprotein with a GlcNac2Man9 chain attached to Asn-75 of each subunit. The sugar chain in SBA lies at the noncanonical interface of the protein, and it is found to interact with the amino acid residues in the adjacent noncanonical interface. These interactions further stabilize SBA with respect to ConA, which is not glycosylated.

The N-terminal to C-terminal motif in protein folding and function

Essentially all proteins known to fold kinetically in a two-state manner have their N- and C-terminal secondary structural elements in contact, and the terminal elements often dock as part of the experimentally measurable initial folding step. Conversely, all N-C no-contact proteins studied so far fold by non-two-state kinetics. By comparison, about half of the single domain proteins in the Protein Data Bank have their N- and C-terminal elements in contact, more than expected on a random probability basis but not nearly enough to account for the bias in protein folding. Possible reasons for this bias relate to the mechanisms for initial protein folding, native state stability, and final turnover.

The Folding Transition State of the Cold Shock Protein is Strongly Polarized

The cold shock protein CspB from Bacillus subtilis consists of a three-stranded (beta1-beta3) and a two stranded (beta4-beta5) sheet, which form a closed beta barrel structure. CspB folds and unfolds rapidly in a two-state reaction, and the unfolded and the folded molecules interconvert with a time constant of 30 ms at the midpoint of the urea-induced transition (at 25 degrees C). The transition state of folding is native-like, as judged by the Tanford betaT value of > or =0.9. By using a mutational approach and Phi value analysis, we find that the folding transition state of CspB is energetically polarized. Despite the high betaT value, most Phi values are low. Values close to 1 were found for only a few residues, particularly in strand beta1 (Lys5, Val6, Lys7, Asn10). The interactions of the Asn10 side-chain with the backbone at positions 12 and 13 define the turn that connects the strands beta1 and beta2. Lys5 and Val6 in beta1 interact with residues in beta4, and their high Phi values indicate that an energetic linkage between beta1 and beta4 and thus between the two sheets exists already in the transition state. We compared our experimental Phi values with theoretical predictions of the folding pathway of cold shock proteins. Several of them suggest that the entire first sheet is formed in the transition state, and some identify the beta1-beta4 pairing as a crucial step in folding. Alternative paths that involve formation of the second sheet and beta3-beta5 pairing reactions were, however, suggested as well. The calculations gave coarse-grained pictures that are limited in resolution to the two sheets of CspB or to the elements of secondary structure. They did not identify the key residues with the high Phi values within these structural elements.

Salt-bridges can stabilize but do not accelerate the folding of the homodimeric coiled-coil peptide GCN4-p1

Double mutant cycle analysis was employed to ascertain the role of intra- and interchain salt-bridges in the folding and stability of the dimeric coiled-coil peptide, GCN4-p1, the 33-residue leucine zipper domain of the transcriptional activator GCN4. Equilibrium circular dichroism studies of the urea-induced unfolding reaction at neutral pH revealed that both types of ionic interactions, localized primarily in the N-terminal portion of the molecule, enhance the stability of the native coiled-coil. By contrast, comparable stopped-flow circular dichroism studies indicate that the salt-bridge interactions, with one possible exception, are not well formed in the transition state for folding. Although the E22Q/R25A double mutant failed to fold, fragmentation studies suggest that the E22/R25 intramolecular salt-bridge may play a critical role in stabilizing C-terminal nascent helices that drive the association reaction. The remaining salt-bridges appear to stabilize the parallel-stranded coiled-coil architecture of GCN4-p1 only after the peptide traverses the rate-limiting, dimeric transition state.

Rapid collapse precedes the fast two-state folding of the cold shock protein

The cold shock protein Bc-Csp folds very rapidly in a reaction that is well described by a kinetic two-state mechanism without intermediates. We measured the shortening of six intra-protein distances during folding by Förster resonance energy transfer (FRET) in combination with stopped-flow experiments. Single tryptophan residues were engineered into the protein as the donors, and single 5-(((acetylamino)ethyl)amino)naphthalene-1-sulfonate (AEDANS) residues were placed as the acceptors at solvent-exposed sites of Bc-Csp. Their R0 value of about 22 A was well suited for following distance changes during the folding of this protein with a high sensitivity. The mutagenesis and the labeling did not alter the refolding kinetics. The changes in energy transfer during folding were monitored by both donor and acceptor emission and reciprocal effects were found. In two cases the donor-acceptor distances were similar in the unfolded and the folded state and, as a consequence, the kinetic changes in energy transfer upon folding were very small. For four donor/acceptor pairs we found that > or =50% of the increase in energy transfer upon folding occurred prior to the rate-limiting step of folding. This reveals that about half of the shortening of the intra-molecular distances upon folding has occurred already before the rate-limiting step and suggests that the fast two-state folding reaction of Bc-Csp is preceded by a very rapid collapse.

Unfolding of the cold shock protein studied with biased molecular dynamics

The cold shock protein from Bacillus caldolyticus is a small beta-barrel protein that folds in a two-state mechanism. For the native protein and for several mutants, a wealth of experimental data are available on stability and folding, so that it is an optimal system to study this process. We compare data from unfolding simulations (trajectories of 5 and up to 12 ns) obtained with a bias potential at room temperature and from unbiased thermal unfolding simulations with experimental data. The unfolding patterns derived from the trajectories starting from different native-like conformations and subject to different unfolding conditions agree. The transition state found in the simulations of unfolding is close to the native structure in agreement with experiment. Moreover, a lower value of the free energy barrier of unfolding was found for the mutant R3E than for the mutant E46A and the native protein, as indicated by experimental data. The first unfolding event involves the three-stranded beta-sheet whose decomposition corresponds to the transition state. In contrast to conclusions drawn from experiments, we found that the two-stranded beta-strand forms the most stable substructure, which decomposes very late in the unfolding process. However, assuming that this structure forms very early in the folding process, our findings would not contradict the experiments but require a different interpretation of them.

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.

Folding rates and low-entropy-loss routes of two-state proteins

We develop a simple model for computing the rates and routes of folding of two-state proteins from the contact maps of their native structures. The model is based on the graph-theoretical concept of effective contact order (ECO). The model predicts that proteins fold by "zipping up" in a sequence of small-loop-closure events, depending on the native chain fold. Using a simple equation, with a few physical rate parameters, we obtain a good correlation with the folding rates of 24 two-state folding proteins. The model rationalizes data from Phi-value analysis that have been interpreted in terms of delocalized or polarized transition states. This model indicates how much of protein folding may take place in parallel, not along a single reaction coordinate or with a single transition state.

Hammond behavior versus ground state effects in protein folding: evidence for narrow free energy barriers and residual structure in unfolded states

Apparent transition state movement upon mutation or changes in solvent conditions is frequently observed in protein folding and is often interpreted in terms of Hammond behavior. This led to the conclusion that barrier regions in protein folding are broad maxima on the free energy landscape. Here, we use the concept of self-interaction and cross-interaction parameters to test experimental data of 21 well-characterized proteins for Hammond behavior. This allows us to characterize the origin of transition state movements along different reaction coordinates. Only one of the 21 proteins shows a small but coherent transition state movement in agreement with the Hammond postulate. In most proteins the structure of the transition state is insensitive to changes in protein stability. The apparent change in the position of the transition state upon mutation, which is frequently observed in phi-value analysis, is in most cases due to ground-state effects caused by structural changes in the unfolded state. This argues for significant residual structure in unfolded polypeptide chains of many proteins. Disruption of these residual interactions by mutation often leads to decreased folding rates, which implies that these interactions are still present in the transition state. The failure to detect Hammond behavior shows that the free energy barriers encountered by a folding polypeptide chain are generally rather narrow and robust maxima for all experimentally explorable reaction coordinates.

Early events in protein folding

Recent advances have significantly increased the time and spectroscopic resolution of protein folding experiments. We can now study the timescale and nature of polypeptide collapse, and how this correlates with secondary and tertiary structure formation. Studies on ultrafast folding proteins and peptides provide experimental benchmarks on a timescale that overlaps directly with that of molecular dynamics simulations. This makes possible direct tests of both simulations and current models of protein folding.

Water contributes actively to the rapid crossing of a protein unfolding barrier

The cold-shock protein CspB folds rapidly in a N <= => U two-state reaction via a transition state that is about 90% native in its interactions with denaturants and water. This suggested that the energy barrier to unfolding is overcome by processes occurring in the protein itself, rather than in the solvent. Nevertheless, CspB unfolding depends on the solvent viscosity. We determined the activation volumes of unfolding and refolding by pressure-jump and high-pressure stopped-flow techniques in the presence of various denaturants. The results obtained by these methods agree well. The activation volume of unfolding is positive (Delta V(++)(NU)=16(+/-4) ml/mol) and virtually independent of the nature and the concentration of the denaturant. We suggest that in the transition state the protein is expanded and water molecules start to invade the hydrophobic core. They have, however, not yet established favorable interactions to compensate for the loss of intra-protein interactions. The activation volume of refolding is positive as well (Delta V(++)(NU)=53(+/-6) ml/mol) and, above 3 M urea, independent of the concentration of the denaturant. At low concentrations of urea or guanidinium thiocyanate, Delta V(++)(UN) decreases significantly, suggesting that compact unfolded forms become populated under these conditions.