Is There or Isn't There? The Case for (and Against) Residual Structure in Chemically Denatured Proteins

Denatured State Conformational Biases in Three-Helix Bundles Containing Divergent Sequences Localize near Turns and Helix Capping Residues

Rhodopseudomonas palustris cytochrome c', a four-helix bundle, and the second ubiquitin-associated domain, UBA(2), a three-helix bundle from the human homologue of yeast Rad23, HHR23A, deviate from random coil behavior under denaturing conditions in a fold-specific manner. The random coil deviations in each of these folds occur near interhelical turns and loops in their tertiary structures. Here, we examine an additional three-helix bundle with an identical fold to UBA(2), but a highly divergent sequence, the first ubiquitin-associated domain, UBA(1), of HHR23A. We use histidine-heme loop formation methods, employing eight single histidine variants, to probe for denatured state conformational bias of a UBA(1) domain fused to the N-terminus of iso-1-cytochrome c (iso-1-Cytc). Guanidine hydrochloride (GuHCl) denaturation shows that the iso-1-Cytc domain unfolds first, followed by the UBA(1) domain. Denatured state (4 and 6 M GuHCl) histidine-heme loop formation studies show that as the size of the histidine-heme loop increases, loop stability decreases, as expected for the Jacobson-Stockmayer relationship. However, loops formed with His35, His31, and His15, of UBA(1), are 0.6-1.1 kcal/mol more stable than expected from the Jacobson-Stockmayer relationship, confirming the importance of deviations of the denatured state from random coil behavior near interhelical turns of helical domains for facilitating folding to the correct topology. For UBA(1) and UBA(2), hydrophobic clusters on either side of the turns partially explain deviations from random coil behavior; however, helix capping also appears to be important.

An "Onion-like" Model of Protein Unfolding: Collective versus Site Specific Approaches

Approximating protein unfolding by an all-or-none cooperative event is a convenient assumption that can provide precious global information on protein stability. It is however quickly emerging that the scenario is far more complex and that global denaturation curves often hide a rich heterogeneity of states that are largely probe dependent. In this review, we revisit the importance of gaining site-specific information on the unfolding process. We focus on nuclear magnetic resonance, as this is the main technique able to provide site-specific information. We review historical and most modern approaches that have allowed an appreciable advancement of the field of protein folding. We also demonstrate how unfolding is a reporter dependent event, suggesting the outmost importance of selecting the reporter carefully.

Structural and Energetic Characterization of the Denatured State from the Perspectives of Peptides, the Coil Library, and Intrinsically Disordered Proteins

The α and polyproline II (PPII) basins are the two most populated regions of the Ramachandran map when constructed from the protein coil library, a widely used denatured state model built from the segments of irregular structure found in the Protein Data Bank. This indicates the α and PPII conformations are dominant components of the ensembles of denatured structures that exist in solution for biological proteins, an observation supported in part by structural studies of short, and thus unfolded, peptides. Although intrinsic conformational propensities have been determined experimentally for the common amino acids in short peptides, and estimated from surveys of the protein coil library, the ability of these intrinsic conformational propensities to quantitatively reproduce structural behavior in intrinsically disordered proteins (IDPs), an increasingly important class of proteins in cell function, has thus far proven elusive to establish. Recently, we demonstrated that the sequence dependence of the mean hydrodynamic size of IDPs in water and the impact of heat on the coil dimensions, provide access to both the sequence dependence and thermodynamic energies that are associated with biases for the α and PPII backbone conformations. Here, we compare results from peptide-based studies of intrinsic conformational propensities and surveys of the protein coil library to those of the sequence-based analysis of heat effects on IDP hydrodynamic size, showing that a common structural and thermodynamic description of the protein denatured state is obtained.

Residual Structure of Unfolded Ubiquitin as Revealed by Hydrogen/Deuterium-Exchange 2D NMR

The characterization of residual structures persistent in unfolded proteins in concentrated denaturant solution is currently an important issue in studies of protein folding because the residual structure present, if any, in the unfolded state may form a folding initiation site and guide the subsequent folding reactions. Here, we studied the hydrogen/deuterium (H/D)-exchange behavior of unfolded human ubiquitin in 6 M guanidinium chloride. We employed a dimethylsulfoxide (DMSO)-quenched H/D-exchange NMR technique with the use of spin desalting columns, which allowed us to perform a quick medium exchange from 6 M guanidinium chloride to a quenching DMSO solution. Based on the backbone resonance assignment of ubiquitin in the DMSO solution, we successfully investigated the H/D-exchange kinetics of 60 identified peptide amide groups in the ubiquitin sequence. Although a majority of these amide groups were not protected, certain amide groups involved in a middle helix (residues 23-34) and an N-terminal β-hairpin (residues 2-16) were significantly protected with a protection factor of 2.1-4.2, indicating that there were residual structures in unfolded ubiquitin and that these amide groups were more than 52% hydrogen bonded in the residual structures. We show that the hydrogen-bonded residual structures in the α-helix and the β-hairpin are formed even in 6 M guanidinium chloride, suggesting that these residual structures may function as a folding initiation site to guide the subsequent folding reactions of ubiquitin.

Towards Controlling the Reactivity of Enzymes in Mechanochemistry: Inert Surfaces Protect β-Glucosidase Activity During Ball Milling

The activity of β-glucosidases-the enzymes responsible for the final step in the enzymatic conversion of cellulose to glucose-can be maintained and manipulated under mechanochemical conditions in the absence of bulk solvent, either through an unexpected stabilization effect of inert surfaces, or by altering the enzymatic equilibrium. The reported results illustrate unique aspects of mechanoenzymatic reactions that are not observable in conventional aqueous solutions. They also represent the first reported strategies to enhance activity and favor either direction of the reaction under mechanochemical conditions.

