ATP-triggered conformational changes delineate substrate-binding and -folding mechanics of the GroEL chaperonin

The chaperonin GroEL assists the folding of nascent or stress-denatured polypeptides by actions of binding and encapsulation. ATP binding initiates a series of conformational changes triggering the association of the cochaperonin GroES, followed by further large movements that eject the substrate polypeptide from hydrophobic binding sites into a GroES-capped, hydrophilic folding chamber. We used cryo-electron microscopy, statistical analysis, and flexible fitting to resolve a set of distinct GroEL-ATP conformations that can be ordered into a trajectory of domain rotation and elevation. The initial conformations are likely to be the ones that capture polypeptide substrate. Then the binding domains extend radially to separate from each other but maintain their binding surfaces facing the cavity, potentially exerting mechanical force upon kinetically trapped, misfolded substrates. The extended conformation also provides a potential docking site for GroES, to trigger the final, 100° domain rotation constituting the "power stroke" that ejects substrate into the folding chamber.     

Free-energy landscape of the villin headpiece in an all-atom force field

We investigate the landscape of the internal free-energy of the 36 amino acid villin headpiece with a modified basin hopping method in the all-atom force field PFF01, which was previously used to predictively fold several helical proteins with atomic resolution. We identify near native conformations of the protein as the global optimum of the force field. More than half of the twenty best simulations started from random initial conditions converge to the folding funnel of the native conformation, but several competing low-energy metastable conformations were observed. From 76,000 independently generated conformations we derived a decoy tree which illustrates the topological structure of the entire low-energy part of the free-energy landscape and characterizes the ensemble of metastable conformations. These emerge as similar in secondary content, but differ in tertiary arrangement.     

Visualization of Protein Folding Funnels in Lattice Models

Protein folding occurs in a very high dimensional phase space with an exponentially large number of states, and according to the energy landscape theory it exhibits a topology resembling a funnel. In this statistical approach, the folding mechanism is unveiled by describing the local minima in an effective one-dimensional representation. Other approaches based on potential energy landscapes address the hierarchical structure of local energy minima through disconnectivity graphs. In this paper, we introduce a metric to describe the distance between any two conformations, which also allows us to go beyond the one-dimensional representation and visualize the folding funnel in 2D and 3D. In this way it is possible to assess the folding process in detail, e.g., by identifying the connectivity between conformations and establishing the paths to reach the native state, in addition to regions where trapping may occur. Unlike the disconnectivity maps method, which is based on the kinetic connections between states, our methodology is based on structural similarities inferred from the new metric. The method was developed in a 27-mer protein lattice model, folded into a 3×3×3 cube. Five sequences were studied and distinct funnels were generated in an analysis restricted to conformations from the transition-state to the native configuration. Consistent with the expected results from the energy landscape theory, folding routes can be visualized to probe different regions of the phase space, as well as determine the difficulty in folding of the distinct sequences. Changes in the landscape due to mutations were visualized, with the comparison between wild and mutated local minima in a single map, which serves to identify different trapping regions. The extension of this approach to more realistic models and its use in combination with other approaches are discussed.     

Funnel hunting in a rough terrain: learning and discriminating native energy funnels

Protein folding and binding is commonly depicted as a search for the minimum energy conformation. Modeling of protein complex structures by RosettaDock often results in a set of low-energy conformations near the native structure. Ensembles of low-energy conformations can appear, however, in other regions, especially when backbone movements occur upon binding. What then characterizes the energy landscape near the correct orientation? We applied a machine learning algorithm to distinguish ensembles of low-energy conformations around the native conformation from other low-energy ensembles. The resulting classifier, FunHunt, identifies the native orientation in 50/52 protein complexes in a test set. The features used by FunHunt teach us about the nature of native interfaces. Remarkably, the energy decrease of trajectories toward near-native orientations is significantly larger than for other orientations. This provides a possible explanation for the stability of association in the native orientation.     

Folding funnels, binding funnels, and protein function.

