Energy barriers, pathways, and dynamics during folding of large, multidomain RNAs

Large, multidomain RNA molecules are generally thought to fold following multiple pathways down rugged landscapes populated with intermediates and traps. A challenge to understanding RNA folding reactions is the complex relationships that exist between the structure of the RNA and its folding landscape. The identification of intermediate species that populate folding landscapes and characterization of elements of their structures are the key components to solving the RNA folding problem. This review explores recent studies that characterize the dominant pathways by which RNA folds, structural and dynamic features of intermediates that populate the folding landscape, and the energy barriers that separate the distinct steps of the folding process.     

Folding simulations of small proteins

Understanding how a protein folds is a long-standing challenge in modern science. We have used an optimized atomistic model (united-residue force field) to simulate folding of small proteins of various structures: HP-36 (alpha protein), protein A (beta), 1fsd (alpha+beta), and betanova (beta). Extensive Monte Carlo folding simulations (ten independent runs with 10(9) Monte Carlo steps at a temperature) starting from non-native conformations are carried out for each protein. In all cases, proteins fold into their native-like conformations at appropriate temperatures, and glassy transitions occur at low temperatures. To investigate early folding trajectories, 200 independent runs with 10(6) Monte Carlo steps are also performed at a fixed temperature for a protein. There are a variety of possible pathways during non-equilibrium early processes (fast process, approximately 10(4) Monte Carlo steps). Finally, these pathways converge to the point unique for each protein. The convergence point of the early folding pathways can be determined only by direct folding simulations. The free energy surface, an equilibrium thermodynamic property, dictates the rest of the folding (slow process, approximately 10(8) Monte Carlo steps).     

Direct Observation of Parallel Folding Pathways Revealed Using a Symmetric Repeat Protein System

Although progress has been made to determine the native fold of a polypeptide from its primary structure, the diversity of pathways that connect the unfolded and folded states has not been adequately explored. Theoretical and computational studies predict that proteins fold through parallel pathways on funneled energy landscapes, although experimental detection of pathway diversity has been challenging. Here, we exploit the high translational symmetry and the direct length variation afforded by linear repeat proteins to directly detect folding through parallel pathways. By comparing folding rates of consensus ankyrin repeat proteins (CARPs), we find a clear increase in folding rates with increasing size and repeat number, although the size of the transition states (estimated from denaturant sensitivity) remains unchanged. The increase in folding rate with chain length, as opposed to a decrease expected from typical models for globular proteins, is a clear demonstration of parallel pathways. This conclusion is not dependent on extensive curve-fitting or structural perturbation of protein structure. By globally fitting a simple parallel-Ising pathway model, we have directly measured nucleation and propagation rates in protein folding, and have quantified the fluxes along each path, providing a detailed energy landscape for folding. This finding of parallel pathways differs from results from kinetic studies of repeat-proteins composed of sequence-variable repeats, where modest repeat-to-repeat energy variation coalesces folding into a single, dominant channel. Thus, for globular proteins, which have much higher variation in local structure and topology, parallel pathways are expected to be the exception rather than the rule.     

Multiscale simulations of protein folding: application to formation of secondary structures

A multiscale simulation method of protein folding is proposed, using atomic representation of protein and solvent, combing genetic algorithms to determine the key protein structures from a global view, with molecular dynamic simulations to reveal the local folding pathways, thus providing an integrated landscape of protein folding. The method is found to be superior to previously investigated global search algorithms or dynamic simulations alone. For secondary structure formation of a selected peptide, RN24, the structures and dynamics produced by this method agree well with corresponding experimental results. Three most populated conformations are observed, including hairpin, β-sheet and α-helix. The energetic barriers separating these three structures are comparable to the kinetic energy of the atoms of the peptide, implying that the transition between these states can be easily triggered by kinetic perturbations, mainly through electrostatic interactions between charged atoms. Transitions between α-helix and β-sheet should jump over at least two energy barriers and may stay in the energetic trap of hairpin. It is proposed that the structure of proteins should be jointly governed by thermodynamic and dynamic factors; free energy is not the exclusive dominant for stability of proteins.     

Folding studies of immunoglobulin-like beta-sandwich proteins suggest that they share a common folding pathway.

