Pathway diversity and concertedness in protein folding: an ab-initio approach

Making use of an ab-initio folding simulator, we generate in vitro pathways leading to the native fold in moderate size single- domain proteins. The assessment of pathway diversity is not biased by any a priori information on the native fold. We focus on two study cases, hyperthermophile variant of protein G domain (1gb4) and ubiquitin (1ubi), with the same topology but different context dependence in their native folds. We demonstrate that a quenching of structural fluctuations is achieved once the proteins find a stationary plateau maximizing the number of highly protected hydrogen bonds. This enables us to identify the folding nucleus and show that folding does not become expeditious until a concerted event takes place generating a topology able to prevent water attack on a maximal number of hydrogen bonds. This result is consistent with the standard nucleation mechanism postulated for two-state folders. Pathway diversity is correlated with the extent of conflict between local structural propensity and large-scale context, rather than with contact order: In highly context-dependent proteins, the success of folding cannot rely on a single fortuitous event in which local propensity is overruled by large-scale effects. We predict mutational Pi values on individual pathways, compute ensemble averages and predict extent of surface burial and percentage of hydrogen bonding on each component of the transition state ensemble, thus deconvoluting individual folding-route contributions to the averaged two-state kinetic picture. Our predicted kinetic isotopic effects find experimental support and lead to further probes. Finally, the molecular redesign potentiality of the method, aimed at increasing folding expediency, is explored.     

Pathway heterogeneity in protein folding

We generate ab initio folding pathways in two single-domain proteins, hyperthermophile variant of protein G domain (1gb4) and ubiquitin (1ubi), both presumed to be two-state folders. Both proteins are endowed with the same topology but, as shown in this work, rely to a different extent on large-scale context to find their native folds. First, we demonstrate a generic feature of two-state folders: A downsizing of structural fluctuations is achieved only when the protein reaches a stationary plateau maximizing the number of highly protected hydrogen bonds. This enables us to identify the folding nucleus and show that folding does not become expeditious until a topology is generated that is able to protect intramolecular hydrogen bonds from water attack. Pathway heterogeneity is shown to be dependent on the extent to which the protein relies on large-scale context to fold, rather than on contact order: Proteins that can only stabilize native secondary structure by packing it against scaffolding hydrophobic moieties are meant to have a heterogeneous transition-state ensemble if they are to become successful folders (otherwise, successful folding would be too fortuitous an event.) We estimate mutational Phi values as ensemble averages and deconvolute individual-route contributions to the averaged two-state kinetic picture. Our results find experimental corroboration in the well-studied chymotrypsin inhibitor (CI2), while leading to verifiable predictions for the other two study cases.     

Time-resolved backbone desolvation and mutational hot spots in folding proteins

A method is presented to identify hot mutational spots and predict the extent of surface burial at the transition state relative to the native fold in two-state folding proteins. The method is based on ab initio simulations of folding histories in which transitions between coarsely defined conformations and pairwise interactions are dependent on the solvent environments created by the chain. The highly conserved mammalian ubiquitin is adopted as a study case to make predictions. The evolution in time of the chain topology suggests a nucleation process with a critical point signaled by a sudden quenching of structural fluctuations. The occurrence of this nucleus is shown to be concurrent with a sudden escalation in the number of three-body correlations whereby hydrophobic units approach residue pairs engaged in amide-carbonyl hydrogen bonding. These correlations determine a pattern designed to structure the surrounding solvent, protecting intramolecular hydrogen bonds from water attack. Such correlations are shown to be required to stabilize the nucleus, with kinetic consequences for the folding process. Those nuclear residues that adopt the dual role of protecting and being protected while engaged in hydrogen bonds are predicted to be the hottest mutational spots. Some such residues are shown not to retain the same protecting role in the native fold. This kinetic treatment of folding nucleation is independently validated vis-a-vis a Phi-value analysis on chymotrypsin inhibitor 2, a protein for which extensive mutational data exists.     

The Folding Trajectory of RNase H Is Dominated by Its Topology and Not Local Stability: A Protein Engineering Study of Variants that Fold via Two-State and Three-State Mechanisms

