Energetic coupling between native-state prolyl isomerization and conformational protein folding

In folded proteins, prolyl peptide bonds are usually thought to be either trans or cis because only one of the isomers can be accommodated in the native folded protein. For the N-terminal domain of the gene-3 protein of the filamentous phage fd (N2 domain), Pro161 resides at the tip of a beta hairpin and was found to be cis in the crystal structure of this protein. Here we show that Pro161 exists in both the cis and the trans conformations in the folded form of the N2 domain. We investigated how conformational folding and prolyl isomerization are coupled in the unfolding and refolding of N2 domain. A combination of single-mixing and double-mixing unfolding and refolding experiments showed that, in unfolded N2 domain, 7% of the molecules contain a cis-Pro161 and 93% of the molecules contain a trans-Pro161. During refolding, the fraction of molecules with a cis-Pro161 increases to 85%. This implies that 10.3 kJ mol(-1) of the folding free energy was used to drive this 75-fold change in the Pro161 cis/trans equilibrium constant during folding. The stabilities of the forms with the cis and the trans isomers of Pro161 and their folding kinetics could be determined separately because their conformational folding is much faster than the prolyl isomerization reactions in the native and the unfolded proteins. The energetic coupling between conformational folding and Pro161 isomerization is already fully established in the transition state of folding, and the two isomeric forms are thus truly native forms. The folding kinetics are well described by a four-species box model, in which the N2 molecules with either isomer of Pro161 can fold to the native state and in which cis/trans isomerization occurs in both the unfolded and the folded proteins.     

Relationship between the native-state hydrogen exchange and folding pathways of a four-helix bundle protein

The hydrogen exchange behavior of a four-helix bundle protein in low concentrations of denaturant reveals some partially unfolded forms that are significantly more stable than the fully unfolded state. Kinetic folding of the protein, however, is apparently two-state in the absence of the accumulation of early folding intermediates. The partially unfolded forms are either as folded as or more folded than the rate-limiting transition state and appear to represent the major intermediates in a folding and unfolding reaction. These results are consistent with the suggestion that partially unfolded intermediates may form after the rate-limiting step for small proteins with apparent two-state folding kinetics.     

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.     

Transient intermediary states with high and low folding probabilities in the apparent two-state folding equilibrium of ACBP at low pH

Measurements of the stability as a function of pH for the acyl-coenzyme A binding protein (ACBP) has shown a significant difference in the pH transition midpoint measured by NMR spectroscopy at pH 3.12 and the transition midpoint measured at pH 2.92 and 2.97 by circular dichroism and by fluorescence spectroscopy, respectively. A similar behavior has not been observed in other proteins. It is suggested that these differences arise because the population of the unfolded molecules still contains significant amounts of native like secondary and tertiary structure. NMR spectroscopy measures the concentration of the two components of the folding unfolding equilibrium individually, whereas circular dichroism and fluorescence measure the concentration of the conformations of the light-absorbing chromophores present in both the folded and the unfolded molecules. In the narrow pH range, nascent structure can be detected as the average amount of secondary structure per unfolded molecule and hydrophobic interactions in the population of unfolded molecules. These structures are not observable immediately by NMR spectroscopy; however, a chemical shift analysis of the peptide backbone (13)C chemical shift indicates strongly the existence of short-lived and transient helical structures at pH 2.3. Magnetization transfer studies have been applied to study the equilibrium between folded and unfolded ACBP near the pH transition point measured by NMR. This study has shown that there are two categories of subpopulations in the population of unfolded ACBP. One for which magnetization can be transferred to the folded form during the folding process, and one for which transfer is not observed. The molecules of the latter population of unfolded protein apparently, do not fold within the time-frame of the magnetization transfer experiment. This result suggests the existence of a subpopulation of the acid-unfolded protein molecules with a high propensity for folding. It is suggested that in this subpopulation, a particular set of native like interactions in the peptide backbone and between side-chains in the peptide chain have to be formed.     

Reversible denaturation of the gene V protein of bacteriophage f1

