The unfolded state of NTL9 is compact in the absence of denaturant

Interest in the unfolded state of proteins has grown with the realization that this state can have considerable structure in the absence of denaturants. Natively unfolded proteins, mutations that unfold proteins under native conditions, and changes in pH that induce unfolding are attractive models for the unfolded state in the absence of denaturant. The unfolded state of the N-terminal domain of ribosomal protein L9 (NTL9) was previously shown to contain significant non-native electrostatic interactions [Cho, J. H., Sato, S., and Raleigh, D. P. (2004) J. Mol. Biol. 338, 827-837]. NTL9 has a mixed alpha-beta structure and folds via a two-state mechanism. We have generated a model of the unfolded state of NTL9 in the absence of denaturant by substitution of an alanine for phenylalanine 5 located in the core of this protein. The CD spectrum of the variant, denoted as F5A, exhibits significantly less structure than the wild type; however, the mean residue ellipticity of F5A at 222 nm (-8200 deg cm(2) dmol(-)(1)) is considerably larger than expected for a fully unfolded protein, indicating that residual secondary structure is populated. F5A also has more residual structure than the urea-unfolded wild type. The stability of F5A is estimated to be at least 1 kcal/mol unfavorable, showing that the unfolded state is populated to 84% or more. NMR pulsed-field gradient measurements yield a hydrodynamic radius of 16.1 A for wild-type NTL9 and 20.8 A for the F5A variant in native buffer. The physiologically relevant unfolded state of wild-type NTL9 is likely to be even more compact than F5A since the mutation should reduce the level of hydrophobic clustering in the unfolded state in the absence of denaturant. The hydrodynamic radius of F5A increases to 25.9 A in 8 M urea, and a value of 23.5 A is obtained for the wild type under similar conditions. The results show that the unfolded state of F5A in the absence of denaturant is more compact and contains more structure than the urea-unfolded form.     

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.     

The folding pathway of barnase: the rate-limiting transition state and a hidden intermediate under native conditions

The nature of the rate-limiting transition state at zero denaturant (TS(1)) and whether there are hidden intermediates are the two major unsolved problems in defining the folding pathway of barnase. In earlier studies, it was shown that TS(1) has small phi values throughout the structure of the protein, suggesting that the transition state has either a defined partially folded secondary structure with all side chains significantly exposed or numerous different partially unfolded structures with similar stability. To distinguish the two possibilities, we studied the effect of Gly mutations on the folding rate of barnase to investigate the secondary structure formation in the transition state. Two mutations in the same region of a beta-strand decreased the folding rate by 20- and 50-fold, respectively, suggesting that the secondary structures in this region are dominantly formed in the rate-limiting transition state. We also performed native-state hydrogen exchange experiments on barnase at pD 5.0 and 25 degrees C and identified a partially unfolded state. The structure of the intermediate was investigated using protein engineering and NMR. The results suggest that the intermediate has an omega loop unfolded. This intermediate is more folded than the rate-limiting transition state previously characterized at high denaturant concentrations (TS(2)). Therefore, it exists after TS(2) in folding. Consistent with this conclusion, the intermediate folds with the same rate and denaturant dependence as the wild-type protein, but unfolds faster with less dependence on the denaturant concentration. These and other results in the literature suggest that barnase folds through partially unfolded intermediates that exist after the rate-limiting step. Such folding behavior is similar to those of cytochrome c and Rd-apocyt b(562). Together, we suggest that other small apparently two-state proteins may also fold through hidden intermediates.     

Comparison between denaturant- and temperature-induced unfolding pathways of protein: a lattice Monte Carlo simulation

Denaturant-induced unfolding of protein is simulated by using a Monte Carlo simulation with a lattice model for protein and denaturant. Following the binding theory for denaturant-induced unfolding, the denaturant molecules are modeled to interact with protein by nearest-neighbor interactions. By analyzing the conformational states on the unfolding pathway of protein, the denaturant-induced unfolding pathway is compared with the temperature-induced unfolding pathway under the same condition; that is, the free energies of unfolding under two different pathways are equal. The two unfoldings show markedly different conformational distributions in unfolded states. From the calculation of the free energy of protein as a function of the number fraction (Q0) of native contacts relative to the total number of contacts, it is found that the free energy of the largely unfolded state corresponding to low Q0 (0.1 < Q0 < 0.5) under temperature-induced unfolding is lower than that under denaturant-induced unfolding, whereas the free energy of the unfolded state close to the native state (Q0 > 0.5) is lower in denaturant-induced unfolding than in temperature-induced unfolding. A comparison of two unfolding pathways reveals that the denaturant-induced unfolding shows a wider conformational distribution than the temperature-induced unfolding, while the temperature-induced unfolding shows a more compact unfolded state than the denaturant-induced unfolding especially in the low Q0 region (0.1 < Q0 < 0.5).     

