Chaperonin chamber accelerates protein folding through passive action of preventing aggregation

The original experiments reconstituting GroEL–GroES-mediated protein folding were carried out under “nonpermissive” conditions, where the chaperonin system was absolutely required and substrate proteins could not achieve the native state if diluted directly from denaturant into solution. Under “permissive” conditions, however, employing lower substrate concentration and lower temperature, some substrate proteins can be refolded both by the chaperonin system and while free in solution. For several of these, the protein refolds more rapidly inside the GroEL–GroES cis chamber than free in solution, suggesting that the chamber may have an active role in assisting protein folding. Here, we observe that the difference is caused by reversible multimolecular association while folding in solution, an avenue of kinetic partitioning that slows the overall rate of renaturation relative to the chaperonin chamber, where such associations cannot occur. For Rubisco, reversible aggregation during folding in solution was observed by gel filtration. For a mutant of maltose-binding protein (DM-MBP), the rate of folding in solution declined with increasing concentration, and the folding reaction produced light scattering. Under solution conditions where chloride was absent, however, light scattering no longer occurred, and DM-MBP folded at the same rate as in the cis cavity. In a further test, dihydrofolate reductase, thermally inactivated in the cis cavity or in solution, was substantially reactivated upon temperature downshift in the cis cavity but not in solution, where aggregation occurred. We conclude that the GroEL–GroES chamber behaves as a passive “Anfinsen cage” whose primary role is to prevent multimolecular association during folding.     

Functional modulation of a protein folding landscape via side-chain distortion

Ultrahigh-resolution (< 1.0 Å) structures have revealed unprecedented and unexpected details of molecular geometry, such as the deformation of aromatic rings from planarity. However, the functional utility of such energetically costly strain is unknown. The 0.83 Å structure of α-lytic protease (αLP) indicated that residues surrounding a conserved Phe side-chain dictate a rotamer which results in a ∼6° distortion along the side-chain, estimated to cost 4 kcal/mol. By contrast, in the closely related protease Streptomyces griseus Protease B (SGPB), the equivalent Phe adopts a different rotamer and is undistorted. Here, we report that the αLP Phe side-chain distortion is both functional and conserved in proteases with large pro regions. Sequence analysis of the αLP serine protease family reveals a bifurcation separating those sequences expected to induce distortion and those that would not, which correlates with the extent of kinetic stability. Structural and folding kinetics analyses of family members suggest that distortion of this side-chain plays a role in increasing kinetic stability within the αLP family members that use a large Pro region. Additionally, structural and kinetic folding studies of mutants demonstrate that strain alters the folding free energy landscape by destabilizing the transition state (TS) relative to the native state (N). Although side-chain distortion comes at a cost of foldability, it suppresses the rate of unfolding, thereby enhancing kinetic stability and increasing protein longevity under harsh extracellular conditions. This ability of a structural distortion to enhance function is unlikely to be unique to αLP family members and may be relevant in other proteins exhibiting side-chain distortions.     

The nature of protein folding pathways

How do proteins fold, and why do they fold in that way? This Perspective integrates earlier and more recent advances over the 50-y history of the protein folding problem, emphasizing unambiguously clear structural information. Experimental results show that, contrary to prior belief, proteins are multistate rather than two-state objects. They are composed of separately cooperative foldon building blocks that can be seen to repeatedly unfold and refold as units even under native conditions. Similarly, foldons are lost as units when proteins are destabilized to produce partially unfolded equilibrium molten globules. In kinetic folding, the inherently cooperative nature of foldons predisposes the thermally driven amino acid-level search to form an initial foldon and subsequent foldons in later assisted searches. The small size of foldon units, ∼20 residues, resolves the Levinthal time-scale search problem. These microscopic-level search processes can be identified with the disordered multitrack search envisioned in the “new view” model for protein folding. Emergent macroscopic foldon–foldon interactions then collectively provide the structural guidance and free energy bias for the ordered addition of foldons in a stepwise pathway that sequentially builds the native protein. These conclusions reconcile the seemingly opposed new view and defined pathway models; the two models account for different stages of the protein folding process. Additionally, these observations answer the “how” and the “why” questions. The protein folding pathway depends on the same foldon units and foldon–foldon interactions that construct the native structure.     

