Protein folding at single-molecule resolution

The protein folding reaction carries great significance for cellular function and hence continues to be the research focus of a large interdisciplinary protein science community. Single-molecule methods are providing new and powerful tools for dissecting the mechanisms of this complex process by virtue of their ability to provide views of protein structure and dynamics without associated ensemble averaging. This review briefly introduces common FRET and force methods, and then explores several areas of protein folding where single-molecule experiments have yielded insights. These include exciting new information about folding landscapes, dynamics, intermediates, unfolded ensembles, intrinsically disordered proteins, assisted folding and biomechanical unfolding. Emerging and future work is expected to include advances in single-molecule techniques aimed at such investigations, and increasing work on more complex systems from both the physics and biology standpoints, including folding and dynamics of systems of interacting proteins and of proteins in cells and organisms. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.     

Membrane proteins enter the fold

Membrane proteins have historically been recalcitrant to biophysical folding studies. However, recent adaptations of methods from the soluble protein folding field have found success in their applications to transmembrane proteins composed of both α-helical and β-barrel conformations. Avoiding aggregation is critical for the success of these experiments. Altogether these studies are leading to discoveries of folding trajectories, foundational stabilizing forces and better-defined endpoints that enable more accurate interpretation of thermodynamic data. Increased information on membrane protein folding in the cell shows that the emerging biophysical principles are largely recapitulated even in the complex biological environment.     

Effects of macromolecular crowding agents on protein folding in vitro and in silico

Proteins fold and function inside cells which are environments very different from that of dilute buffer solutions most often used in traditional experiments. The crowded milieu results in excluded-volume effects, increased bulk viscosity and amplified chances for inter-molecular interactions. These environmental factors have not been accounted for in most mechanistic studies of protein folding executed during the last decades. The question thus arises as to how these effects—present when polypeptides normally fold in vivo—modulate protein biophysics. To address excluded volume effects, we use synthetic macromolecular crowding agents, which take up significant volume but do not interact with proteins, in combination with strategically selected proteins and a range of equilibrium and time-resolved biophysical (spectroscopic and computational) methods. In this review, we describe key observations on macromolecular crowding effects on protein stability, folding and structure drawn from combined in vitro and in silico studies. As expected based on Minton’s early predictions, many proteins (apoflavodoxin, VlsE, cytochrome c, and S16) became more thermodynamically stable (magnitude depends inversely on protein stability in buffer) and, unexpectedly, for apoflavodoxin and VlsE, the folded states changed both secondary structure content and, for VlsE, overall shape in the presence of macromolecular crowding. For apoflavodoxin and cytochrome c, which have complex kinetic folding mechanisms, excluded volume effects made the folding energy landscapes smoother (i.e., less misfolding and/or kinetic heterogeneity) than in buffer.     

Rapid mixing methods for exploring the kinetics of protein folding

Information on the time-dependence of molecular species is critical for elucidating reaction mechanisms in chemistry and biology. Rapid flow experiments involving turbulent mixing of two or more solutions continue to be the main source of kinetic information on protein folding and other biochemical processes, such as ligand binding and enzymatic reactions. Recent advances in mixer design and detection methods have opened a new window for exploring conformational changes in proteins on the microsecond time scale. These developments have been especially important for exploring early stages of protein folding.     

Protein folding and binding in confined spaces and in crowded solutions

Simple theoretical models are presented to illustrate the effects of spatial confinement and macromolecular crowding on the equilibria and rates of protein folding and binding. Confinement is expected to significantly stabilize the folded state, but for crowding only a marginal effect on protein stability is expected. In confinement the unfolded chain is restricted to a cage but in crowding the unfolded chain may explore different interstitial voids. Because confinement and crowding eliminate the more expanded conformations of the unfolded state, folding from the compact unfolded state is expected to speed up. Crowding will shift the binding equilibrium of proteins toward the bound state. The significant slowing down in protein diffusion by crowding, perhaps beneficial for chaperonin action, could result in a decrease in protein binding rates.     

Desolvation is a likely origin of robust enthalpic barriers to protein folding

