Gatekeepers in the ribosomal protein s6: thermodynamics, kinetics, and folding pathways revealed by a minimalist protein model

We investigate the effect of structural gatekeepers on the folding of the ribosomal protein S6. Folding thermodynamics and early refolding kinetics are studied for this system utilizing computer simulations of a minimalist protein model. When gatekeepers are eliminated, the thermodynamic signature of a folding intermediate emerges, and a marked decrease in folding efficiency is observed. We explain the prerequisites that determine the "strength" of a given gatekeeper. The investigated gatekeepers are found to have distinct functions, and to guide the folding and time-dependent packing of non-overlapping secondary structure elements in the protein. Gatekeepers avoid kinetic traps during folding by favoring the formation of "productive topologies" on the way to the native state. The trends in folding rates in the presence/absence of gatekeepers observed for our minimalist model of S6 are in very good agreement with experimental data on this protein.     

Lattice simulations of aggregation funnels for protein folding.

A computer model of protein aggregation competing with productive folding is proposed. Our model adapts techniques from lattice Monte Carlo studies of protein folding to the problem of aggregation. However, rather than starting with a single string of residues, we allow independently folding strings to undergo collisions and consider their interactions in different orientations. We first present some background into the nature and significance of protein aggregation and the use of lattice Monte Carlo simulations in understanding other aspects of protein folding. The results of a series of simulation experiments involving simple versions of the model illustrate the importance of considering aggregation in simulations of protein folding and provide some preliminary understanding of the characteristics of the model. Finally, we discuss the value of the model in general and of our particular design decisions and experiments. We conclude that computer simulation techniques developed to study protein folding can provide insights into protein aggregation, and that a better understanding of aggregation may in turn provide new insights into and constraints on the more general protein folding problem.     

Conformational detours during folding of a collapsed state

The protein S6 is a useful model to probe the role of partially folded states in the folding process. In the absence of salt, S6 folds from the denatured state D to the native state N without detectable intermediates. High concentrations of sodium sulfate induce the accumulation of a collapsed state C, which is off the direct folding route. However, the mutation VA85 enables S6 to fold from C directly to N through the transition state TS(C). According to the denaturant dependence of this reaction, TS(C) and C are equally compact, but the data are difficult to deconvolute. Therefore, I have measured the heat capacities (DeltaC(p)) for the D-->C and C-->TS(C) transitions. The DeltaC(p)-values suggest that C needs to increase its surface area in order to fold directly to N. This underlines that it is a misfolded state that can only fold by at least partial unfolding. In contrast to the C-state formed by S6 wildtype, the VA85 C-state is just as compact as the native state, and this may be a prerequisite for direct folding. Individual "gatekeeper" residues may thus play a disproportionately large role in guiding proteins through different folding pathways.     

Protein aggregation and amyloidosis: confusion of the kinds?

Recent years have witnessed major advances in our understanding of the structural basis of protein aggregation on several fronts. Firstly, high-resolution structural information that remained elusive for many years was provided by a series of studies of amyloid fibers using NMR, X-ray crystallography and electron microscopy, thereby confirming earlier models based on lower resolution observations. Secondly, studies of the sequence determinants of protein aggregation culminated in the development of computer algorithms that predict aggregation-prone sequences with good accuracy, allowing the design of mutations that reduce aggregation. Thirdly, based on the first results from such predictions and on statistical analysis of naturally occurring aggregating sequences, a picture is emerging in which aggregation-prone sequences are capped by gatekeeper residues that oppose aggregation. In addition to their aggregation-opposing function, it seems that gatekeeper residues are also important in determining chaperone selectivity for strongly aggregating regions. Finally, recent computational and experimental work shows that preventing aggregation does not necessarily mean that amyloid formation is prevented and vice versa. Thus, although aggregation and amyloidosis correlate to a certain extent, they are different processes and should be treated as such.     

Optimal combination of theory and experiment for the characterization of the protein folding landscape of S6: how far can a minimalist model go?

