Where the complex things are: single molecule and ensemble spectroscopic investigations of protein folding dynamics

Pathway regulation mechanism revealed by cotranslational folding of villin headpiece subdomain HP35

Cotranslational folding is one of the most important features of protein folding in vivo. Although many studies have shown that the folding pathways of cotranslational folding are different from free folding in vitro, the detailed mechanism of folding dynamics is lacking. Here we combine all-atom molecular simulations with an ideal ribosome tunnel model to investigate the cotranslational folding of villin headpiece subdomain HP35. By comparing the folding dynamics between cotranslational folding and free folding, we found that cotranslational folding tends to fold along the pathway that is easier to fold into native state in the latter. In addition, the roles of the ribosome tunnel and sequential folding are analyzed separately. Our results show that the ribosome can prevent the untimely folding of the C segment of HP35 to reduce the non-native interactions, while the translation speed can regulate the amounts of native and non-native interactions and the balance between them. Overall, these results give insights into the general mechanisms of cotranslational protein folding.

Dramatic Shape Changes Occur as Cytochrome c Folds

Extensive experimental studies on the folding of cytochrome c (Cyt c) make this small protein an ideal target for atomic detailed simulations for the purposes of quantitatively characterizing the structural transitions and the associated time scales for folding to the native state from an ensemble of unfolded states. We use previously generated atomically detailed folding trajectories by the stochastic difference equation in length to calculate the time-dependent changes in the small-angle X-ray scattering (SAXS) profiles. Excellent agreement is obtained between experiments and simulations for the time-dependent SAXS spectra, allowing us to identify the structures of the folding intermediates, which shows that Cyt c reaches the native state by a sequential folding mechanism. Using the ensembles of structures along the folding pathways, we show that compaction and the sphericity of Cyt c change dramatically from the prolate ellipsoid shape in the unfolded state to the spherical native state. Our data, which are in unprecedented quantitative agreement with all aspects of time-resolved SAXS experiments, show that hydrophobic collapse and amide group protection coincide on the 100 microseconds time scale, which is in accordance with ultrafast hydrogen/deuterium exchange studies. Based on these results, we propose that compaction of polypeptide chains, accompanied by dramatic shape changes, is a universal characteristic of globular proteins, regardless of the underlying folding mechanism.

Dynamics of proteins in solution

The dynamics of proteins in solution includes a variety of processes, such as backbone and side-chain fluctuations, interdomain motions, as well as global rotational and translational (i.e. center of mass) diffusion. Since protein dynamics is related to protein function and essential transport processes, a detailed mechanistic understanding and monitoring of protein dynamics in solution is highly desirable. The hierarchical character of protein dynamics requires experimental tools addressing a broad range of time- and length scales. We discuss how different techniques contribute to a comprehensive picture of protein dynamics, and focus in particular on results from neutron spectroscopy. We outline the underlying principles and review available instrumentation as well as related analysis frameworks.

Hypothesis: structural heterogeneity of the unfolded proteins originating from the coupling of the local clusters and the long-range distance distribution

We propose a hypothesis that explains two apparently contradicting observations for the heterogeneity of the unfolded proteins. First, the line confocal method of the single-molecule Förster resonance energy transfer (sm-FRET) spectroscopy revealed that the unfolded proteins possess broad peaks in the FRET efficiency plot, implying the significant heterogeneity that lasts longer than milliseconds. Second, the fluorescence correlation method demonstrated that the unfolded proteins fluctuate in the time scale shorter than 100 ns. To formulate the hypothesis, we first summarize the recent consensus for the structure and dynamics of the unfolded proteins. We next discuss the conventional method of the sm-FRET spectroscopy and its limitations for the analysis of the unfolded proteins, followed by the advantages of the line confocal method that revealed the heterogeneity. Finally, we propose that the structural heterogeneity formed by the local clustering of hydrophobic residues modulates the distribution of the long-range distance between the labeled chromophores, resulting in the broadening of the peak in the FRET efficiency plot. A clustering of hydrophobic residues around the chromophore might further contribute to the broadening. The proposed clusters are important for the understanding of protein folding mechanism.

Microsecond resolved single-molecule FRET time series measurements based on the line confocal optical system combined with hybrid photodetectors