Impact of Heat on Coil Hydrodynamic Size Yields the Energetics of Denatured State Conformational Bias

Conformational equilibria in the protein denatured state have key roles regulating folding, stability, and function. The extent of conformational bias in the protein denatured state under folding conditions, however, has thus far proven elusive to quantify, particularly with regard to its sequence dependence and energetic character. To better understand the structural preferences of the denatured state, we analyzed both the sequence dependence to the mean hydrodynamic size of disordered proteins in water and the impact of heat on the coil dimensions, showing that the sequence dependence and thermodynamic energies associated with intrinsic biases for the α and polyproline II (PPII) backbone conformations can be obtained. Experiments that evaluate how the hydrodynamic size changes with compositional changes in the protein reveal amino acid specific preferences for PPII that are in good quantitative agreement with calorimetry-measured values from unfolded peptides and those inferred by survey of the protein coil library. At temperatures above 25 °C, the denatured state follows the predictions of a PPII-dominant ensemble. Heat effects on coil hydrodynamic size indicate the α bias is comparable to the PPII bias at cold temperatures. Though historically thought to give poor resolution to structural details, the hydrodynamic size of the unfolded state is found to be an effective reporter on the extent of the biases for the α and PPII backbone conformations.

Sequence Reversal Prevents Chain Collapse and Yields Heat-Sensitive Intrinsic Disorder

Sequence patterns of charge, hydrophobicity, hydrogen bonding, and other amino acid physicochemical properties contribute to mechanisms of protein folding, but how sequence composition and patterns influence the conformational dynamics of the denatured state ensemble is not fully understood. To investigate structure-sequence relationships in the denatured state, we reversed the sequence of staphylococcal nuclease and characterized its structure, thermodynamic character, and hydrodynamic radius using circular dichroism spectroscopy, dynamic light scattering, analytical ultracentrifugation, and size-exclusion chromatography as a function of temperature. The macromolecular size of "Retro-nuclease" is highly expanded in solution with characteristics similar to biological intrinsically disordered proteins. In contradistinction to a disordered state, Retro-nuclease exhibits a broad sigmoid transition of its hydrodynamic dimensions as temperature is increased, indicating a thermodynamically controlled compaction. Counterintuitively, the magnitude of these temperature-induced hydrodynamic changes exceed that observed from thermal denaturation of folded unaltered staphylococcal nuclease. Undetectable by calorimetry and intrinsic tryptophan fluorescence, the lack of heat capacity or fluorescence changes throughout the thermal transition indicate canonical hydrophobic collapse did not drive the Retro-nuclease structural transitions. Temperature-dependent circular dichroism spectroscopy performed on Retro-nuclease and computer simulations correlate to temperature sensitivity in the intrinsic sampling of backbone conformations for polyproline II and α-helix. The experimental results indicate a role for sequence direction in mediating the collapse of the polypeptide chain, whereas the simulation trends illustrate the generality of the observed heat effects on disordered protein structure.

Thermally versus Chemically Denatured Protein States

Protein unfolding thermodynamic parameters are conventionally extracted from equilibrium thermal and chemical denaturation experiments. Despite decades of work, the degree of structure and the compactness of denatured states populated in these experiments are still open questions. Here, building on previous works, we show that thermally and chemically denatured protein states are distinct from the viewpoint of far-ultraviolet circular dichroism experiments that report on the local conformational status of peptide bonds. The differences identified are independent of protein length, structural class, or experimental conditions, highlighting the presence of two sub-ensembles within the denatured states. The sub-ensembles, U_T and U_D for thermally induced and denaturant-induced unfolded states, respectively, can exclusively exchange populations as a function of temperature at high chemical denaturant concentrations. Empirical analysis suggests that chemically denatured states are ∼50% more expanded than the thermally denatured chains of the same protein. Our observations hint that the strength of protein backbone-backbone interactions and identity versus backbone-solvent/co-solvent interactions determine the conformational distributions. We discuss the implications for protein folding mechanisms, the heterogeneity in relaxation rates, and folding diffusion coefficients.

Biophysical and structural characterization of the small heat shock protein HspA from Thermosynechococcus vulcanus in 2 M urea

Small heat shock proteins (sHSPs) belong to the superfamily of molecular chaperones. They prevent aggregation of partially unfolded or misfolded client proteins, providing protection to organisms under stress conditions. Here, we report the biophysical and structural characterization of a small heat shock protein (HspA) from a thermophilic cyanobacterium Thermosynechococcus vulcanus in the presence of 2 M urea. HspA has been shown to be important for the protection of Photosystem II and the Phycobilisome antenna complex at elevated temperatures. Heterologously expressed HspA requires the presence of 1-2 M urea to maintain its solubility at concentrations required for most characterization methods. Spectroscopic studies reveal the presence of the β-sheet structure and intactness of the tertiary fold in HspA. In vitro assays show that the HspA maintains chaperone-like activity in protecting soluble proteins from thermal aggregation. Chromatography and electron microscopy show that the HspA exists as a mixture of oligomeric forms in the presence of 2 M urea. HspA was successfully crystallized only in the presence of 2 M urea. The crystal structure of HspA shows urea-induced loss of about 30% of the secondary structure without major alteration in the tertiary structure of the protein. The electron density maps reveal changes in the hydrogen bonding network which we attribute to the presence of urea. The crystal structure of HspA demonstrates a mixture of both direct interactions between urea and protein functionalities and interactions between urea and the surrounding solvent that indirectly affect the protein, which are in accordance with previously published studies.

Investigating the Molecular Basis of the Aggregation Propensity of the Pathological D76N Mutant of Beta-2 Microglobulin: Role of the Denatured State

