Structural basis for the self-chaperoning function of an RNA collapsed state

Prior to folding to a native functional structure, many large RNAs form conformationally collapsed states. Formation of the near-native collapsed state for the bI5 group I intron RNA plays an obligatory role in self-chaperoning assembly with its CBP2 protein cofactor by preventing formation of stable, misassembled complexes. We show that the collapsed state is essential because CBP2 assembles indiscriminately with the bI5 RNA in any folding state to form long-lived complexes. The most stable protein interaction site in the expanded state-CBP2 complex overlaps, but is not identical to, the native site. Folding to the collapsed state circumvents two distinct misassembly events: inhibitory binding by multiple equivalents of CBP2 and formation of bridged complexes in which CBP2 straddles cognate and noncognate RNAs. Strikingly, protein-bound sites in the expanded state RNA complex are almost the inverse of native RNA-RNA and RNA-protein interactions, indicating that folding to the collapsed state significantly reduces the fraction of RNA surfaces accessible for misassembly. The self-chaperoning function for the bI5 collapsed state is likely to be conserved in other ribonucleoproteins where a protein cofactor binds tightly at a simple RNA substructure or has an RNA binding surface composed of multiple functional sites.     

Folding pathway dependence on energetic frustration and interaction heterogeneity for a three-dimensional hydrophobic protein model

Monte Carlo simulations of a hydrophobic protein model of 40 monomers in the cubic lattice are used to explore the effect of energetic frustration and interaction heterogeneity on its folding pathway. The folding pathway is described by the dependence of relevant conformational averages on an appropriate reaction coordinate, pfold, defined as the probability for a given conformation to reach the native structure before unfolding. We compare the energetically frustrated and heterogeneous hydrophobic potential, according to which individual monomers have a higher or lower tendency to form contacts unspecifically depending on their hydrophobicities, to an unfrustrated homogeneous Go-type potential with uniformly attractive native interactions and neutral non-native interactions (called Go1 in this study), and to an unfrustrated heterogeneous potential with neutral non-native interactions and native interactions having the same energy as the hydrophobic potential (called Go2 in this study). Folding kinetics are slowed down dramatically when energetic frustration increases, as expected and previously observed in a two-dimensional model. Contrary to our previous results in two dimensions, however, it appears that the folding pathway and transition state ensemble can be significantly dependent on the energy function used to stabilize the native structure. The sequence of events along the reaction coordinate, or the order along this coordinate in which different regions of the native conformation become structured, turns out to be similar for the hydrophobic and Go2 potentials, but with analogous events tending to occur at lower pfold values in the first case. In particular, the transition state obtained from the ensemble around pfold = 0.5 is more structured for the hydrophobic potential. For Go1, not only the transition state ensemble but the order of events itself is modified, suggesting that interaction heterogeneity, in addition to energetic frustration, can have significant effects on the folding mechanism, most likely by modifying the probability of different contacts in the unfolded state, the starting point for the folding reaction. Although based on a simple model, these results provide interesting insight into how sequence-dependent switching between folding pathways might occur in real proteins.     

Topological and energetic factors: what determines the structural details of the transition state ensemble and "en-route" intermediates for protein folding? An investigation for small globular proteins

Recent experimental results suggest that the native fold, or topology, plays a primary role in determining the structure of the transition state ensemble, at least for small, fast-folding proteins. To investigate the extent of the topological control of the folding process, we studied the folding of simplified models of five small globular proteins constructed using a Go-like potential to retain the information about the native structures but drastically reduce the energetic frustration and energetic heterogeneity among residue-residue native interactions. By comparing the structure of the transition state ensemble (experimentally determined by Phi-values) and of the intermediates with those obtained using our models, we show that these energetically unfrustrated models can reproduce the global experimentally known features of the transition state ensembles and "en-route" intermediates, at least for the analyzed proteins. This result clearly indicates that, as long as the protein sequence is sufficiently minimally frustrated, topology plays a central role in determining the folding mechanism.     

Unfolding of the cold shock protein studied with biased molecular dynamics