Experimental data from global analyses of temperature (T) and denaturant dependence of the folding rates of small proteins led to an intrinsic enthalpic folding barrier hypothesis: to a good approximation, the T-dependence of folding rate under constant native stability conditions is Arrhenius. Furthermore, for a given protein, the slope of isostability folding rate versus 1/T is essentially independent of native stability. This hypothesis implies a simple relationship between chevron and Eyring plots of folding that is easily discernible when both sets of rates are expressed as functions of native stability. Using experimental data in the literature, we verify the predicted chevron-Eyring relationship for 14 proteins and determine their intrinsic enthalpic folding barriers, which vary approximately from 15 kcal/mol to 40 kcal/mol for different proteins. These enthalpic barriers do not appear to correlate with folding rates, but they exhibit correlation with equilibrium unfolding enthalpy at room temperature. Intrinsic enthalpic barriers with similarly high magnitudes apply as well to at least two cases of peptide-peptide and peptide-protein association, suggesting that these barriers are a hallmark of certain general and fundamental kinetic processes during folding and binding. Using a class of explicit-chain C(alpha) protein models with constant elementary enthalpic desolvation barriers between C(alpha) positions, we show that small microscopic pairwise desolvation barriers, which are a direct consequence of the particulate nature of water, can act cooperatively to give rise to a significant overall enthalpic barrier to folding. This theoretical finding provides a physical rationalization for the high intrinsic enthalpic barriers in protein folding energetics. Ramifications of entropy-enthalpy compensation in hydrophobic association for the height of enthalpic desolvation barrier are discussed.     

Computational design of proteins stereochemically optimized in size, stability, and folding speed

Artificial proteins potentially barrier-free in the folding kinetics are approached computationally under the guidance of protein-folding theories. The smallest and fastest folding globular protein triple-helix-bundle (THB) is so modified as to minimize or eliminate its presumed barriers in folding speed. As the barriers may reside in the ordering of either secondary or tertiary structure, the elements of both secondary and tertiary structure in the protein are targeted for prenucleation with suitable stereochemically constrained amino acid residues. The required elements of topology and sequence for the THB are optimized independently; first the topology is optimized with simulated annealing in polypeptides of highly simplified alphabet; next, the sequence in side chains is optimized using the standard inverse design methods. The resultant three best-adapted THBs, variable in topology and distinctive in sequences, are assessed by comparing them with a few benchmark proteins. The results of mainly molecular dynamics (MD) comparisons, undertaken in explicit water at different temperatures, show that the designed sequences are favorably placed against the chosen benchmarks as THB proteins potentially thermostable in the native folds. Folding simulation experiments with MD establish that the designed sequences are rapid in the folding of individual helices, but not in the evolution of tertiary structure; energetic cum topological frustrations remain but could be the artifacts of the starting conformations that were chosen in the THBs in the folding simulations. Overall, a practical high-throughput approach for de novo protein design has been developed that may have fruitful application for any type of tertiary structure.     

A framework for describing topological frustration in models of protein folding

In a natively folded protein of moderate or larger size, the protein backbone may weave through itself in complex ways, raising questions about what sequence of events might have to occur in order for the protein to reach its native configuration from the unfolded state. A mathematical framework is presented here for describing the notion of a topological folding barrier, which occurs when a protein chain must pass through a hole or opening, formed by other regions of the protein structure. Different folding pathways encounter different numbers of such barriers and therefore different degrees of frustration. A dynamic programming algorithm finds the optimal theoretical folding path and minimal degree of frustration for a protein based on its natively folded configuration. Calculations over a database of protein structures provide insights into questions such as whether the path of minimal frustration might tend to favor folding from one or from many sites of folding nucleation, or whether proteins favor folding around the N terminus, thereby providing support for the hypothesis that proteins fold co-translationally. The computational methods are applied to a multi-disulfide bonded protein, with computational findings that are consistent with the experimentally observed folding pathway. Attention is drawn to certain complex protein folds for which the computational method suggests there may be a preferred site of nucleation or where folding is likely to proceed through a relatively well-defined pathway or intermediate. The computational analyses lead to testable models for protein folding.     

The protein folding 'speed limit'

How fast can a protein possibly fold? This question has stimulated experimentalists to seek fast folding proteins and to engineer them to fold even faster. Proteins folding at or near the speed limit are prime candidates for all-atom molecular dynamics simulations. They may also have no free energy barrier, allowing the direct observation of intermediate structures on the pathways from the unfolded to the folded state. Both experimental and theoretical approaches predict a speed limit of approximately N/100micros for a generic N-residue single-domain protein, with alpha proteins folding faster than beta or alphabeta. The predicted limits suggest that most known ultrafast folding proteins can be engineered to fold more than ten times faster.     

Kinetics of unfolding and folding from amide hydrogen exchange in native ubiquitin