Beta-2 microglobulin (β2m) is a protein responsible for a pathologic condition, known as dialysis-related amyloidosis (DRA), caused by its aggregation and subsequent amyloid formation. A naturally occurring mutation of β2m, D76N, presents a higher amyloidogenic propensity compared to the wild type counterpart. Since the three-dimensional structure of the protein is essentially unaffected by the mutation, the increased aggregation propensity of D76N has been generally ascribed to its lower thermodynamic stability and increased dynamics. In this study we compare the equilibrium unfolding and the aggregation propensity of wild type β2m and D76N variant at different experimental conditions. Our data revealed a surprising effect of the D76N mutation in the residual structure of the denatured state, which appears less compact than that of the wild type protein. A careful investigation of the structural malleability of the denatured state of wild type β2m and D76N pinpoint a clear role of the denatured state in triggering the amyloidogenic propensity of the protein. The experimental results are discussed in the light of the previous work on β2m and its role in disease.

Protein unfolding mechanisms and their effects on folding experiments

In this review, I discuss the various methods researchers use to unfold proteins in the lab in order to understand protein folding both in vitro and in vivo. The four main techniques, chemical-, heat-, pressure- and force-denaturation, produce distinctly different unfolded conformational ensembles. Recent measurements have revealed different folding kinetics from different unfolding mechanisms. Thus, comparing these distinct unfolded ensembles sheds light on the underlying free energy landscape of folding.

The Folding Pathway of the KIX Domain

The KIX domain is an 89-residues globular domain with an important role in mediating protein-protein interactions. The presence of two distinct binding sites in such a small domain makes KIX a suitable candidate to investigate the effect of the potentially divergent demands between folding and function. Here, we report an extensive mutational analysis of the folding pathway of the KIX domain, based on 30 site-directed mutants, which allow us to assess the structures of both the transition and denatured states. Data reveal that, while the transition state presents mostly native-like interactions, the denatured state is somewhat misfolded. We mapped some of the non-native contacts in the denatured state using a second round of mutagenesis, based on double mutant cycles on 15 double mutants. Interestingly, such a misfolding arises from non-native interactions involving the residues critical for the function of the protein. The results described in this work appear to highlight the diverging demands between folding and function that may lead to misfolding, which may be observed in the early stages of folding.

Conformations of a Metastable SH3 Domain Characterized by smFRET and an Excluded-Volume Polymer Model

Conformational states of the metastable drkN SH3 domain were characterized using single-molecule fluorescence techniques. Under nondenaturing conditions, two Förster resonance energy transfer (FRET) populations were observed that corresponded to a folded and an unfolded state. FRET-estimated radii of gyration and hydrodynamic radii estimated by fluorescence correlation spectroscopy of the two coexisting conformations are in agreement with previous ensemble x-ray scattering and NMR measurements. Surprisingly, when exposed to high concentrations of urea and GdmCl denaturants, the protein still exhibits two distinct FRET populations. The dominant conformation is expanded, showing a low FRET efficiency, consistent with the expected behavior of a random chain with excluded volume. However, approximately one-third of the drkN SH3 conformations showed high, nearly 100%, FRET efficiency, which is shown to correspond to denaturation-induced looped conformations that remain stable on a timescale of at least 100 μs. These loops may contain interconverting conformations that are more globally collapsed, hairpin-like, or circular, giving rise to the observed heterogeneous broadening of this population. Although the underlying mechanism of chain looping remains elusive, FRET experiments in formamide and dimethyl sulfoxide suggest that interactions between hydrophobic groups in the distal regions may play a significant role in the formation of the looped state.

Insights into Unfolded Proteins from the Intrinsic ϕ/ψ Propensities of the AAXAA Host-Guest Series

Various host-guest peptide series are used by experimentalists as reference conformational states. One such use is as a baseline for random-coil NMR chemical shifts. Comparison to this random-coil baseline, through secondary chemical shifts, is used to infer protein secondary structure. The use of these random-coil data sets rests on the perception that the reference chemical shifts arise from states where there is little or no conformational bias. However, there is growing evidence that the conformational composition of natively and nonnatively unfolded proteins fail to approach anything that can be construed as random coil. Here, we use molecular dynamics simulations of an alanine-based host-guest peptide series (AAXAA) as a model of unfolded and denatured states to examine the intrinsic propensities of the amino acids. We produced ensembles that are in good agreement with the experimental NMR chemical shifts and confirm that the sampling of the 20 natural amino acids in this peptide series is be far from random. Preferences toward certain regions of conformational space were both present and dependent upon the environment when compared under conditions typically used to denature proteins, i.e., thermal and chemical denaturation. Moreover, the simulations allowed us to examine the conformational makeup of the underlying ensembles giving rise to the ensemble-averaged chemical shifts. We present these data as an intrinsic backbone propensity library that forms part of our Structural Library of Intrinsic Residue Propensities to inform model building, to aid in interpretation of experiment, and for structure prediction of natively and nonnatively unfolded states.

β-Structure within the Denatured State of the Helical Protein Domain BBL

Protein denatured states are the origin of both healthy and toxic conformational species. Denatured states of ultrafast folding proteins are of interest in mechanistic studies because they are energetically close to the kinetic bottleneck of folding. However, their transient nature makes them elusive to experiment. Here, we generated the denatured state of the helical domain BBL that is poised to fold in microseconds by a single-point mutation and combined circular dichroism spectroscopy, single-molecule fluorescence fluctuation analysis, and computer simulation to characterize its structure and dynamics. Circular dichroism showed a largely unfolded ensemble with marginal helix but significant β-sheet content. Main-chain structure and dynamics were unaffected by side-chain interactions that stabilize the native state, as revealed by site-directed mutagenesis and nanosecond loop closure kinetics probed by fluorescence correlation spectroscopy. Replica-exchange and constant-temperature molecular dynamics simulations showed a highly collapsed, hydrogen-bonded denatured state containing turn and β-sheet structure and few nucleating helices in an otherwise unfolded ensemble. An irregular β-hairpin element that connects helices in the native fold was poised to be formed. The surprising observation of β-structure in regions that form helices in the native state is reconciled by a generic low-energy pathway from the northwest quadrant of Ramachandran space to the helical basin present under folding conditions, proposed recently. Our results show that, indeed, rapid nucleation of helix emanates from β-structure formed early within a collapsed ensemble of unfolded conformers.