Redrawing the Ramachandran plot after inclusion of hydrogen-bonding constraints

A protein backbone has two degrees of conformational freedom per residue, described by its φ,ψ-angles. Accordingly, the energy landscape of a blocked peptide unit can be mapped in two dimensions, as shown by Ramachandran, Sasisekharan, and Ramakrishnan almost half a century ago. With atoms approximated as hard spheres, the eponymous Ramachandran plot demonstrated that steric clashes alone eliminate ¾ of φ,ψ-space, a result that has guided all subsequent work. Here, we show that adding hydrogen-bonding constraints to these steric criteria eliminates another substantial region of φ,ψ-space for a blocked peptide; for conformers within this region, an amide hydrogen is solvent-inaccessible, depriving it of a hydrogen-bonding partner. Yet, this “forbidden” region is well populated in folded proteins, which can provide longer-range intramolecular hydrogen-bond partners for these otherwise unsatisfied polar groups. Consequently, conformational space expands under folding conditions, a paradigm-shifting realization that prompts an experimentally verifiable conjecture about likely folding pathways.     

Transiently populated intermediate functions as a branching point of the FF domain folding pathway

Studies of protein folding and the intermediates that are formed along the folding pathway provide valuable insights into the process by which an unfolded ensemble forms a functional native conformation. However, because intermediates on folding pathways can serve as initiation points of aggregation (implicated in a number of diseases), their characterization assumes an even greater importance. Establishing the role of such intermediates in folding, misfolding, and aggregation remains a major challenge due to their often low populations and short lifetimes. We recently used NMR relaxation dispersion methods and computational techniques to determine an atomic resolution structure of the folding intermediate of a small protein module—the FF domain—with an equilibrium population of 2–3% and a millisecond lifetime, 25 °C. Based on this structure a variant FF domain has been designed in which the native state is selectively destabilized by removing the carboxyl-terminal helix in the native structure to produce a highly populated structural mimic of the intermediate state. Here, we show via solution NMR studies of the designed mimic that the mimic forms distinct conformers corresponding to monomeric and dimeric (K_d = 0.2 mM) forms of the protein. The conformers exchange on the seconds timescale with a monomer association rate of 1.1·10⁴ M^-1 s^-1 and with a region responsible for dimerization localized to the amino-terminal residues of the FF domain. This study establishes the FF domain intermediate as a central player in both folding and misfolding pathways and illustrates how incomplete folding can lead to the formation of higher-order structures.     

Direct observation of multiple misfolding pathways in a single prion protein molecule

Protein misfolding is a ubiquitous phenomenon associated with a wide range of diseases. Single-molecule approaches offer a powerful tool for deciphering the mechanisms of misfolding by measuring the conformational fluctuations of a protein with high sensitivity. We applied single-molecule force spectroscopy to observe directly the misfolding of the prion protein PrP, a protein notable for having an infectious misfolded state that is able to propagate by recruiting natively folded PrP. By measuring folding trajectories of single PrP molecules held under tension in a high-resolution optical trap, we found that the native folding pathway involves only two states, without evidence for partially folded intermediates that have been proposed to mediate misfolding. Instead, frequent but fleeting transitions were observed into off-pathway intermediates. Three different misfolding pathways were detected, all starting from the unfolded state. Remarkably, the misfolding rate was even higher than the rate for native folding. A mutant PrP with higher aggregation propensity showed increased occupancy of some of the misfolded states, suggesting these states may act as intermediates during aggregation. These measurements of individual misfolding trajectories demonstrate the power of single-molecule approaches for characterizing misfolding directly by mapping out nonnative folding pathways.     

Direct observation of ultrafast folding and denatured state dynamics in single protein molecules

