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

Single-molecule spectroscopy reveals chaperone-mediated expansion of substrate protein

Molecular chaperones are an essential part of the machinery that avoids protein aggregation and misfolding in vivo. However, understanding the molecular basis of how chaperones prevent such undesirable interactions requires the conformational changes within substrate proteins to be probed during chaperone action. Here we use single-molecule fluorescence spectroscopy to investigate how the DnaJ-DnaK chaperone system alters the conformational distribution of the denatured substrate protein rhodanese. We find that in a first step the ATP-independent binding of DnaJ to denatured rhodanese results in a compact denatured ensemble of the substrate protein. The following ATP-dependent binding of multiple DnaK molecules, however, leads to a surprisingly large expansion of denatured rhodanese. Molecular simulations indicate that hard-core repulsion between the multiple DnaK molecules provides the underlying mechanism for disrupting even strong interactions within the substrate protein and preparing it for processing by downstream chaperone systems.

Single-molecule spectroscopy of the unexpected collapse of an unfolded protein at low pH

The dimensions of intrinsically disordered and unfolded proteins critically depend on the solution conditions, such as temperature, pH, ionic strength, and osmolyte or denarurant concentration. However, a quantitative understanding of how the complex combination of chain-chain and chain-solvent interactions is affected by the solvent is still missing. Here, we take a step towards this goal by investigating the combined effect of pH and denaturants on the dimensions of an unfolded protein. We use single-molecule fluorescence spectroscopy to extract the dimensions of unfolded cold shock protein (CspTm) in mixtures of the denaturants urea and guanidinium chloride (GdmCl) at neutral and acidic pH. Surprisingly, even though a change in pH from 7 to 2.9 increases the net charge of CspTm from -3.8 to +10.2, the radius of gyration of the chain is very similar under both conditions, indicating that protonation of acidic side chains at low pH results in additional hydrophobic interactions. We use a simple shared binding site model that describes the joint effect of urea and GdmCl, together with polyampholyte theory and an ion cloud model that includes the chemical free energy of counterion interactions and side chain protonation, to quantify this effect.

Temperature-dependent solvation modulates the dimensions of disordered proteins

For disordered proteins, the dimensions of the chain are an important property that is sensitive to environmental conditions. We have used single-molecule Förster resonance energy transfer to probe the temperature-induced chain collapse of five unfolded or intrinsically disordered proteins. Because this behavior is sensitive to the details of intrachain and chain-solvent interactions, the collapse allows us to probe the physical interactions governing the dimensions of disordered proteins. We find that each of the proteins undergoes a collapse with increasing temperature, with the most hydrophobic one, λ-repressor, undergoing a reexpansion at the highest temperatures. Although such a collapse might be expected due to the temperature dependence of the classical "hydrophobic effect," remarkably we find that the largest collapse occurs for the most hydrophilic, charged sequences. Using a combination of theory and simulation, we show that this result can be rationalized in terms of the temperature-dependent solvation free energies of the constituent amino acids, with the solvation properties of the most hydrophilic residues playing a large part in determining the collapse.

Microsecond barrier-limited chain collapse observed by time-resolved FRET and SAXS

It is generally held that random-coil polypeptide chains undergo a barrier-less continuous collapse when the solvent conditions are changed to favor the fully folded native conformation. We test this hypothesis by probing intramolecular distance distributions during folding in one of the paradigms of folding reactions, that of cytochrome c. The Trp59-to-heme distance was probed by time-resolved Förster resonance energy transfer in the microsecond time range of refolding. Contrary to expectation, a state with a Trp59-heme distance close to that of the guanidinium hydrochloride (GdnHCl) denatured state is present after ~27 μs of folding. A concomitant decrease in the population of this state and an increase in the population of a compact high-FRET (Förster resonance energy transfer) state (efficiency>90%) show that the collapse is barrier limited. Small-angle X-ray scattering (SAXS) measurements over a similar time range show that the radius of gyration under native favoring conditions is comparable to that of the GdnHCl denatured unfolded state. An independent comprehensive global thermodynamic analysis reveals that marginally stable partially folded structures are also present in the nominally unfolded GdnHCl denatured state. These observations suggest that specifically collapsed intermediate structures with low stability in rapid equilibrium with the unfolded state may contribute to the apparent chain contraction observed in previous fluorescence studies using steady-state detection. In the absence of significant dynamic averaging of marginally stable partially folded states and with the use of probes sensitive to distance distributions, barrier-limited chain contraction is observed upon transfer of the GdnHCl denatured state ensemble to native-like conditions.