Proteins can sample a variety of partially folded conformations during the transition between the unfolded and native states. Some proteins never significantly populate these high-energy states and fold by an apparently two-state process. However, many proteins populate detectable, partially folded forms during the folding process. The role of such intermediates is a matter of considerable debate. A single amino acid change can convert Escherichia coli ribonuclease H from a three-state folder that populates a kinetic intermediate to one that folds in an apparent two-state fashion. We have compared the folding trajectories of the three-state RNase H and the two-state RNase H, proteins with the same native-state topology but altered regional stability, using a protein engineering approach. Our data suggest that both versions of RNase H fold through a similar trajectory with similar high-energy conformations. Mutations in the core and the periphery of the protein affect similar aspects of folding for both variants, suggesting a common trajectory with folding of the core region followed by the folding of the periphery. Our results suggest that formation of specific partially folded conformations may be a general feature of protein folding that can promote, rather than hinder, efficient folding.     

Unification of the folding mechanisms of non-two-state and two-state proteins

We have collected the kinetic folding data for non-two-state and two-state globular proteins reported in the literature, and investigated the relationships between the folding kinetics and the native three-dimensional structure of these proteins. The rate constants of formation of both the intermediate and the native state of non-two-state folders were found to be significantly correlated with protein chain length and native backbone topology, which is represented by the absolute contact order and sequence-distant native pairs. The folding rate of two-state folders, which is known to be correlated with the native backbone topology, apparently does not correlate significantly with protein chain length. On the basis of a comparison of the folding rates of the non-two-state and two-state folders, it was found that they are similarly dependent on the parameters that reflect the native backbone topology. This suggests that the mechanisms behind non-two-state and two-state folding are essentially identical. The present results lead us to propose a unified mechanism of protein folding, in which folding occurs in a hierarchical manner, reflecting the hierarchy of the native three-dimensional structure, as embodied in the case of non-two-state folding with an accumulation of the intermediate. Apparently, two-state folding is merely a simplified version of hierarchical folding caused either by an alteration in the rate-limiting step of folding or by destabilization of the intermediate.     

Structural analysis of the rate-limiting transition states in the folding of Im7 and Im9: similarities and differences in the folding of homologous proteins

The bacterial immunity proteins Im7 and Im9 fold with mechanisms of different kinetic complexity. Whilst Im9 folds in a two-state transition at pH 7.0 and 10 degrees C, Im7 populates an on-pathway intermediate under these conditions. In order to assess the role of sequence versus topology in the folding of these proteins, and to analyse the effect of populating an intermediate on the landscape for folding, we have determined the conformational properties of the rate-limiting transition state for Im9 folding/unfolding using Phi(F)-value analysis and have compared the results with similar data obtained previously for Im7. The data show that the rate-limiting transition states for Im9 and Im7 folding/unfolding are similar: both are compact (beta(T)=0.94 and 0.89, respectively) and contain three of the four native helices docked around a specific hydrophobic core. Significant differences are observed, however, in the magnitude of the Phi(F)-values obtained for the two proteins. Of the 20 residues studied in both proteins, ten have Phi(F)-values in Im7 that exceed those in Im9 by more than 0.2, and of these five differ by more than 0.4. The data suggest that the population of an intermediate in Im7 results in folding via a transition state ensemble that is conformationally restricted relative to that of Im9. The data are consistent with the view that topology is an important determinant of folding. Importantly, however, they also demonstrate that while the folding transition state may be conserved in homologous proteins that fold with two and three-state kinetics, the population of an intermediate can have a significant effect on the breadth of the transition state ensemble.     

Monomer topology defines folding speed of heptamer

Small monomeric proteins often fold in apparent two-state processes with folding speeds dictated by their native-state topology. Here we test, for the first time, the influence of monomer topology on the folding speed of an oligomeric protein: the heptameric cochaperonin protein 10 (cpn10), which in the native state has seven beta-barrel subunits noncovalently assembled through beta-strand pairing. Cpn10 is a particularly useful model because equilibrium-unfolding experiments have revealed that the denatured state in urea is that of a nonnative heptamer. Surprisingly, refolding of the nonnative cpn10 heptamer is a simple two-state kinetic process with a folding-rate constant in water (2.1 sec(-1); pH 7.0, 20 degrees C) that is in excellent agreement with the prediction based on the native-state topology of the cpn10 monomer. Thus, the monomers appear to fold as independent units, with a speed that correlates with topology, although the C and N termini are trapped in beta-strand pairing with neighboring subunits. In contrast, refolding of unfolded cpn10 monomers is dominated by a slow association step.     

First principles prediction of protein folding rates