BACKGROUND: Are folding pathways conserved in protein families? To test this explicitly and ask to what extent structure specifies folding pathways requires comparison of proteins with a common fold. Our strategy is to choose members of a highly diverse protein family with no conservation of function and little or no sequence identity, but with structures that are essentially the same. The immunoglobulin-like fold is one of the most common structural families, and is subdivided into superfamilies with no detectable evolutionary or functional relationship. RESULTS: We compared the folding of a number of immunoglobulin-like proteins that have a common structural core and found a strong correlation between folding rate and stability. The results suggest that the folding pathways of these immunoglobulin-like proteins share common features. CONCLUSIONS: This study is the first to compare the folding of structurally related proteins that are members of different superfamilies. The most likely explanation for the results is that interactions that are important in defining the structure of immunoglobulin-like proteins are also used to guide folding.     

The family feud: do proteins with similar structures fold via the same pathway?

Theoretical and experimental studies of protein folding have suggested that the topology of the native state may be the most important factor determining the folding pathway of a protein, independent of its specific amino acid sequence. To test this concept, many experimental studies have been carried out with the aim of comparing the folding pathways of proteins that possess similar tertiary structures, but divergent sequences. Many of these studies focus on quantitative comparisons of folding transition state structures, as determined by Phi(f) value analysis of folding kinetic data. In some of these studies, folding transition state structures are found to be highly conserved, whereas in others they are not. We conclude that folds displaying more conserved transition state structures may have the most restricted number of possible folding pathways and that folds displaying low transition state structural conservation possess many potential pathways for reaching the native state.     

Autonomous subdomains in protein folding.

Proteolytic dissection of native trp repressor and horse heart cytochrome c has been used to infer some of the steps in the folding pathways of the intact proteins. For both proteins, small fragments are capable of undergoing spontaneous noncovalent association to form subdomains with native-like secondary and/or tertiary structural features, suggesting that dissection/reassembly may be a general method to gain insight into the structures of folding intermediates. The importance of this approach is its simplicity and potential applicability to studying the folding pathways of a wide range of proteins. The proteases report on the structure and dynamics of the native state, circumventing the need for prior knowledge of the structures of folding intermediates. The observation that small fragments of proteins can associated noncovalently suggests that protein folding can be viewed as an intramolecular "recognition" process. The results imply that substantial information about protein structure and folding is encoded at the level of subdomains, and that chain connectivity has only a minor role in determining the fold.     

Parallel folding pathways in the SH3 domain protein

The transition-state ensemble (TSE) is the set of protein conformations with an equal probability to fold or unfold. Its characterization is crucial for an understanding of the folding process. We determined the TSE of the src-SH3 domain protein by using extensive molecular dynamics simulations of the Go model and computing the folding probability of a generated set of TSE candidate conformations. We found that the TSE possesses a well-defined hydrophobic core with variable enveloping structures resulting from the superposition of three parallel folding pathways. The most preferred pathway agrees with the experimentally determined TSE, while the two least preferred pathways differ significantly. The knowledge of the different pathways allows us to design the interactions between amino acids that guide the protein to fold through the least preferred pathway. This particular design is akin to a circular permutation of the protein. The finding motivates the hypothesis that the different experimentally observed TSEs in homologous proteins and circular permutants may represent potentially available pathways to the wild-type protein.     

Parallel protein-unfolding pathways revealed and mapped

Theoretical studies of protein folding suggest that multiple folding pathways should exist, but there is little experimental evidence to support this. Here we demonstrate changes in the flux between different transition states on parallel folding pathways, resulting in unprecedented upward curvature in the denaturant-dependent unfolding kinetics of a beta-sandwich protein. As denaturant concentration increases, the highly compact transition state of one pathway becomes destabilized and the dominant flux of protein molecules shifts toward another pathway with a less structured transition state. Furthermore, point mutations alter the relative accessibility of the pathways, allowing the structure of two transition states on separate, direct folding pathways to be mapped by systematic Phi-value analysis. It has been suggested that pathways with diffuse rather than localized transition states are evolutionarily selected to prevent misfolding, and indeed we find that the transition state favored at high concentrations of denaturant is more polarized than the physiologically relevant one.     

Importance of sequence specificity for predicting protein folding pathways: perturbed Gaussian chain model