Physicochemical bases for protein folding, dynamics, and protein-ligand binding

Proteins are essential parts of living organisms and participate in virtually every process within cells. As the genomic sequences for increasing number of organisms are completed, research into how proteins can perform such a variety of functions has become much more intensive because the value of the genomic sequences relies on the accuracy of understanding the encoded gene products. Although the static three-dimensional structures of many proteins are known, the functions of proteins are ultimately governed by their dynamic characteristics, including the folding process, conformational fluctuations, molecular motions, and protein-ligand interactions. In this review, the physicochemical principles underlying these dynamic processes are discussed in depth based on the free energy landscape (FEL) theory. Questions of why and how proteins fold into their native conformational states, why proteins are inherently dynamic, and how their dynamic personalities govern protein functions are answered. This paper will contribute to the understanding of structure-function relationship of proteins in the post-genome era of life science research.

Fast closure of long loops at the initiation of the folding transition of globular proteins studied by time-resolved FRET-based methods

The protein folding problem would be considered “solved” when it will be possible to “read genes”, i.e., to predict the native fold of proteins, their dynamics, and the mechanism of fast folding based solely on sequence data. The long-term goal should be the creation of an algorithm that would simulate the stepwise mechanism of folding, which constrains the conformational space and in which random search for stable interactions is possible. Here, we focus attention on the initial phases of the folding transition starting with the compact disordered collapsed ensemble, in search of the initial sub-domain structural biases that direct the otherwise stochastic dynamics of the backbone. Our studies are designed to test the “loop hypothesis”, which suggests that fast closure of long loop structures by non-local interactions between clusters of mainly non-polar residues is an essential conformational step at the initiation of the folding transition of globular proteins. We developed and applied experimental methods based on time-resolved resonance excitation energy transfer (trFRET) measurements combined with fast mixing methods and studied the initial phases of the folding of

Comparing protein folding in vitro and in vivo: foldability meets the fitness challenge

In this review, we compare and contrast current knowledge about in vitro and in vivo protein folding. Major advances in understanding fundamental principles underlying protein folding in optimized in vitro conditions have yielded detailed physicochemical principles of folding landscapes for small, single domain proteins. In addition, there has been increased research focusing on the key features of protein folding in the cell that differentiate it from in vitro folding, such as co-translational folding, chaperone-facilitated folding, and folding in crowded conditions with many weak interactions. Yet these two research areas have not been bridged effectively in research carried out to date. This review points to gaps between the two that are ripe for future research. Moreover, we emphasize the biological selection pressures that impact protein folding in vivo and how fitness drives the evolution of protein sequences in ways that may place foldability in tension with other requirements on a given protein. We suggest that viewing the physicochemical process of protein folding through the lens of evolution will unveil new insights and pose novel challenges about in-cell folding landscapes.

Computational and theoretical methods for protein folding

A computational approach is essential whenever the complexity of the process under study is such that direct theoretical or experimental approaches are not viable. This is the case for protein folding, for which a significant amount of data are being collected. This paper reports on the essential role of in silico methods and the unprecedented interplay of computational and theoretical approaches, which is a defining point of the interdisciplinary investigations of the protein folding process. Besides giving an overview of the available computational methods and tools, we argue that computation plays not merely an ancillary role but has a more constructive function in that computational work may precede theory and experiments. More precisely, computation can provide the primary conceptual clues to inspire subsequent theoretical and experimental work even in a case where no preexisting evidence or theoretical frameworks are available. This is cogently manifested in the application of machine learning methods to come to grips with the folding dynamics. These close relationships suggested complementing the review of computational methods within the appropriate theoretical context to provide a self-contained outlook of the basic concepts that have converged into a unified description of folding and have grown in a synergic relationship with their computational counterpart. Finally, the advantages and limitations of current computational methodologies are discussed to show how the smart analysis of large amounts of data and the development of more effective algorithms can improve our understanding of protein folding.

Folding of a large protein at high structural resolution