Thermodynamics and Kinetics of Single-Chain Monellin Folding with Structural Insights into Specific Collapse in the Denatured State Ensemble

Proteins, which behave as random coils in high denaturant concentrations undergo collapse transition similar to polymers on denaturant dilution. We study collapse in the denatured ensemble of single-chain monellin (MNEI) using a coarse-grained protein model and molecular dynamics simulations. The model is validated by quantitatively comparing the computed guanidinium chloride and pH-dependent thermodynamic properties of MNEI folding with the experiments. The computed properties such as the fraction of the protein in the folded state and radius of gyration (R_g) as function of [GuHCl] are in good agreement with the experiments. The folded state of MNEI is destabilized with an increase in pH due to the deprotonation of the residues Glu24 and Cys42. On decreasing [GuHCl], the protein in the unfolded ensemble showed specific compaction. The R_g of the protein decreased steadily with [GuHCl] dilution due to increase in the number of native contacts in all the secondary structural elements present in the protein. MNEI folding kinetics is complex with multiple folding pathways and transiently stable intermediates are populated in these pathways. In strong stabilizing conditions, the protein in the unfolded ensemble showed transition to a more compact unfolded state where R_g decreased by ≈ 17% due to the formation of specific native contacts in the protein. The intermediate populated in the dominant MNEI folding pathway satisfies the structural features of the dry molten globule inferred from experiments.     

An effective solvent theory connecting the underlying mechanisms of osmolytes and denaturants for protein stability

An all-atom Gō model of Trp-cage protein is simulated using discontinuous molecular dynamics in an explicit minimal solvent, using a single, contact-based interaction energy between protein and solvent particles. An effective denaturant or osmolyte solution can be constructed by making the interaction energy attractive or repulsive. A statistical mechanical equivalence is demonstrated between this effective solvent model and models in which proteins are immersed in solutions consisting of water and osmolytes or denaturants. Analysis of these studies yields the following conclusions: 1), Osmolytes impart extra stability to the protein by reducing the entropy of the unfolded state. 2), Unfolded states in the presence of osmolyte are more collapsed than in water. 3), The folding transition in osmolyte solutions tends to be less cooperative than in water, as determined by the ratio of van 't Hoff to calorimetric enthalpy changes. The decrease in cooperativity arises from an increase in native structure in the unfolded state, and thus a lower thermodynamic barrier at the transition midpoint. 4), Weak denaturants were observed to destabilize small proteins not by lowering the unfolded enthalpy, but primarily by swelling the unfolded state and raising its entropy. However, adding a strong denaturant destabilizes proteins enthalpically. 5), The folding transition in denaturant-containing solutions is more cooperative than in water. 6), Transfer to a concentrated osmolyte solution with purely hard-sphere steric repulsion significantly stabilizes the protein, due to excluded volume interactions not present in the canonical Tanford transfer model. 7), Although a solution with hard-sphere interactions adds a solvation barrier to native contacts, the folding is nevertheless less cooperative for reasons 1–3 above, because a hard-sphere solvent acts as a protecting osmolyte.     

Hydrophobic Core Mutations Associated with Cataract Development in Mice Destabilize Human ��D-Crystallin

The human eye lens is composed of fiber cells packed with crystallins up to 450 mg/ml. Human γD-crystallin (HγD-Crys) is a monomeric, two-domain protein of the lens central nucleus. Both domains of this long lived protein have double Greek key β-sheet folds with well packed hydrophobic cores. Three mutations resulting in amino acid substitutions in the γ-crystallin buried cores (two in the N-terminal domain (N-td) and one in the C-terminal domain (C-td)) cause early onset cataract in mice, presumably an aggregated state of the mutant crystallins. It has not been possible to identify the aggregating precursor within lens tissues. To compare in vivo cataract-forming phenotypes with in vitro unfolding and aggregation of γ-crystallins, mouse mutant substitutions were introduced into HγD-Crys. The mutant proteins L5S, V75D, and I90F were expressed and purified from Escherichia coli. WT HγD-Crys unfolds in vitro through a three-state pathway, exhibiting an intermediate with the N-td unfolded and the C-td native-like. L5S and V75D in the N-td also displayed three-state unfolding transitions, with the first transition, unfolding of the N-td, shifted to significantly lower denaturant concentrations. I90F destabilized the C-td, shifting the overall unfolding transition to lower denaturant concentrations. During thermal denaturation, the mutant proteins exhibited lowered thermal stability compared with WT. Kinetic unfolding experiments showed that the N-tds of L5S and V75D unfolded faster than WT. I90F was globally destabilized and unfolded more rapidly. These results support models of cataract formation in which generation of partially unfolded species are precursors to the aggregated cataractous states responsible for light scattering.     

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.     