Single-molecule fluorescence resonance energy transfer (smFRET) experiments are extremely useful in studying protein folding but are generally limited to time scales of greater than ≈100 μs and distances greater than ≈2 nm. We used single-molecule fluorescence quenching by photoinduced electron transfer, detecting short-range events, in combination with fluorescence correlation spectroscopy (PET-FCS) to investigate folding dynamics of the small binding domain BBL with nanosecond time resolution. The kinetics of folding appeared as a 10-μs decay in the autocorrelation function, resulting from stochastic fluctuations between denatured and native conformations of individual molecules. The observed rate constants were probe independent and in excellent agreement with values derived from conventional temperature-jump (T-jump) measurements. A submicrosecond relaxation was detected in PET-FCS data that reported on the kinetics of intrachain contact formation within the thermally denatured state. We engineered a mutant of BBL that was denatured under the reaction conditions that favored folding of the parent wild type (“D_phys”). D_phys had the same kinetic signature as the thermally denatured state and revealed segmental diffusion with a time constant of intrachain contact formation of 500 ns. This time constant was more than 10 times faster than folding and in the range estimated to be the “speed limit” of folding. D_phys exhibited significant deviations from a random coil. The solvent viscosity and temperature dependence of intrachain diffusion showed that chain motions were slaved by the presence of intramolecular interactions. PET-FCS in combination with protein engineering is a powerful approach to study the early events and mechanism of ultrafast protein folding.     

Fibrillation precursor of superoxide dismutase 1 revealed by gradual tuning of the protein-folding equilibrium

Although superoxide dismutase 1 (SOD1) stands out as a relatively soluble protein in vitro, it can be made to fibrillate by mechanical agitation. The mechanism of this fibrillation process is yet poorly understood, but attains considerable interest due to SOD1’s involvement in the neurodegenerative disease amyotrophic lateral sclerosis (ALS). In this study, we map out the apoSOD1 fibrillation process from how it competes with the global folding events at increasing concentrations of urea: We determine how the fibrillation lag time (τ_lag) and maximum growth rate (ν_max) depend on gradual titration of the folding equilibrium, from the native to the unfolded state. The results show that the agitation-induced fibrillation of apoSOD1 uses globally unfolded precursors and relies on fragmentation-assisted growth. Mutational screening and fibrillation m-values (∂ log τ_lag/∂[urea] and ∂ log ν_max/∂[urea]) indicate moreover that the fibrillation pathway proceeds via a diffusely bound transient complex that responds to the global physiochemical properties of the SOD1 sequence. Fibrillation of apoSOD1, as it bifurcates from the denatured ensemble, seems thus mechanistically analogous to that of disordered peptides, save the competing folding transition to the native state. Finally, we examine by comparison with in vivo data to what extent this mode of fibrillation, originating from selective amplification of mechanically brittle aggregates by sample agitation, captures the mechanism of pathological SOD1 aggregation in ALS.     

Multimolecule test-tube simulations of protein unfolding and aggregation

Molecular dynamics simulations of protein folding or unfolding, unlike most in vitro experimental methods, are performed on a single molecule. The effects of neighboring molecules on the unfolding/folding pathway are largely ignored experimentally and simply not modeled computationally. Here, we present two all-atom, explicit solvent molecular dynamics simulations of 32 copies of the Engrailed homeodomain (EnHD), an ultrafast-folding and -unfolding protein for which the folding/unfolding pathway is well-characterized. These multimolecule simulations, in comparison with single-molecule simulations and experimental data, show that intermolecular interactions have little effect on the folding/unfolding pathway. EnHD unfolded by the same mechanism whether it was simulated in only water or also in the presence of other EnHD molecules. It populated the same native state, transition state, and folding intermediate in both simulation systems, and was in good agreement with experimental data available for each of the three states. Unfolding was slowed slightly by interactions with neighboring proteins, which were mostly hydrophobic in nature and ultimately caused the proteins to aggregate. Protein–water hydrogen bonds were also replaced with protein–protein hydrogen bonds, additionally contributing to aggregation. Despite the increase in protein–protein interactions, the protein aggregates formed in simulation did not do so at the total exclusion of water. These simulations support the use of single-molecule techniques to study protein unfolding and also provide insight into the types of interactions that occur as proteins aggregate at high temperature at an atomic level.     

Exploring the folding energy landscape of a series of designed consensus tetratricopeptide repeat proteins