The cold shock protein from Bacillus caldolyticus is a small beta-barrel protein that folds in a two-state mechanism. For the native protein and for several mutants, a wealth of experimental data are available on stability and folding, so that it is an optimal system to study this process. We compare data from unfolding simulations (trajectories of 5 and up to 12 ns) obtained with a bias potential at room temperature and from unbiased thermal unfolding simulations with experimental data. The unfolding patterns derived from the trajectories starting from different native-like conformations and subject to different unfolding conditions agree. The transition state found in the simulations of unfolding is close to the native structure in agreement with experiment. Moreover, a lower value of the free energy barrier of unfolding was found for the mutant R3E than for the mutant E46A and the native protein, as indicated by experimental data. The first unfolding event involves the three-stranded beta-sheet whose decomposition corresponds to the transition state. In contrast to conclusions drawn from experiments, we found that the two-stranded beta-strand forms the most stable substructure, which decomposes very late in the unfolding process. However, assuming that this structure forms very early in the folding process, our findings would not contradict the experiments but require a different interpretation of them.     

Kinetics, thermodynamics and evolution of non-native interactions in a protein folding nucleus

A lattice model with side chains was used to investigate protein folding with computer simulations. In this model, we rigorously demonstrate the existence of a specific folding nucleus. This nucleus contains specific interactions not present in the native state that, when weakened, slow folding but do not change protein stability. Such a decoupling of folding kinetics from thermodynamics has been observed experimentally for real proteins. From our results, we conclude that specific non-native interactions in the transition state would give rise to straight phi-values that are negative or larger than unity. Furthermore, we demonstrate that residue Ile 34 in src SH3, which has been shown to be kinetically, but not thermodynamically, important, is universally conserved in proteins with the SH3 fold. This is a clear example of evolution optimizing the folding rate of a protein independent of its stability and function.     

A new protein folding recognition potential function

On the study of protein inverse folding problem, one goal is to find simple and efficient potential to evaluate the compatibility between structure and a given sequence. We present here a novo empirical mean force potential to address the importance of electrostatic interactions in protein inverse folding study. It is based on protein main chain polar fraction and constructed in a way similar with Sippl's from a database of 64 known independent three-dimensional protein structures. This potential was applied to recognize the protein native conformations among a conformation pool. Calculated results show that this potential is powerful in picking out native conformations, in addition it can also find structure similarity between proteins with low sequence similarity. The success of this new potential clearly shows the importance of electrostatic factors in protein inverse folding studies. © 1995 Wiley-Liss, Inc.     

Discrimination of the native from misfolded protein models with an energy function including implicit solvation.

An essential requirement for theoretical protein structure prediction is an energy function that can discriminate the native from non-native protein conformations. To date most of the energy functions used for this purpose have been extracted from a statistical analysis of the protein structure database, without explicit reference to the physical interactions responsible for protein stability. The use of the statistical functions has been supported by the widespread belief that they are superior for such discrimination to physics-based energy functions. An effective energy function which combined the CHARMM vacuum potential with a Gaussian model for the solvation free energy is tested for its ability to discriminate the native structure of a protein from misfolded conformations; the results are compared with those obtained with the vacuum CHARMM potential. The test is performed on several sets of misfolded structures prepared by others, including sets of about 650 good decoys for six proteins, as well as on misfolded structures of chymotrypsin inhibitor 2. The vacuum CHARMM potential is successful in most cases when energy minimized conformations are considered, but fails when applied to structures relaxed by molecular dynamics. With the effective energy function the native state is always more stable than grossly misfolded conformations both in energy minimized and molecular dynamics-relaxed structures. The present results suggest that molecular mechanics (physics-based) energy functions, complemented by a simple model for the solvation free energy, should be tested for use in the inverse folding problem, and supports their use in studies of the effective energy surface of proteins in solution. Moreover, the study suggests that the belief in the superiority of statistical functions for these purposes may be ill founded.     

Analyzing the effect of homogeneous frustration in protein folding