Kinetic folding of the large two-domain maltose binding protein (MBP; 370 residues) was studied at high structural resolution by an advanced hydrogen-exchange pulse-labeling mass-spectrometry method (HX MS). Dilution into folding conditions initiates a fast molecular collapse into a polyglobular conformation (<20 ms), determined by various methods including small angle X-ray scattering. The compaction produces a structurally heterogeneous state with widespread low-level HX protection and spectroscopic signals that match the equilibrium melting posttransition-state baseline. In a much slower step (7-s time constant), all of the MBP molecules, although initially heterogeneously structured, form the same distinct helix plus sheet folding intermediate with the same time constant. The intermediate is composed of segments that are distant in the MBP sequence but adjacent in the native protein where they close the longest residue-to-residue contact. Segments that are most HX protected in the early molecular collapse do not contribute to the initial intermediate, whereas the segments that do participate are among the less protected. The 7-s intermediate persists through the rest of the folding process. It contains the sites of three previously reported destabilizing mutations that greatly slow folding. These results indicate that the intermediate is an obligatory step on the MBP folding pathway. MBP then folds to the native state on a longer time scale (~100 s), suggestively in more than one step, the first of which forms structure adjacent to the 7-s intermediate. These results add a large protein to the list of proteins known to fold through distinct native-like intermediates in distinct pathways.

Slow proton transfer coupled to unfolding explains the puzzling results of single-molecule experiments on BBL, a paradigmatic downhill folding protein

A battery of thermodynamic, kinetic, and structural approaches has indicated that the small α-helical protein BBL folds-unfolds via the one-state downhill scenario. Yet, single-molecule fluorescence spectroscopy offers a more conflicting view. Single-molecule experiments at pH 6 show a unique half-unfolded conformational ensemble at mid denaturation, whereas other experiments performed at higher pH show a bimodal distribution, as expected for two-state folding. Here we use thermodynamic and laser T-jump kinetic experiments combined with theoretical modeling to investigate the pH dependence of BBL stability, folding kinetics and mechanism within the pH 6-11 range. We find that BBL unfolding is tightly coupled to the protonation of one of its residues with an apparent pKa of ~ 7. Therefore, in chemical denaturation experiments around neutral pH BBL unfolds gradually, and also converts in binary fashion to the protonated species. Moreover, under the single-molecule experimental conditions (denaturant midpoint and 279 K), we observe that proton transfer is much slower than the ~ 15 microseconds folding-unfolding kinetics of BBL. The relaxation kinetics is distinctly biphasic, and the overall relaxation time (i.e. 0.2-0.5 ms) becomes controlled by the proton transfer step. We then show that a simple theoretical model of protein folding coupled to proton transfer explains quantitatively all these results as well as the two sets of single-molecule experiments, including their more puzzling features. Interestingly, this analysis suggests that BBL unfolds following a one-state downhill folding mechanism at all conditions. Accordingly, the source of the bimodal distributions observed during denaturation at pH 7-8 is the splitting of the unique conformational ensemble of BBL onto two slowly inter-converting protonation species. Both, the unprotonated and protonated species unfold gradually (one-state downhill), but they exhibit different degree of unfolding at any given condition because the native structure is less stable for the protonated form.

Submillisecond protein folding events monitored by rapid mixing and mass spectrometry-based oxidative labeling

Kinetic measurements can provide insights into protein folding mechanisms. However, the initial (submillisecond) stages of folding still represent a formidable analytical challenge. A number of ultrarapid triggering techniques have been available for some time, but coupling of these techniques with detection methods that are capable of providing detailed structural information has proven to be difficult. The current work addresses this issue by combining submillisecond mixing with laser-induced oxidative labeling. Apomyoglobin (aMb) serves as a model system for our measurements. Exposure of the protein to a brief pulse of hydroxyl radical (·OH) at different time points during folding introduces covalent modifications at solvent accessible side chains. The extent of labeling is monitored using mass spectrometry-based peptide mapping, providing spatially resolved measurements of changes in solvent accessibility. The submillisecond mixer used here improves the time resolution by a factor of 50 compared to earlier ·OH labeling experiments from our laboratory. Data obtained in this way indicate that early aMb folding events are driven by both local and sequence-remote docking of hydrophobic side chains. Assembly of a partially formed A(E)G(H) scaffold after 0.2 ms is followed by stepwise consolidation that ultimately yields the native state. Major conformational changes go to completion within 0.1 s. The technique introduced here is capable of providing in-depth structural information on very short time scales that have thus far been dominated by low resolution (global) spectroscopic probes. By employing submillisecond mixing in conjunction with slower mixing techniques, it is possible to observe complete folding pathways, from fractions of a millisecond all the way to minutes.

