The unfolded state of the C-terminal domain of the ribosomal protein L9 contains both native and non-native structure

Residual Structure in the Denatured State of the Fast-Folding UBA(1) Domain from the Human DNA Excision Repair Protein HHR23A

The structure of the first ubiquitin-associated domain from HHR23A, UBA(1), was determined by X-ray crystallography at a 1.60 Å resolution, and its stability, folding kinetics, and residual structure under denaturing conditions have been investigated. The concentration dependence of thermal denaturation and size-exclusion chromatography indicate that UBA(1) is monomeric. Guanidine hydrochloride (GdnHCl) denaturation experiments reveal that the unfolding free energy, ΔGu°'(H2O), of UBA(1) is 2.4 kcal mol-1. Stopped-flow folding kinetics indicates sub-millisecond folding with only proline isomerization phases detectable at 25 °C. The full folding kinetics are observable at 4 °C, yielding a folding rate constant, kf, in the absence of a denaturant of 13,000 s-1 and a Tanford β-value of 0.80, consistent with a compact transition state. Evaluation of the secondary structure via circular dichroism shows that the residual helical structure in the denatured state is replaced by polyproline II structure as the GdnHCl concentration increases. Analysis of NMR secondary chemical shifts for backbone 15NH, 13CO, and 13Cα atoms between 4 and 7 M GdnHCl shows three islands of residual helical secondary structure that align in sequence with the three native-state helices. Extrapolation of the NMR data to 0 M GdnHCl demonstrates that helical structure would populate to 17-33% in the denatured state under folding conditions. Comparison with NMR data for a peptide corresponding to helix 1 indicates that this helix is stabilized by transient tertiary interactions in the denatured state of UBA(1). The high helical content in the denatured state, which is enhanced by transient tertiary interactions, suggests a diffusion-collision folding mechanism.

Helix-Coil Transition at a Glycine Following a Nascent α-Helix: A Synergetic Guidance Mechanism for Helix Growth

A detailed understanding of forces guiding the rapid folding of a polypeptide from an apparently random coil state to an ordered α-helical structure following the rate-limiting preorganization of the initial three residue backbones into helical conformation is imperative to comprehending and regulating protein folding and for the rational design of biological mimetics. However, several details of this process are still unknown. First, although the helix-coil transition was proposed to originate at the residue level (J. Chem. Phys. 1959, 31, 526-535; J. Chem. Phys. 1961, 34, 1963-1974), all helix-folding studies have only established it between time-averaged bulk states of a long-lived helix and several transiently populated random coils, along the whole helix model sequence. Second, the predominant thermodynamic forces driving either this two-state transition or the faster helix growth following helix nucleation are still unclear. Third, the conformational space of the random coil state is not well-defined unlike its corresponding α-helix. Here we investigate the restrictions placed on the conformational space of a Gly residue backbone, as a result of it immediately succeeding a nascent α-helical turn. Analyses of the temperature-dependent 1D-, 2D-NMR, FT-IR, and CD spectra and GROMACS MD simulation trajectory of a Gly residue backbone following a model α-helical turn, which is artificially rigidified by a covalent hydrogen bond surrogate, reveal that: (i) the α-helical turn guides the ϕ torsion of the Gly exclusively into either a predominantly populated entropically favored α-helical (α-ϕ) state or a scarcely populated random coil (RC-ϕ) state; (ii) the α-ϕ state of Gly in turn favors the stability of the preceding α-helical turn, while the RC-ϕ state disrupts it, revealing an entropy-driven synergetic guidance for helix growth in the residue following helix nucleation. The applicability of a current synergetic guidance mechanism to explain rapid helix growth in folded and unfolded states of proteins and helical peptides is discussed.

Modeling of Disordered Protein Structures Using Monte Carlo Simulations and Knowledge-Based Statistical Force Fields

The description of protein disordered states is important for understanding protein folding mechanisms and their functions. In this short review, we briefly describe a simulation approach to modeling protein interactions, which involve disordered peptide partners or intrinsically disordered protein regions, and unfolded states of globular proteins. It is based on the CABS coarse-grained protein model that uses a Monte Carlo (MC) sampling scheme and a knowledge-based statistical force field. We review several case studies showing that description of protein disordered states resulting from CABS simulations is consistent with experimental data. The case studies comprise investigations of protein⁻peptide binding and protein folding processes. The CABS model has been recently made available as the simulation engine of multiscale modeling tools enabling studies of protein⁻peptide docking and protein flexibility. Those tools offer customization of the modeling process, driving the conformational search using distance restraints, reconstruction of selected models to all-atom resolution, and simulation of large protein systems in a reasonable computational time. Therefore, CABS can be combined in integrative modeling pipelines incorporating experimental data and other modeling tools of various resolution.