Cooperative Unfolding of Residual Structure in Heat Denatured Proteins by Urea and Guanidinium Chloride

The denatured states of proteins have always attracted our attention due to the fact that the denatured state is the only experimentally achievable state of a protein, which can be taken as initial reference state for considering the in vitro folding and defining the native protein stability. It is known that heat and guanidinium chloride (GdmCl) give structurally different states of RNase-A, lysozyme, α-chymotrypsinogen A and α-lactalbumin. On the contrary, differential scanning calorimetric (DSC) and isothermal titration calorimetric measurements, reported in the literature, led to the conclusion that heat denatured and GdmCl denatured states are thermodynamically and structurally identical. In order to resolve this controversy, we have measured changes in the far-UV CD (circular dichroism) of these heat-denatured proteins on the addition of different concentrations of GdmCl. The observed sigmoidal curve of each protein was analyzed for Gibbs free energy change in the absence of the denaturant (ΔG0X→D) associated with the process heat denatured (X) state ↔ GdmCl denatured (D) state. To confirm that this thermodynamic property represents the property of the protein alone and is not a manifestation of salvation effect, we measured urea-induced denaturation curves of these heat denatured proteins under the same experimental condition in which GdmCl-induced denaturation was carried out. In this paper we report that (a) heat denatured proteins contain secondary structure, and GdmCl (or urea) induces a cooperative transition between X and D states, (b) for each protein at a given pH and temperature, thermodynamic cycle connects quantities, ΔG0N→X (native (N) state ↔ X state), ΔG0X→D and ΔG0N→D (N state ↔ D state), and

Diffuse binding of Zn(2+) to the denatured ensemble of Cu/Zn superoxide dismutase 1

•

Nearly complete backbone assignments for monomeric SOD1 were obtained.

•

Monomeric SOD1 is unstructured in 9 M urea.

•

Zn²⁺-binding to monomeric SOD1 occurs through diffuse coordination to most His residues.

•

The binding geometry of Zn²⁺ is different than in the native, folded SOD1.

The stability and structural properties of the metalloprotein superoxide dismutase 1 (SOD1) are found to depend critically on metal ions. Native SOD1 monomers coordinate one structural Zn²⁺ and one redox-active Cu^2+/1+ to the active site. To do this, the Zn²⁺ ions need to interact with the SOD1 protein on the denatured side of the folding barrier, prior to the formation of the folding nucleus. In this study, we have examined at residue level the nature of this early Zn²⁺ binding by NMR studies on the urea denatured-state of SOD1. Nearly complete backbone chemical shift assignments were obtained in 9 M urea at physiological pH, conditions at which NMR studies are scarce. Our results demonstrate that SOD1 is predominantly unstructured under these conditions. Chemical-shift changes upon Zn²⁺ titration show that denatured SOD1 retains a significant affinity to Zn²⁺ ions, even in 9 M urea. However, the Zn²⁺ interactions are not limited to the native metal-binding ligands in the two binding sites, but are seen for all His residues. Moreover, the native Cu^2+/1+ ligand H46 seems not to bind as well as the other His residues, while the nearby non-native H43 does bind, indicating that the binding geometry is relaxed. The result suggests that the Zn²⁺-binding observed to catalyze folding of SOD1 in physiological buffer is initiated by diffuse, non-specific coordination to the coil, which subsequently funnels by ligand exchange into the native coordination geometry of the folded monomer. Altogether, this diffuse binding is a result with fundamental implications for folding of metalloproteins in general.

The kinetics of folding of frataxin

The role of the denatured state in protein folding represents a key issue for the proper evaluation of folding kinetics and mechanisms. The yeast ortholog of the human frataxin, a mitochondrial protein essential for iron homeostasis and responsible for Friedreich's ataxia, has been shown to undergo cold denaturation above 0 °C, in the absence of chemical denaturants. This interesting property provides the unique opportunity to explore experimentally the molecular mechanism of both the hot and cold denaturation. In this work, we present the characterization of the temperature and urea dependence of the folding kinetics of yeast frataxin, and show that while at neutral pH and in the absence of a denaturant a simple two-state model may satisfactorily describe the temperature dependence of the folding and unfolding rate constants, the results obtained in urea over a wide range of pH reveal an intriguing complexity, suggesting that folding of frataxin involves a broad smooth free energy barrier.

Protein Denaturation with Guanidinium: A 2D-IR Study

Guanidinium (Gdm⁺) is a widely used denaturant, but it is still largely unknown how it operates at the molecular level. In particular, the effect of guanidinium on the different types of secondary structure motifs of proteins is at present not clear. Here, we use two-dimensional infrared spectroscopy (2D-IR) to investigate changes in the secondary structure of two proteins with mainly α-helical or β-sheet content upon addition of Gdm-¹³C¹⁵N₃·Cl. We find that upon denaturation, the β-sheet protein shows a complete loss of β-sheet structure, whereas the α-helical protein maintains most of its secondary structure. These results suggest that Gdm⁺ disrupts β-sheets much more efficiently than α-helices, possibly because in the former, hydrophobic interactions are more important and the number of dangling hydrogen bonds is larger.