Repeat proteins contain short, tandem arrays of simple structural motifs (20−40 aa). These stack together to form nonglobular structures that are stabilized by short-range interactions from residues close in primary sequence. Unlike globular proteins, they have few, if any, long-range nonlocal stabilizing interactions. One ubiquitous repeat is the tetratricopeptide motif (TPR), a 34-aa helix-turn-helix motif. In this article we describe the folding kinetics of a series of 7 designed TPR proteins that are assembled from arraying identical designed consensus repeats (CTPRan). These range from the smallest 2-repeat protein to a large 10-repeat protein (≈350 aa). In particular, we describe how the energy landscape changes with the addition of repeat units. The data reveal that although the CTPRa proteins have low local frustration, their highly symmetric, modular native structure is reflected in their multistate kinetics of unfolding and folding. Moreover, although the initial folding of all CTPRan proteins involves a nucleus with similar solvent accessibility, their subsequent folding to the native structure depends directly on repeat number. This corresponds to an increasingly complex landscape that culminates in CTPRa10 populating a misfolded, off-pathway intermediate. These results extend our current understanding of the malleable folding pathways of repeat proteins and highlight the consequences of adding identical repeats to the energy landscape.     

Folding domain B of protein A on a dynamically partitioned free energy landscape

The B domain of staphylococcal protein A (BdpA) is a small helical protein that has been studied intensively in kinetics experiments and detailed computer simulations that include explicit water. The simulations indicate that BdpA needs to reorganize in crossing the transition barrier to facilitate folding its C-terminal helix (H3) onto the nucleus formed from helices H1 and H2. This process suggests frustration between two partially ordered forms of the protein, but recent φ value measurements indicate that the transition structure is relatively constant over a broad range of temperatures. Here we develop a simplistic model to investigate the folding transition in which properties of the free energy landscape can be quantitatively compared with experimental data. The model is a continuation of the Muñoz–Eaton model to include the intermittency of contacts between structured parts of the protein, and the results compare variations in the landscape with denaturant and temperature to φ value measurements and chevron plots of the kinetic rates. The topography of the model landscape (in particular, the feature of frustration) is consistent with detailed simulations even though variations in the φ values are close to measured values. The transition barrier is smaller than indicated by the chevron data, but it agrees in order of magnitude with a similar α-carbon type of model. Discrepancies with the chevron plots are investigated from the point of view of solvent effects, and an approach is suggested to account for solvent participation in the model.     

Heat shock protein 70 kDa chaperone/DnaJ cochaperone complex employs an unusual dynamic interface

The heat shock protein 70 kDa (Hsp70)/DnaJ/nucleotide exchange factor system assists in intracellular protein (re)folding. Using solution NMR, we obtained a three-dimensional structure for a 75-kDa Hsp70–DnaJ complex in the ADP state, loaded with substrate peptide. We establish that the J domain (residues 1–70) binds with its positively charged helix II to a negatively charged loop in the Hsp70 nucleotide-binding domain. The complex shows an unusual “tethered” binding mode which is stoichiometric and saturable, but which has a dynamic interface. The complex represents part of a triple complex of Hsp70 and DnaJ both bound to substrate protein. Mutagenesis data indicate that the interface is also of relevance for the interaction of Hsp70 and DnaJ in the ATP state. The solution complex is completely different from a crystal structure of a disulfide-linked complex of homologous proteins [Jiang, et al. (2007) Mol Cell 28:422–433].     

Measuring hydrogen exchange rates in invisible protein excited states