Single-molecule FRET study of denaturant induced unfolding of RNase H

Single-molecule fluorescence (F?rster) resonance energy transfer (FRET) experiments were performed on surface-immobilized RNase H molecules as a function of the concentration of the chemical denaturant guanidinium chloride (GdmCl). For comparison, we measured ensemble FRET on RNase H solutions. The single-molecule approach allowed us to study FRET distributions of the subpopulation of unfolded molecules without interference from the folded population. The unfolded ensemble experienced a continuous shift of the FRET efficiency distribution with increasing concentration of GdmCl, indicating a heterogeneous population of expanding, unfolded polypeptide chains. We have analyzed the behavior of the unfolded state quantitatively with a model in which the unfolded state is described by a continuum of substates, with the free energy of each substate linearly coupled to its m-value, the proportionality coefficient between free energy and denaturant activity. By fitting this model to the data, we have derived energetic and structural parameters that describe the unfolded state ensemble. Specifically, we have found that the average size of the unfolded state increases from 23-38 A between 0 and 6 M denaturant. Excellent agreement was achieved between the fitted model and our FRET measurements, and with previously published nuclear magnetic resonance and small-angle X-ray scattering data.     

Small-angle X-ray scattering and single-molecule FRET spectroscopy produce highly divergent views of the low-denaturant unfolded state

The results of more than a dozen single-molecule F?rster resonance energy transfer (smFRET) experiments suggest that chemically unfolded polypeptides invariably collapse from an expanded random coil to more compact dimensions as the denaturant concentration is reduced. In sharp contrast, small-angle X-ray scattering (SAXS) studies suggest that, at least for single-domain proteins at non-zero denaturant concentrations, such compaction may be rare. Here, we explore this discrepancy by studying protein L, a protein previously studied by SAXS (at 5?°C), which suggested fixed unfolded-state dimensions from 1.4 to 5?M guanidine hydrochloride (GuHCl), and by smFRET (at 25?°C), which suggested that, in contrast, the chain contracts by 15-30% over this same denaturant range. Repeating the earlier SAXS study under the same conditions employed in the smFRET studies, we observe little, if any, evidence that the unfolded state of protein L contracts as the concentration of GuHCl is reduced. For example, scattering profiles (and thus the shape and dimensions) collected within ～4?ms after dilution to as low as 0.67?M GuHCl are effectively indistinguishable from those observed at equilibrium at higher denaturant. Our results thus argue that the disagreement between SAXS and smFRET is statistically significant and that the experimental evidence in favor of obligate polypeptide collapse at low denaturant cannot be considered conclusive yet.     

PII structure in the model peptides for unfolded proteins: studies on ubiquitin fragments and several alanine-rich peptides containing QQQ, SSS, FFF, and VVV

A great deal of attention has been paid lately to the structures in unfolded proteins due to the recent discovery of many biologically functional but natively unfolded proteins and the far-reaching implications of order in unfolded states for protein folding. Recently, studies on oligo-Ala, oligo-Lys, oligo-Asp, and oligo-Glu, as well as oligo-Pro, have indicated that the left-handed polyproline II (PII) is the major local structure in these short peptides. Here, we show by NMR and CD studies that ubiquitin fragments, model unfolded peptides composed of nonrepeating amino acids, and four alanine-rich peptides containing QQQ, SSS, FFF, and VVV sequences are all present in aqueous solution predominantly in the extended PII or beta conformation. The results from this and related studies indicate that PII might be a major backbone conformation in unfolded proteins. The presence of defined local backbone structure in unfolded proteins is inconsistent with predictions from random coil models.     

Preorganized secondary structure as an important determinant of fast protein folding