Experimental studies have demonstrated that many small, single-domain proteins fold via simple two-state kinetics. We present a first principles approach for predicting these experimentally determined folding rates. Our approach is based on a nucleation-condensation folding mechanism, where the rate-limiting step is a random, diffusive search for the native tertiary topology. To estimate the rates of folding for various proteins via this mechanism, we first determine the probability of randomly sampling a conformation with the native fold topology. Next, we convert these probabilities into folding rates by estimating the rate that a protein samples different topologies during diffusive folding. This topology-sampling rate is calculated using the Einstein diffusion equation in conjunction with an experimentally determined intra-protein diffusion constant. We have applied our prediction method to the 21 topologically distinct small proteins for which two-state rate data is available. For the 18 beta-sheet and mixed alpha-beta native proteins, we predict folding rates within an average factor of 4, even though the experimental rates vary by a factor of approximately 4 x 10(4). Interestingly, the experimental folding rates for the three four-helix bundle proteins are significantly underestimated by this approach, suggesting that proteins with significant helical content may fold by a faster, alternative mechanism. This method can be applied to any protein for which the structure is known and hence can be used to predict the folding rates of many proteins prior to experiment.     

Evidence concerning rate-limiting steps in protein folding from the effects of trifluoroethanol

The refolding kinetics of 13 proteins have been studied in the presence of 2,2,2-trifluoroethanol (TFE). Low concentrations of TFE increased the folding rates of all the proteins, whereas higher concentrations have the opposite effect. The extent of deceleration of folding correlates closely with similar effects of guanidine hydrochloride and can be related to the burial of accessible surface area during folding. For those proteins folding in a two-state manner, the extent of acceleration of folding correlates closely with the number of local backbone hydrogen bonds in the native structure. For those proteins that fold in a multistate manner, however, the extent of acceleration is much smaller than that predicted from the data for two-state proteins. These results support the concept that for two-state proteins the search for native-like contacts is a key aspect of the folding reaction, whereas the rate-determining steps for folding of multistate proteins are associated with the reorganization of stable structure within a collapsed state or with the search for native-like interactions within less structured regions.     

Probing possible downhill folding: native contact topology likely places a significant constraint on the folding cooperativity of proteins with approximately 40 residues

Experiments point to appreciable variations in folding cooperativity among natural proteins with approximately 40 residues, indicating that the behaviors of these proteins are valuable for delineating the contributing factors to cooperative folding. To explore the role of native topology in a protein's propensity to fold cooperatively and how native topology might constrain the degree of cooperativity achievable by a given set of physical interactions, we compared folding/unfolding kinetics simulated using three classes of native-centric C^α chain models with different interaction schemes. The approach was applied to two homologous 45-residue fragments from the peripheral subunit-binding domain family and a 39-residue fragment of the N-terminal domain of ribosomal protein L9. Free-energy profiles as functions of native contact number were computed to assess the heights of thermodynamic barriers to folding. In addition, chevron plots of folding/unfolding rates were constructed as functions of native stability to facilitate comparison with available experimental data. Although common Gō-like models with pairwise Lennard-Jones-type interactions generally fold less cooperatively than real proteins, the rank ordering of cooperativity predicted by these models is consistent with experiment for the proteins investigated, showing increasing folding cooperativity with increasing nonlocality of a protein's native contacts. Models that account for water-expulsion (desolvation) barriers and models with many-body (nonadditive) interactions generally entail higher degrees of folding cooperativity indicated by more linear model chevron plots, but the rank ordering of cooperativity remains unchanged. A robust, experimentally valid rank ordering of model folding cooperativity independent of the multiple native-centric interaction schemes tested here argues that native topology places significant constraints on how cooperatively a protein can fold.     

Protein unfolding rates correlate as strongly as folding rates with native structure

Although the folding rates of proteins have been studied extensively, both experimentally and theoretically, and many native state topological parameters have been proposed to correlate with or predict these rates, unfolding rates have received much less attention. Moreover, unfolding rates have generally been thought either to not relate to native topology in the same manner as folding rates, perhaps depending on different topological parameters, or to be more difficult to predict. Using a dataset of 108 proteins including two-state and multistate folders, we find that both unfolding and folding rates correlate strongly, and comparably well, with well-established measures of native topology, the absolute contact order and the long range order, with correlation coefficient values of 0.75 or higher. In addition, compared to folding rates, the absolute values of unfolding rates vary more strongly with native topology, have a larger range of values, and correlate better with thermodynamic stability. Similar trends are observed for subsets of different protein structural classes. Taken together, these results suggest that choosing a scaffold for protein engineering may require a compromise between a simple topology that will fold sufficiently quickly but also unfold quickly, and a complex topology that will unfold slowly and hence have kinetic stability, but fold slowly. These observations, together with the established role of kinetic stability in determining resistance to thermal and chemical denaturation as well as proteases, have important implications for understanding fundamental aspects of protein unfolding and folding and for protein engineering and design.     