Hydrogen exchange rates have become a valuable probe for studying the relationship between dynamics and structure and for dissecting the mechanism by which proteins fold to their native conformation. Typically measured rates correspond to averages over all protein states from which hydrogen exchange can occur. Here we describe a new NMR experiment based on chemical exchange saturation transfer that provides an avenue for obtaining uncontaminated, per-residue amide hydrogen exchange rates for interconverting native and invisible states so long as they can be separated on the basis of distinct ¹⁵N chemical shifts. The approach is applied to the folding reaction of the Fyn SH3 domain that exchanges between a highly populated, NMR-visible native state and a conformationally excited, NMR-invisible state, corresponding to the unfolded ensemble. Excellent agreement between experimentally derived hydrogen exchange rates of the excited state at a pair of pHs is obtained, taking into account the expected dependence of exchange on pH. Extracted rates for the unfolded ensemble have been used to test hydrogen exchange predictions based on the primary protein sequence that are used in many analyses of solvent exchange rates, with a Pearson correlation coefficient of 0.84 obtained.The energy landscape of a protein is a multidimensional surface composed of many local minima in addition to the global minimum that is the native conformation (1, 2). An understanding of the relation between protein structure, dynamics, and function is predicated, therefore, on an analysis of the various conformational states that populate the minima on the landscape. This requires detailed structural and dynamics studies of each of the conformers and quantification of their relative energies as well as the kinetics of exchange between them. Biophysical techniques such as X-ray diffraction and NMR spectroscopy are available for obtaining detailed structural information on the molecules populating the lowest-energy regions of the landscape, providing insight into the structure–function paradigm for a great number of proteins. However, it is becoming increasingly well understood that studies of the ground states of proteins are not sufficient. Additional states that can be sparsely populated and transiently formed, referred to as excited states in what follows, are often important for processes that include molecular recognition, ligand binding, enzyme catalysis, and folding (3–9). Detailed studies of such excited states are, however, challenging because they are not “visible” to standard biophysical methods and as a consequence atomic resolution information is lacking.Recent developments in solution NMR spectroscopy are changing this paradigm by providing an avenue for quantifying excited states at a level of detail that has typically only been possible in studies of highly populated native protein conformers (10). Backbone ¹H, ¹³C, and ¹⁵N chemical shifts for invisible protein states can now be measured so long as they are populated at levels of 0.5% or higher and exchange with an NMR visible state with rates on the order of approximately 100 to several thousand per second (10–12). These chemical shifts have, in turn, been used along with database computational approaches to generate atomic resolution models of excited states (7, 13, 14), providing detailed insights into a range of biochemical processes. It is now possible in some cases to measure scalar and residual dipolar couplings that, in turn, provide additional structural insights (15). Further advances have led to the measurement of side-chain ¹³C and ¹H chemical shifts in invisible states (16–18) that can be sensitive to hydrophobic contacts in these rare conformers. Experiments for quantifying pico- to nanosecond time-scale side-chain dynamics have also emerged (19), showing in some cases large differences in motion between ground and excited states that directly relate to function (20).Although the tool kit of NMR experiments for studying rare protein conformers is expanding, it remains significantly smaller than that for highly populated states, and there is a continuing need for the development of further methodologies. Notably absent from the tool kit is an approach for measuring hydrogen exchange rates in excited-state protein conformers. Ever since the pioneering work of Linderstrøm-Lang in the 1950s it has been recognized that the rates of exchange of amide hydrogens with solvent protons are important site-specific parameters of protein structure and dynamics (21). More recently, amide hydrogen exchange measurements have been used to investigate regions of local protein disorder, folding/unfolding processes, hydrogen bonding, allostery, and ligand binding (22–24). NMR spectroscopy has emerged as a primary tool for quantifying amide hydrogen exchange rates, k_H-EX, on a per-residue basis. However, because the overall transfer of hydrogens from solvent is measured, the extracted rates often contain contributions from all of the exchange-accessible states of the protein (25–27). This is the case for H–D exchange experiments where time-dependent intensity changes of NMR signals are measured to extract k_H-EX (22), as well as experiments that perturb the water signal and quantify the perturbation at exchangeable amide sites (26). Similar bulk measures of k_H-EX are obtained using other methods as well (28), from which insight into a particular excited state can often only be inferred. Herein we describe an approach by which site-specific hydrogen exchange rates of amides in individual states along the energy landscape can be obtained for a protein system at equilibrium and under native conditions. An application to the folding reaction of the G48A Fyn SH3 domain is provided, where k_H-EX values of the native conformation (ground state) and an excited state (unfolded ensemble) are measured. Because exchange rates are quantified by measurement of intensities of cross-peaks in NMR spectra of the ground state, the method is sensitive to k_H-EX values as large as several hundred per second for the excited state. The measurement of “pure” k_H-EX values provides a powerful new approach for understanding the dynamical properties of excited protein states.     

