Protein folding by motion planning

We investigate a novel approach for studying protein folding that has evolved from robotics motion planning techniques called probabilistic roadmap methods (PRMs). Our focus is to study issues related to the folding process, such as the formation of secondary and tertiary structures, assuming we know the native fold. A feature of our PRM-based framework is that the large sets of folding pathways in the roadmaps it produces, in just a few hours on a desktop PC, provide global information about the protein's energy landscape. This is an advantage over other simulation methods such as molecular dynamics or Monte Carlo methods which require more computation and produce only a single trajectory in each run. In our initial studies, we obtained encouraging results for several small proteins. In this paper, we investigate more sophisticated techniques for analyzing the folding pathways in our roadmaps. In addition to more formally revalidating our previous results, we present a case study showing that our technique captures known folding differences between the structurally similar proteins G and L.     

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.     

Folding path of P5abc RNA involves direct coupling of secondary and tertiary structures

Folding mechanisms in which secondary structures are stabilized through the formation of tertiary interactions are well documented in protein folding but challenge the folding hierarchy normally assumed for RNA. However, it is increasingly clear that RNA could fold by a similar mechanism. P5abc, a small independently folding tertiary domain of the Tetrahymena thermophila group I ribozyme, is known to fold by a secondary structure rearrangement involving helix P5c. However, the extent of this rearrangement and the precise stage of folding that triggers it are unknown. We use experiments and simulations to show that the P5c helix switches to the native secondary structure late in the folding pathway and is directly coupled to the formation of tertiary interactions in the A-rich bulge. P5c mutations show that the switch in P5c is not rate-determining and suggest that non-native interactions in P5c aid folding rather than impede it. Our study illustrates that despite significant differences in the building blocks of proteins and RNA, there may be common ways in which they self-assemble.     

Early collapse is not an obligate step in protein folding

The dimensions and secondary structure content of two proteins which fold in a two-state manner are measured within milliseconds of denaturant dilution using synchrotron-based, stopped-flow small-angle X-ray scattering and far-UV circular dichroism spectroscopy. Even upon a jump to strongly native conditions, neither ubiquitin nor common-type acylphosphatase contract prior to the major folding event. Circular dichroism and fluorescence indicate that negligible amounts of secondary and tertiary structures form in the burst phase. Thus, for these two denatured states, collapse and secondary structure formation are not energetically downhill processes even under aqueous, low-denaturant conditions. In addition, water appears to be as good a solvent as that with high concentrations of denaturant, when considering the over-all dimensions of the denatured state. However, the removal of denaturant does subtly alter the distribution of backbone dihedral phi,psi angles, most likely resulting in a shift from the polyproline II region to the helical region of the Ramachandran map. We consider the thermodynamic origins of these behaviors along with implications for folding mechanisms and computer simulations thereof.     

Is there a unifying mechanism for protein folding?

Proteins appear to fold by diverse pathways, but variations of a simple mechanism - nucleation-condensation - describe the overall features of folding of most domains. In general, secondary structure is inherently unstable and its stability is enhanced by tertiary interactions. Consequently, an extensive interplay of secondary and tertiary interactions determines the transition-state for folding, which is structurally similar to the native state, being formed in a general collapse (condensation) around a diffuse nucleus. As the propensity for stable secondary structure increases, folding becomes more hierarchical and eventually follows a framework mechanism where the transition state is assembled from pre-formed secondary structural elements.     

The denatured state dictates the topology of two proteins with almost identical sequence but different native structure and function

The protein folding problem is often studied by comparing the mechanisms of proteins sharing the same structure but different sequence. The recent design of the two proteins G_A88 and G_B88, displaying different structures and functions while sharing 88% sequence identity (49 out of 56 amino acids), allows the unique opportunity for a complementary approach. At which stage of its folding pathway does a protein commit to a given topology? Which residues are crucial in directing folding mechanisms to a given structure? By using a combination of biophysical and computational techniques, we have characterized the folding of both G_A88 and G_B88. We show that, contrary to expectation, G_B88, characterized by a native α+β fold, displays in the denatured state a content of native-like helical structure greater than G_A88, which is all-α in its native state. Both experiments and simulations indicate that such residual structure may be tuned by changing pH. Thus, despite the high sequence identity, the folding pathways for these two proteins appear to diverge as early as in the denatured state. Our results suggest a mechanism whereby protein topology is committed very early along the folding pathway, being imprinted in the residual structure of the denatured state.     