The energy landscape theory has been an invaluable theoretical framework in the understanding of biological processes such as protein folding, oligomerization, and functional transitions. According to the theory, the energy landscape of protein folding is funneled toward the native state, a conformational state that is consistent with the principle of minimal frustration. It has been accepted that real proteins are selected through natural evolution, satisfying the minimum frustration criterion. However, there is evidence that a low degree of frustration accelerates folding. We examined the interplay between topological and energetic protein frustration. We employed a C_α structure‐based model for simulations with a controlled nonspecific energetic frustration added to the potential energy function. Thermodynamics and kinetics of a group of 19 proteins are completely characterized as a function of increasing level of energetic frustration. We observed two well‐separated groups of proteins: one group where a little frustration enhances folding rates to an optimal value and another where any energetic frustration slows down folding. Protein energetic frustration regimes and their mechanisms are explained by the role of non‐native contact interactions in different folding scenarios. These findings strongly correlate with the protein free‐energy folding barrier and the absolute contact order parameters. These computational results are corroborated by principal component analysis and partial least square techniques. One simple theoretical model is proposed as a useful tool for experimentalists to predict the limits of improvements in real proteins.Proteins 2013; 81:1727–1737. © 2013 Wiley Periodicals, Inc.     

Exploring the Levinthal limit in protein folding

According to the thermodynamic hypothesis, the native state of proteins is uniquely defined by their amino acid sequence. On the other hand, according to Levinthal, the native state is just a local minimum of the free energy and a given amino acid sequence, in the same thermodynamic conditions, can assume many, very different structures that are as thermodynamically stable as the native state. This is the Levinthal limit explored in this work. Using computer simulations, we compare the interactions that stabilize the native state of four different proteins with those that stabilize three non-native states of each protein and find that the nature of the interactions is very similar for all such 16 conformers. Furthermore, an enhancement of the degree of fluctuation of the non-native conformers can be explained by an insufficient relaxation to their local free energy minimum. These results favor Levinthal’s hypothesis that protein folding is a kinetic non-equilibrium process.     

Selective stabilization of a partially unfolded protein by a metabolite

When proteins fold in vivo, the intermediates that exist transiently on their folding pathways are exposed to the potential interactions with a plethora of metabolites within the cell. However, these potential interactions are commonly ignored. Here, we report a case in which a ubiquitous metabolite interacts selectively with a nonnative conformation of a protein and facilitates protein folding and unfolding process. From our previous proteomics study, we have discovered that Escherichia coli glyceraldehyde-3-phosphate dehydrogenase (GAPDH), which is not known to bind ATP under native conditions, is apparently destabilized in the presence of a physiological concentration of ATP. To decipher the origin of this surprising effect, we investigated the thermodynamics and kinetics of folding and unfolding of GAPDH in the presence of ATP. Equilibrium unfolding of the protein in urea showed that a partially unfolded equilibrium intermediate accumulates in the presence of ATP. This intermediate has a quaternary structure distinct from the native protein. Also, ATP significantly accelerates the unfolding of GAPDH by selectively stabilizing a transition state that is distinct from the native state of the protein. Moreover, ATP also significantly accelerates the folding of GAPDH. These results demonstrate that ATP interacts specifically with a partially unfolded form of GAPDH and affects the kinetics of folding and unfolding of this protein. This unusual effect of ATP on the folding of GAPDH implies that endogenous metabolites may facilitate protein folding in vivo by interacting with partially unfolded intermediates.     

Statistical mechanics of protein folding by cluster distance geometry

This is our second type of model for protein folding where the configurational parameters and the effective potential energy function are chosen in such a way that all conformations are described and the canonical partition function can be evaluated analytically. Structure is described in terms of distances between pairs of sequentially contiguous blocks of eight residues, and all possible conformations are grouped into 71 subsets in terms of bounds on these distances. The energy is taken to be a sum of pairwise interactions between such blocks. The 210 energy parameters were adjusted so that the native folds of 32 small proteins are favored in free energy over the denatured state. We then found 146 proteins having negligible sequence similarity to any of the training proteins, yet the free energy of the respective correct native states were favored over the denatured state.     

Improved Gō-like models demonstrate the robustness of protein folding mechanisms towards non-native interactions