Amide hydrogen (NH) exchange is one of the few experimental techniques with the potential for determining the thermodynamics and kinetics of conformational motions at nearly every residue in native proteins. Quantitative interpretation of NH exchange in terms of molecular motions relies on a simple two-state kinetic model: at any given slowly exchanging NH, a closed or exchange-incompetent conformation is in equilibrium with an open or exchange-competent conformation. Previous studies have demonstrated the accuracy of this model in measuring conformational equilibria by comparing exchange data with the thermodynamics of protein unfolding. We report here a test of the accuracy of the model in determining the kinetics of conformational changes in native proteins. The kinetics of folding and unfolding for ubiquitin have been measured by conventional methods and compared with those derived from a comprehensive analysis of the pH dependence of exchange in native ubiquitin. Rate constants for folding and unfolding from these two very different types of experiments show good agreement. The simple model for NH exchange thus appears to be a robust framework for obtaining quantitative information about molecular motions in native proteins.     

Threading a peptide through a peptide: protein loops, rotaxanes, and knots

Proteins adopt complex folds in nature that typically avoid conformations that are knotted or “threaded” through closed loops. Is this the result of fundamental barriers to folding, or have proteins simply evolved to avoid threaded conformations? Organic synthesis has been used in supramolecular chemistry to install topological links in small molecules. By following these principles, we now show that it is possible to assemble a topologically linked protein complex by threading a linear protein through a cyclic protein to form a [2]pseudo‐rotaxane. Subsequent ring closure using native chemical ligation cyclizes the linear protein, forming a [2]heterocatenane. Although the kinetics of protein threading are slower than the folding kinetics of the native protein, threading appears to be a highly efficient process.     

Thermodynamics, kinetics, and salt dependence of folding of YopM, a large leucine-rich repeat protein

Small globular proteins have many contacts between residues that are distant in primary sequence. These contacts create a complex network between sequence-distant segments of secondary structure, which may be expected to promote the cooperative folding of globular proteins. Although repeat proteins, which are composed of tandem modular units, lack sequence-distant contacts, several of considerable length have been shown to undergo cooperative two-state folding. To explore the limits of cooperativity in repeat proteins, we have studied the unfolding of YopM, a leucine-rich repeat (LRR) protein of over 400 residues. Despite its large size and modular architecture (15 repeats), YopM equilibrium unfolding is highly cooperative, and shows a very strong dependence on the concentration of urea. In contrast, kinetic studies of YopM folding indicate a mechanism that includes one or more transient intermediates. The urea dependence of the folding and unfolding rates suggests a relatively small transition state ensemble. As with the urea dependence, we have found an extreme dependence of the free energy of unfolding on the concentration of salt. This salt dependence likely results from general screening of a large number of unfavorable columbic interactions in the folded state, rather than from specific cation binding.     

Expression, folding, and characterization of small proteins with increasing disulfide complexity by a pT7-7-derived phagemid

The expression, folding, and characterization of a series of small proteins with increasingly complex disulfide bond patterns were characterized. A phagemid was prepared from the pT7-7 plasmid to facilitate mutagenic studies with these proteins. cDNAs coding for bovine, rat, and human prolactin; human growth hormone; and bovine alpha-lactalbumin were amplified by PCR using primers that inserted restriction sites at the 5' and 3' ends and reduced the coding sequence to the mature methionyl protein with bacterially preferred codons in the 5' region. The expressed proteins were folded and oxidized by methods that allowed disulfide bond formation to occur either during or following folding. The effectiveness of the folding procedures was determined for each protein by electrophoresis, absorption spectroscopy, and functional studies. The redox conditions required for folding functional proteins varied as the number of disulfide bonds per unit molecular weight increased. Human growth hormone, 22 kDa; human prolactin, 23 kDa; and bovine prolactin, 23 kDa, contain two, three, and three disulfides, respectively, and are folded correctly by air oxidation performed during renaturation under alkaline conditions. Proper disulfide bond formation of rat prolactin, 23 kDa, containing three disulfide bonds required the addition of a reducing agent at the initiation of renaturation. Bovine alpha-lactalbumin, 14 kDa with four disulfide bonds, required complete renaturation prior to the removal of a reducing agent. SDS-gel electrophoresis under nonreducing conditions provided information regarding the proper folding of these proteins. The absorption of 250-nm light by disulfide bonds also provided information regarding the proper folding of rat prolactin and bovine alpha-lactalbumin.     

Membrane proteins shape up: understanding in vitro folding