The guanidine hydrochloride (GuHCl)-induced denaturation of the gene V protein of bacteriophage f1 has been studied, using the chemical reactivity of a cysteine residue that is buried in the folded protein and the circular dichroism (CD) at 211 and 229 nm as measures of the fraction of polypeptide chains in the folded form. It is found that this dimeric protein unfolds in a single cooperative transition from a folded dimer to two unfolded monomers. A folded, monomeric form of the gene V protein was not detected at equilibrium. The kinetics of unfolding of the gene V protein in 3 M GuHCl and the refolding in 2 M GuHCl are also consistent with a transition between a folded dimer and two unfolded monomers. The GuHCl concentration dependence of the rates of folding and unfolding suggests that the transition state for folding is near the folded conformation.     

NMR and CD spectroscopy show that imino acid restriction of the unfolded state leads to efficient folding

Protein folding is determined by molecular features in the unfolded state, as well as the native folded structure. In the unfolded state, imino acids both restrict conformational space and present cis-trans isomerization barriers to folding. Because of its high proline and hydroxyproline content, the collagen triple-helix offers an opportunity to characterize the impact of imino acids on the unfolded state and folding kinetics. Here, NMR and CD spectroscopy are used to characterize the role of imino acids in a triple-helical peptide, T1-892, which contains an 18-residue sequence from type I collagen and a C-terminal (Gly-Pro-Hyp)(4) domain. The replacement of Pro or Hyp by an Ala in the (Gly-Pro-Hyp)(4) region significantly decreases the folding rate at low but not high concentrations, consistent with less efficient nucleation. To understand the molecular basis of the decreased folding rate, changes in the unfolded as well as the folded states of the peptides were characterized. While the trimer states of the peptides are all similar, NMR dynamics studies show monomers with all trans (Gly-Pro-Hyp)(4) are less flexible than monomers containing Pro --> Ala or Hyp --> Ala substitutions. Nucleation requires all trans bonds in the (Gly-Pro-Hyp)(4) domain and the constrained monomer state of the all trans nucleation domain in T1-892 increases its competency to initiate triple-helix formation and illustrates the impact of the unfolded state on folding kinetics.     

Dependence of Protein Folding Stability and Dynamics on the Density and Composition of Macromolecular Crowders

We investigate the effect of macromolecular crowding on protein folding, using purely repulsive crowding particles and a self-organizing polymer model of protein folding. We find that the variation in folding stability with crowder size for typical α-, β-, and α/β-proteins is well described by an adaptation of the scaled particle theory. The native state, the transition state, and the unfolded protein are treated as effective hard spheres, with the folded and transition state radii independent of the size and concentration of the crowders. Remarkably, we find that, as the effective unfolded state radius is very weakly dependent on the crowder concentration, it can also be approximated by a single size. The same model predicts the effect of crowding on the folding barrier and therefore refolding rates with no adjustable parameters. A simple extension of the scaled-particle theory model, assuming additivity, can also describe the behavior of mixtures of crowding particles.     

Recombination of S-peptide with S-protein during folding of ribonuclease S. II. Kinetic characterization of a stable folding intermediate shown by S-protein at pH 1.7.

At pH 1.7 S-peptide dissociates from S-protein but S-protein remains partly folded below 30 °C. A folded form of S-protein, labeled I₃, is detected and measured by its ability to combine rapidly with S-peptide at pH 6.8 and then to form native ribonuclease S. The second-order combination reaction (k = 0.7 × 10⁶m^?1s^?1 at 20 °C) can be monitored either by tyrosine absorbance or fluorescence emission; the subsequent first-order folding reaction (half-time, 68 ms; 20 °C) is monitored by 2′CMP ² binding. Combination with S-peptide and folding to form native RNAase S is considerably slower for both classes of unfolded S-protein (see preceding paper).I₃ shows a thermal folding transition at pH 1.7: it is completely unfolded above 32 °C and reaches a limiting low-temperature value of 65% below 10 °C. The 35% S-protein remaining at 10 °C is unfolded as judged by its refolding behavior in forming native RNAase S at pH 6.8. The folding transition of S-protein at pH 1.7 is a broad, multi-state transition. This is shown both by the large fraction of unfolded S-protein remaining at low temperatures and by the large differences between the folding transition curves monitored by I₃ and by tyrosine absorbance.The fact that S-protein remains partly folded after dissociation of S-peptide at pH 1.7 but not at pH 6.8 may be explained by two earlier observations. (1) Native RNAase A is stable in the temperature range of the S-protein folding transition at pH 1.7, and (2) the binding constant of S-protein for S-peptide falls steadily as the pH is lowered, by more than four orders of magnitude between pH 8.3 and pH 2.7, at 0 °C. The following explanation is suggested for why folding intermediates are observed easily in the transition of S-protein but not of RNAase A. The S-protein transition is shifted to lower temperatures, where folding intermediates should be more stable: consequently, intermediates in the folding of RNAase A which do not involve the S-peptide moiety and which are populated to almost detectable levels can be observed at the lower temperatures of the S-protein transition.     