Folding the proteome

Protein folding is an essential prerequisite for protein function and hence cell function. Kinetic and thermodynamic studies of small proteins that refold reversibly were essential for developing our current understanding of the fundamentals of protein folding mechanisms. However, we still lack sufficient understanding to accurately predict protein structures from sequences, or the effects of disease-causing mutations. To date, model proteins selected for folding studies represent only a small fraction of the complexity of the proteome and are unlikely to exhibit the breadth of folding mechanisms used in vivo. We are in urgent need of new methods - both theoretical and experimental - that can quantify the folding behavior of a truly broad set of proteins under in vivo conditions. Such a shift in focus will provide a more comprehensive framework from which to understand the connections between protein folding, the molecular basis of disease, and cell function and evolution.

The folding pathway of a functionally competent C-terminal domain of nucleophosmin: protein stability and denatured state residual structure

Nucleophosmin (NPM1) is a nucleolar protein implicated in ribosome biogenesis, centrosome duplication and cell cycle control; the NPM1 gene is the most frequent target for mutations in Acute Myeloid Leukemia. Mutations map to the C-terminal domain of the protein and cause its unfolding, loss of DNA binding properties and aberrant cellular localization. Here we investigate the folding pathway and denatured state properties of a NPM1 C-terminal domain construct encompassing the last 70 residues in the reference sequence. This construct is more stable than the previously characterized domain, which consisted of the last 53 residues. Data reveal that, similarly to what was discovered for the shorter construct, also the 70-residue construct of NPM1 displays a detectable residual structure in its denatured state. The higher stability of the latter domain allows us to conclude that the denatured state is robust to changes in solvent composition and that it consists of a discrete state in equilibrium with the expanded fully unfolded conformation. This observation, which might appear as a technicality, is in fact of general importance for the understanding of the folding of proteins. The implications of our results are discussed in the context of previous works on single domain helical proteins.

Secondary structure, a missing component of sequence-based minimotif definitions

Minimotifs are short contiguous segments of proteins that have a known biological function. The hundreds of thousands of minimotifs discovered thus far are an important part of the theoretical understanding of the specificity of protein-protein interactions, posttranslational modifications, and signal transduction that occur in cells. However, a longstanding problem is that the different abstractions of the sequence definitions do not accurately capture the specificity, despite decades of effort by many labs. We present evidence that structure is an essential component of minimotif specificity, yet is not used in minimotif definitions. Our analysis of several known minimotifs as case studies, analysis of occurrences of minimotifs in structured and disordered regions of proteins, and review of the literature support a new model for minimotif definitions that includes sequence, structure, and function.

Measurement of energy landscape roughness of folded and unfolded proteins

The dynamics of protein conformational changes, from protein folding to smaller changes, such as those involved in ligand binding, are governed by the properties of the conformational energy landscape. Different techniques have been used to follow the motion of a protein over this landscape and thus quantify its properties. However, these techniques often are limited to short timescales and low-energy conformations. Here, we describe a general approach that overcomes these limitations. Starting from a nonnative conformation held by an aromatic disulfide bond, we use time-resolved spectroscopy to observe nonequilibrium backbone dynamics over nine orders of magnitude in time, from picoseconds to milliseconds, after photolysis of the disulfide bond. We find that the reencounter probability of residues that initially are in close contact decreases with time following an unusual power law that persists over the full time range and is independent of the primary sequence. Model simulations show that this power law arises from subdiffusional motion, indicating a wide distribution of trapping times in local minima of the energy landscape, and enable us to quantify the roughness of the energy landscape (4-5 k(B)T). Surprisingly, even under denaturing conditions, the energy landscape remains highly rugged with deep traps (>20 k(B)T) that result from multiple nonnative interactions and are sufficient for trapping on the millisecond timescale. Finally, we suggest that the subdiffusional motion of the protein backbone found here may promote rapid folding of proteins with low contact order by enhancing contact formation between nearby residues.

NMR insights into folding and self-association of Plasmodium falciparum P2

The eukaryotic 60S-ribosomal stalk is composed of acidic ribosomal proteins (P1 and P2) and neutral protein P0, which are thought to be associated as a pentameric structure, [2P1, 2P2, P0]. Plasmodium falciparum P2 (PfP2) appears to play additional non-ribosomal functions associated with its tendency for homo-oligomerization. Recombinant bacterially expressed PfP2 protein also undergoes self-association, as shown by SDS-PAGE analysis and light scattering studies. Secondary structure prediction algorithms predict the native PfP2 protein to be largely helical and this is corroborated by circular dichroism investigation. The ¹H-¹⁵N HSQC spectrum of native P2 showed only 43 cross peaks compared to the expected 138. The observed peaks were found to belong to the C-terminal region, suggesting that this segment is flexible and solvent exposed. In 9 M urea denaturing conditions the chain exhibited mostly non-native β structural propensity. ¹⁵N Relaxation data for the denatured state indicated substantial variation in ms-µs time scale motion along the chain. Average area buried upon folding (AABUF) calculations on the monomer enabled identification of hydrophobic patches along the sequence. Interestingly, the segments of slower motion in the denatured state coincided with these hydrophobic patches, suggesting that in the denatured state the monomeric chain undergoes transient hydrophobic collapse. The implications of these results for the folding mechanism and self-association of PfP2 are discussed.

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.

Stable folding core in the folding transition state of an alpha-helical integral membrane protein

