Alessandro Borgia

Single-Molecule Studies of Protein Folding

Although protein-folding studies began several decades ago, it is only recently that the tools to analyze protein folding at the single-molecule level have been developed. Advances in single-molecule fluorescence and force spectroscopy techniques allow investigation of the folding and dynamics of single protein molecules, both at equilibrium and as they fold and unfold. The experiments are far from simple, however, both in execution and in interpretation of the results. In this review, we discuss some of the highlights of the work so far and concentrate on cases where comparisons with the classical experiments can be made. We conclude that, although there have been relatively few startling insights from single-molecule studies, the rapid progress that has been made suggests that these experiments have significant potential to advance our understanding of protein folding. In particular, new techniques offer the possibility to explore regions of the energy landscape that are inaccessible to classical ensemble measurements and, perhaps, to observe rare events undetectable by other means.

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.

Experimental evidence for a frustrated energy landscape in a three-helix-bundle protein family

Energy landscape theory is a powerful tool for understanding the structure and dynamics of complex molecular systems, in particular biological macromolecules. The primary sequence of a protein defines its free-energy landscape and thus determines the folding pathway and the rate constants of folding and unfolding, as well as the protein’s native structure. Theory has shown that roughness in the energy landscape will lead to slower folding, but derivation of detailed experimental descriptions of this landscape is challenging. Simple folding models show that folding is significantly influenced by chain entropy; proteins in which the contacts are local fold quickly, owing to the low entropy cost of forming stabilizing, native contacts during folding. For some protein families, stability is also a determinant of folding rate constants. Where these simple metrics fail to predict folding behaviour, it is probable that there are features in the energy landscape that are unusual. Such general observations cannot explain the folding behaviour of the R15, R16 and R17 domains of α-spectrin. R15 folds ∼3,000 times faster than its homologues, although they have similar structures, stabilities and, as far as can be determined, transition-state stabilities. Here we show that landscape roughness (internal friction) is responsible for the slower folding and unfolding of R16 and R17. We use chimaeric domains to demonstrate that this internal friction is a property of the core, and suggest that frustration in the landscape of the slow-folding spectrin domains may be due to misdocking of the long helices during folding. Theoretical studies have suggested that rugged landscapes will result in slower folding; here we show experimentally that such a phenomenon directly influences the folding kinetics of a ‘normal’ protein, that is, one with a significant energy barrier that folds on a relatively slow, millisecond–second, timescale.

Single-molecule fluorescence reveals sequence-specific misfolding in multidomain proteins.

A large range of debilitating medical conditions is linked to protein misfolding, which may compete with productive folding particularly in proteins containing multiple domains. Seventy-five per cent of the eukaryotic proteome consists of multidomain proteins, yet it is not understood how interdomain misfolding is avoided. It has been proposed that maintaining low sequence identity between covalently linked domains is a mechanism to avoid misfolding. Here we use single-molecule Förster resonance energy transfer to detect and quantify rare misfolding events in tandem immunoglobulin domains from the I band of titin under native conditions. About 5.5 per cent of molecules with identical domains misfold during refolding in vitro and form an unexpectedly stable state with an unfolding half-time of several days. Tandem arrays of immunoglobulin-like domains in humans show significantly lower sequence identity between neighbouring domains than between non-adjacent domains. In particular, the sequence identity of neighbouring domains has been found to be preferentially below 40 per cent. We observe no misfolding for a tandem of naturally neighbouring domains with low sequence identity (24 per cent), whereas misfolding occurs between domains that are 42 per cent identical. Coarse-grained molecular simulations predict the formation of domain-swapped structures that are in excellent agreement with the observed transfer efficiency of the misfolded species. We infer that the interactions underlying misfolding are very specific and result in a sequence-specific domain-swapping mechanism. Diversifying the sequence between neighbouring domains seems to be a successful evolutionary strategy to avoid misfolding in multidomain proteins.

Localizing internal friction along the reaction coordinate of protein folding by combining ensemble and single-molecule fluorescence spectroscopy