The detailed characterization of the overall free energy landscape associated with the folding process of a protein is the ultimate goal in protein folding studies. Modern experimental techniques provide accurate thermodynamic and kinetic measurements on restricted regions of a protein landscape. Although simplified protein models can access larger regions of the landscape, they are oftentimes built on assumptions and approximations that affect the accuracy of the results. We present a new methodology that allows to combine the complementary strengths of theory and experiment for a more complete characterization of a protein folding landscape. We prove that this new procedure allows a simplified protein model to reproduce remarkably well (correlation coefficient > 0.9) all experimental data available on free energies differences upon single mutations for S6 ribosomal protein and two circular permutants. Our results confirm and quantify the hypothesis, recently formulated on the basis of experimental data, that the folding landscape of protein S6 is strongly affected by an atypical distribution of contact energies.     

Are residues in a protein folding nucleus evolutionarily conserved?

Protein is the working molecule of the cell, and evolution is the hallmark of life. It is important to understand how protein folding and evolution influence each other. Several studies correlating experimental measurement of residue participation in folding nucleus and sequence conservation have reached different conclusions. These studies are based on assessment of sequence conservation at folding nucleus sites using entropy or relative entropy measurement derived from multiple sequence alignment. Here we report analysis of conservation of folding nucleus using an evolutionary model alternative to entropy-based approaches. We employ a continuous time Markov model of codon substitution to distinguish mutation fixed by evolution and mutation fixed by chance. This model takes into account bias in codon frequency, bias-favoring transition over transversion, as well as explicit phylogenetic information. We measure selection pressure using the ratio omega of synonymous versus non-synonymous substitution at individual residue site. The omega-values are estimated using the PAML method, a maximum-likelihood estimator. Our results show that there is little correlation between the extent of kinetic participation in protein folding nucleus as measured by experimental phi-value and selection pressure as measured by omega-value. In addition, two randomization tests failed to show that folding nucleus residues are significantly more conserved than the whole protein, or the median omega value of all residues in the protein. These results suggest that at the level of codon substitution, there is no indication that folding nucleus residues are significantly more conserved than other residues. We further reconstruct candidate ancestral residues of the folding nucleus and suggest possible test tube mutation studies for testing folding behavior of ancient folding nucleus.     

Comparison between long-range interactions and contact order in determining the folding rate of two-state proteins: application of long-range order to folding rate prediction.

The contact order is believed to be an important factor for understanding protein folding mechanisms. In our earlier work, we have shown that the long-range interactions play a vital role in protein folding. In this work, we analyzed the contribution of long-range contacts to determine the folding rate of two-state proteins. We found that the residues that are close in space and are separated by at least ten to 15 residues in sequence are important determinants of folding rates, suggesting the presence of a folding nucleus at an interval of approximately 25 residues. A novel parameter "long-range order" has been proposed to predict protein folding rates. This parameter shows as good a relationship with the folding rate of two-state proteins as contact order. Further, we examined the minimum limit of residue separation to determine the long-range contacts for different structural classes. We observed an excellent correlation between long-range order and folding rate for all classes of globular proteins. We suggest that in mixed-class proteins, a larger number of residues can serve as folding nuclei compared to all-alpha and all-beta proteins. A simple statistical method has been developed to predict the folding rates of two-state proteins using the long-range order that produces an agreement with experimental results that is better or comparable to other methods in the literature.     

In vivo aggregation of bovine beta-lactoglobulin is affected by Cys at position 121

Bovine beta-lactoglobulin (BLG) has been widely used as a model system to study protein folding and aggregation and for biotechnology applications. Native BLG contains two disulfide bonds and one free cysteine at position 121. This free thiol group has been shown to be responsible for the irreversibility of BLG denaturation in vitro, but nothing is known about its relevance during protein folding inside the cell. Here, we report the expression of soluble wild type recombinant BGL in Escherichia coli cells at about 109 mg rBLG/g wet weight cells and a comparison between the aggregation of wt BLG and its variant C121S upon intracellular expression. We show that in E. coli C121SBLG is more prone to aggregation than the wild type protein and that their different behavior depends on the oxidation of disulfide bonds. Our results underline the key contribution of the unpaired cysteine residue during the oxidative folding pathway and indicate BLG as a useful tool for the study of protein aggregation in vivo.     