Atomic-level structure characterization of an ultrafast folding mini-protein denatured state

Atomic-level analyses of non-native protein ensembles constitute an important aspect of protein folding studies to reach a more complete understanding of how proteins attain their native form exhibiting biological activity. Previously, formation of hydrophobic clusters in the 6 M urea-denatured state of an ultrafast folding mini-protein known as TC5b from both photo-CIDNP NOE transfer studies and FCS measurements was observed. Here, we elucidate the structural properties of this mini-protein denatured in 6 M urea performing (15)N NMR relaxation studies together with a thorough NOE analysis. Even though our results demonstrate that no elements of secondary structure persist in the denatured state, the heterogeneous distribution of R(2) rate constants together with observing pronounced heteronuclear NOEs along the peptide backbone reveals specific regions of urea-denatured TC5b exhibiting a high degree of structural rigidity more frequently observed for native proteins. The data are complemented with studies on two TC5b point mutants to verify the importance of hydrophobic interactions for fast folding. Our results corroborate earlier findings of a hydrophobic cluster present in urea-denatured TC5b comprising both native and non-native contacts underscoring their importance for ultra rapid folding. The data assist in finding ways of interpreting the effects of pre-existing native and/or non-native interactions on the ultrafast folding of proteins; a fact, which might have to be considered when defining the starting conditions for molecular dynamics simulation studies of protein folding.     

Multiple parallel-pathway folding of proline-free Staphylococcal nuclease

When a protein exhibits complex kinetics of refolding, we often ascribe the complexity to slow isomerization events in the denatured protein, such as cis/trans isomerization of peptidyl prolyl bonds. Does the complex folding kinetics arise only from this well-known reason? Here, we have investigated the refolding of a proline-free variant of staphylococcal nuclease by stopped-flow, double-jump techniques, to examine the folding reactions without the slow prolyl isomerizations. As a result, the protein folds into the native state along at least two accessible parallel pathways, starting from a macroscopically single denatured-state ensemble. The presence of intermediates on the individual folding pathways has revealed the existence of multiple parallel pathways, and is characterized by multi-exponential folding kinetics with a lag phase. Therefore, a "single" amino acid sequence can fold along the multiple parallel pathways. This observation in staphylococcal nuclease suggests that the multiple folding may be more general than we have expected, because the multiple parallel-pathway folding cannot be excluded from proteins that show simpler kinetics.     

New pathways in folding of the Tetrahymena group I RNA enzyme.

Previous studies have shown that the earliest detectable step in folding of the Tetrahymena ribozyme is tertiary structure formation of the peripheral element P5abc. This, along with other results, has suggested that P5abc may serve as a scaffold upon which additional tertiary structure is built. Herein we use the onset of oligonucleotide cleavage activity as a readout for native state formation and investigate the effect of P5abc on the rate of folding to the native structure. Despite the early folding of P5abc, its removal to give the E delta P5abc variant decreases the rate of attainment of an active structure less than fivefold (20-100 mM Mg2+, 15-50 degrees C). Furthermore, P5abc added in trans is able to bind the folded E delta P5abc ribozyme and promote oligonucleotide cleavage at least tenfold more rapidly than folding of the wild-type ribozyme, indicating that E delta P5abc does not have to first unfold before productively binding P5abc to form the true native state. This suggests that a state with the overall tertiary structure formed but with P5abc unfolded represents a viable on-pathway intermediate for the wild-type ribozyme. These results provide strong evidence for the existence of two pathways to the native state: in one pathway P5abc forms tertiary structure first, and in another it forms late. The pathway in which P5abc forms first is favored because P5abc can fold quickly and because its tertiary structure is stable in the absence of additional structured elements, not because P5abc formation is required for subsequent folding steps. In the course of these experiments, we also found that most of the ribozyme population does not reach the native state directly under standard conditions in vitro, but instead forms an inactive structure that is stable for hours. Finally, the fraction that does fold to the native state folds with a single rate constant of 1 min-1, suggesting that there are no significantly populated "fast-track" pathways that reach the native state directly by avoiding slow folding steps.     