An experimental investigation of conformational fluctuations in proteins G and L

The B1 domains of streptococcal proteins G and L are structurally similar, but they have different sequences and they fold differently. We have measured their NMR spectra at variable temperature using a range of concentrations of denaturant. Many residues have curved amide proton temperature dependence, indicating that they significantly populate alternative, locally unfolded conformations. The results, therefore, provide a view of the locations of low-lying, locally unfolded conformations. They indicate approximately 4-6 local minima for each protein, all within ca. 2.5 kcal/mol of the native state, implying a locally rough energy landscape. Comparison with folding data for these proteins shows that folding involves most molecules traversing a similar path, once a transition state containing a beta hairpin has been formed, thereby defining a well-populated pathway down the folding funnel. The hairpin that directs the folding pathway differs for the two proteins and remains the most stable part of the folded protein.     

Characterization of the hydrodynamic properties of the folding transition state of an SH3 domain by magnetization transfer NMR spectroscopy

Protein folding kinetic data have been obtained for the marginally stable N-terminal Src homology 3 domain of the Drosophila protein drk (drkN SH3) in an investigation of the hydrodynamic properties of its folding transition state. Due to the presence of NMR resonances of both folded and unfolded states at equilibrium, kinetic data can be derived from NMR magnetization transfer techniques under equilibrium conditions. Kinetic analysis as a function of urea (less than approximately 1 M) and glycerol enables determination of alpha values, measures of the energetic sensitivity of the transition state to the perturbation relative to the end states of the protein folding reaction (the folded and unfolded states). Both end states have previously been studied experimentally by NMR spectroscopic and other biophysical methods in great detail and under nondenaturing conditions. Combining these results with the kinetic folding data obtained here, we can characterize the folding transition state without requiring empirical models for the unfolded state structure. We are thus able to give a reliable measure of the solvent-accessible surface area of the transition state of the drkN SH3 domain (4730 +/- 360 A(2)) based on urea titration data. Glycerol titration data give similar results and additionally demonstrate that folding of this SH3 domain is dependent on solvent viscosity, which is indicative of at least partial hydration of the transition state. Because SH3 domains appear to fold by a common folding mechanism, the data presented here provide valuable insight into the transition states of the drkN and other SH3 domains.     

Acid- and pressure-induced (un)folding of yeast glutathione reductase: competition between protein oligomerization and aggregation

Glutathione reductase (GR) is a homodimeric flavoenzyme involved in cellular defense against oxidative stress. In the present study, we have used a combination of acidic pH and hydrostatic pressure to investigate the (un)folding transition of yeast GR. Our results indicate that at pH 2 a distinct partially folded state is stabilized, as judged by intrinsic fluorescence, bis ANS binding and circular dichroism (CD) analysis. Further characterization of this partially folded state by size exclusion chromatography revealed that it corresponds to expanded GR monomers. CD analysis at pH 2 showed a significant loss of secondary structure. The partially folded GR monomers stabilized at pH 2 were fully and reversibly unfolded using hydrostatic pressure (up to 3.5kbar) as a thermodynamic perturbant. By contrast, return to physiological pH after exposure to acidic pH led to a competing reaction between refolding dimerization and aggregation of GR. These results support the notion that a partially folded intermediate state is not only critical for folding of GR but also appears to be a seed for protein aggregation.     

Kinetic analysis of folding and unfolding the 56 amino acid IgG-binding domain of streptococcal protein G.