Towards a consistent modeling of protein thermodynamic and kinetic cooperativity: how applicable is the transition state picture to folding and unfolding?

To what extent do general features of folding/unfolding kinetics of small globular proteins follow from their thermodynamic properties? To address this question, we investigate a new simplified protein chain model that embodies a cooperative interplay between local conformational preferences and hydrophobic burial. The present four-helix-bundle 55mer model exhibits protein-like calorimetric two-state cooperativity. It rationalizes native-state hydrogen exchange observations. Our analysis indicates that a coherent, self-consistent physical account of both the thermodynamic and kinetic properties of the model leads naturally to the concept of a native state ensemble that encompasses considerable conformational fluctuations. Such a multiple-conformation native state is seen to involve conformational states similar to those revealed by native-state hydrogen exchange. Many of these conformational states are predicted to lie below native baselines commonly used in interpreting calorimetric data. Folding and unfolding kinetics are studied under a range of intrachain interaction strengths as in experimental chevron plots. Kinetically determined transition midpoints match well with their thermodynamic counterparts. Kinetic relaxations are found to be essentially single-exponential over an extended range of model interaction strengths. This includes the entire unfolding regime and a significant part of a folding regime with a chevron rollover, as has been observed for real proteins that fold with non-two-state kinetics. The transition state picture of protein folding and unfolding is evaluated by comparing thermodynamic free energy profiles with actual kinetic rates. These analyses suggest that some chevron rollovers may arise from an internal frictional effect that increasingly impedes chain motions with more native conditions, rather than being caused by discrete deadtime folding intermediates or shifts of the transition state peak as previously posited.     

Protein topology affects the appearance of intermediates during the folding of proteins with a flavodoxin-like fold

The topology of a native protein influences the rate with which it is formed, but does topology affect the appearance of folding intermediates and their specific role in kinetic folding as well? This question is addressed by comparing the folding data recently obtained on apoflavodoxin from Azotobacter vinelandii with those available on all three other alpha-beta parallel proteins the kinetic folding mechanism of which has been studied, i.e. Anabaena apoflavodoxin, Fusarium solani pisi cutinase and CheY. Two kinetic folding intermediates, one on-pathway and the other off-pathway, seem to be present during the folding of proteins with an alpha-beta parallel, also called flavodoxin-like, topology. The on-pathway intermediate lies on a direct route from the unfolded to the native state of the protein involved. The off-pathway intermediate needs to unfold to allow the production of native protein. Available simulation data of the folding of CheY show the involvement of two intermediates with characteristics that resemble those of the two intermediates experimentally observed. Apparently, protein topology governs the appearance and kinetic roles of protein folding intermediates during the folding of proteins that have a flavodoxin-like fold.     

Protein folding and the organization of the protein topology universe

The mechanism by which proteins fold to their native states has been the focus of intense research in recent years. The rate-limiting event in the folding reaction is the formation of a conformation in a set known as the transition-state ensemble. The structural features present within such ensembles have now been analysed for a series of proteins using data from a combination of biochemical and biophysical experiments together with computer-simulation methods. These studies show that the topology of the transition state is determined by a set of interactions involving a small number of key residues and, in addition, that the topology of the transition state is closer to that of the native state than to that of any other fold in the protein universe. Here, we review the evidence for these conclusions and suggest a molecular mechanism that rationalizes these findings by presenting a view of protein folds that is based on the topological features of the polypeptide backbone, rather than the conventional view that depends on the arrangement of different types of secondary-structure elements. By linking the folding process to the organization of the protein structure universe, we propose an explanation for the overwhelming importance of topology in the transition states for protein folding.     

Characterizing the unstructured intermediates in oxidative folding