Defining the structural features of a transition state is important in understanding a folding reaction. Here, we use Φ-value and double mutant analyses to probe the folding transition state of the membrane protein bacteriorhodopsin. We focus on the final C-terminal helix, helix G, of this seven transmembrane helical protein. Φ-values could be derived for 12 amino acid residues in helix G, most of which have low or intermediate values, suggesting that native structure is disrupted at these amino acid positions in the transition state. Notably, a cluster of residues between E204 and M209 all have Φ-values close to zero. Disruption of helix G is further confirmed by a low Φ-value of 0.2 between residues T170 on helix F and S226 on helix G, suggesting the absence of a native hydrogen bond between helices F and G. Φ-values for paired mutations involved in four interhelical hydrogen bonds revealed that all but one of these bonds is absent in the transition state. The unstructured helix G contrasts with Φ-values along helix B that are generally high, implying native structure in helix B in the transition state. Thus helix B seems to constitute part of a stable folding nucleus while the consolidation of helix G is a relatively late folding event. Polarization of secondary structure correlates with sequence position, with a structured helix B near the N terminus contrasting with an unstructured C-terminal helix G.

Diverse pathways of oxidative folding of disulfide proteins: underlying causes and folding models

The pathway of oxidative folding of disulfide proteins exhibits a high degree of diversity, which is manifested mainly by distinct structural heterogeneity and diverse rearrangement pathways of folding intermediates. During the past two decades, the scope of this diversity has widened through studies of more than 30 disulfide-rich proteins by various laboratories. A more comprehensive landscape of the mechanism of protein oxidative folding has emerged. This review will cover three themes. (1) Elaboration of the scope of diversity of disulfide folding pathways, including the two opposite extreme models, represented by bovine pancreatic trypsin inhibitor (BPTI) and hirudin. (2) Demonstration of experimental evidence accounting for the underlying mechanism of the folding diversity. (3) Discussion of the convergence between the extreme models of oxidative folding and models of conventional conformational folding (framework model, hydrophobic collapse model).

Manifestations of native topology in the denatured state ensemble of Rhodopseudomonas palustris cytochrome c'

To provide insight into the role of local sequence in the nonrandom coil behavior of the denatured state, we have extended our measurements of histidine-heme loop formation equilibria for cytochrome c' to 6 M guanidine hydrochloride. We observe that there is some reduction in the scatter about the best fit line of loop stability versus loop size data in 6 M versus 3 M guanidine hydrochloride, but the scatter is not eliminated. The scaling exponent, ν(3), of 2.5 ± 0.2 is also similar to that found previously in 3 M guanidine hydrochloride (2.6 ± 0.3). Rates of histidine-heme loop breakage in the denatured state of cytochrome c' show that some histidine-heme loops are significantly more persistent than others at both 3 and 6 M guanidine hydrochloride. Rates of histidine-heme loop formation more closely approximate random coil behavior. This observation indicates that heterogeneity in the denatured state ensemble results mainly from contact persistence. When mapped onto the structure of cytochrome c', the histidine-heme loops with slow breakage rates coincide with chain reversals between helices 1 and 2 and between helices 2 and 3. Molecular dynamics simulations of the unfolding of cytochrome c' at 498 K show that these reverse turns persist in the unfolded state. Thus, these portions of the primary structure of cytochrome c' set up the topology of cytochrome c' in the denatured state, predisposing the protein to fold efficiently to its native structure.

Contribution of long-range interactions to the secondary structure of an unfolded globin

This work explores the effect of long-range tertiary contacts on the distribution of residual secondary structure in the unfolded state of an alpha-helical protein. N-terminal fragments of increasing length, in conjunction with multidimensional nuclear magnetic resonance, were employed. A protein representative of the ubiquitous globin fold was chosen as the model system. We found that, while most of the detectable alpha-helical population in the unfolded ensemble does not depend on the presence of the C-terminal region (corresponding to the native G and H helices), specific N-to-C long-range contacts between the H and A-B-C regions enhance the helical secondary structure content of the N terminus (A-B-C regions). The simple approach introduced here, based on the evaluation of N-terminal polypeptide fragments of increasing length, is of general applicability to identify the influence of long-range interactions in unfolded proteins.

The folding of single domain proteins--have we reached a consensus?

Rather than stressing the most recent advances in the field, this review highlights the fundamental topics where disagreement remains and where adequate experimental data are lacking. These topics include properties of the denatured state and the role of residual structure, the nature of the fundamental steps and barriers, the extent of pathway heterogeneity and non-native interactions, recent comparisons between theory and experiment, and finally, dynamical properties of the folding reaction.

Unfolded-state dynamics and structure of protein L characterized by simulation and experiment

While several experimental techniques now exist for characterizing protein unfolded states, all-atom simulation of unfolded states has been challenging due to the long time scales and conformational sampling required. We address this problem by using a combination of accelerated calculations on graphics processor units and distributed computing to simulate tens of thousands of molecular dynamics trajectories each up to approximately 10 mus (for a total aggregate simulation time of 127 ms). We used this approach in conjunction with Trp-Cys contact quenching experiments to characterize the unfolded structure and dynamics of protein L. We employed a polymer theory method to make quantitative comparisons between high-temperature simulated and chemically denatured experimental ensembles and find that reaction-limited quenching rates calculated from simulation agree remarkably well with experiment. In both experiment and simulation, we find that unfolded-state intramolecular diffusion rates are very slow compared to highly denatured chains and that a single-residue mutation can significantly alter unfolded-state dynamics and structure. This work suggests a view of the unfolded state in which surprisingly low diffusion rates could limit folding and opens the door for all-atom molecular simulation to be a useful predictive tool for characterizing protein unfolded states along with experiments that directly measure intramolecular diffusion.

The hydrodynamic and conformational properties of denatured proteins in dilute solutions