Single-molecule (sm) fluorescence time series measurements based on the line confocal optical system are a powerful strategy for the investigation of the structure, dynamics, and heterogeneity of biological macromolecules. This method enables the detection of more than several thousands of fluorescence photons per millisecond from single fluorophores, implying that the potential time resolution for measurements of the fluorescence resonance energy transfer (FRET) efficiency is 10 μs. However, the necessity of using imaging photodetectors in the method limits the time resolution in the FRET efficiency measurements to approximately 100 μs. In this investigation, a new photodetector called a hybrid photodetector (HPD) was incorporated into the line confocal system to improve the time resolution without sacrificing the length of the time series detection. Among several settings examined, the system based on a slit width of 10 μm and a high-speed counting device made the best of the features of the line confocal optical system and the HPD. This method achieved a time resolution of 10 μs and an observation time of approximately 5 ms in the sm-FRET time series measurements. The developed device was used for the native state of the B domain of protein A.

Unified understanding of folding and binding mechanisms of globular and intrinsically disordered proteins

Extensive experimental and theoretical studies have advanced our understanding of the mechanisms of folding and binding of globular proteins, and coupled folding and binding of intrinsically disordered proteins (IDPs). The forces responsible for conformational changes and binding are common in both proteins; however, these mechanisms have been separately discussed. Here, we attempt to integrate the mechanisms of coupled folding and binding of IDPs, folding of small and multi-subdomain proteins, folding of multimeric proteins, and ligand binding of globular proteins in terms of conformational selection and induced-fit mechanisms as well as the nucleation–condensation mechanism that is intermediate between them. Accumulating evidence has shown that both the rate of conformational change and apparent rate of binding between interacting elements can determine reaction mechanisms. Coupled folding and binding of IDPs occurs mainly by induced-fit because of the slow folding in the free form, while ligand binding of globular proteins occurs mainly by conformational selection because of rapid conformational change. Protein folding can be regarded as the binding of intramolecular segments accompanied by secondary structure formation. Multi-subdomain proteins fold mainly by the induced-fit (hydrophobic collapse) mechanism, as the connection of interacting segments enhances the binding (compaction) rate. Fewer hydrophobic residues in small proteins reduce the intramolecular binding rate, resulting in the nucleation–condensation mechanism. Thus, the folding and binding of globular proteins and IDPs obey the same general principle, suggesting that the coarse-grained, statistical mechanical model of protein folding is promising for a unified theoretical description of all mechanisms.

Protein folding transition path times from single molecule FRET

The transition path is the tiny segment of a single molecule trajectory when the free energy barrier between states is crossed and for protein folding contains all of the information about the self-assembly mechanism. As a first step toward obtaining structural information during the transition path from experiments, single molecule FRET spectroscopy has been used to determine average transition path times from a photon-by-photon analysis of fluorescence trajectories. These results, obtained for several different proteins, have already provided new and demanding tests that support both the accuracy of all-atom molecular dynamics simulations and the basic postulates of energy landscape theory of protein folding.

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.

Can Local Probes Go Global? A Joint Experiment-Simulation Analysis of λ(6-85) Folding

The process of protein folding is known to involve global motions in a cooperative affair; the structure of most of the protein sequences is gained or lost over a narrow range of temperature, denaturant, or pressure perturbations. At the same time, recent simulations and experiments reveal a complex structural landscape with a rich set of local motions and conformational changes. We couple experimental kinetic and thermodynamic measurements with specifically tailored analysis of simulation data to isolate local versus global folding probes. We find that local probes exhibit lower melting temperatures, smaller surface area changes, and faster kinetics compared to global ones. We also see that certain local probes of folding match the global behavior more closely than others. Our work highlights the importance of using multiple probes to fully characterize protein folding dynamics by theory and experiment.

Fast Folding Dynamics of an Intermediate State in RNase H Measured by Single-Molecule FRET

We have studied the folding kinetics of the core intermediate (I) state of RNase H by using a combination of single-molecule FRET (smFRET) and hidden Markov model analysis. To measure fast dynamics in thermal equilibrium as a function of the concentration of the denaturant GdmCl, a special FRET labeled variant, RNase H 60-113, which is sensitive to folding of the protein core, was immobilized on PEGylated surfaces. Conformational transitions between the unfolded (U) state and the I state could be described by a two-state model within our experimental time resolution, with millisecond mean residence times. The I state population was always a minority species in the entire accessible range of denaturant concentrations. By introducing the measured free energy differences between the U and I states as constraints in global fits of the GdmCl dependence of FRET histograms of a differently labeled RNase H variant (RNase H 3-135), we were able to reveal the free energy differences and, thus, population ratios of all three macroscopic state ensembles, U, I and F (folded state) as a function of denaturant concentration.