On the origins of the weak folding cooperativity of a designed ββα ultrafast protein FSD-1

FSD-1, a designed small ultrafast folder with a ββα fold, has been actively studied in the last few years as a model system for studying protein folding mechanisms and for testing of the accuracy of computational models. The suitability of this protein to describe the folding of naturally occurring α/β proteins has recently been challenged based on the observation that the melting transition is very broad, with ill-resolved baselines. Using molecular dynamics simulations with the AMBER protein force field (ff96) coupled with the implicit solvent model (IGB = 5), we shed new light into the nature of this transition and resolve the experimental controversies. We show that the melting transition corresponds to the melting of the protein as a whole, and not solely to the helix-coil transition. The breadth of the folding transition arises from the spread in the melting temperatures (from ∼325 K to ∼302 K) of the individual transitions: formation of the hydrophobic core, β-hairpin and tertiary fold, with the helix formed earlier. Our simulations initiated from an extended chain accurately predict the native structure, provide a reasonable estimate of the transition barrier height, and explicitly demonstrate the existence of multiple pathways and multiple transition states for folding. Our exhaustive sampling enables us to assess the quality of the Amber ff96/igb5 combination and reveals that while this force field can predict the correct native fold, it nonetheless overstabilizes the α-helix portion of the protein (Tm = ∼387K) as well as the denatured structures.     

Residue-level NMR view of the urea-driven equilibrium folding transition of SUMO-1 (1-97): native preferences do not increase monotonously

SUMO-1 (1-97) is a crucial protein in the machinery of post-translational modifications. We observed by circular dichroism and fluorescence spectroscopy that urea-induced unfolding of this protein is a complex process with the possibility of occurrence of detectable intermediates along the way. The tertiary structure is completely lost around approximately 4.5 M urea with a transition mid-point at 2.53 M urea, while the secondary structure unfolding seems to show two transitions, with mid-points at 2.42 M and 5.69 M urea. We have elucidated by systematic urea titration, the equilibrium residue level structural and dynamics changes along the entire folding/unfolding transition by multidimensional NMR. With urea dilution, the protein is seen to progressively lose most of the broad beta-domain structural preferences present at 8 M urea, acquire some helical propensities at 5 M urea, and lose some of them again on further dilution of urea. Between 3 M and 2 M urea, the protein starts afresh to acquire native structural features. These observations are contrary to the conventional notion that proteins fold with monotonously increasing native-type preferences. For folding below approximately 3 M urea, the region around the alpha1 helix appears to be a potential folding initiation site. The folding seems to start with a collapse into native-like topologies, at least in parts, and is followed by formation of secondary and tertiary structure, perhaps by cooperative rearrangements. The motional characteristics of the protein show sequence-dependent variation as the concentration of urea is progressively reduced. At the sub-nanosecond level, the features are extremely unusual for denatured states, and only certain segments corresponding to the flexible regions in the native protein display these motions at the different concentrations of urea.     

Fast in silico protein folding by introduction of alternating hydrogen bond potentials

We accelerate protein folding in all-atom molecular dynamics simulations by introducing alternating hydrogen bond potentials as a supplement to the force field. The alternating hydrogen bond potentials result in accelerated hydrogen bond reordering, which leads to rapid formation of secondary structure elements. The method does not require knowledge of the native state but generates the potentials based on the development of the tertiary structure in the simulation. In protein folding, the formation of secondary structure elements, especially α-helix and β-sheet, is very important, and we show that our method can fold both efficiently and with great speed.     

Prediction of protein folding pathways.