The use of simple theoretical models has provided a considerable contribution to our present understanding of the means by which proteins adopt their native fold from the plethora of available unfolded states. A common assumption in building computationally tractable models has been the neglect of stabilizing non-native interactions in the class of models described as "Gō-like." The focus of this study is the characterization of the folding of a number of proteins via a Gō-like model, which aims to map a maximal amount of information reflecting the protein sequence onto a "minimalist" skeleton. This model is shown to contain sufficient information to reproduce the folding transition states of a number of proteins, including topologically analogous proteins that fold via different transition states. Remarkably, these models also demonstrate consistency with the general features of folding transition states thought to be stabilized by non-native interactions. This suggests that native interactions are the primary determinant of most protein folding transition states, and that non-native interactions lead only to local structural perturbations. A prediction is also included for an asymmetrical folding transition state of bacteriophage lambda protein W, which has yet to be subjected to experimental characterization.     

Amphiphilic ��-Helical Potential: A Putative Folding Motif Adding Few Constraints to Protein Evolution

Evidence from a number of studies indicates that protein folding is dictated not only by factors stabilizing the native state, but also by potentially independent factors that create folding pathways. How natural selection might cope simultaneously with two independent factors was addressed in this study within the framework of the "Lim-model" of protein folding, which postulates that the early stages of folding of all globular proteins, regardless of their native structure, are directed at least in part by potential to form amphiphilic α-helices. For this purpose, the amphiphilic α-helical potential in randomly ordered amino acid sequences and the conservation in phylogeny of amphiphilic α-helical potential within various proteins were assessed. These analyses revealed that amphiphilic α-helical potential is a common occurrence in random sequences, and that the presence of amphiphilic α-helical potential is present but not conserved in phylogeny within a given protein. The results suggest that the rapid formation of molten globules and the variable behavior of those globules depending on the protein may be a fundamental property of polymers of naturally occurring amino acids more so than a trait that must be derived or maintained by natural selection. Further, the results point toward the utility of randomly occurring process in protein function and evolution, and suggest that the formation of efficient pathways that determine early processes in protein folding, unlike the formation of stable, native protein structure, does not present a substantial hurdle during the evolution of amino acid sequences.     

Mutational analysis of acylphosphatase suggests the importance of topology and contact order in protein folding.

Muscle acylphosphatase (AcP) is a small protein that folds very slowly with two-state behavior. The conformational stability and the rates of folding and unfolding have been determined for a number of mutants of AcP in order to characterize the structure of the folding transition state. The results show that the transition state is an expanded version of the native protein, where most of the native interactions are partially established. The transition state of AcP turns out to be remarkably similar in structure to that of the activation domain of procarboxypeptidase A2 (ADA2h), a protein having the same overall topology but sharing only 13% sequence identity with AcP. This suggests that transition states are conserved between proteins with the same native fold. Comparison of the rates of folding of AcP and four other proteins with the same topology, including ADA2h, supports the concept that the average distance in sequence between interacting residues (that is, the contact order) is an important determinant of the rate of protein folding.     

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.     

The importance of explicit chain representation in protein folding models: an examination of Ising-like models

A class of models that represents a protein chain as a sequence of "folded" and "unfolded" residues has recently been used to correlate rates and mechanisms of protein folding with the protein native structure. In order to better understand the conditions under which these "Ising-like" models apply, we compare results from this model to those obtained from an off-lattice model which uses the same potential function. We find that Ising-like models by construction impose folding via a highly sequential nucleation-condensation mechanism, which in turn leads to more rugged energy landscapes, fewer "pathways" to the native state, and in the specific case examined here, the cold shock protein A from Escherichia coli, a qualitative difference in the most likely order of events in folding.     

Intrinsically disordered protein from a pathogenic mesophile Mycobacterium tuberculosis adopts structured conformation at high temperature