Molecular dynamics simulations of certain mutant peptide models from staphylococcal nuclease reveal that initial hydrophobic collapse associated with turn propensity drives β-hairpin folding

An important nucleation event during the folding of staphylococcal nuclease involves the formation of a β-hairpin by the sequence (21) DTVKLMYKGQPMTFR(35) . Earlier studies show that the turn sequence 'YKGQP' has an important role in the folding of this β-hairpin. To understand the active or passive nature of the turn sequence 'YKGQP' in the folding of the aforementioned β-hairpin sequence, we studied glycine mutant peptides Ac-(2) DTVKLMYGGQPMTFR(16) -NMe (K9G:15), Ac-(2) DTVKLMYKGGPMTFR(16) -NMe (Q11G:15), Ac-(2) DTVKLMYGGGPMTFR(16) -NMe (K9G/Q11G:15), and Ac-(2) DTVKLMGGGGGMTFR(16) -NMe (penta-G:15) by using molecular dynamics simulations, starting with two different unfolded states, polyproline II and extended conformational forms. Further, 5mer mutant turn peptides Ac-(2) YGGQP(6) -NMe (K3G:5), Ac-(2) YKGGP(6) -NMe (Q5G:5), Ac-(2) YGGGP(6) -NMe (K3G/Q5G:5), and Ac-(2) GGGGG(6) -NMe (penta-G:5) were also studied individually. Our results show that an initial hydrophobic collapse and loop closure occurs in all 15mer mutants, but only K9G:15 mutant forms a stable native-like β-hairpin. In the other 15mer mutants, the hydrophobic collapsed state would not proceed to β-hairpin formation. Of the different simulations performed for the penta-G:15 mutant, in only one simulation a nonnative β-hairpin conformation is sampled with highly flexible loop region ((8) GGGGG(12) ), which has no specific conformational preference as a 5mer. While the sequence 'YGGQP' in the K3G:5 simulation shows relatively higher β-turn propensity, the presence of this sequence in K9G:15 peptide seems to be driving the β-hairpin formation. Thus, these results seem to suggest that for the formation of a stable β-hairpin, the initial hydrophobic collapse is to be assisted by a turn propensity. Initial hydrophobic collapse alone is not sufficient to guide β-hairpin formation.

Quantifying the topography of the intrinsic energy landscape of flexible biomolecular recognition

Biomolecular functions are determined by their interactions with other molecules. Biomolecular recognition is often flexible and associated with large conformational changes involving both binding and folding. However, the global and physical understanding for the process is still challenging. Here, we quantified the intrinsic energy landscapes of flexible biomolecular recognition in terms of binding-folding dynamics for 15 homodimers by exploring the underlying density of states, using a structure-based model both with and without considering energetic roughness. By quantifying three individual effective intrinsic energy landscapes (one for interfacial binding, two for monomeric folding), the association mechanisms for flexible recognition of 15 homodimers can be classified into two-state cooperative "coupled binding-folding" and three-state noncooperative "folding prior to binding" scenarios. We found that the association mechanism of flexible biomolecular recognition relies on the interplay between the underlying effective intrinsic binding and folding energy landscapes. By quantifying the whole global intrinsic binding-folding energy landscapes, we found strong correlations between the landscape topography measure Λ (dimensionless ratio of energy gap versus roughness modulated by the configurational entropy) and the ratio of the thermodynamic stable temperature versus trapping temperature, as well as between Λ and binding kinetics. Therefore, the global energy landscape topography determines the binding-folding thermodynamics and kinetics, crucial for the feasibility and efficiency of realizing biomolecular function. We also found "U-shape" temperature-dependent kinetic behavior and a dynamical cross-over temperature for dividing exponential and nonexponential kinetics for two-state homodimers. Our study provides a unique way to bridge the gap between theory and experiments.

Stepwise protein folding at near amino acid resolution by hydrogen exchange and mass spectrometry