Pressure-Temperature Analysis of the Stability of the CTL9 Domain Reveals Hidden Intermediates

The observation of two-state unfolding for many small single-domain proteins by denaturants has led to speculation that protein sequences may have evolved to limit the population of partially folded states that could be detrimental to fitness. How such strong cooperativity arises from a multitude of individual interactions is not well understood. Here, we investigate the stability and folding cooperativity of the C-terminal domain of the ribosomal protein L9 in the pressure-temperature plane using site-specific NMR. In contrast to apparent cooperative unfolding detected with denaturant-induced and thermal-induced unfolding experiments and stopped-flow refolding studies at ambient pressure, NMR-detected pressure unfolding revealed significant deviation from two-state behavior, with a core region that was selectively destabilized by increasing temperature. Comparison of pressure-dependent NMR signals from both the folded and unfolded states revealed the population of at least one invisible excited state at atmospheric pressure. The core destabilizing cavity-creating I98A mutation apparently increased the cooperativity of the loss of folded-state peak intensity while also increasing the population of this invisible excited state present at atmospheric pressure. These observations highlight how local stability is subtly modulated by sequence to tune protein conformational landscapes and illustrate the ability of pressure- and temperature-dependent studies to reveal otherwise hidden states.

Helical Propensity Affects the Conformational Properties of the Denatured State of Cytochrome c'

Changing the helical propensity of a polypeptide sequence might be expected to affect the conformational properties of the denatured state of a protein. To test this hypothesis, alanines at positions 83 and 87 near the center of helix 3 of cytochrome c' from Rhodopseudomonas palustris were mutated to serine to decrease the stability of this helix. A set of 13 single histidine variants in the A83S/A87S background were prepared to permit assessment of the conformational properties of the denatured state using histidine-loop formation in 3 M guanidine hydrochloride. The data are compared with previous histidine-heme loop formation data for wild-type cytochrome c'. As expected, destabilization of helix 3 decreases the global stabilities of the histidine variants in the A83S/A87S background relative to the wild-type background. Loop stability versus loop size data yields a scaling exponent of 2.1 ± 0.2, similar to the value of 2.3 ± 0.2 obtained for wild-type cytochrome c'. However, the stabilities of all histidine-heme loops, which contain the helix 3 sequence segment, are increased in the A83S/A87S background compared to the wild-type background. Rate constants for histidine-heme loop breakage are similar for the wild-type and A83S/A87S variants. However, for histidine-heme loops that contain the helix 3 sequence segment, the rate constants for loop formation increase in the A83S/A87S background compared to the wild-type background. Thus, residual helical structure appears to stiffen the polypeptide chain slowing loop formation in the denatured state. The implications of these results for protein folding mechanisms are discussed.

Effect of an Imposed Contact on Secondary Structure in the Denatured State of Yeast Iso-1-cytochrome c

There is considerable evidence that long-range interactions stabilize residual protein structure under denaturing conditions. However, evaluation of the effect of a specific contact on structure in the denatured state has been difficult. Iso-1-cytochrome c variants with a Lys54 → His mutation form a particularly stable His-heme loop in the denatured state, suggestive of loop-induced residual structure. We have used multidimensional nuclear magnetic resonance methods to assign ¹H and ¹⁵N backbone amide and ¹³C backbone and side chain chemical shifts in the denatured state of iso-1-cytochrome c carrying the Lys54 → His mutation in 3 and 6 M guanidine hydrochloride and at both pH 6.4, where the His54-heme loop is formed, and pH 3.6, where the His54-heme loop is broken. Using the secondary structure propensity score, with the 6 M guanidine hydrochloride chemical shift data as a random coil reference state for data collected in 3 M guanidine hydrochloride, we found residual helical structure in the denatured state for the 60s helix and the C-terminal helix, but not in the N-terminal helix in the presence or absence of the His54-heme loop. Non-native helical structure is observed in two regions that form Ω-loops in the native state. There is more residual helical structure in the C-terminal helix at pH 6.4 when the loop is formed. Loop formation also appears to stabilize helical structure near His54, consistent with induction of helical structure observed when His-heme bonds form in heme-peptide model systems. The results are discussed in the context of the folding mechanism of cytochrome c.

High Pressure ZZ-Exchange NMR Reveals Key Features of Protein Folding Transition States