Recent experimental and theoretical studies suggest that rates and pathways of protein folding are largely decided by topology of the native structures, at least for small proteins. However, some exceptions are known; for example, protein L and protein G have the same topology, but exhibit different characteristics of the TSE. Thus, folding pathways of some proteins are critically affected by detailed information on amino acid sequences. To investigate the sequence specificity, we calculate folding pathways of 20 small proteins using the perturbed Gaussian chain model developed by Portman et al. (Phys Rev Lett 1998;81:5237-5240; J Chem Phys 2001;114:5069-5081). Characteristics of the TSE predicted by the model are in good agreement with experimental phi-value data for many proteins at coarse-grained level. Especially, estimation of folding TSE for protein G and protein L based on both topology and additional sequence information are consistent with experimental phi-value data. With only topology information, however, the model predicts the TSE of protein G incorrectly. Moreover, the model that uses topology and sequence information describes free energy profiles of two-state and three-state folders consistently with experiment, whereas the topology only model predicts free energy profiles of some proteins incorrectly. This indicates that sequence specificity also has critical roles in determining the folding pathways for some proteins.     

Searching for multiple folding pathways of a nearly symmetrical protein: temperature dependent phi-value analysis of the B domain of protein A

The B domain of protein A (BdpA) is a popular paradigm for simulating protein folding pathways. The discrepancies between so many simulations and subsequent experimental testing may be attributable to the protein being highly symmetrical: changing experimental conditions could perturb the subtle interplay between the effects of symmetry in the native structure and the effects of asymmetry from specific interactions in a given sequence. If the protein folds via multiple pathways, perturbations, such as temperature, denaturant concentration, and mutation, should change the flux of micro pathways, leading to changes in the bulk properties of the transition state. We tested this hypothesis by conducting a Phi-analysis of BdpA as a function of temperature from 25.0 degrees C to 60.0 degrees C. The Phi-values had no significant dependence on temperature and the values at 55.0 degrees C (denaturing conditions) are very similar to those at 25.0 degrees C (folding conditions), indicating the structure of the transition state does not significantly change although the experimental conditions are considerably altered. The results suggest that BdpA folds via a single dominant folding pathway.     

Compact dimension of denatured states of staphylococcal nuclease