Published data on the characterization of unfolded proteins in dilute solutions in aqueous guanidine hydrochloride are analyzed to show that the data are not fit by either the random flight or wormlike chain models for linear chains. The analysis includes data on the intrinsic viscosity, root-mean-square radius of gyration, from small-angle X-ray scattering, and hydrodynamic radius, from the translational diffusion coefficient. It is concluded that residual structure consistent with that deduced from nuclear magnetic resonance on these solutions can explain the dilute solution results in a consistent manner through the presence of ring structures, which otherwise have an essentially flexible coil conformation. The ring structures could be in a state of continual flux and rearrangement. Calculation of the radius of gyration for the random-flight model gives a similar reduction of this measure for chains joined at their endpoints, or those containing loop with two dangling ends, each one-fourth the total length of the chain. This relative insensitivity to the details of the ring structure is taken to support the behavior observed across a range of proteins.

Deciphering the folding transition state structure and denatured state properties of nucleophosmin C-terminal domain

Nucleophosmin (NPM1), one of the most abundant nucleolar proteins, is a frequent target of oncogenic mutations in acute myeloid leukaemia (AML). Mutation-induced changes at the C-terminal domain of NPM1 (Cter-NPM1) compromise its stability and cause the aberrant translocation of NPM1 to the cytosol. Hence, this protein represents a suitable candidate to investigate the relations between folding and disease. Since Cter-NPM1 folds via a compact denatured state, stabilization of the folded state of the mutated variants demands detailed structural information on both the native and denatured states. Here, we present the characterization of the complete folding pathway of Cter-NPM1 and provide molecular details for both the transition and the denatured states. The structure of the transition state was assessed by Phi-value analysis, whereas residual structure in the denatured state was mapped by evaluating the effect of mutations as modulated by conditions promoting denatured state compaction. Data reveal that folding of Cter-NPM1 proceeds via an extended nucleus and that the denatured state retains significant malleable structure at the interface between the second and third helices. Our observations constitute the essential prerequisite for structure-based drug-design studies, aimed at identifying molecules that may rescue pathological NPM1 mutants by stabilizing the native-like state.

Modulating native-like residual structure in the fully denatured state of photoactive yellow protein affects its refolding

Residual structure in the fully unfolded state is a key element for understanding protein folding. We show that the residual structure in fully denatured photoactive yellow protein (PYP) is affected by isomerization of its p-coumaric acid (pCA) chromophore. The exposure of total surface area and hydrophobic surface area upon unfolding was quantified by denaturant m values and heat capacity changes (DeltaC(p)), respectively. The exposure of the buried surface area upon the unfolding of the acid-denatured state of PYP containing trans-pCA is approximately 20% smaller than that of the native state. In contrast, for the partially unfolded pB photocycle intermediate containing cis-pCA, unfolding-induced exposure of the surface area is not decreased. These results show that pCA photoisomerization reduces residual structure in the fully unfolded state. Thus, residual structure in the fully unfolded state of PYP is under direct experimental control by photoexcitation. The sensitivity of the unfolded state to pCA isomerization provides a novel criterion that residual structure in the unfolded state of PYP is native-like, involving native-like protein-chromophore interactions. A largely untested prediction is that native-like residual structure facilitates the conformational search during folding. In the case of PYP, refolding from the less disordered fully unfolded state containing trans-pCA indeed is substantially accelerated. The burial of hydrophobic surface area in the fully unfolded state suggests that a significant part of the hydrophobic collapse process already has occurred in the denatured state.

Mathematics, thermodynamics, and modeling to address ten common misconceptions about protein structure, folding, and stability

To fully understand the roles proteins play in cellular processes, students need to grasp complex ideas about protein structure, folding, and stability. Our current understanding of these topics is based on mathematical models and experimental data. However, protein structure, folding, and stability are often introduced as descriptive, qualitative phenomena in undergraduate classes. In the process of learning about these topics, students often form incorrect ideas. For example, by learning about protein folding in the context of protein synthesis, students may come to an incorrect conclusion that once synthesized on the ribosome, a protein spends its entire cellular life time in its fully folded native confirmation. This is clearly not true; proteins are dynamic structures that undergo both local fluctuations and global unfolding events. To prevent and address such misconceptions, basic concepts of protein science can be introduced in the context of simple mathematical models and hands-on explorations of publicly available data sets. Ten common misconceptions about proteins are presented, along with suggestions for using equations, models, sequence, structure, and thermodynamic data to help students gain a deeper understanding of basic concepts relating to protein structure, folding, and stability.

Structural heterogeneity of 6 M GdmCl-denatured proteins: implications for the mechanism of protein folding

An in vitro experiment with protein folding is typically initiated with 6 M GdmCl-denatured proteins, which are generally considered fully unfolded. However, studies conducted by various laboratories have shown that many 6 M GdmCl-denatured proteins are structurally heterogeneous and still retain nativelike residual structures. The extent of conformational heterogeneity of the 6 M GdmCl-denatured protein has significant implications for the folding landscape as well as the interpretation of the observed early stage folding mechanism. Using the method of disulfide scrambling, we are able to gain rough insight into the diverse structural properties of 6 M GdmCl-denatured proteins. It demonstrates that most 6 M GdmCl-denatured proteins are approximately fully denatured, but partially unfolded. Most of them comprise diverse conformational isomers. We review here the cumulative evidence obtained from various laboratories and also provide experimental data obtained in our laboratory.

Paramagnetic relaxation enhancements in unfolded proteins: theory and application to drkN SH3 domain