Understanding protein folding mechanisms and their sequence dependence requires the determination of residue-specific apparent kinetic rate constants for the folding and unfolding reactions. Conventional two-dimensional NMR, such as HSQC experiments, can provide residue-specific information for proteins. However, folding is generally too fast for such experiments. ZZ-exchange NMR spectroscopy allows determination of folding and unfolding rates on much faster time scales, yet even this regime is not fast enough for many protein folding reactions. The application of high hydrostatic pressure slows folding by orders of magnitude due to positive activation volumes for the folding reaction. We combined high pressure perturbation with ZZ-exchange spectroscopy on two autonomously folding protein domains derived from the ribosomal protein, L9. We obtained residue-specific apparent rates at 2500 bar for the N-terminal domain of L9 (NTL9), and rates at atmospheric pressure for a mutant of the C-terminal domain (CTL9) from pressure dependent ZZ-exchange measurements. Our results revealed that NTL9 folding is almost perfectly two-state, while small deviations from two-state behavior were observed for CTL9. Both domains exhibited large positive activation volumes for folding. The volumetric properties of these domains reveal that their transition states contain most of the internal solvent excluded voids that are found in the hydrophobic cores of the respective native states. These results demonstrate that by coupling it with high pressure, ZZ-exchange can be extended to investigate a large number of protein conformational transitions.

How cooperative are protein folding and unfolding transitions?

A thermodynamically and kinetically simple picture of protein folding envisages only two states, native (N) and unfolded (U), separated by a single activation free energy barrier, and interconverting by cooperative two-state transitions. The folding/unfolding transitions of many proteins occur, however, in multiple discrete steps associated with the formation of intermediates, which is indicative of reduced cooperativity. Furthermore, much advancement in experimental and computational approaches has demonstrated entirely non-cooperative (gradual) transitions via a continuum of states and a multitude of small energetic barriers between the N and U states of some proteins. These findings have been instrumental towards providing a structural rationale for cooperative versus noncooperative transitions, based on the coupling between interaction networks in proteins. The cooperativity inherent in a folding/unfolding reaction appears to be context dependent, and can be tuned via experimental conditions which change the stabilities of N and U. The evolution of cooperativity in protein folding transitions is linked closely to the evolution of function as well as the aggregation propensity of the protein. A large activation energy barrier in a fully cooperative transition can provide the kinetic control required to prevent the accumulation of partially unfolded forms, which may promote aggregation. Nevertheless, increasing evidence for barrier-less "downhill" folding, as well as for continuous "uphill" unfolding transitions, indicate that gradual non-cooperative processes may be ubiquitous features on the free energy landscape of protein folding.

Quantitative assessments of the distinct contributions of polypeptide backbone amides versus side chain groups to chain expansion via chemical denaturation

In aqueous solutions with high concentrations of chemical denaturants such as urea and guanidinium chloride (GdmCl) proteins expand to populate heterogeneous conformational ensembles. These denaturing environments are thought to be good solvents for generic protein sequences because properties of conformational distributions align with those of canonical random coils. Previous studies showed that water is a poor solvent for polypeptide backbones, and therefore, backbones form collapsed globular structures in aqueous solvents. Here, we ask if polypeptide backbones can intrinsically undergo the requisite chain expansion in aqueous solutions with high concentrations of urea and GdmCl. We answer this question using a combination of molecular dynamics simulations and fluorescence correlation spectroscopy. We find that the degree of backbone expansion is minimal in aqueous solutions with high concentrations of denaturants. Instead, polypeptide backbones sample conformations that are denaturant-specific mixtures of coils and globules, with a persistent preference for globules. Therefore, typical denaturing environments cannot be classified as good solvents for polypeptide backbones. How then do generic protein sequences expand in denaturing environments? To answer this question, we investigated the effects of side chains using simulations of two archetypal sequences with amino acid compositions that are mixtures of charged, hydrophobic, and polar groups. We find that side chains lower the effective concentration of backbone amides in water leading to an intrinsic expansion of polypeptide backbones in the absence of denaturants. Additional dilution of the effective concentration of backbone amides is achieved through preferential interactions with denaturants. These effects lead to conformational statistics in denaturing environments that are congruent with those of canonical random coils. Our results highlight the role of side chain-mediated interactions as determinants of the conformational properties of unfolded states in water and in influencing chain expansion upon denaturation.

Advances in turbulent mixing techniques to study microsecond protein folding reactions