The kinetic folding of ribonuclease H was studied by hydrogen exchange (HX) pulse labeling with analysis by an advanced fragment separation mass spectrometry technology. The results show that folding proceeds through distinct intermediates in a stepwise pathway that sequentially incorporates cooperative native-like structural elements to build the native protein. Each step is seen as a concerted transition of one or more segments from an HX-unprotected to an HX-protected state. Deconvolution of the data to near amino acid resolution shows that each step corresponds to the folding of a secondary structural element of the native protein, termed a "foldon." Each folded segment is retained through subsequent steps of foldon addition, revealing a stepwise buildup of the native structure via a single dominant pathway. Analysis of the pertinent literature suggests that this model is consistent with experimental results for many proteins and some current theoretical results. Two biophysical principles appear to dictate this behavior. The principle of cooperativity determines the central role of native-like foldon units. An interaction principle termed "sequential stabilization" based on native-like interfoldon interactions orders the pathway.

Extremely Slow Dynamics of an Abiotic Helical Assembly: Unusual Relevance to the Secondary Structure of Proteins

Serendipitously, we found that isoleucine methylester functionalized perylenediimide 1 undergoes an extremely slow supramolecular helical assembly over a day's time. Surprisingly, heating led to irreversible chiral denaturation. However, reversible helical assembly could be achieved only in the presence of nondenatured aggregates of 1, which act as seeds. The intriguing functional relevance deduced from 1 was employed to draw parallels with the secondary structure of proteins, envisaging its plausible implications.

Synthesis and evaluation of fluorescent cap analogues for mRNA labelling

We describe the synthesis and properties of five dinucleotide fluorescent cap analogues labelled at the ribose of the 7-methylguanosine moiety with either anthraniloyl (Ant) or N-methylanthraniloyl (Mant), which have been designed for the preparation of fluorescent mRNAs via transcription in vitro. Two of the analogues bear a methylene modification in the triphosphate bridge, providing resistance against either the Dcp2 or DcpS decapping enzymes. All these compounds were prepared by ZnCl₂-mediated coupling of a nucleotide P-imidazolide with a fluorescently labelled mononucleotide. To evaluate the utility of these compounds for studying interactions with cap-binding proteins and cap-related cellular processes, both biological and spectroscopic features of those compounds were determined. The results indicate acceptable quantum yields of fluorescence, pH independence, environmental sensitivity, and photostability. The cap analogues are incorporated by RNA polymerase into mRNA transcripts that are efficiently translated in vitro. Transcripts containing fluorescent caps but unmodified in the triphosphate chain are hydrolysed by Dcp2 whereas those containing a α-β methylene modification are resistant. Model studies exploiting sensitivity of Mant to changes of local environment demonstrated utility of the synthesized compounds for studying cap-related proteins.

Polymer scaling laws of unfolded and intrinsically disordered proteins quantified with single-molecule spectroscopy

The dimensions of unfolded and intrinsically disordered proteins are highly dependent on their amino acid composition and solution conditions, especially salt and denaturant concentration. However, the quantitative implications of this behavior have remained unclear, largely because the effective theta-state, the central reference point for the underlying polymer collapse transition, has eluded experimental determination. Here, we used single-molecule fluorescence spectroscopy and two-focus correlation spectroscopy to determine the theta points for six different proteins. While the scaling exponents of all proteins converge to 0.62 ± 0.03 at high denaturant concentrations, as expected for a polymer in good solvent, the scaling regime in water strongly depends on sequence composition. The resulting average scaling exponent of 0.46 ± 0.05 for the four foldable protein sequences in our study suggests that the aqueous cellular milieu is close to effective theta conditions for unfolded proteins. In contrast, two intrinsically disordered proteins do not reach the Θ-point under any of our solvent conditions, which may reflect the optimization of their expanded state for the interactions with cellular partners. Sequence analyses based on our results imply that foldable sequences with more compact unfolded states are a more recent result of protein evolution.

Early hydrophobic collapse of α₁-antitrypsin facilitates formation of a metastable state: insights from oxidative labeling and mass spectrometry