Compared to eukaryotes, the occurrence of "intrinsically disordered" or "natively unfolded" proteins in prokaryotes has not been explored extensively. Here, we report the occurrence of an intrinsically disordered protein from the mesophilic human pathogen Mycobacterium tuberculosis. The Histidine-tagged recombinant Rv3221c biotin-binding protein is intrinsically disordered at ambient and physiological growth temperatures as revealed by circular dichroism and Fourier transform infrared (FTIR) spectroscopic studies. However, an increase in temperature induces a transition from disordered to structured state with a folding temperature of approximately 53 degrees C. Addition of a structure inducing solvent trifluoroethanol (TFE) causes the protein to fold at lower temperatures suggesting that TFE fosters hydrophobic interactions, which drives protein folding. Differential Scanning Calorimetry studies revealed that folding is endothermic and the transition from a disordered to structured state is continuous (higher-order), implying existence of intermediates during folding process. Secondary structure analysis revealed that the protein has propensity to form beta-sheets. This is in conformity with FTIR spectrum that showed an absorption peak at wave number of 1636 cm(-1), indicative of disordered beta-sheet conformation in the native state. These data suggest that although Rv3221c may be disordered under ambient or optimal growth temperature conditions, it has the potential to fold into ordered structure at high temperature driven by increased hydrophobic interactions. In contrast to the generally known behavior of other intrinsically disordered proteins folding at high temperature, Rv3221c does not appear to oligomerize or aggregate as revealed through numerous experiments including Congo red binding, Thioflavin T-binding, turbidity measurements, and examining molar ellipticity as a function of protein concentration. The amino acid composition of Rv3221c reveals that it has 24% charged and 54.9% hydrophobic amino acid residues. In this respect, this protein, although belonging to the class of intrinsically disordered proteins, has distinct features. The intrinsically disordered state and the biotin-binding feature of this protein suggest that it may participate in many biochemical processes requiring biotin as a cofactor and adopt suitable conformations upon binding other folded targets.     

Three topologically equivalent core residues affect the transition state ensemble in a protein folding reaction

The single domain protein, interleukin-1beta, is representative of a distinct class of proteins characterized by their beta-trefoil topology. Each subdomain of this structural class is composed of a beta beta beta loop beta (betabetabetaLbeta) motif comprised of approximately 50 residues and gives the protein a pseudo- 3-fold axis of symmetry. A common feature of proteins in this topological family appears to be that they are slow folders, which reach the native state on the order of tens to 100s of seconds. Sequence analysis of interleukin-1beta indicates that three phenylalanine residues located at positions 42, 101, and 146 are well conserved, separated by approximately 50 residues in the primary sequence, located in similar positions in the pseudo-symmetric units of the trefoil, and are juxtaposed to one another in conformational space. These residues surround the hydrophobic cavity and "pin" the hairpin triplet cap to the core beta-barrel. To determine if cap-barrel interactions are involved in maintaining the structural stability and cooperativity or in controlling the slow formation of the native state, we performed a series of mutational studies. The results indicate that interleukin-1beta tolerates large increases in side-chain volume at these three topologically conserved sites with little effect on stability, while the kinetics show significant differences in both the unfolding and refolding rates. Taken together, our results indicate that these conserved core residues are essential contacts in the transition-state ensemble for folding.     

Modulating protein folding rates in vivo and in vitro by side-chain interactions between the parallel beta strands of green fluorescent protein

We have identified pairs of residues across the two parallel beta strands of green fluorescent protein that facilitate native strand register of the surface-exposed beta barrel. After constructing a suitable host environment around two guest residues, minimizing interactions of the guest residues with surrounding side-chains yet maintaining the wild-type protein structure and the chromophore environment, we introduced a library of cross-strand pairings by cassette mutagenesis. Colonies of Escherichia coli transformed with the library differ in intracellular fluorescence. Most of the fluorescent pairs have predominantly charged and polar guest site residues. The magnitude and the rate of fluorescence acquisition in vivo from transformed E. coli cells varies among the mutants despite comparable levels of protein expression. Spectroscopic measurements of purified mutants show that the native protein structure is maintained. Kinetic studies using purified protein with fully matured chromophores demonstrate that the mutants span a 10-fold range in folding rates with undetectable differences in unfolding rates. Thus, green fluorescent protein provides an ideal system for monitoring determinants of in vivo protein folding. Cross-strand pairings affect both protein stability and folding kinetics by favoring the formation of native strand register preferentially to non-native strand alignments.