The 56 amino acid B domain of protein G (GB) is a stable globular folding unit with no disulfide cross-links. The physical properties of GB offer extraordinary flexibility for evaluating the energetics of the folding reaction. The protein is monomeric and very soluble in both folded and unfolded forms. The folding reaction has been previously examined by differential scanning calorimetry (Alexander et al., 1992) and found to exhibit two-state unfolding behavior over a wide pH range with an unfolding transition near 90 degrees C (GB1) at neutral pH. Here, the kinetics of folding and unfolding two naturally occurring versions of GB have been measured using stopped-flow mixing methods and analyzed according to transition-state theory. GB contains no prolines, and the kinetics of folding and unfolding can be fit to a single, first-order rate constant over the temperature range of 5-35 degrees C. The major thermodynamic changes going from the unfolded state to the transition state are (1) a large decrease in heat capacity (delta Cp), indicating that the transition state is compact and solvent inaccessible relative to the unfolded state; (2) a large loss of entropy; and (3) a small increase in enthalpy. The most surprising feature of the folding of GB compared to that of previously studied proteins is that its folding approximates a rapid diffusion controlled process with little increase in enthalpy going from the unfolded to the transition state.     

How, when and why proteins collapse: the relation to folding

Unfolded proteins under strongly denaturing conditions are highly expanded. However, when the conditions are more close to native, an unfolded protein may collapse to a compact globular structure distinct from the folded state. This transition is akin to the coil-globule transition of homopolymers. Single-molecule FRET experiments have been particularly conducive in revealing the collapsed state under conditions of coexistence with the folded state. The collapse can be even more readily observed in natively unfolded proteins. Time-resolved studies, using FRET and small-angle scattering, have shown that the collapse transition is a very fast event, probably occurring on the submicrosecond time scale. The forces driving collapse are likely to involve both hydrophobic and backbone interactions. The loss of configurational entropy during collapse makes the unfolded state less stable compared to the folded state, thus facilitating folding.     

Salt-Dependent Folding Energy Landscape of RNA Three-Way Junction

RNAs are highly negatively charged chain molecules. Metal ions play a crucial role in RNA folding stability and conformational changes. In this work, we employ the recently developed tightly bound ion (TBI) model, which accounts for the correlation between ions and the fluctuation of ion distributions, to investigate the ion-dependent free energy landscape for the three-way RNA junction in a 16S rRNA domain. The predicted electrostatic free energy landscape suggests that 1), ion-mediated electrostatic interactions cause an ensemble of unfolded conformations narrowly populated around the maximally extended structure; and 2), Mg²⁺ ion-induced correlation effects help bring the helices to the folded state. Nonelectrostatic interactions, such as noncanonical interactions within the junctions and between junctions and helix stems, might further limit the conformational diversity of the unfolded state, resulting in a more ordered unfolded state than the one predicted from the electrostatic effect. Moreover, the folded state is predominantly stabilized by the coaxial stacking force. The TBI-predicted folding stability agrees well with the experimental measurements for the different Na⁺ and Mg²⁺ ion concentrations. For Mg²⁺ solutions, the TBI model, which accounts for the Mg²⁺ ion correlation effect, gives more improved predictions than the Poisson-Boltzmann theory, which tends to underestimate the role of Mg²⁺ in stabilizing the folded structure. Detailed control tests indicate that the dominant ion correlation effect comes from the charge-charge Coulombic correlation rather than the size (excluded volume) correlation between the ions. Furthermore, the model gives quantitative predictions for the ion size effect in the folding energy landscape and folding cooperativity.     

Two-dimensional 1H nuclear magnetic resonance study of the (5-55) single-disulphide folding intermediate of bovine pancreatic trypsin inhibitor

An analogue of the bovine pancreatic trypsin inhibitor (BPTI) folding intermediate that contains only the disulphide bond between Cys5 and Cys55 has been prepared in Escherichia coli by protein engineering methods, with the other four Cys residues replaced by Ser. Two-dimensional 1H nuclear magnetic resonance studies of the analogue have resulted in essentially complete resonance assignments of the folded form of the protein. The folded protein has a compact conformation that is structurally very similar to that of native BPTI, although there are subtle differences and the folded conformation is not very stable. Approximately half of the protein molecules are unfolded at 3 degrees C, and this proportion increases at higher temperatures. The folded and unfolded conformations are in slow exchange. The conformational properties of the analogue can explain many aspects of the kinetic role that the normal (5-55) intermediate plays in the folding of BPTI.     

Molecular basis of cooperativity in protein folding IV. Core: A general cooperative folding model