A recently developed method is used here to characterize some of the folding intermediates, and the oxidative folding processes, of RNase A. This method is based on the ability of trans-[Pt(en)(2)Cl(2)](2+) to oxidize cysteine residues to form disulfide bonds faster than the disulfide bonds can be rearranged by reshuffling or reduction. Variations of this method have enabled us to address three issues. (i) How the nature of the residual structure and/or conformational order that is present, or develops, during the initial stages of folding can be elucidated. It is shown here that there is a 10-fold increase in the propensity of the unfolded reduced forms of RNase A to form the native set of disulfides directly, compared to the propensity under strongly denaturing conditions (4-6 M GdnHCl). Thus, the unfolded reduced forms of RNase A are not statistical coils with a more condensed form than in the GdnHCl-denatured state; rather, it is suggested that reduced RNase A has a little bias toward a native topology. (ii) The structural characterization of oxidative folding intermediates in terms of disulfide pairing is demonstrated; specifically, a lower-limit estimate is made of the percentage of native disulfide-containing molecules in the two-disulfide ensemble of RNase A. (iii) The critical role of structured intermediate species in determining the oxidative folding pathways of proteins was shown previously. Here, we demonstrate that the presence of a structured intermediate in the oxidative folding of proteins can be revealed by this method.     

Acidic conditions stabilise intermediates populated during the folding of Im7 and Im9

The helical bacterial immunity proteins Im7 and Im9 have been shown to fold via kinetic mechanisms of differing complexity, despite having 60 % sequence identity. At pH 7.0 and 10 degrees C, Im7 folds in a three-state mechanism involving an on-pathway intermediate, while Im9 folds in an apparent two-state transition. In order to examine the folding mechanisms of these proteins in more detail, the folding kinetics of both Im7 and Im9 (at 10 degrees C in 0.4 M sodium sulphate) have been examined as a function of pH. Kinetic modelling of the folding and unfolding data for Im7 between pH 5.0 and 8.0 shows that the on-pathway intermediate is stabilised by more acidic conditions, whilst the native state is destabilised. The opposing effect of pH on the stability of these states results in a significant population of the intermediate at equilibrium at pH 6.0 and below. At pH 7.0, the folding and unfolding kinetics for Im9 can be fitted adequately by a two-state model, in accord with previous results. However, under acidic conditions there is a clear change of slope in the plot of the logarithm of the folding rate constant versus denaturant concentration, consistent with the population of one or more intermediate(s) early during folding. The kinetic data for Im9 at these pH values can be fitted to a three-state model, where the intermediate ensemble is stabilised and the native state destabilised as the pH is reduced, rationalising previous results that showed that an intermediate is not observed experimentally at pH 7.0. The data suggest that intermediate formation is a general step in immunity protein folding and demonstrate that it is necessary to explore a wide range of refolding conditions in order to show that intermediates do not form in the folding of other small, single-domain proteins.     

Critical nucleation size in the folding of small apparently two-state proteins

For apparently two-state proteins, we found that the size (number of folded residues) of a transition state is mostly encoded by the topology, defined by total contact distance (TCD) of the native state, and correlates with its folding rate. This is demonstrated by using a simple procedure to reduce the native structures of the 41 two-state proteins with native TCD as a constraint, and is further supported by analyzing the results of eight proteins from protein engineering studies. These results support the hypothesis that the major rate-limiting process in the folding of small apparently two-state proteins is the search for a critical number of residues with the topology close to that of the native state.     

Understanding the Role of Three-Dimensional Topology in Determining the?Folding Intermediates of Group I Introns

Many RNA molecules exert their biological function only after folding to unique three-dimensional structures. For long, noncoding RNA molecules, the complexity of finding the native topology can be a major impediment to correct folding to the biologically active structure. An RNA molecule may fold to a near-native structure but not be able to continue to the correct structure due to a topological barrier such as crossed strands or incorrectly stacked helices. Achieving the native conformation thus requires unfolding and refolding, resulting in a long-lived intermediate. We investigate the role of topology in the folding of two phylogenetically related catalytic group I introns, the Twort and Azoarcus group I ribozymes. The kinetic models describing the Mg²⁺-mediated folding of these ribozymes were previously determined by time-resolved hydroxyl (⋅OH) radical footprinting. Two intermediates formed by parallel intermediates were resolved for each RNA. These data and analytical ultracentrifugation compaction analyses are used herein to constrain coarse-grained models of these folding intermediates as we investigate the role of nonnative topology in dictating the lifetime of the intermediates. Starting from an ensemble of unfolded conformations, we folded the RNA molecules by progressively adding native constraints to subdomains of the RNA defined by the ⋅OH time-progress curves to simulate folding through the different kinetic pathways. We find that nonnative topologies (arrangement of helices) occur frequently in the folding simulations despite using only native constraints to drive the reaction, and that the initial conformation, rather than the folding pathway, is the major determinant of whether the RNA adopts nonnative topology during folding. From these analyses we conclude that biases in the initial conformation likely determine the relative flux through parallel RNA folding pathways.