Folding funnels have been the focus of considerable attention during the last few years. These have mostly been discussed in the general context of the theory of protein folding. Here we extend the utility of the concept of folding funnels, relating them to biological mechanisms and function. In particular, here we describe the shape of the funnels in light of protein synthesis and folding; flexibility, conformational diversity, and binding mechanisms; and the associated binding funnels, illustrating the multiple routes and the range of complexed conformers. Specifically, the walls of the folding funnels, their crevices, and bumps are related to the complexity of protein folding, and hence to sequential vs. nonsequential folding. Whereas the former is more frequently observed in eukaryotic proteins, where the rate of protein synthesis is slower, the latter is more frequent in prokaryotes, with faster translation rates. The bottoms of the funnels reflect the extent of the flexibility of the proteins. Rugged floors imply a range of conformational isomers, which may be close on the energy landscape. Rather than undergoing an induced fit binding mechanism, the conformational ensembles around the rugged bottoms argue that the conformers, which are most complementary to the ligand, will bind to it with the equilibrium shifting in their favor. Furthermore, depending on the extent of the ruggedness, or of the smoothness with only a few minima, we may infer nonspecific, broad range vs. specific binding. In particular, folding and binding are similar processes, with similar underlying principles. Hence, the shape of the folding funnel of the monomer enables making reasonable guesses regarding the shape of the corresponding binding funnel. Proteins having a broad range of binding, such as proteolytic enzymes or relatively nonspecific endonucleases, may be expected to have not only rugged floors in their folding funnels, but their binding funnels will also behave similarly, with a range of complexed conformations. Hence, knowledge of the shape of the folding funnels is biologically very useful. The converse also holds: If kinetic and thermodynamic data are available, hints regarding the role of the protein and its binding selectivity may be obtained. Thus, the utility of the concept of the funnel carries over to the origin of the protein and to its function.     

Noninteracting local-structure model of folding and unfolding transition in globular proteins. I. Formulation

A statistical-mechanical model (a noninteracting local structure model) of folding and unfolding transition in globular proteins is described and a formulation is given to calculate the partition function. The process of transition is discussed in this model within the framework of equilibrium statistical mechanics. In order to clarify the range of applicability of such an approach, the characteristics of the folding and unfolding transition in globular proteins are analyzed from the statistical-physical point of view. A theoretical advantage is pointed out in studying folding and unfolding processes taking place as conformational fluctuations in individual protein molecules under macroscopic equilibrium at the melting temperature. In this case, paths of folding and unfolding are shown to be identical in the statistical sense. A key to the noninteracting local structure model lies in the concept of local structures and the assumption of the absence of interactions between local structures. A local structure is defined as a continuous section of the chain which takes the same or similar local conformation as in the native conformation. The assumption of the absence of inter-actions between local structures endows the model with the remarkable character that its partition function can be calculated exactly; thereby the equilibrium population of various conformations along the folding and unfolding paths can be discussed only by a knowledge of the folded native conformation.     

Protein folding in the cell: reshaping the folding funnel

Models of protein folding have historically focused on a subset of 'well-behaved' proteins that can be successfully refolded from denaturants in vitro. Energy landscapes, including folding funnel 'cartoons', describe the largely uncomplicated folding of these isolated chains at infinite dilution. However, the frequent failure of many polypeptides to fold to their native state requires more comprehensive models of folding to accommodate the crucial role of interactions between partially folded intermediates. By incorporating additional deep minima, which reflect off-pathway interchain interactions, the folding funnel concept can be extended to describe the behavior of a more diverse set of proteins under more physiologically relevant conditions. In particular, the effects of ribosomes (translation), molecular chaperones and other aspects of the cellular environment on early chain conformations can be included to account for the folding behavior of polypeptide chains in cells.     

Relaxation folding and the principle of the minimum rate of energy dissipation for conformational motions in a viscous medium

A numerical simulation of the folding of a model polymer chain of 50 units with valence bonds of a fixed length and fixed valence angle values has been performed using the strong friction approximation. The rate of energy dissipation in the system has been analyzed for conformational motions along a trajectory determined by the equations of mechanics and the trajectories characterized by random and variable deviations from the mechanical path. The validity of the principle of the minimum average rate of the energy dissipation for the conformational relaxation of a macromolecule in a viscous medium has been demonstrated. A profile of the relaxation energy funnel for the folding of a macromolecular chain has been constructed. Slow and rapid stages of folding could be distinguished in the energy funnel profile; the final state was separated from the nearest conformations of the folded chain by an energy gap.     

Circular dichroic investigation of the native and non-native conformational states of the growth factor receptor-binding protein 2 N-terminal src homology domain 3: effect of binding to a proline-rich peptide from guanine nucleotide exchange factor