Fluorescence and circular dichroism stopped-flow have been widely used to determine the kinetics of protein folding including folding rates and possible folding pathways. Yet, these measurements are not able to provide spatial information of protein folding/unfolding. Especially, conformations of denatured states cannot be elaborated in detail. In this study, we apply the method of fluorescence energy transfer with a stopped-flow technique to study global structural changes of the staphylococcal nuclease (SNase) mutant K45C, where lysine 45 is replaced by cysteine, during folding and unfolding. By labeling the thiol group of cysteine with TNB (5,5'-dithiobis-2-nitrobenzoic acid) as an energy acceptor and the tryptophan at position 140 as a donor, distance changes between the acceptor and the donor during folding and unfolding are measured from the efficiency of energy transfer. Results indicate that the denatured states of SNase are highly compact regardless of how the denatured states (pH-induced or GdmCl-induced) are induced. The range of distance changes between two probes is between 25.6 and 25.4 A while it is 20.4 A for the native state. Furthermore, the folding process consists of three kinetic phases while the unfolding process is a single phase. These observations agree with our previous sequential model: N(0) left arrow over right arrow D(1) left arrow over right arrow D(2) left arrow over right arrow D(3) (Chen et al., J Mol Biol 1991;220:771-778). The efficiency of protein folding may be attributed to initiating the folding process from these compact denatured structures.     

Increasing stability reduces conformational heterogeneity in a protein folding intermediate ensemble

A multi-site, time-resolved fluorescence resonance energy transfer methodology has been used to study structural heterogeneity in a late folding intermediate ensemble, IL, of the small protein barstar. Four different intra-molecular distances have been measured within the structural components of IL. The IL ensemble is shown to consist of different sub-populations of molecules, in each of which one or more of the four distances are native-like and the remaining distances are unfolded-like. In very stable conditions that favor formation of IL, all four distances are native-like in most molecules. In less stable conditions, one or more distances are unfolded-like. As stability is decreased, the proportion of molecules with unfolded-like distances increases. Thus, the results show that protein folding intermediates are ensembles of different structural forms, and they demonstrate that conformational entropy increases as structures become less stable. These observations provide direct experimental evidence in support of a basic tenet of energy landscape theory for protein folding, that available conformational space, as represented by structural heterogeneity in IL, becomes restricted as the stability is increased. The results also vindicate an important prediction of energy landscape theory, that different folding pathways may become dominant under different folding conditions. In more stable folding conditions, uniformly native-like compactness is achieved during folding to IL, whereas in less stable conditions, uniformly native-like compactness is achieved only later during the folding of IL to N.     

Deconvoluting Protein (Un)folding Structural Ensembles Using X-Ray Scattering,Nuclear Magnetic Resonance Spectroscopy and Molecular Dynamics Simulation

The folding and unfolding of protein domains is an apparently cooperative process, but transient intermediates have been detected in some cases. Such (un)folding intermediates are challenging to investigate structurally as they are typically not long-lived and their role in the (un)folding reaction has often been questioned. One of the most well studied (un)folding pathways is that of Drosophila melanogaster Engrailed homeodomain (EnHD): this 61-residue protein forms a three helix bundle in the native state and folds via a helical intermediate. Here we used molecular dynamics simulations to derive sample conformations of EnHD in the native, intermediate, and unfolded states and selected the relevant structural clusters by comparing to small/wide angle X-ray scattering data at four different temperatures. The results are corroborated using residual dipolar couplings determined by NMR spectroscopy. Our results agree well with the previously proposed (un)folding pathway. However, they also suggest that the fully unfolded state is present at a low fraction throughout the investigated temperature interval, and that the (un)folding intermediate is highly populated at the thermal midpoint in line with the view that this intermediate can be regarded to be the denatured state under physiological conditions. Further, the combination of ensemble structural techniques with MD allows for determination of structures and populations of multiple interconverting structures in solution.     

Effect of solvent conditions upon refolding pathways and intermediates for a simple lattice protein

Folding pathways and intermediates for a two-dimensional lattice protein have been investigated via computer simulation at various denaturant concentrations. The protein is represented as a chain of 8 hydrophobic (H) and 12 polar (P) beads on a square lattice sequenced in such a way that the native state is a compact hydrophobic core surrounded by a shell of polar beads. Two nonbonded H beads are said to attract each other with a potential of mean force of strength ϵ. Increasing |ϵ/kT| mimics decreasing the denaturant concentration in the solution. Dynamic Monte Carlo simulations have been performed in order to investigate the folding transition and the folding pathways. Sharp folding—unfolding transitions are observed and the folding process proceeds along well-defined pathways that are populated by partially folded intermediates. The folding pathways as well as the populations of the intermediates are strongly dependent upon the denaturant concentration. Generally, intermediates containing long open stretches of H beads are more populated at high denaturant concentration, whereas compact intermediates containing a substantial number of hydrophobic contacts are more populated at low denaturant concentrations. The folding process is also observed to be cooperative in nature in that the chain does not start folding until a key fold in the middle section of the chain is formed correctly. © 1997 John Wiley & Sons, Inc. Biopoly 42: 399–409, 1997     

Polymer Uncrossing and Knotting in Protein Folding,and Their Role in Minimal Folding Pathways

We introduce a method for calculating the extent to which chain non-crossing is important in the most efficient, optimal trajectories or pathways for a protein to fold. This involves recording all unphysical crossing events of a ghost chain, and calculating the minimal uncrossing cost that would have been required to avoid such events. A depth-first tree search algorithm is applied to find minimal transformations to fold , , , and knotted proteins. In all cases, the extra uncrossing/non-crossing distance is a small fraction of the total distance travelled by a ghost chain. Different structural classes may be distinguished by the amount of extra uncrossing distance, and the effectiveness of such discrimination is compared with other order parameters. It was seen that non-crossing distance over chain length provided the best discrimination between structural and kinetic classes. The scaling of non-crossing distance with chain length implies an inevitable crossover to entanglement-dominated folding mechanisms for sufficiently long chains. We further quantify the minimal folding pathways by collecting the sequence of uncrossing moves, which generally involve leg, loop, and elbow-like uncrossing moves, and rendering the collection of these moves over the unfolded ensemble as a multiple-transformation “alignment”. The consensus minimal pathway is constructed and shown schematically for representative cases of an , , and knotted protein. An overlap parameter is defined between pathways; we find that proteins have minimal overlap indicating diverse folding pathways, knotted proteins are highly constrained to follow a dominant pathway, and proteins are somewhere in between. Thus we have shown how topological chain constraints can induce dominant pathway mechanisms in protein folding.     

Dynamic Monte Carlo simulations of globular protein folding/unfolding pathways. II. Alpha-helical motifs

Dynamic Monte Carlo simulations of the folding pathways of alpha-helical protein motifs have been undertaken in the context of a diamond lattice model of globular proteins. The first question addressed in the nature of the assembly process of an alpha-helical hairpin. While the hairpin could, in principle, be formed via the diffusion-collision-adhesion of isolated performed helices, this is not the dominant mechanism of assembly found in the simulations. Rather, the helices that form native hairpins are constructed on-site, with folding initiating at or near the turn in almost all cases. Next, the folding/unfolding pathways of four-helix bundles having tight bends and one and two long loops in the native state are explored. Once again, an on-site construction mechanism of folding obtains, with a hairpin forming first, followed by the formation of a three-helix bundle, and finally the fourth helix of the native bundle assembles. Unfolding is essentially the reverse of folding. A simplified analytic theory is developed that reproduces the equilibrium folding transitions obtained from the simulations remarkably well and, for the dominant folding pathway, correctly identifies the intermediates seen in the simulations. The analytic theory provides the free energy along the reaction co-ordinate and identifies the transition state for all three motifs as being quite close to the native state, with three of the four helices assembled, and approximately one turn of the fourth helix in place. The transition state is separated from the native conformation by a free-energy barrier of mainly energetic origin and from the denatured state by a barrier of mainly entropic origin. The general features of the folding pathway seen in all variants of the model four-helix bundles are similar to those observed in the folding of beta-barrel, Greek key proteins; this suggests that many of the qualitative aspects of folding are invariant to the particular native state topology and secondary structure.     

Predicting Secondary Structural Folding Kinetics for Nucleic Acids

We report a new computational approach to the prediction of RNA secondary structure folding kinetics. In this approach, each elementary kinetic step is represented as the transformation between two secondary structures that differ by a helix. Based on the free energy landscape analysis, we identify three types of dominant pathways and the rate constants for the kinetic steps: 1), formation; 2), disruption of a helix stem; and 3), helix formation with concomitant partial melting of a competing (incompatible) helix. The third pathway, termed the tunneling pathway, is the low-barrier dominant pathway for the conversion between two incompatible helices. Comparisons with experimental data indicate that this new method is quite reliable in predicting the kinetics for RNA secondary structural folding and structural rearrangements. The approach presented here may provide a robust first step for further systematic development of a predictive theory for the folding kinetics for large RNAs.     