Recent experimental and computational advances in the protein folding arena have shown that the readout of the one-dimensional sequence information into three-dimensional structure begins within the first few microseconds of folding. The initiation of refolding reactions has been achieved by several means, including temperature jumps, flash photolysis, pressure jumps, and rapid mixing methods. One of the most commonly used means of initiating refolding of chemically denatured proteins is by turbulent flow mixing with refolding dilution buffer, where greater than 99% mixing efficiency has been achieved within 10's of microseconds. Successful interfacing of turbulent flow mixers with complementary detection methods, including time-resolved Fluorescence Spectroscopy (trFL), Förster Resonance Energy Transfer, Circular Dichroism, Small-Angle X-ray Scattering, Hydrogen Exchange followed by Mass Spectrometry and Nuclear Magnetic Resonance Spectroscopy, Infrared Spectroscopy (IR), and Fourier Transform IR Spectroscopy, has made this technique very attractive for monitoring various aspects of structure formation during folding. Although continuous-flow (CF) mixing devices interfaced with trFL detection have a dead time of only 30 µs, burst phases have been detected in this time scale during folding of peptides and of large proteins (e.g., CheY and TIM barrels). Furthermore, a major limitation of the CF mixing technique has been the requirement of large quantities of sample. In this brief communication, we will discuss the recent flurry of activity in micromachining and microfluidics, guided by computational simulations, which are likely to lead to dramatic improvements in time resolution and sample consumption for CF mixers over the next few years.

Residual structure in unfolded proteins

The denatured state ensemble (DSE) of unfolded proteins, once considered to be well-modeled by an energetically featureless random coil, is now well-known to contain flickering elements of residual structure. The position and nature of DSE residual structure may provide clues toward deciphering the protein folding code. This review focuses on recent advances in our understanding of the nature of DSE collapse under folding conditions, the quantification of the stability of residual structure in the DSE, the determination of the location and types of residues involved in thermodynamically significant residual structure and advances in detection of long-range interactions in the DSE.

Counting peptide-water hydrogen bonds in unfolded proteins

It is often assumed that the peptide backbone forms a substantial number of additional hydrogen bonds when a protein unfolds. We challenge that assumption in this article. Early surveys of hydrogen bonding in proteins of known structure typically found that most, but not all, backbone polar groups are satisfied, either by intramolecular partners or by water. When the protein is folded, these groups form approximately two hydrogen bonds per peptide unit, one donor or acceptor for each carbonyl oxygen or amide hydrogen, respectively. But when unfolded, the backbone chain is often believed to form three hydrogen bonds per peptide unit, one partner for each oxygen lone pair or amide hydrogen. This assumption is based on the properties of small model compounds, like N-methylacetamide, or simply accepted as self-evident fact. If valid, a chain of N residues would have approximately 2N backbone hydrogen bonds when folded but 3N backbone hydrogen bonds when unfolded, a sufficient difference to overshadow any uncertainties involved in calculating these per-residue averages. Here, we use exhaustive conformational sampling to monitor the number of H-bonds in a statistically adequate population of blocked polyalanyl-six-mers as the solvent quality ranges from good to poor. Solvent quality is represented by a scalar parameter used to Boltzmann-weight the population energy. Recent experimental studies show that a repeating (Gly-Ser) polypeptide undergoes a denaturant-induced expansion accompanied by breaking intramolecular peptide H-bonds. Results from our simulations augment this experimental finding by showing that the number of H-bonds is approximately conserved during such expansion⇋compaction transitions.

Minireview: structural insights into early folding events using continuous-flow time-resolved small-angle X-ray scattering

Small-angle X-ray scattering (SAXS) is a powerful method for obtaining quantitative structural information on the size and shape of proteins, and it is increasingly used in kinetic studies of folding and association reactions. In this minireview, we discuss recent developments in using SAXS to obtain structural information on the unfolded ensemble and early folding intermediates of proteins using continuous-flow mixing devices. Interfacing of these micromachined devices to SAXS beamlines has allowed access to the microsecond time regime. The experimental constraints in implementation of turbulence and laminar flow-based mixers with SAXS detection and a comparison of the two approaches are presented. Current improvements and future prospects of microsecond time-resolved SAXS and the synergy with ab initio structure prediction and molecular dynamics simulations are discussed.

Taming the complexity of protein folding

Protein folding is an important problem in structural biology with significant medical implications, particularly for misfolding disorders like Alzheimer's disease. Solving the folding problem will ultimately require a combination of theory and experiment, with theoretical models providing a comprehensive view of folding and experiments grounding these models in reality. Here we review progress towards this goal over the past decade, with an emphasis on recent theoretical advances that are empowering chemically detailed models of folding and the new results these technologies are providing. In particular, we discuss new insights made possible by Markov state models (MSMs), including the role of non-native contacts and the hub-like character of protein folded states.

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.

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

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