The folding and unfolding kinetics of the B-domain of staphylococcal protein A, a small three-helix bundle protein, were probed by NMR. The lineshape of a single histidine resonance was fit as a function of denaturant to give folding and unfolding rate constants. The B-domain folds extremely rapidly in a two-state manner, with a folding rate constant of 120,000 s-1, making it one of the fastest-folding proteins known. Diffusion-collision theory predicts folding and unfolding rate constants that are in good agreement with the experimental values. The apparent rate constant as a function of denaturant ('chevron plot') is predicted within an order of magnitude. Our results are consistent with a model whereby fast-folding proteins utilize a diffusion-collision mechanism, with the preorganization of one or more elements of secondary structure in the unfolded protein.     

Reversible unfolding of FtsZ cell division proteins from archaea and bacteria. Comparison with eukaryotic tubulin folding and assembly

The stability, refolding, and assembly properties of FtsZ cell division proteins from Methanococcus jannaschii and Escherichia coli have been investigated. Their guanidinium chloride unfolding has been studied by circular dichroism spectroscopy. FtsZ from E. coli and tubulin released the bound guanine nucleotide, coinciding with an initial unfolding stage at low denaturant concentrations, followed by unfolding of the apoprotein. FtsZ from M. jannaschii released its nucleotide without any detectable secondary structural change. It unfolded in an apparently two-state transition at larger denaturant concentrations. Isolated FtsZ polypeptide chains were capable of spontaneous refolding and GTP-dependent assembly. The homologous eukaryotic tubulin monomers misfold in solution, but fold within the cytosolic chaperonin CCT. Analysis of the extensive tubulin loop insertions in the FtsZ/tubulin common core and of the intermolecular contacts in model microtubules and tubulin-CCT complexes shows a loop insertion present at every element of lateral protofilament contact and at every contact of tubulin with CCT (except at loop T7). The polymers formed by purified FtsZ have a distinct limited protofilament association in comparison with microtubules. We propose that the loop insertions of tubulin and its CCT-assisted folding coevolved with the lateral association interfaces responsible for extended two-dimensional polymerization into microtubule polymers.     

Unfolding the unique c-type heme protein, Chlamydomonas reinhardtii cytochrome f

We have studied the unfolding reaction of cytochrome f from the green alga Chlamydomonas reinhardtii. Cytochrome f is different from all other c-type heme proteins in that it is a large, two-domain protein with predominantly beta-sheet structure. Moreover, the sixth axial ligand to the heme-iron is unique in cytochrome f: it is provided by the N-terminal alpha-amino group. Unfolding of oxidized and reduced cytochrome f by guanidine hydrochloride (GuHCl) was monitored by far-UV circular dichroism (CD), Soret absorption, and tyrosine emission: the same unfolding curves were obtained regardless of method. Neither oxidized nor reduced unfolded cytochrome f can be refolded at neutral pH. At pH 3.5 refolding takes place (upon dilution to lower denaturant concentrations or by electron injection to the unfolded, oxidized form), although the reaction is extremely slow. Reduced cytochrome f appears much more resistant towards denaturant perturbation than the oxidized form (in pH range 7-3.5). The heme in unfolded cytochrome f remains low-spin to pH 4 but turns high-spin at pH 3.5 (presumably due to protonation of the N-terminal amino group). Our results suggest that the unfolding process for cytochrome f is complex, involving kinetically trapped intermediates not resolvable by spectroscopy.     

Direct observation of microscopic reversibility in single-molecule protein folding

Both folded and unfolded conformations should be observed for a protein at its melting temperature (T(m)), where DeltaG between these states is zero. In an all-atom molecular dynamics simulation of chymotrypsin inhibitor 2 (CI2) at its experimental T(m), the protein rapidly loses its low-temperature native structure; it then unfolds before refolding to a stable, native-like conformation. The initial unfolding follows the unfolding pathway described previously for higher-temperature simulations: the hydrophobic core is disrupted, the beta-sheet pulls apart and the alpha-helix unravels. The unfolded state reached under these conditions maintains a kernel of structure in the form of a non-native hydrophobic cluster. Refolding simply reverses this path, the side-chain interactions shift, the helix refolds, and the native packing and hydrogen bonds are recovered. The end result of this refolding is not the initial crystal structure; it contains the proper topology and the majority of the native contacts, but the structure is expanded and the contacts are long. We believe this to be the native state at elevated temperature, and the change in volume and contact lengths is consistent with experimental studies of other native proteins at elevated temperature and the chemical denaturant equivalent of T(m).     

Protein denaturant binding polynomials