Transition states for folding of circular-permuted proteins

To explore the role of entropy and chain connectivity in protein folding, a particularly interesting scheme, namely, the circular permutation, has been used. Recently, experimental observations showed that there are large differences in the folding mechanisms between the wild-type proteins and their circular permutants. These differences are strongly related to the change in the intrachain connectivity. Some results obtained by molecular dynamics simulations also showed a good agreement with the experimental findings. Here, we use a topology-based free-energy functional method to study the role of the chain connectivity in folding by comparing features of transition states of the wild-type proteins with those of their circular permutants. We concentrate our study on 3 small globular proteins, namely, the alpha-spectrin SH3 domain (SH3), the chymotrypsin inhibitor 2 (CI2), and the ribosomal protein S6, and obtain exciting results that are consistent with the available experimental and simulation results. A heterogeneity of the interaction energies between contacts for protein CI2 and for protein S6 is also introduced, which characterizes the strong interactions between contacts with long loops, as speculated from experiments for protein S6. The comparison between the folding nucleus of the wild-type proteins and those of their circular permutants indicates that chain connectivity affects remarkably the shapes of the energy profiles and thus the folding mechanism. Further comparisons between our theoretical calculated phi(th) values and the experimental observed phi(exp) values for the 3 proteins and their permutants show that our results are in good agreement with experimental ones and that correlations between them are high. These indicate that the free-energy functional method really provides a way to analyze the folding behavior of the circular-permuted proteins and therefore the folding mechanism of the wild-type proteins.     

Folding of Cu, Zn superoxide dismutase and familial amyotrophic lateral sclerosis

Cu, Zn superoxide dismutase (SOD1) has been implicated in the familial form of the neurodegenerative disease amyotrophic lateral sclerosis (ALS). It has been suggested that mutant mediated SOD1 misfolding/aggregation is an integral part of the pathology of ALS. We study the folding thermodynamics and kinetics of SOD1 using a hybrid molecular dynamics approach. We reproduce the experimentally observed SOD1 folding thermodynamics and find that the residues which contribute the most to SOD1 thermal stability are also crucial for apparent two-state folding kinetics. Surprisingly, we find that these residues are located on the surface of the protein and not in the hydrophobic core. Mutations in some of the identified residues are found in patients with the disease. We argue that the identified residues may play an important role in aggregation. To further characterize the folding of SOD1, we study the role of cysteine residues in folding and find that non-native disulfide bond formation may significantly alter SOD1 folding dynamics and aggregation propensity.     

Effect of denaturant and protein concentrations upon protein refolding and aggregation: a simple lattice model.