SH3 (src homology domain 3) domains are small protein modules that interact with proline-rich peptides. The structure of the N-terminal SH3 domain from growth factor receptor-binding protein 2 (Grb2), an adapter protein in the intracellular signaling pathway to Ras, was investigated by circular dichroic (CD) spectroscopy. The compact native beta-barrel conformation, previously elucidated by NMR spectroscopy, was largely predominant at pH = 4.8, in the absence of salt. From the structural changes induced by varying pH, ionic strength, temperature, or hydrophobicity of the environment, evidence for the existence of distinct nonnative conformations was obtained in the far- and near-UV domains. Along the free energy scale, these appear to distribute into two conformational ensembles, depending on the extent of structural and thermodynamic differences compared to the native conformation. The first ensemble consists of non-native conformations with a nativelike secondary structure, and the second is composed of partially unfolded conformations having short alpha-helical fragments or turnlike motifs in their nonnative secondary structure. Most of the observed nonnative conformations exist in mild or nondenaturing conditions. They probably have distinct compactness of their inner structure, depending on the strength of nonlocal interactions, but only the native all-beta conformation possesses a condensed protein exterior, appropriate for the binding to the VPPPVPPRRR decapeptide from Sos. Upon binding, the native conformation undergoes a local tertiary structure change in a hydrophobic pocket at the binding site. This is accompanied by the PP-II helix folding of the proline-rich peptide. Interestingly, in the near-UV domain, a significant change in the spectral contribution of an aromatic exciton was observed, thus allowing quantitative tracking of the binding process.     

Structured disorder and conformational selection

Traditionally, molecular disorder has been viewed as local or global instability. Molecules or regions displaying disorder have been considered inherently unstructured. The term has been routinely applied to cases for which no atomic coordinates can be derived from crystallized molecules. Yet, even when it appears that the molecules are disordered, prevailing conformations exist, with population times higher than those of all alternate conformations. Disordered molecules are the outcome of rugged energy landscapes away from the native state around the bottom of the funnel. Ruggedness has a biological function, creating a distribution of structured conformers that bind via conformational selection, driving association and multimolecular complex formation, whether chain-linked in folding or unlinked in binding. We classify disordered molecules into two types. The first type possesses a hydrophobic core. Here, even if the native conformation is unstable, it still has a large enough population time, enabling its experimental detection. In the second type, no such hydrophobic core exists. Hence, the native conformations of molecules belonging to this category have shorter population times, hindering their experimental detection. Although there is a continuum of distribution of hydrophobic cores in proteins, an empirical, statistically based hydrophobicity function may be used as a guideline for distinguishing the two disordered molecule types. Furthermore, the two types relate to steps in the protein folding reaction. With respect to protein design, this leads us to propose that engineering-optimized specific electrostatic interactions to avoid electrostatic repulsion would reduce the type I disordered state, driving the molten globule (MG) --> native (N) state. In contrast, for overcoming the type II disordered state, in addition to specific interactions, a stronger hydrophobic core is also indicated, leading to the denatured --> MG --> N state.     

A size dependent folding contour for cytochrome C

The paper describes an experimental construct of the folding route of the heme protein cytochrome-C. The construct highlights a slowing down near the nose of the folding funnel caused by the multiplicity of the energy traps near the native conformation created as a result of complex heme-peptide interaction. Interestingly the hydrodynamic size, the size heterogeneity and peroxidase activity serve as a triple measure of the distance of this near equilibrium departure from native conformation. Accordingly, the folding process is marked with a gradual and reversible reduction of mean hydrodynamic size, size heterogeneity and peroxidase activity (higher in unfolded state). The Dynamic Light Scattering based straightforward illustration of hydrodynamic size variation may serve as a model to slow folding observed in case of heme proteins, the heme itself serving as a natural facilitator for the native peptide conformation.     

Biochemical characterization of different conformational states of the Sf9 cell-purified p53His175 mutant protein