The biologically active conformation of α₁-antitrypsin (α₁AT) and other serine protease inhibitors represents a metastable state, characterized by an exposed reactive center loop (RCL) that acts as bait for the target enzyme. The protein can also adopt an inactive "latent" conformation that has the RCL inserted as a central strand in β-sheet A. This latent form is thermodynamically more stable than the active conformation. Nonetheless, folding of α₁AT consistently yields the active state. The reasons that the metastable form is kinetically preferred remain controversial. The current work demonstrates that a carefully orchestrated folding mechanism prevents RCL insertion into sheet A. Temporal changes in solvent accessibility during folding are monitored using pulsed oxidative labeling and mass spectrometry. The data obtained in this way complement recent hydrogen/deuterium exchange results. Those hydrogen/deuterium exchange measurements revealed that securing of the RCL by hydrogen bonding of the first β-strand in sheet C is one factor that favors formation of the active conformation. The oxidative labeling data presented here reveal that this anchoring is preceded by the formation of hydrophobic contacts in a confined region of the protein. This partial collapse sequesters the RCL insertion site early on and is therefore instrumental in steering α₁AT towards its active conformation. RCL anchoring by hydrogen bonding starts to contribute at a later stage. Together, these two factors ensure that formation of the active conformation is kinetically favored. This work demonstrates how the use of complementary labeling techniques can provide insights into the mechanisms of protracted folding reactions.

Counteracting chemical chaperone effects on the single-molecule α-synuclein structural landscape

Protein structure and function depend on a close interplay between intrinsic folding energy landscapes and the chemistry of the protein environment. Osmolytes are small-molecule compounds that can act as chemical chaperones by altering the environment in a cellular context. Despite their importance, detailed studies on the role of these chemical chaperones in modulating structure and dimensions of intrinsically disordered proteins have been limited. Here, we used single-molecule Förster resonance energy transfer to test the counteraction hypothesis of counterbalancing effects between the protecting osmolyte trimethylamine-N-oxide (TMAO) and denaturing osmolyte urea for the case of α-synuclein, a Parkinson's disease-linked protein whose monomer exhibits significant disorder. The single-molecule experiments, which avoid complications from protein aggregation, do not exhibit clear solvent-induced cooperative protein transitions for these osmolytes, unlike results from previous studies on globular proteins. Our data demonstrate the ability of TMAO and urea to shift α-synuclein structures towards either more compact or expanded average dimensions. Strikingly, the experiments directly reveal that a 21 [urea][TMAO] ratio has a net neutral effect on the protein's dimensions, a result that holds regardless of the absolute osmolyte concentrations. Our findings shed light on a surprisingly simple aspect of the interplay between urea and TMAO on α-synuclein in the context of intrinsically disordered proteins, with potential implications for the biological roles of such chemical chaperones. The results also highlight the strengths of single-molecule experiments in directly probing the chemical physics of protein structure and disorder in more chemically complex environments.

Slow unfolded-state structuring in Acyl-CoA binding protein folding revealed by simulation and experiment

Protein folding is a fundamental process in biology, key to understanding many human diseases. Experimentally, proteins often appear to fold via simple two- or three-state mechanisms involving mainly native-state interactions, yet recent network models built from atomistic simulations of small proteins suggest the existence of many possible metastable states and folding pathways. We reconcile these two pictures in a combined experimental and simulation study of acyl-coenzyme A binding protein (ACBP), a two-state folder (folding time ~10 ms) exhibiting residual unfolded-state structure, and a putative early folding intermediate. Using single-molecule FRET in conjunction with side-chain mutagenesis, we first demonstrate that the denatured state of ACBP at near-zero denaturant is unusually compact and enriched in long-range structure that can be perturbed by discrete hydrophobic core mutations. We then employ ultrafast laminar-flow mixing experiments to study the folding kinetics of ACBP on the microsecond time scale. These studies, along with Trp-Cys quenching measurements of unfolded-state dynamics, suggest that unfolded-state structure forms on a surprisingly slow (~100 μs) time scale, and that sequence mutations strikingly perturb both time-resolved and equilibrium smFRET measurements in a similar way. A Markov state model (MSM) of the ACBP folding reaction, constructed from over 30 ms of molecular dynamics trajectory data, predicts a complex network of metastable stables, residual unfolded-state structure, and kinetics consistent with experiment but no well-defined intermediate preceding the main folding barrier. Taken together, these experimental and simulation results suggest that the previously characterized fast kinetic phase is not due to formation of a barrier-limited intermediate but rather to a more heterogeneous and slow acquisition of unfolded-state structure.

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.