We present a study of the competition between protein refolding and aggregation for simple lattice model proteins. The effect of solvent conditions (i.e., the denaturant concentration and the protein concentration) on the folding and aggregation behavior of a system of simple, two-dimensional lattice protein molecules has been investigated via (dynamic Monte Carlo simulations. The population profiles and aggregation propensities of the nine most populated intermediate configurations exhibit a complex dependence on the solution conditions that can be understood by considering the competition between intra- and interchain interactions. Some of these configurations are not even seen in isolated chain simulations; they are observed to be highly aggregation prone and are stabilized primarily by the aggregation reaction in multiple-chain systems. Aggregation arises from the association of partially folded intermediates rather than from the association of denatured random-coil states. The aggregation reaction dominates over the folding reaction at high protein concentration and low denaturant concentration, resulting in low refolding yields at those conditions. However, optimum folding conditions exist at which the refolding yield is a maximum, in agreement with some experimental observations.     

Transition-states in protein folding kinetics: the structural interpretation of Phi values

Phi values are experimental measures of the effects of mutations on the folding kinetics of a protein. A central question is what structural information Phi values give about the transition-state of folding. Traditionally, a Phi value is interpreted as representing the "nativeness" of a mutated residue in the transition-state. However, this interpretation is often problematic. We present here a better structural interpretation of Phi values for mutations within a given helix. Our interpretation is based on a simple physical model that distinguishes between secondary and tertiary free energy contributions of helical residues. From a linear fit of the model to experimental data, we obtain two structural parameters: the extent of helix formation in the transition-state, and the nativeness of tertiary interactions in the transition-state. We apply the model to all proteins with well-characterized helices for which more than 10 Phi values are available: protein A, CI2, and protein L. The model is simple to apply to experimental data, captures nonclassical Phi values <0 or >1 in these helices, and explains how different mutations at a given site can lead to different Phi values.     

Folding behavior of chaperonin-mediated substrate protein

Chaperonin-mediated protein folding is complex. There have been diverse results on folding behavior, and the chaperonin molecules have been investigated as enhancing or retarding the folding rate. To understand the diversity of chaperonin-mediated protein folding, we report a study based on simulations using a simplified Gō-type model. By considering effects of affinity between the substrate protein and the chaperonin wall and spatial confinement of the chaperonin cavity, we study the thermodynamics and kinetics of folding of an unfrustrated substrate protein encapsulated in a chaperonin cavity. The affinity makes the hydrophobic residues of the protein bind to the chaperonin wall, and a strong (or weak) affinity results in a large (or small) effect of binding. Compared with the folding in bulk, the folding in chaperonin cavity with different strengths of affinity shows two kinds of behaviors: one with less dependence on the affinity but more reliance on the spatial confinement effect and the other relying strongly on the affinity. It is found that the enhancement or retardation of the folding rate depends on the competition between the spatial confinement and the affinity due to the chaperonin cavity, and a strong affinity produces a slow folding while a weak affinity induces a fast folding. The crossover between two kinds of folding behaviors happens in the case that the favorable effect of confinement is balanced by the unfavorable effect of the affinity, and a critical affinity strength is roughly defined. By analyzing the contacts formed between the residues of the protein and the chaperonin wall and between the residues of the protein themselves, the role of the affinity in the folding processes is studied. The binding of the residues with the chaperonin wall reduces the formation of both native contacts and nonnative contact or mis-contacts, providing a loose structure for further folding after allosteric change of the chaperonin cavity. In addition, 15 single-site-mutated mutants are simulated in order to test the validity of our model and to investigate the importance of affinity. Inspiringly, our results of the folding rates have a good correlation with those obtained from experiments. The folding rates are inversely correlated with the strength of the binding interactions, i.e., the weaker the binding, the faster the folding. We also find that the inner hydrophobic residues have larger effects on the folding kinetics than those of the exterior hydrophobic residues. We suggest that, besides the confinement effect, the affinity acts as another important factor to affect the folding of the substrate proteins in chaperonin systems, providing an understanding of the folding mechanism of the molecular chaperonin systems.     

Modulation of S6 fibrillation by unfolding rates and gatekeeper residues

We present a protein engineering analysis of the fibrillation of a protein from a thermophilic organism, the 101 residue S6 from Thermus thermophilus. When agitated, S6 fibrillates at pH 2.0 in 0.4 M NaCl. Under these solvent conditions, S6 has native-like secondary structure and also unfolds and refolds cooperatively. However, its tertiary structure appears to be more plastic than at neutral pH, and some regions of the protein may be partially unstructured. At 42 degrees C, there is a lag phase of several days after which fibrillation takes place over several hours. Data from the fibrillation behaviour of a comprehensive series of single and double mutants of S6 suggests that several factors control the onset of fibrillation. Firstly, there appears to be a contiguous region of "gatekeeper" residues that inhibit fibrillation, since their truncation significantly reduces the duration of the lag phase. This region overlaps extensively with the partially unstructured region of the protein, suggesting that residues with enhanced flexibility and solvent-accessibility are important for the initiation of fibrillation. Secondly, longer lag phases correlate with faster rates of unfolding. We interpret this to mean that kinetic stability also controls fibrillation but in the sense that the quasi-native state, rather than the denatured state, is the species that participates in nucleation. This implies that fibrillation can also occur from a quasi-native state as opposed to an ensemble of highly fluctuating structures, and highlights the delicate balance between flexibility and structure required to form organized assemblies of polypeptide chains.     

Accounting for protein-solvent contacts facilitates design of nonaggregating lattice proteins

The folding specificity of proteins can be simulated using simplified structural models and knowledge-based pair-potentials. However, when the same models are used to simulate systems that contain many proteins, large aggregates tend to form. In other words, these models cannot account for the fact that folded, globular proteins are soluble. Here we show that knowledge-based pair-potentials, which include explicitly calculated energy terms between the solvent and each amino acid, enable the simulation of proteins that are much less aggregation-prone in the folded state. Our analysis clarifies why including a solvent term improves the foldability. The aggregation for potentials without water is due to the unrealistically attractive interactions between polar residues, causing artificial clustering. When a water-based potential is used instead, polar residues prefer to interact with water; this leads to designed protein surfaces rich in polar residues and well-defined hydrophobic cores, as observed in real protein structures. We developed a simple knowledge-based method to calculate interactions between the solvent and amino acids. The method provides a starting point for modeling the folding and aggregation of soluble proteins. Analysis of our simple model suggests that inclusion of these solvent terms may also improve off-lattice potentials for protein simulation, design, and structure prediction.     

Derivation of a solubility condition for proteins from an analysis of the competition between folding and aggregation

Failure in maintaining protein solubility in vivo impairs protein homeostasis and results in protein misfolding and aggregation, which are often associated with severe neurodegenerative and systemic disorders that include Alzheimer's and Parkinson's diseases and type II diabetes. In this work we formulate a model of the competition between folding and aggregation, and derive a condition on the solubility of proteins in terms of the stability of their folded states, their aggregation propensities and their degradation rates. From our model, the bistability between folding and aggregation emerges as an intrinsic aspect of protein homeostasis. The analysis of the conditions that determine such a bistability provides a rationalization of the recently observed relationship between the cellular abundance and the aggregation propensity of proteins. We then discuss how the solubility condition that we derive can help rationalise the correlation that has been reported between evolutionary rates and expression levels or proteins, as well as in vivo protein solubility and expression level measurements, and recently elucidated trends of proteome evolution.     

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.     

Simulating the minimum core for hydrophobic collapse in globular proteins.

To investigate the nature of hydrophobic collapse considered to be the driving force in protein folding, we have simulated aqueous solutions of two model hydrophobic solutes, methane and isobutylene. Using a novel methodology for determining contacts, we can precisely follow hydrophobic aggregation as it proceeds through three stages: dispersed, transition, and collapsed. Theoretical modeling of the cluster formation observed by simulation indicates that this aggregation is cooperative and that the simulations favor the formation of a single cluster midway through the transition stage. This defines a minimum solute hydrophobic core volume. We compare this with protein hydrophobic core volumes determined from solved crystal structures. Our analysis shows that the solute core volume roughly estimates the minimum core size required for independent hydrophobic stabilization of a protein and defines a limiting concentration of nonpolar residues that can cause hydrophobic collapse. These results suggest that the physical forces driving aggregation of hydrophobic molecules in water is indeed responsible for protein folding.     

Conserved signature proposed for folding in the lipocalin superfamily

We systematically identify a group of evolutionarily conserved residues proposed for folding in a model beta-barrel superfamily, the lipocalins. The nature of conservation at the structural level is defined and we show that the conserved residues are involved in a network of interactions that form the core of the fold. Exploratory kinetic studies are conducted with a model superfamily member, human serum retinol-binding protein, to examine their role. The present results, coupled with key experimental studies conducted with another lipocalin beta-lactoglobulin, suggest that the evolutionarily conserved regions fold on a faster folding time-scale than the non-conserved regions.