Theory, simulations and experimental results have suggested an important role of internal friction in the kinetics of protein folding. Recent experiments on spectrin domains provided the first evidence for a pronounced contribution of internal friction in proteins that fold on the millisecond timescale. However, it has remained unclear how this contribution is distributed along the reaction and what influence it has on the folding dynamics. Here we use a combination of single-molecule Förster resonance energy transfer, nanosecond fluorescence correlation spectroscopy, microfluidic mixing and denaturant- and viscosity-dependent protein-folding kinetics to probe internal friction in the unfolded state and at the early and late transition states of slow- and fast-folding spectrin domains. We find that the internal friction affecting the folding rates of spectrin domains is highly localized to the early transition state, suggesting an important role of rather specific interactions in the rate-limiting conformational changes.

Quantifying heterogeneity and conformational dynamics from single molecule FRET of diffusing molecules: recurrence analysis of single particles (RASP)

Single molecule Förster resonance energy transfer (FRET) experiments are a versatile method for investigating the conformational distributions and dynamics of biological macromolecules. In a common type of experiment, the fluorescence bursts from individual molecules freely diffusing in solution are detected as they pass through the observation volume of a confocal microscope. Correlation analysis of the fluorescence bursts shows that under typical experimental conditions, for time scales up to several tens of milliseconds, the probability for a molecule to return to the confocal volume is greater than the probability of a new molecule being detected. Here we present RASP (recurrence analysis of single particles), a method that is based on this recurrence behavior and allows us to significantly extend the information that can be extracted from single molecule FRET experiments. The number and peak shapes of subpopulations within the sample can be identified essentially in a model-free way by constructing recurrence FRET efficiency histograms. These are obtained by first selecting photon bursts from a small transfer efficiency range (initial bursts), and then building the FRET efficiency histogram only from bursts detected within a short time (the recurrence interval) after the initial bursts. Systematic variation of the recurrence interval allows the kinetics of interconversion between subpopulations to be determined on time scales from ~50 μs up to ~100 ms from equilibrium measurements. We demonstrate the applicability of the method on measurements of several peptides and proteins with different degrees of conformational heterogeneity and folding dynamics. The concepts presented here can be extended to other observables available from single molecule fluorescence experiments.

Consistent View of Polypeptide Chain Expansion in Chemical Denaturants from Multiple Experimental Methods.

There has been a long-standing controversy regarding the effect of chemical denaturants on the dimensions of unfolded and intrinsically disordered proteins: A wide range of experimental techniques suggest that polypeptide chains expand with increasing denaturant concentration, but several studies using small-angle X-ray scattering (SAXS) have reported no such increase of the radius of gyration (Rg). This inconsistency challenges our current understanding of the mechanism of chemical denaturants, which are widely employed to investigate protein folding and stability. Here, we use a combination of single-molecule Förster resonance energy transfer (FRET), SAXS, dynamic light scattering (DLS), and two-focus fluorescence correlation spectroscopy (2f-FCS) to characterize the denaturant dependence of the unfolded state of the spectrin domain R17 and the intrinsically disordered protein ACTR in two different denaturants. Standard analysis of the primary data clearly indicates an expansion of the unfolded state with increasing denaturant concentration irrespective of the protein, denaturant, or experimental method used. This is the first case in which SAXS and FRET have yielded even qualitatively consistent results regarding expansion in denaturant when applied to the same proteins. To more directly illustrate this self-consistency, we used both SAXS and FRET data in a Bayesian procedure to refine structural ensembles representative of the observed unfolded state. This analysis demonstrates that both of these experimental probes are compatible with a common ensemble of protein configurations for each denaturant concentration. Furthermore, the resulting ensembles reproduce the trend of increasing hydrodynamic radius with denaturant concentration obtained by 2f-FCS and DLS. We were thus able to reconcile the results from all four experimental techniques quantitatively, to obtain a comprehensive structural picture of denaturant-induced unfolded state expansion, and to identify the most likely sources of earlier discrepancies.

Probing the Action of Chemical Denaturant on an Intrinsically Disordered Protein by Simulation and Experiment.