Quantifying internal friction in unfolded and intrinsically disordered proteins with single-molecule spectroscopy

Internal friction, which reflects the “roughness” of the energy landscape, plays an important role for proteins by modulating the dynamics of their folding and other conformational changes. However, the experimental quantification of internal friction and its contribution to folding dynamics has remained challenging. Here we use the combination of single-molecule Förster resonance energy transfer, nanosecond fluorescence correlation spectroscopy, and microfluidic mixing to determine the reconfiguration times of unfolded proteins and investigate the mechanisms of internal friction contributing to their dynamics. Using concepts from polymer dynamics, we determine internal friction with three complementary, largely independent, and consistent approaches as an additive contribution to the reconfiguration time of the unfolded state. We find that the magnitude of internal friction correlates with the compactness of the unfolded protein: its contribution dominates the reconfiguration time of approximately 100 ns of the compact unfolded state of a small cold shock protein under native conditions, but decreases for more expanded chains, and approaches zero both at high denaturant concentrations and in intrinsically disordered proteins that are expanded due to intramolecular charge repulsion. Our results suggest that internal friction in the unfolded state will be particularly relevant for the kinetics of proteins that fold in the microsecond range or faster. The low internal friction in expanded intrinsically disordered proteins may have implications for the dynamics of their interactions with cellular binding partners.     

Mapping backbone and side-chain interactions in the transition state of a coupled protein folding and binding reaction

Understanding the mechanism of protein folding requires a detailed knowledge of the structural properties of the barriers separating unfolded from native conformations. The S-peptide from ribonuclease S forms its α-helical structure only upon binding to the folded S-protein. We characterized the transition state for this binding-induced folding reaction at high resolution by determining the effect of site-specific backbone thioxylation and side-chain modifications on the kinetics and thermodynamics of the reaction, which allows us to monitor formation of backbone hydrogen bonds and side-chain interactions in the transition state. The experiments reveal that α-helical structure in the S-peptide is absent in the transition state of binding. Recognition between the unfolded S-peptide and the S-protein is mediated by loosely packed hydrophobic side-chain interactions in two well defined regions on the S-peptide. Close packing and helix formation occurs rapidly after binding. Introducing hydrophobic residues at positions outside the recognition region can drastically slow down association.     

The dual-basin landscape in GFP folding

Recent experimental studies suggest that the mature GFP has an unconventional landscape composed of an early folding event with a typical funneled landscape, followed by a very slow search and rearrangement step into the locked, active chromophore-containing structure. As we have shown previously, the substantial difference in time scales is what generates the observed hysteresis in thermodynamic folding. The interconversion between locked and the soft folding structures at intermediate denaturant concentrations is so slow that it is not observed under the typical experimental observation time. Simulations of a coarse-grained model were used to describe the fast folding event as well as identify native-like intermediates on energy landscapes enroute to the fluorescent native fold. Interestingly, these simulations reveal structural features of the slow dynamic transition to chromophore activation. Experimental evidence presented here shows that the trapped, native-like intermediate has structural heterogeneity in residues previously linked to chromophore formation. We propose that the final step of GFP folding is a “locking” mechanism leading to chromophore formation and high stability. The combination of previous experimental work and current simulation work is explained in the context of a dual-basin folding mechanism described above.     

Golden triangle for folding rates of globular proteins

The ability of protein chains to spontaneously form their spatial structures is a long-standing puzzle in molecular biology. Experimentally measured rates of spontaneous folding of single-domain globular proteins range from microseconds to hours: the difference (11 orders of magnitude) is akin to the difference between the life span of a mosquito and the age of the universe. Here, we show that physical theory with biological constraints outlines a “golden triangle” limiting the possible range of folding rates for single-domain globular proteins of various size and stability, and that the experimentally measured folding rates fall within this narrow triangle built without any adjustable parameters, filling it almost completely. In addition, the golden triangle predicts the maximal size of protein domains that fold under solely thermodynamic (rather than kinetic) control. It also predicts the maximal allowed size of the “foldable” protein domains, and the size of domains found in known protein structures is in a good agreement with this limit.