The cold denatured state of the C-terminal domain of protein L9 is compact and contains both native and non-native structure

Cold denaturation is a general property of globular proteins, and the process provides insight into the origins of the cooperativity of protein folding and the nature of partially folded states. Unfortunately, studies of protein cold denaturation have been hindered by the fact that the cold denatured state is normally difficult to access experimentally. Special conditions such as addition of high concentrations of denaturant, encapsulation into reverse micelles, the formation of emulsified solutions, high pressure, or extremes of pH have been applied, but these can perturb the unfolded state of proteins. The cold denatured state of the C-terminal domain of the ribosomal protein L9 can be populated under native-like conditions by taking advantage of a destabilizing point mutation which leads to cold denaturation at temperatures above 0 degrees C. This state is in slow exchange with the native state on the NMR time scale. Virtually complete backbone (15)N, (13)C, and (1)H as well as side-chain (13)C(beta) and (1)H(beta) chemical shift assignments were obtained for the cold denatured state at pH 5.7, 12 degrees C. Chemical shift analysis, backbone N-H residual dipolar couplings, amide proton NOEs, and R(2) relaxation rates all indicate that the cold denatured state of CTL9 (the C-terminal domain of the ribosomal protein L9) not only contains significant native-like secondary structure but also non-native structure. The regions corresponding to the two native alpha-helices show a strong tendency to populate helical Phi and Psi angles. The segment which connects alpha-helix 2 and beta-strand 2 (residues 107-124) in the native state exhibits a significant preference to form non-native helical structure in the cold denatured state. The structure observed in the cold denatured state of the I98A mutant is similar to that observed in the pH 3.8 unfolded state of wild type CTL9 at 25 degrees C, suggesting that it is a robust feature of the denatured state ensemble of this protein. The implications for protein folding and for studies of cold denatured states are discussed.

Cooperative formation of native-like tertiary contacts in the ensemble of unfolded states of a four-helix protein

In studies of the ensembles of unfolded structures of a four-helix bundle protein, we have detected the presence of potential precursors of native tertiary structures. These observations were based on the perturbation of NMR chemical shifts of the protein backbone atoms by single site mutations. Some mutations change the chemical shifts of residues remote from the site of mutation indicating the presence of an interaction between the mutated and the remote residues, suggesting that the formation of helix segments and helix-helix interactions is cooperative. We can begin to track down the folding mechanism of this protein using only experimental data by combining the information available for the rate limiting structure formation during the folding process with measurements of the site specific hydrogen bond formation in the burst phase, and with the existence prior to the folding reaction of tertiary structures in the ensemble of otherwise unfolded structures observed in the present study.

Protein folded states are kinetic hubs

Understanding molecular kinetics, and particularly protein folding, is a classic grand challenge in molecular biophysics. Network models, such as Markov state models (MSMs), are one potential solution to this problem. MSMs have recently yielded quantitative agreement with experimentally derived structures and folding rates for specific systems, leaving them positioned to potentially provide a deeper understanding of molecular kinetics that can lead to experimentally testable hypotheses. Here we use existing MSMs for the villin headpiece and NTL9, which were constructed from atomistic simulations, to accomplish this goal. In addition, we provide simpler, humanly comprehensible networks that capture the essence of molecular kinetics and reproduce qualitative phenomena like the apparent two-state folding often seen in experiments. Together, these models show that protein dynamics are dominated by stochastic jumps between numerous metastable states and that proteins have heterogeneous unfolded states (many unfolded basins that interconvert more rapidly with the native state than with one another) yet often still appear two-state. Most importantly, we find that protein native states are hubs that can be reached quickly from any other state. However, metastability and a web of nonnative states slow the average folding rate. Experimental tests for these findings and their implications for other fields, like protein design, are also discussed.

Urea denatured state ensembles contain extensive secondary structure that is increased in hydrophobic proteins

The goal of this article is to gain a better understanding of the denatured state ensemble (DSE) of proteins through an experimental and computational study of their denaturation by urea. Proteins unfold to different extents in urea and the most hydrophobic proteins have the most compact DSE and contain almost as much secondary structure as folded proteins. Proteins that unfold to the greatest extent near pH 7 still contain substantial amounts of secondary structure. At low pH, the DSE expands due to charge-charge interactions and when the net charge per residue is high, most of the secondary structure is disrupted. The proteins in the DSE appear to contain substantial amounts of polyproline II conformation at high urea concentrations. In all cases considered, including staph nuclease, the extent of unfolding by urea can be accounted for using the data and approach developed in the laboratory of Wayne Bolen (Auton et al., Proc Natl Acad Sci 2007; 104:15317-15323).