Chemical denaturants are the most commonly used agents for unfolding proteins and are thought to act by better solvating the unfolded state. Improved solvation is expected to lead to an expansion of unfolded chains with increasing denaturant concentration, providing a sensitive probe of the denaturant action. However, experiments have so far yielded qualitatively different results concerning the effects of chemical denaturation. Studies using Förster resonance energy transfer (FRET) and other methods found an increase in radius of gyration with denaturant concentration, but most small-angle X-ray scattering (SAXS) studies found no change. This discrepancy therefore challenges our understanding of denaturation mechanism and more generally the accuracy of these experiments as applied to unfolded or disordered proteins. Here, we use all-atom molecular simulations to investigate the effect of urea and guanidinium chloride on the structure of the intrinsically disordered protein ACTR, which can be studied by experiment over a wide range of denaturant concentration. Using unbiased molecular simulations with a carefully calibrated denaturant model, we find that the protein chain indeed swells with increasing denaturant concentration. This is due to the favorable association of urea or guanidinium chloride with the backbone of all residues and with the side-chains of almost all residues, with denaturant-water transfer free energies inferred from this association in reasonable accord with experimental estimates. Interactions of the denaturants with the backbone are dominated by hydrogen bonding, while interactions with side-chains include other contributions. By computing FRET efficiencies and SAXS intensities at each denaturant concentration, we show that the simulation trajectories are in accord with both experiments on this protein, demonstrating that there is no fundamental inconsistency between the two types of experiment. Agreement with experiment also supports the picture of chemical denaturation described in our simulations, driven by weak association of denaturant with the protein. Our simulations support some assumptions needed for each experiment to accurately reflect changes in protein size, namely, that the commonly used FRET chromophores do not qualitatively alter the results and that possible effects such as preferential solvent partitioning into the interior of the chain do not interfere with the determination of radius of gyration from the SAXS experiments.

Unveiling a Hidden Folding Intermediate in c-Type Cytochromes by Protein Engineering

Several investigators have highlighted a correlation between the basic features of the folding process of a protein and its topology, which dictates the folding pathway. Within this conceptual framework we proposed that different members of the cytochrome c (cyt c) family share the same folding mechanism, involving a consensus partially structured state. Pseudomonas aeruginosa cyt c551 (Pa cyt c551) folds via an apparent two-state mechanism through a high energy intermediate. Here we present kinetic evidence demonstrating that it is possible to switch its folding mechanism from two to three state, stabilizing the high energy intermediate by rational mutagenesis. Characterization of the folding kinetics of one single-site mutant of the Pa cyt c551 (Phe7 to Ala) indeed reveals an additional refolding phase and a fast unfolding process which are explained by the accumulation of a partially folded species. Further kinetic analysis highlights the presence of two parallel processes both leading to the native state, suggesting that the above mentioned species is a non obligatory on-pathway intermediate. Determination of the crystallographic structure of F7A shows the presence of an extended internal cavity, which hosts three “bound” water molecules and a H-bond in the N-terminal helix, which is shorter than in the wild type protein. These two features allow us to propose a detailed structural interpretation for the stabilization of the native and especially the intermediate states induced by a single crucial mutation. These results show how protein engineering, x-ray crystallography and state-of-the-art kinetics concur to unveil a folding intermediate and the structural determinants of its stability.

Transient misfolding dominates multidomain protein folding

Neighbouring domains of multidomain proteins with homologous tandem repeats have divergent sequences, probably as a result of evolutionary pressure to avoid misfolding and aggregation, particularly at the high cellular protein concentrations. Here we combine microfluidic-mixing single-molecule kinetics, ensemble experiments and molecular simulations to investigate how misfolding between the immunoglobulin-like domains of titin is prevented. Surprisingly, we find that during refolding of tandem repeats, independent of sequence identity, more than half of all molecules transiently form a wide range of misfolded conformations. Simulations suggest that a large fraction of these misfolds resemble an intramolecular amyloid-like state reported in computational studies. However, for naturally occurring neighbours with low sequence identity, these transient misfolds disappear much more rapidly than for identical neighbours. We thus propose that evolutionary sequence divergence between domains is required to suppress the population of long-lived, potentially harmful misfolded states, whereas large populations of transient misfolded states appear to be tolerated.