Recent 1H nuclear magnetic resonance (n.m.r.) hydrogen exchange experiments on five different proteins have delineated the secondary structures formed in trapped, partially folded intermediates. The early forming structural elements are identifiable through a technique described in this work to predict folding pathways. The method assumes that the sequential selection of structural fragments such as alpha-helices and beta-strands involved in the folding process is founded upon the maximal burial of solvent accessible surface from both the formation of internal structure and substructure association. The substructural elements were defined objectively by major changes in main-chain direction. The predicted folding pathways are in complete correspondence with the n.m.r. results in that the formed structural fragments found in the folding intermediates are those predicted earliest in the pathways. The technique was also applied to proteins of known tertiary structure and with fold similar to one of the five proteins examined by 1H n.m.r. The pathways for these structures also showed general consistency with the n.m.r. observations, suggesting conservation of a secondary structural framework or molten globule about which folding nucleates and proceeds.     

Using motion planning to map protein folding landscapes and analyze folding kinetics of known native structures.

We investigate a novel approach for studying the kinetics of protein folding. Our framework has evolved from robotics motion planning techniques called probabilistic roadmap methods (PRMs) that have been applied in many diverse fields with great success. In our previous work, we presented our PRM-based technique and obtained encouraging results studying protein folding pathways for several small proteins. In this paper, we describe how our motion planning framework can be used to study protein folding kinetics. In particular, we present a refined version of our PRM-based framework and describe how it can be used to produce potential energy landscapes, free energy landscapes, and many folding pathways all from a single roadmap which is computed in a few hours on a desktop PC. Results are presented for 14 proteins. Our ability to produce large sets of unrelated folding pathways may potentially provide crucial insight into some aspects of folding kinetics, such as proteins that exhibit both two-state and three-state kinetics that are not captured by other theoretical techniques.     

Folding and oxidation of the antibody domain C(H)3

The non-covalent homodimer formed by the C-terminal domains of the IgG1 heavy chains (C(H)3) is the simplest naturally occurring model system for studying immunoglobulin folding and assembly. In the native state, the intrachain disulfide bridge, which connects a three-stranded and a four-stranded beta-sheet is buried in the hydrophobic core of the protein. Here, we show that the disulfide bridge is not required for folding and association, since the reduced C(H)3 domain folds to a dimer with defined secondary and tertiary structure. However, the thermodynamic stability of the reduced C(H)3 dimer is much lower than that of the oxidized state. This allows the formation of disulfide bonds either concomitant with folding (starting from the reduced, denatured state) or after folding (starting from the reduced dimer). The analysis of the two processes revealed that, under all conditions investigated, one of the cysteine residues, Cys 86, reacts preferentially with oxidized glutathione to a mixed disulfide that subsequently interacts with the less-reactive second thiol group of the intra-molecular disulfide bond. For folded C(H)3, the second step in the oxidation process is slow. In contrast, starting from the unfolded and reduced protein, the oxidation reaction is faster. However, the overall folding reaction of C(H)3 during oxidative folding is a slow process. Especially, dimerization is slow, compared to the association starting from the denatured oxidized state. This deceleration may be due to misfolded conformations trapped by the disulfide bridge.     

Cooperative folding of the isolated alpha-helical domain of hen egg-white lysozyme.

Proteins in the alpha-lactalbumin and c-type lysozyme family have been studied extensively as model systems in protein folding. Early formation of the alpha-helical domain is observed in both alpha-lactalbumin and c-type lysozyme; however, the details of the kinetic folding pathways are significantly different. The major folding intermediate of hen egg-white lysozyme has a cooperatively formed tertiary structure, whereas the intermediate of alpha-lactalbumin exhibits the characteristics of a molten globule. In this study, we have designed and constructed an isolated alpha-helical domain of hen egg-white lysozyme, called Lyso-alpha, as a model of the lysozyme folding intermediate that is stable at equilibrium. Disulfide-exchange studies show that under native conditions, the cysteine residues in Lyso-alpha prefer to form the same set of disulfide bonds as in the alpha-helical domain of full-length lysozyme. Under denaturing conditions, formation of the nearest-neighbor disulfide bonds is strongly preferred. In contrast to the isolated alpha-helical domain of alpha-lactalbumin, Lyso-alpha with two native disulfide bonds exhibits a well-defined tertiary structure, as indicated by cooperative thermal unfolding and a well-dispersed NMR spectrum. Thus, the determinants for formation of the cooperative side-chain interactions are located mainly in the alpha-helical domain. Our studies suggest that the difference in kinetic folding pathways between alpha-lactalbumin and lysozyme can be explained by the difference in packing density between secondary structural elements and support the hypothesis that the structured regions in a protein folding intermediate may correspond to regions that can fold independently.     