The cooperative nature of the protein folding process is independent of the characteristic fold and the specific secondary structure attributes of a globular protein. A general folding/unfolding model should, therefore, be based upon structural features that transcend the peculiarities of α-helices, β-sheets, and other structural motifs found in proteins. The studies presented in this paper suggest that a single structural characteristic common to all globular proteins is essential for cooperative folding. The formation of a partly folded state from the native state results in the exposure to solvent of two distinct regions: (1) the portions of the protein that are unfolded; and (2) the “complementary surfaces,” located in the regions of the protein that remain folded. The cooperative character of the folding/unfolding transition is determined largely by the energetics of exposing complementary surface regions to the solvent. By definition, complementary regions are present only in partly folded states; they are absent from the native and unfolded states. An unfavorable free energy lowers the probability of partly folded states and increases the cooperativity of the transition. In this paper we present a mathematical formulation of this behavior and develop a general cooperative folding/unfolding model, termed the “complementary region” (CORE) model. This model successfully reproduces the main properties of folding/unfolding transitions without limiting the number of partly folded states accessible to the protein, thereby permitting a systematic examination of the structural and solvent conditions under which intermediates become populated. It is shown that the CORE model predicts two-state folding/unfolding behavior, even though the two-state character is not assumed in the model. © 1993 Wiley-Liss, Inc.     

Cunning simplicity of a hierarchical folding

A hierarchic scheme of protein folding does not solve the Levinthal paradox since it cannot provide a simultaneous explanation for major features observed for protein folding: (i) folding within non-astronomical time, (ii) independence of the native structure on large variations in the folding rates of given protein under different conditions, and (iii) co-existence, in a visible quantity, of only the native and the unfolded molecules during folding of moderate size (single-domain) proteins. On the contrary, a nucleation mechanism can account for all these major features simultaneously and resolves the Levinthal paradox.     

A framework for describing topological frustration in models of protein folding

In a natively folded protein of moderate or larger size, the protein backbone may weave through itself in complex ways, raising questions about what sequence of events might have to occur in order for the protein to reach its native configuration from the unfolded state. A mathematical framework is presented here for describing the notion of a topological folding barrier, which occurs when a protein chain must pass through a hole or opening, formed by other regions of the protein structure. Different folding pathways encounter different numbers of such barriers and therefore different degrees of frustration. A dynamic programming algorithm finds the optimal theoretical folding path and minimal degree of frustration for a protein based on its natively folded configuration. Calculations over a database of protein structures provide insights into questions such as whether the path of minimal frustration might tend to favor folding from one or from many sites of folding nucleation, or whether proteins favor folding around the N terminus, thereby providing support for the hypothesis that proteins fold co-translationally. The computational methods are applied to a multi-disulfide bonded protein, with computational findings that are consistent with the experimentally observed folding pathway. Attention is drawn to certain complex protein folds for which the computational method suggests there may be a preferred site of nucleation or where folding is likely to proceed through a relatively well-defined pathway or intermediate. The computational analyses lead to testable models for protein folding.     

Folding and unfolding of gammaTIM monomers and dimers

Kinetic simulations of the folding and unfolding of triosephosphate isomerase (TIM) from yeast were conducted using a single monomer gammaTIM polypeptide chain that folds as a monomer and two gammaTIM chains that fold to the native dimer structure. The basic protein model used was a minimalist Gō model using the native structure to determine attractive energies in the protein chain. For each simulation type--monomer unfolding, monomer refolding, dimer unfolding, and dimer refolding--thirty simulations were conducted, successfully capturing each reaction in full. Analysis of the simulations demonstrates four main conclusions. First, all four simulation types have a similar "folding order", i.e., they have similar structures in intermediate stages of folding between the unfolded and folded state. Second, despite this similarity, different intermediate stages are more or less populated in the four different simulations, with 1), no intermediates populated in monomer unfolding; 2), two intermediates populated with beta(2)-beta(4) and beta(1)-beta(5) regions folded in monomer refolding; 3), two intermediates populated with beta(2)-beta(3) and beta(2)-beta(4) regions folded in dimer unfolding; and 4), two intermediates populated with beta(1)-beta(5) and beta(1)-beta(5) + beta(6) + beta(7) + beta(8) regions folded in dimer refolding. Third, simulations demonstrate that dimer binding and unbinding can occur early in the folding process before complete monomer-chain folding. Fourth, excellent agreement is found between the simulations and MPAX (misincorporation proton alkyl exchange) experiments. In total, this agreement demonstrates that the computational Gō model is accurate for gammaTIM and that the energy landscape of gammaTIM appears funneled to 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.