The nature of protein folding pathways: The classical versus the new view

Summary Pulsed hydrogen exchange and other studies of the kinetic refolding pathways of several small proteins have established that folding intermediates with native-like secondary structures are well populated, but these studies have also shown that the folding kinetics are not well synchronized. Older studies of the kinetics of formation of the native protein, monitored by optical probes, indicate that the folding kinetics should be synchronized. The model commonly used in these studies is the simple sequential model, which postulates a unique folding pathway with defined and sequential intermediates. Theories of the folding process and Monte Carlo simulations of folding suggest that neither the folding pathway nor the set of folding intermediates is unique, and that folding intermediates accumulate because of kinetic traps caused by partial misfolding. Recent experiments with cytochrome c lend support to this new view of folding pathways. These different views of the folding process are discussed. Misfolding and consequent slowing down of the folding process as a result of cis-trans isomerization about prolyl peptide bonds in the unfolded protein are well known; isomerization occurs before refolding is initiated. The occurrence of equilibrium intermediates on the kinetic folding pathways of some proteins, such as -lactalbumin and apomyoglobin, argues that these intermediates are not caused by kinetic traps but rather are stable intermediates under certain conditions, and this conclusion is consistent with a sequential model of folding. Folding reactions with successive kinetic intermediates, in which late intermediates are more highly folded than early intermediates, indicate that folding is hierarchical. New experiments that test the predictions of the classical and the new views are needed.