In this study, we expressed and purified the p53 mutant encoded by the His175 allele (p53His175) in a baculovirus expression system in order to study the folding and the DNA binding activity of the protein. A two-site ELISA revealed that purified p53His175 protein preferentially displayed a PAb1620 conformation, which appeared to be not sufficient to interact specifically with DNA. The cryptic DNA binding activity of this mutant was then investigated by electrophoretic mobility shift assay in the presence of anti-p53 antibodies, and shown to be refractory to significant activation by PAb421 (a potent allosteric activator of wild-type p53's DNA binding activity). Nevertheless, p53His175 DNA binding was regulated by antibodies targeting the N-terminal region of the protein. Furthermore, while the protein preferentially displayed a PAb1620 conformation, our data suggested the existence of an equilibrium between at least two folding states of the protein (PAb1620 and PAb240 conformations). A model rationalizing the conformation, antibody-interacting ability and DNA binding regulation potential of p53His175 is presented.     

Characterization of cytochrome c folding intermediates induced by sucrose and phosphate

Spectroscopic and calorimetric investigations of the folding of denatured cytochrome c in the presence of phosphate ion and sugar were carried out to understand subtle differences in the nature of induced conformation and folding energy landscape. Altered conformations of cyt c induced by sucrose and phosphate, with same absorbance wavelength maxima, exhibit lack of tertiary interactions in segment 70-85 and similar α-helical content. However, compactness, the exposure of the heme to solvent and the secondary structure content in the two conformations are different. Although downhill folding was observed for both conformations, extent of cooperativity is higher in case of phosphate-induced conformation.     

Cotranslational folding promotes beta-helix formation and avoids aggregation in vivo

Newly synthesized proteins must form their native structures in the crowded environment of the cell, while avoiding non-native conformations that can lead to aggregation. Yet, remarkably little is known about the progressive folding of polypeptide chains during chain synthesis by the ribosome or of the influence of this folding environment on productive folding in vivo. P22 tailspike is a homotrimeric protein that is prone to aggregation via misfolding of its central β-helix domain in vitro. We have produced stalled ribosome:tailspike nascent chain complexes of four fixed lengths in vivo, in order to assess cotranslational folding of newly synthesized tailspike chains as a function of chain length. Partially synthesized, ribosome-bound nascent tailspike chains populate stable conformations with some native-state structural features even prior to the appearance of the entire β-helix domain, regardless of the presence of the chaperone trigger factor, yet these conformations are distinct from the conformations of released, refolded tailspike truncations. These results suggest that organization of the aggregation-prone β-helix domain occurs cotranslationally, prior to chain release, to a conformation that is distinct from the accessible energy minimum conformation for the truncated free chain in solution.     

Induced fit of RNA on binding the L7Ae protein to the kink-turn motif

The kink-turn is a widespread motif in RNA consisting of a three-nucleotide bulge flanked on one side by consecutive A3G mismatches. Important examples are found in the ribosome, U4 RNA, and in snoRNAs involved in RNA modification. The motif is a common protein binding site, and the RNA has been found to adopt a tightly kinked conformation in crystal structures. However, in free solution there is a dynamic exchange between kinked and extended conformations, with the equilibrium driven toward the kinked form by the addition of metal ions. Here we used fluorescence resonance energy transfer (FRET) to show that the L7Ae protein of Archaeoglobus fulgidus binds to RNA containing a kink-turn with nanomolar affinity, and induces folding into the tightly kinked conformation even in the absence of metal ions. Thus this RNA may act as a relatively flexible hinge during RNA folding, until fixed into its ultimate kinked structure by the binding of L7 or related protein.     

Induced folding in RNA recognition by Arabidopsis thaliana DCL1

DCL1 is the ribonuclease that carries out miRNA biogenesis in plants. The enzyme has two tandem double stranded RNA binding domains (dsRBDs) in its C-terminus. Here we show that the first of these domains binds precursor RNA fragments when isolated and cooperates with the second domain in the recognition of substrate RNA. Remarkably, despite showing RNA binding activity, this domain is intrinsically disordered. We found that it acquires a folded conformation when bound to its substrate, being the first report of a complete dsRBD folding upon binding. The free unfolded form shows tendency to adopt folded conformations, and goes through an unfolded bound state prior to the folding event. The significance of these results is discussed by comparison with the behavior of other dsRBDs.     

An evolutionary strategy for all-atom folding of the 60-amino-acid bacterial ribosomal protein l20

We have investigated an evolutionary algorithm for de novo all-atom folding of the bacterial ribosomal protein L20. We report results of two simulations that converge to near-native conformations of this 60-amino-acid, four-helix protein. We observe a steady increase of "native content" in both simulated ensembles and a large number of near-native conformations in their final populations. We argue that these structures represent a significant fraction of the low-energy metastable conformations, which characterize the folding funnel of this protein. These data validate our all-atom free-energy force field PFF01 for tertiary structure prediction of a previously inaccessible structural family of proteins. We also compare folding simulations of the evolutionary algorithm with the basin-hopping technique for the Trp-cage protein. We find that the evolutionary algorithm generates a dynamic memory in the simulated population, which leads to faster overall convergence.     

Ubiquitin folds via a flip-twist-lock mechanism

To perform specific functional activities, the majority of proteins should fold into their distinct three-dimensional conformations. However, the biologically active conformation of a protein is generally found to be marginally stable than the other conformations that the chain can adopt. How a protein finds its native conformation from its post-synthesis unfolded structure in a complex conformational landscape is the unsolved question that still drives the protein folding community. Here, we report the folding mechanism of a globular protein, ubiquitin, from its chemically denatured state using all-atom molecular dynamics simulations. From the kinetic analysis of the simulated trajectories we show that the folding process can be described by the hydrophobic collapse mechanism, initiated by the “dewetting transition”, and subsequently assisted by the origination of an N-terminal folding nucleus, and finally supported by a native salt-bridge interaction between K11 and E34. We show that ubiquitin folds via an intermediate. Finally, we confirm the presence of “biological water” and explain its role to the folding process.