Site-directed spin labeling in combination with paramagnetic relaxation enhancement (PRE) measurements is one of the most promising techniques for studying unfolded proteins. Since the pioneering work of Gillespie and Shortle (J Mol Biol 1997;268:158), PRE data from unfolded proteins have been interpreted using the theory that was originally developed for rotational spin relaxation. At the same time, it can be readily recognized that the relative motion of the paramagnetic tag attached to the peptide chain and the reporter spin such as (1)H(N) is best described as a translation. With this notion in mind, we developed a number of models for the PRE effect in unfolded proteins: (i) mutual diffusion of the two tethered spheres, (ii) mutual diffusion of the two tethered spheres subject to a harmonic potential, (iii) mutual diffusion of the two tethered spheres subject to a simulated mean-force potential (Smoluchowski equation); (iv) explicit-atom molecular dynamics simulation. The new models were used to predict the dependences of the PRE rates on the (1)H(N) residue number and static magnetic field strength; the results are appreciably different from the Gillespie-Shortle model. At the same time, the Gillespie-Shortle approach is expected to be generally adequate if the goal is to reconstruct the distance distributions between (1)H(N) spins and the paramagnetic center (provided that the characteristic correlation time is known with a reasonable accuracy). The theory has been tested by measuring the PRE rates in three spin-labeled mutants of the drkN SH3 domain in 2M guanidinium chloride. Two modifications introduced into the measurement scheme-using a reference compound to calibrate the signals from the two samples (oxidized and reduced) and using peak volumes instead of intensities to determine the PRE rates-lead to a substantial improvement in the quality of data. The PRE data from the denatured drkN SH3 are mostly consistent with the model of moderately expanded random-coil protein, although part of the data point toward a more compact structure (local hydrophobic cluster). At the same time, the radius of gyration reported by Choy et al. (J Mol Biol 2002;316:101) suggests that the protein is highly expanded. This seemingly contradictory evidence can be reconciled if one assumes that denatured drkN SH3 forms a conformational ensemble that is dominated by extended conformations, yet also contains compact (collapsed) species. Such behavior is apparently more complex than predicted by the model of a random-coil protein in good solvent/poor solvent.

Non-sequence-specific interactions can account for the compaction of proteins unfolded under "native" conditions

Proteins unfolded by high concentrations of chemical denaturants adopt expanded, largely structure-free ensembles of conformations that are well approximated as random coils. In contrast, globular proteins unfolded under less denaturing conditions (via mutations, or transiently unfolded after a rapid jump to native conditions) and molten globules (arising due to mutations or cosolvents) are often compact. Here we explore the origins of this compaction using a truncated equilibrium-unfolded variant of the 57-residue FynSH3 domain. As monitored by far-UV circular dichroism, NMR spectroscopy, and hydrogen-exchange kinetics, CDelta4 (a 4-residue carboxy-terminal deletion variant of FynSH3) appears to be largely unfolded even in the absence of denaturant. Nevertheless, CDelta4 is quite compact under these conditions, with a hydrodynamic radius only slightly larger than that of the native protein. In order to understand the origins of this molten-globule-like compaction, we have characterized a random sequence polypeptide of identical amino acid composition to CDelta4. Notably, we find that the hydrodynamic radius of this random sequence polypeptide also approaches that of the native protein. Thus, while native-like interactions may contribute to the formation of compact "unfolded" states, it appears that non-sequence-specific monomer-monomer interactions can also account for the dramatic compaction observed for molten globules and the "physiological" unfolded state.

Deconstructing the native state: energy landscapes, function, and dynamics of globular proteins

Proteins are highly complex molecules with features exquisitely selected by nature to carry out essential biological functions. Physical chemistry and polymer physics provide us with the tools needed to make sense of this complexity. Upon translation, many proteins fold to a thermodynamically stable form known as the native state. The native state is not static, but consists of a hierarchy of conformations, that are continuously explored through dynamics. In this review we provide a brief introduction to some of the core concepts required in the discussion of the protein native dynamics using energy landscapes ideas. We first discuss recent works which have challenged the structure-function paradigm by demonstrating function in disordered proteins. Next we examine the hierarchical organization in the energy landscapes using atomistic molecular dynamics simulations and principal component analysis. In particular, the role of direct and water-mediated contacts in sculpting the landscape is elaborated. Another approach to studying the native state ensemble is based on choosing high-resolution order parameters for computing one- or two-dimensional free energy surfaces. We demonstrate that 2D free energy surfaces provide rich thermodynamic and kinetic information about the native state ensemble. Brownian dynamics simulations on such a surface indicate that protein conformational dynamics is weakly activated. Finally, we briefly discuss implicit and coarse-grained protein models and emphasize the solvent role in determining native state structure and dynamics.

Thermodynamics of loop formation in the denatured state of rhodopseudomonas palustris cytochrome c': scaling exponents and the reconciliation problem

The observation that denatured proteins yield scaling exponents, nu, consistent with random-coil behavior and yet can also have pockets of residual or nonrandom structure has been termed the "reconciliation problem". To provide greater insight into the denatured state of a foldable sequence, we have measured histidine-heme loop formation equilibria in the denatured state of a class II c-type cytochrome, cytochrome c' from Rhodopseudomonas palustris. We have prepared a series of variants that provide His-heme loop stabilities, pK(loop)(His), for loop sizes ranging from 10 to 111 residues at intervals of 7 to 11 residues along the sequence of the protein. We observe a scaling exponent for loop formation, nu(3), of 2.5+/-0.3. Theoretical values for nu(3) range from 1.8 to 2.4; thus, the observed nu(3) is consistent with random-coil behavior. However, in contrast to data for loop formation as a function of loop size obtained with peptides of homogeneous sequence, we observe considerable scatter about the linear dependence of loop stability on loop size. Thus, foldable sequences behave very differently from homogeneous peptide sequences. The observed scatter suggests that there is considerable variation in the conformational properties along the backbone of a foldable sequence, consistent with alternating compact and extended regions. With regard to the reconciliation problem, it is evident that a scaling exponent consistent with a random coil is necessary but not sufficient to demonstrate random-coil behavior.