We show how moments of the denaturant binding distribution function can be extracted from experimental data on the denaturation of a protein as a function of the concentration of denaturant and how in turn these moments can be used to construct the denaturant binding distribution function. This approach is similar to our recent work on using the maximum-entropy method to construct ligand-binding distributions from moments obtained from titration curves for nucleic acids and proteins. As an example we take literature data on the denaturation of ferro- and ferricytochrome c by guanidine hydrochloride and from it construct the denaturant binding polynomial and binding distribution function for the unfolded protein.     

Unfolding rates of barstar determined in native and low denaturant conditions indicate the presence of intermediates

The folding and unfolding rates of the small protein, barstar, have been monitored using stopped-flow measurements of intrinsic tryptophan fluorescence at 25 degrees C, pH 8.5, and have been compared over a wide range of urea and guanidine hydrochloride (GdnHCl) concentrations. When the logarithms of the rates of folding from urea and from GdnHCl unfolded forms are extrapolated linearly with denaturant concentration, the same rate is obtained for folding in zero denaturant. Similar linear extrapolations of rates of unfolding in urea and GdnHCl yield, however, different unfolding rates in zero denaturant, indicating that such linear extrapolations are not valid. It has been difficult, for any protein, to determine unfolding rates under nativelike conditions in direct kinetic experiments. Using a novel strategy of coupling the reactivity of a buried cysteine residue with 5,5'-dithiobis(2-nitrobenzoic acid) (DTNB) to the unfolding reaction of barstar, the global unfolding and refolding rates have now been determined in low denaturant concentrations. The logarithms of unfolding rates obtained at low urea and GdnHCl concentrations show a markedly nonlinear dependence on denaturant concentration and converge to the same unfolding rate in the absence of denaturant. It is shown that the native protein can sample the fully unfolded conformation even in the absence of denaturant. The observed nonlinear dependences of the logarithms of the refolding and unfolding rates observed for both denaturants are shown to be due to the presence of (un)folding intermediates and not due to movements in the position of the transition state with a change in denaturant concentration.     

Rate of intrachain contact formation in an unfolded protein: temperature and denaturant effects

We have measured the effect of temperature and denaturant concentration on the rate of intrachain diffusion in an unfolded protein. After photodissociating a ligand from the heme iron of unfolded horse cytochrome c, we use transient optical absorption spectroscopy to measure the time scale of the diffusive motions that bring the heme, located at His18, into contact with its native ligand, Met80. Measuring the rate at which this 62 residue intrachain loop forms under both folding and unfolding conditions, we find a significant effect of denaturant on the chain dynamics. The diffusion of the chain accelerates as denaturant concentration decreases, with the contact formation rate approaching a value near approximately 6x10(5) s(-1) in the absence of denaturant. This result agrees well with an extrapolation from recent loop formation measurements in short synthetic peptides. The temperature dependence of the rate of contact formation indicates an Arrhenius activation barrier, Ea approximately 20 kJ/mol, at high denaturant concentrations, comparable to what is expected from solvent viscosity effects alone. Although Ea increases by several kBT as denaturant concentration decreases, the overall rate of diffusion nevertheless increases. These results indicate that inter-residue energetic interactions do not control conformational diffusion in unfolded states, even under folding conditions.     

Denaturant-induced Expansion and Compaction of a Multi-domain Protein: IgG

It is generally believed that unfolded or denatured proteins show random-coil statistics and hence their radius of gyration simply scales with solvent quality (or concentration of denaturant). Indeed, nearly all proteins studied thus far have been shown to undergo a gradual and continuous expansion with increasing concentration of denaturant. Here, we use fluorescence correlation spectroscopy (FCS) to show that while protein A, a multi-domain and predominantly helical protein, expands gradually and continuously with increasing concentration of guanidine hydrochloride (GdnHCl), the F(ab′)₂ fragment of goat anti-rabbit antibody IgG, a multi-subunit all β-sheet protein does not show such continuous expansion behavior. Instead, it first expands and then contracts with increasing concentration of GdnHCl. Even more striking is the fact that the hydrodynamic radius of the most expanded F(ab′)₂ ensemble, observed at 3-4 M GdnHCl, is ∼ 3.6 times that of the native protein. Further FCS measurements involving urea and NaCl show that the unusually expanded F(ab′)₂ conformations might be due to electrostatic repulsions. Taken together, these results suggest that specific interactions need to be considered while assessing the conformational and statistical properties of unfolded proteins, particularly under conditions of low solvent quality.