Protein folding includes oligomerization - examples from the endoplasmic reticulum and cytosol

A correct three-dimensional structure is a prerequisite for protein functionality, and therefore for life. Thus, it is not surprising that our cells are packed with proteins that assist protein folding, the process in which the native three-dimensional structure is formed. In general, plasma membrane and secreted proteins, as well as those residing in compartments along the endocytic and exocytic pathways, fold and oligomerize in the endoplasmic reticulum. The proteins residing in the endoplasmic reticulum are specialized in the folding of this subset of proteins, which renders this compartment a protein-folding factory. This review focuses on protein folding in the endoplasmic reticulum, and discusses the challenge of oligomer formation in the endoplasmic reticulum as well as the cytosol.     

Folding kinetics of the outer membrane proteins OmpA and FomA into phospholipid bilayers

The folding mechanism of outer membrane proteins (OMPs) of Gram-negative bacteria into lipid bilayers has been studied using OmpA of E. coli and FomA of F. nucleatum as examples. Both, OmpA and FomA are soluble in unfolded form in urea and insert and fold into phospholipid bilayers upon strong dilution of the denaturant urea. OmpA is a structural protein and forms a small ion channel, composed of an 8-stranded transmembrane beta-barrel domain. FomA is a voltage-dependent porin, predicted to form a 14 stranded beta-barrel. Both OMPs fold into a range of model membranes of very different phospholipid compositions. Three membrane-bound folding intermediates of OmpA were discovered in folding studies with dioleoylphosphatidylcholine bilayers that demonstrated a highly synchronized mechanism of secondary and tertiary structure formation of beta-barrel membrane proteins. A study on FomA folding into lipid bilayers indicated the presence of parallel folding pathways for OMPs with larger transmembrane beta-barrels.     

Protein folding and misfolding: mechanism and principles

Two fundamentally different views of how proteins fold are now being debated. Do proteins fold through multiple unpredictable routes directed only by the energetically downhill nature of the folding landscape or do they fold through specific intermediates in a defined pathway that systematically puts predetermined pieces of the target native protein into place? It has now become possible to determine the structure of protein folding intermediates, evaluate their equilibrium and kinetic parameters, and establish their pathway relationships. Results obtained for many proteins have serendipitously revealed a new dimension of protein structure. Cooperative structural units of the native protein, called foldons, unfold and refold repeatedly even under native conditions. Much evidence obtained by hydrogen exchange and other methods now indicates that cooperative foldon units and not individual amino acids account for the unit steps in protein folding pathways. The formation of foldons and their ordered pathway assembly systematically puts native-like foldon building blocks into place, guided by a sequential stabilization mechanism in which prior native-like structure templates the formation of incoming foldons with complementary structure. Thus the same propensities and interactions that specify the final native state, encoded in the amino-acid sequence of every protein, determine the pathway for getting there. Experimental observations that have been interpreted differently, in terms of multiple independent pathways, appear to be due to chance misfolding errors that cause different population fractions to block at different pathway points, populate different pathway intermediates, and fold at different rates. This paper summarizes the experimental basis for these three determining principles and their consequences. Cooperative native-like foldon units and the sequential stabilization process together generate predetermined stepwise pathways. Optional misfolding errors are responsible for 3-state and heterogeneous kinetic folding.