Progress in structural biology has begun to reveal the precise architecture of integral membrane proteins. However, the manner in which these complex structures are achieved remains unclear. Recent developments are starting to shed light on the unfolding and folding of a small but growing number of membrane proteins. Mechanistic details derived from kinetic and thermodynamic experiments now enable comparison of the folding of different membrane proteins and their water-soluble cousins. This work also has important implications for other structural and functional studies of membrane proteins in vitro.     

Kinetics and mechanism of the folding of cytochrome c.

The reversible folding of cytochrome c in urea at pH 4.0 was investigated by repetitive pressure perturbation kinetics and by equilibrium spectroscopic methods. Two folding reactions were observed in the 1 ms to 10 s time range. The rates and amplitudes of these reactions depend on urea concentration in a complex manner, which is different for each process. The absorbance spectra of the kinetic amplitudes of the two reactions also differ from each other. A model with a three-state mechanism can quantitatively account for all of the kinetic and equilibrium data, and it enables us to determine the rate constants and volume changes of the two steps. If a rapid protonation step is added to the mechanism, the analysis can be extended to calculate the pH dependence of the rate and amplitude of the faster folding step. This pH dependence is in excellent agreement with previously published data [Tsong, T. Y. (1977) J. Biol. Chem. 252, 8778-8780]. Kinetic experiments in the 695-nm band show clearly that the axial ligand methionine-80 is involved in the slow folding process and the other axial ligand, histidine-18, is involved in the fast process. Additional experiments with a cyanogen bromide fragment of the protein, and fluorescence detection of the folding kinetics of the intact protein, support an interpretation of the model in terms of known structural elements of cytochrome c. This work provides new information about the mechanism of the folding of cytochrome c, resolves conflicts in earlier interpretations, and demonstrates the applicability of the repetitive pressure perturbation kinetics method to protein folding.     

Absence of a stable intermediate on the folding pathway of protein A.

The B-domain of protein A has one of the simplest protein topologies, a three-helix bundle. Its folding has been studied as a model for elementary steps in the folding of larger proteins. Earlier studies suggested that folding might occur by way of a helical hairpin intermediate. Equilibrium hydrogen exchange measurements indicate that the C-terminal helical hairpin could be a potential folding intermediate. Kinetic refolding experiments were performed using stopped-flow circular dichroism and NMR hydrogen-deuterium exchange pulse labeling. Folding of the entire molecule is essentially complete within the 6 ms dead time of the quench-flow apparatus, indicating that the intermediate, if formed, progresses rapidly to the final folded state. Site-directed mutagenesis of the isoleucine residue at position 16 was used to generate a variant protein containing tryptophan (the 116 W mutant). The formation of the putative folding intermediate was expected to be favored in this mutant at the expense of the native folded form, due to predicted unfavorable steric interactions of the bulky tryptophan side chain in the folded state. The 116 W mutant refolds completely within the dead time of a stopped-flow fluorescence experiment. No partly folded intermediate could be detected by either kinetic or equilibrium measurements. Studies of peptide fragments suggest that the protein A sequence has an intrinsic propensity to form a helix II/helix III hairpin. However, its stability appears to be marginal (of the order of 1/2 kT) and it could not be an obligatory intermediate on a defined folding pathway. These results explicitly demonstrate that the protein A B domain folds extremely rapidly by an apparent two-state mechanism without formation of stable partly folded intermediates. Similar mechanisms may also be involved in the rapid folding of subdomains of larger proteins to form the compact molten globule intermediates that often accumulate during the folding process.     

Desolvation Barrier Effects Are a Likely Contributor to the Remarkable Diversity in the Folding Rates of Small Proteins

The variation in folding rate among single-domain natural proteins is tremendous, but common models with explicit representations of the protein chain are either demonstrably insufficient or unclear as to their capability for rationalizing the experimental diversity in folding rates. In view of the critical role of water exclusion in cooperative folding, we apply native-centric, coarse-grained chain modeling with elementary desolvation barriers to investigate solvation effects on folding rates. For a set of 13 proteins, folding rates simulated with desolvation barriers cover ∼ 4.6 orders of magnitude, spanning a range essentially identical to that observed experimentally. In contrast, folding rates simulated without desolvation barriers cover only ∼ 2.2 orders of magnitude. Following a Hammond-like trend, the folding transition-state ensemble (TSE) of a protein model with desolvation barriers generally has a higher average number of native contacts and is structurally more specific, that is, less diffused, than the TSE of the corresponding model without desolvation barriers. Folding is generally significantly slower in models with desolvation barriers because of their higher overall macroscopic folding barriers as well as slower conformational diffusion speeds in the TSE that are ≈ 1/50 times those in models without desolvation barriers. Nonetheless, the average root-mean-square deviation between the TSE and the native conformation is often similar in the two modeling approaches, a finding suggestive of a more robust structural requirement for the folding rate-limiting step. The increased folding rate diversity in models with desolvation barriers originates from the tendency of these microscopic barriers to cause more heightening of the overall macroscopic folding free-energy barriers for proteins with more nonlocal native contacts than those with fewer such contacts. Thus, the enhancement of folding cooperativity by solvation effects is seen as positively correlated with a protein's native topological complexity.     

Folding and stability of sweet protein single-chain monellin. An insight to protein engineering

Engineered single-chain monellin (SCM) proteins were constructed by recombinant technology without disrupting the topology and sweet activity of native protein. Data from 8-anilinonaphthalene-1-sulfonic acid fluorescence, size-exclusion chromatography, and heteronuclear NMR strongly suggest the presence of a folding intermediate at 1.5 m GdnHCl for SCM protein. The structural feature of the folding intermediate from NMR data reveals that the secondary structures became mostly unstable, and protein experiences a dynamic equilibrium between native and unfolded state. All backbone amide protons exchange within 10 min, which imply that no stable hydrogen bonds exist in the secondary structural regions in the folding intermediate. From equilibrium unfolding and mutagenesis studies, the unfolding transition midpoints of mutant proteins gradually shifted toward lower denaturant concentration, indicating stability reductions of mutant proteins. Our results suggest that stability and folding pathways of SCM proteins could be regulated by a combined study of spectroscopy and mutagenesis, and these studies will provide useful information for understanding the folding kinetics of novel engineered proteins.     

Designing a High Throughput Refolding Array Using a Combination of the GroEL Chaperonin and Osmolytes

Although GroE chaperonins and osmolytes had been used separately as protein folding aids, combining these two methods provides a considerable advantage for folding proteins that cannot fold with either osmolytes or chaperonins alone. This technique rapidly identifies superior folding solution conditions for a broad array of proteins that are difficult or impossible to fold by other methods. While testing the broad applicability of this technique, we have discovered that osmolytes greatly simplify the chaperonin reaction by eliminating the requirement for the co-chaperonin GroES which is normally involved in encapsulating folding proteins within the GroEL–GroES cavity. Therefore, combinations of soluble or immobilized GroEL, osmolytes and ATP or even ADP are sufficient to refold the test proteins. The first step in the chaperonin/osmolyte process is to form a stable long-lived chaperonin–substrate protein complex in the absence of nucleotide. In the second step, different osmolyte solutions are added along with nucleotides, thus forming a ‘folding array’ to identify superior folding conditions. The stable chaperonin–substrate protein complex can be concentrated or immobilized prior to osmolyte addition. This procedure prevents-off pathway aggregation during folding/refolding reactions and more importantly allows one to refold proteins at concentrations (~mg/ml) that are substantially higher than the critical aggregation concentration for given protein. This technique can be used for successful refolding of proteins from purified inclusion bodies. Recently, other investigators have used our chaperonin/osmolyte method to demonstrate that a mutant protein that misfolds in human disease can be rescued by GroEL/osmolyte system. Soluble or immobilized GroEL can be easily removed from the released folded protein using simple separation techniques. The method allows for isolation of folded monomeric or oligomeric proteins in quantities sufficient for X-ray crystallography or NMR structural determinations.     

Structure- and sequence-based function prediction for non-homologous proteins

The structural genomics projects have been accumulating an increasing number of protein structures, many of which remain functionally unknown. In parallel effort to experimental methods, computational methods are expected to make a significant contribution for functional elucidation of such proteins. However, conventional computational methods that transfer functions from homologous proteins do not help much for these uncharacterized protein structures because they do not have apparent structural or sequence similarity with the known proteins. Here, we briefly review two avenues of computational function prediction methods, i.e. structure-based methods and sequence-based methods. The focus is on our recent developments of local structure-based and sequence-based methods, which can effectively extract function information from distantly related proteins. Two structure-based methods, Pocket-Surfer and Patch-Surfer, identify similar known ligand binding sites for pocket regions in a query protein without using global protein fold similarity information. Two sequence-based methods, protein function prediction and extended similarity group, make use of weakly similar sequences that are conventionally discarded in homology based function annotation. Combined together with experimental methods we hope that computational methods will make leading contribution in functional elucidation of the protein structures.