A simple model of chaperonin-mediated protein folding

Protein chain collapse modulation and folding stimulation by GroEL-ES

The collapse of polypeptides is thought important to protein folding, aggregation, intrinsic disorder, and phase separation. However, whether polypeptide collapse is modulated in cells to control protein states is unclear. Here, using integrated protein manipulation and imaging, we show that the chaperonin GroEL-ES can accelerate the folding of proteins by strengthening their collapse. GroEL induces contractile forces in substrate chains, which draws them into the cavity and triggers a general compaction and discrete folding transitions, even for slow-folding proteins. This collapse enhancement is strongest in the nucleotide-bound states of GroEL and is aided by GroES binding to the cavity rim and by the amphiphilic C-terminal tails at the cavity bottom. Collapse modulation is distinct from other proposed GroEL-ES folding acceleration mechanisms, including steric confinement and misfold unfolding. Given the prevalence of collapse throughout the proteome, we conjecture that collapse modulation is more generally relevant within the protein quality control machinery.

The Protein Folding Problem: The Role of Theory

The protein folding problem was first articulated as question of how order arose from disorder in proteins: How did the various native structures of proteins arise from interatomic driving forces encoded within their amino acid sequences, and how did they fold so fast? These matters have now been largely resolved by theory and statistical mechanics combined with experiments. There are general principles. Chain randomness is overcome by solvation-based codes. And in the needle-in-a-haystack metaphor, native states are found efficiently because protein haystacks (conformational ensembles) are funnel-shaped. Order-disorder theory has now grown to encompass a large swath of protein physical science across biology.

Understanding the role of hydrophobic patches in protein disaggregation

Protein folding is a very complex process and, so far, the mechanism of folding still intrigues the research community. Despite a large conformational space available (O(10⁴⁷) for a 100 amino acid residue), most proteins fold into their native state within a very short time. While small proteins fold relatively fast (a few microseconds) large globular proteins may take as long as several milliseconds to fold. During the folding process, the protein synthesized in the ribosome is exposed to the crowded environment of the cell and is easily prone to misfolding and aggregation due to interactions with other proteins or biomacromolecules present within the cell. These large proteins, therefore, rely on chaperones for their folding and repair. Chaperones are known to have hydrophobic patchy domains that play a crucial role in shielding the protein against misfolding and disaggregation of aggregated proteins. In the current article, Monte Carlo simulations carried out in the framework of the hydrophobic-polar (H-P) lattice model indicate that hydrophobic patchy domains drastically reduce the inter-protein interactions and are efficient in disaggregating proteins. The effectiveness of the disaggregation depends on the size and distribution of these patches on the surface and also on the strength of the interaction between the protein and the surface. Further, our results indicate that when the patch is complementary to the exposed hydrophobic patch of the protein, protein disaggregation is accompanied by stabilization of the protein even relative to its bulk behavior due to favorable protein-surface interactions. We believe that these findings shed light on the role of the class of chaperones known as heat shock proteins (Hsps) on protein disaggregation and refolding.

Energy landscape remodeling mechanism of Hsp70-chaperone-accelerated protein folding

Hsp70 chaperone is one of the key protein machines responsible for the quality control of protein production in cells. Facilitating in vivo protein folding by counteracting misfolding and aggregation is the essence of its biological function. Although the allosteric cycle during its functional actions has been well characterized both experimentally and computationally, the mechanism by which Hsp70 assists protein folding is still not fully understood. In this work, we studied the Hsp70-mediated folding of model proteins with rugged energy landscape by using molecular simulations. Different from the canonical scenario of Hsp70 functioning, which assumes that folding of substrate proteins occurs spontaneously after releasing from chaperones, our results showed that the substrate protein remains in contacts with the chaperone during its folding process. The direct chaperone-substrate interactions in the open conformation of Hsp70 tend to shield the substrate sites prone to form non-native contacts, which therefore avoids the frustrated folding pathway, leading to a higher folding rate and less probability of misfolding. Our results suggest that in addition to the unfoldase and holdase functions widely addressed in previous studies, Hsp70 can facilitate the folding of its substrate proteins by remodeling the folding energy landscape and directing the folding processes, demonstrating the foldase scenario. These findings add new, to our knowledge, insights into the general molecular mechanisms of chaperone-mediated protein folding.

The concept of protein folding/unfolding and its impacts on human health

Proteins have evolved in specific 3D structures and play different functions in cells and determine various reactions and pathways. The newly synthesized amino acid chains once depart ribosome must crumple into three-dimensional structures so can be biologically active. This process of protein that makes a functional molecule is called protein folding. The protein folding is both a biological and a physicochemical process that depends on the sequence of it. In fact, this process occurs more complicated and in some cases and in exposure to some molecules like glucose (glycation), mistaken folding leads to amyloid structures and fatal disorders called conformational diseases. Such conditions are detected by the quality control system of the cell and these abnormal proteins undergo renovation or degradation. This scenario takes place by the chaperones, chaperonins, and Ubiquitin-proteasome complex. Understanding of protein folding mechanisms from different views including experimental and computational approaches has revealed some intermediate ensembles such as molten globule and has been subjected to biophysical and molecular biology attempts to know more about prevalent conformational diseases.

Hydrophobic confinement modulates thermal stability and assists knotting in the folding of tangled proteins

There is growing support for the idea that the in vivo folding process of knotted proteins is assisted by chaperonins, but the mechanism of chaperonin assisted folding remains elusive. Here, we conduct extensive Monte Carlo simulations of lattice and off-lattice models to explore the effects of confinement and hydrophobic intermolecular interactions with the chaperonin cage in the folding and knotting processes. We find that moderate to high protein-cavity interactions (which are likely to be established in the beginning of the chaperonin working cycle) cause an energetic destabilization of the protein that overcomes the entropic stabilization driven by excluded volume, and leads to a decrease of the melting temperature relative to bulk conditions. Moreover, mild-to-moderate hydrophobic interactions with the cavity (which would be established later in the cycle) lead to a significant enhancement of knotting probability in relation to bulk conditions while simultaneously moderating the effect of steric confinement in the enhancement of thermal stability. Our results thus indicate that the chaperonin may be able to assist knotting without simultaneously thermally stabilizing potential misfolded states to a point that would hamper productive folding thus compromising its functional role.

Functional principles and regulation of molecular chaperones

To be able to perform their biological function, a protein needs to be correctly folded into its three dimensional structure. The protein folding process is spontaneous and does not require the input of energy. However, in the crowded cellular environment where there is high risk of inter-molecular interactions that may lead to protein molecules sticking to each other, hence forming aggregates, protein folding is assisted. Cells have evolved robust machinery called molecular chaperones to deal with the protein folding problem and to maintain proteins in their functional state. Molecular chaperones promote efficient folding of newly synthesized proteins, prevent their aggregation and ensure protein homeostasis in cells. There are different classes of molecular chaperones functioning in a complex interplay. In this review, we discuss the principal characteristics of different classes of molecular chaperones, their structure-function relationships, their mode of regulation and their involvement in human disorders.

Steric confinement and enhanced local flexibility assist knotting in simple models of protein folding

The chaperonin complex GroEL-GroES is able to accelerate the folding process of knotted proteins considerably. However, the folding mechanism inside the chaperonin cage is elusive. Here we use a combination of lattice and off-lattice Monte Carlo simulations of simple Gō models to study the effect of physical confinement and local flexibility on the folding process of protein model systems embedding a trefoil knot in their native structure. This study predicts that steric confinement plays a specific role in the folding of knotted proteins by increasing the knotting probability for very high degrees of confinement. This effect is observed for protein MJ0366 even above the melting temperature for confinement sizes compatible with the size of the GroEL/GroES chaperonin cage. An enhanced local flexibility produces the same qualitative effects on the folding process. In particular, we observe that knotting probability increases up to 40% in the transition state of protein MJ0366 when flexibility is enhanced. This is underlined by a structural change in the transition state, which becomes devoid of helical content. No relation between the knotting mechanism and flexibility was found in the context of the off-lattice model adopted in this work.

Chloroplast Chaperonin: An Intricate Protein Folding Machine for Photosynthesis

Group I chaperonins are large cylindrical-shaped nano-machines that function as a central hub in the protein quality control system in the bacterial cytosol, mitochondria and chloroplasts. In chloroplasts, proteins newly synthesized by chloroplast ribosomes, unfolded by diverse stresses, or translocated from the cytosol run the risk of aberrant folding and aggregation. The chloroplast chaperonin system assists these proteins in folding into their native states. A widely known protein folded by chloroplast chaperonin is the large subunit of ribulose 1,5-bisphosphate carboxylase/oxygenase (Rubisco), an enzyme responsible for the fixation of inorganic CO₂ into organic carbohydrates during photosynthesis. Chloroplast chaperonin was initially identified as a Rubisco-binding protein. All photosynthetic eucaryotes genomes encode multiple chaperonin genes which can be divided into α and β subtypes. Unlike the homo-oligomeric chaperonins from bacteria and mitochondria, chloroplast chaperonins are more complex and exists as intricate hetero-oligomers containing both subtypes. The Group I chaperonin requires proper interaction with a detachable lid-like co-chaperonin in the presence of ATP and Mg²⁺ for substrate encapsulation and conformational transition. Besides the typical Cpn10-like co-chaperonin, a unique co-chaperonin consisting of two tandem Cpn10-like domains joined head-to-tail exists in chloroplasts. Since chloroplasts were proposed as sensors to various environmental stresses, this diversified chloroplast chaperonin system has the potential to adapt to complex conditions by accommodating specific substrates or through regulation at both the transcriptional and post-translational levels. In this review, we discuss recent progress on the unique structure and function of the chloroplast chaperonin system based on model organisms Chlamydomonas reinhardtii and Arabidopsis thaliana. Knowledge of the chloroplast chaperonin system may ultimately lead to successful reconstitution of eukaryotic Rubisco in vitro.

Protein Folding Mediated by Trigger Factor and Hsp70: New Insights from Single-Molecule Approaches

Chaperones assist in protein folding, but what this common phrase means in concrete terms has remained surprisingly poorly understood. We can readily measure chaperone binding to unfolded proteins, but how they bind and affect proteins along folding trajectories has remained obscure. Here we review recent efforts by our labs and others that are beginning to pry into this issue, with a focus on the chaperones trigger factor and Hsp70. Single-molecule methods are central, as they allow the stepwise process of folding to be followed directly. First results have already revealed contrasts with long-standing paradigms: rather than acting only "early" by stabilizing unfolded chain segments, these chaperones can bind and stabilize partially folded structures as they grow to their native state. The findings suggest a fundamental redefinition of the protein folding problem and a more extensive functional repertoire of chaperones than previously assumed.

Decoding Structural Properties of a Partially Unfolded Protein Substrate: En Route to Chaperone Binding

Many proteins comprising of complex topologies require molecular chaperones to achieve their unique three-dimensional folded structure. The E.coli chaperone, GroEL binds with a large number of unfolded and partially folded proteins, to facilitate proper folding and prevent misfolding and aggregation. Although the major structural components of GroEL are well defined, scaffolds of the non-native substrates that determine chaperone-mediated folding have been difficult to recognize. Here we performed all-atomistic and replica-exchange molecular dynamics simulations to dissect non-native ensemble of an obligate GroEL folder, DapA. Thermodynamics analyses of unfolding simulations revealed populated intermediates with distinct structural characteristics. We found that surface exposed hydrophobic patches are significantly increased, primarily contributed from native and non-native β-sheet elements. We validate the structural properties of these conformers using experimental data, including circular dichroism (CD), 1-anilinonaphthalene-8-sulfonic acid (ANS) binding measurements and previously reported hydrogen-deutrium exchange coupled to mass spectrometry (HDX-MS). Further, we constructed network graphs to elucidate long-range intra-protein connectivity of native and intermediate topologies, demonstrating regions that serve as central “hubs”. Overall, our results implicate that genomic variations (or mutations) in the distinct regions of protein structures might disrupt these topological signatures disabling chaperone-mediated folding, leading to formation of aggregates.

Several non-native proteins require molecular chaperones for proper folding. Many unfolded proteins if not folded accurately, become causal factors in various types of misfolding or aggregation induced diseases such as Alzheimer′s, Huntington′s and several other neurodegenerative disorders. However, structural information of non-folded proteins especially chaperone-dependent proteins is difficult to probe experimentally due to their inherent aggregation propensities. In this work, we study DapA protein, which exhibits obligate requirement on GroEL chaperonin machinery for its folding. We use molecular dynamics simulations to reveal populated intermediate structures of DapA in atomic details. The most plausible intermediate was found to be in agreement with recently reported hydrogen-exchange experimental data. Significant increase in surface exposed hydrophobicity was observed in intermediates compared to native, which was further validated using ANS binding experiments. We also constructed network model of these intermediates that provides remarkable insights into stable hubs (or important residues) underlying diverse states of unfolded proteins. In summary, our work provides a molecular picture of an unfolded protein that is en route to chaperone binding, and these underlying structural properties might act as a molecular signal for their productive folding.

Sequence determines degree of knottedness in a coarse-grained protein model

Knots are abundant in globular homopolymers but rare in globular proteins. To shed new light on this long-standing conundrum, we study the influence of sequence on the formation of knots in proteins under native conditions within the framework of the hydrophobic-polar lattice protein model. By employing large-scale Wang-Landau simulations combined with suitable Monte Carlo trial moves we show that even though knots are still abundant on average, sequence introduces large variability in the degree of self-entanglements. Moreover, we are able to design sequences which are either almost always or almost never knotted. Our findings serve as proof of concept that the introduction of just one additional degree of freedom per monomer (in our case sequence) facilitates evolution towards a protein universe in which knots are rare.

On the characterization and software implementation of general protein lattice models

models of proteins have been widely used as a practical means to computationally investigate general properties of the system. In lattice models any sterically feasible conformation is represented as a self-avoiding walk on a lattice, and residue types are limited in number. So far, only two- or three-dimensional lattices have been used. The inspection of the neighborhood of alpha carbons in the core of real proteins reveals that also lattices with higher coordination numbers, possibly in higher dimensional spaces, can be adopted. In this paper, a new general parametric lattice model for simplified protein conformations is proposed and investigated. It is shown how the supporting software can be consistently designed to let algorithms that operate on protein structures be implemented in a lattice-agnostic way. The necessary theoretical foundations are developed and organically presented, pinpointing the role of the concept of main directions in lattice-agnostic model handling. Subsequently, the model features across dimensions and lattice types are explored in tests performed on benchmark protein sequences, using a Python implementation. Simulations give insights on the use of square and triangular lattices in a range of dimensions. The trend of potential minimum for sequences of different lengths, varying the lattice dimension, is uncovered. Moreover, an extensive quantitative characterization of the usage of the so-called "move types" is reported for the first time. The proposed general framework for the development of lattice models is simple yet complete, and an object-oriented architecture can be proficiently employed for the supporting software, by designing ad-hoc classes. The proposed framework represents a new general viewpoint that potentially subsumes a number of solutions previously studied. The adoption of the described model pushes to look at protein structure issues from a more general and essential perspective, making computational investigations over simplified models more straightforward as well.

Role of nonspecific interactions in molecular chaperones through model-based bioinformatics

Molecular chaperones are large proteins or protein complexes from which many proteins require assistance in order to fold. One unique property of molecular chaperones is the cavity they provide in which proteins fold. The interior surface residues which make up the cavities of molecular chaperone complexes from different organisms has recently been identified, including the well-studied GroEL-GroES chaperonin complex found in Escherichia coli. It was found that the interior of these protein complexes is significantly different than other protein surfaces and that the residues found on the protein surface are able to resist protein adsorption when immobilized on a surface. Yet it remains unknown if these residues passively resist protein binding inside GroEL-GroEs (as demonstrated by experiments that created synthetic mimics of the interior cavity) or if the interior also actively stabilizes protein folding. To answer this question, we have extended entropic models of substrate protein folding inside GroEL-GroES to include interaction energies between substrate proteins and the GroEL-GroES chaperone complex. This model was tested on a set of 528 proteins and the results qualitatively match experimental observations. The interior residues were found to strongly discourage the exposure of any hydrophobic residues, providing an enhanced hydrophobic effect inside the cavity that actively influences protein folding. This work provides both a mechanism for active protein stabilization in GroEL-GroES and a model that matches contemporary understanding of the chaperone protein.

Mimicking the action of folding chaperones by Hamiltonian replica-exchange molecular dynamics simulations: application in the refinement of de novo models

The efficiency of using a variant of Hamiltonian replica-exchange molecular dynamics (Chaperone H-replica-exchange molecular dynamics [CH-REMD]) for the refinement of protein structural models generated de novo is investigated. In CH-REMD, the interaction between the protein and its environment, specifically, the electrostatic interaction between the protein and the solvating water, is varied leading to cycles of partial unfolding and refolding mimicking some aspects of folding chaperones. In 10 of the 15 cases examined, the CH-REMD approach sampled structures in which the root-mean-square deviation (RMSD) of secondary structure elements (SSE-RMSD) with respect to the experimental structure was more than 1.0 Å lower than the initial de novo model. In 14 of the 15 cases, the improvement was more than 0.5 Å. The ability of three different statistical potentials to identify near-native conformations was also examined. Little correlation between the SSE-RMSD of the sampled structures with respect to the experimental structure and any of the scoring functions tested was found. The most effective scoring function tested was the DFIRE potential. Using the DFIRE potential, the SSE-RMSD of the best scoring structures was on average 0.3 Å lower than the initial model. Overall the work demonstrates that targeted enhanced-sampling techniques such as CH-REMD can lead to the systematic refinement of protein structural models generated de novo but that improved potentials for the identification of near-native structures are still needed.

Assisted peptide folding by surface pattern recognition

Natively disordered proteins belong to a unique class of biomolecules whose function is related to their flexibility and their ability to adopt desired conformations upon binding to substrates. In some cases these proteins can bind multiple partners, which can lead to distinct structures and promiscuity in functions. In other words, the capacity to recognize molecular patterns on the substrate is often essential for the folding and function of intrinsically disordered proteins. Biomolecular pattern recognition is extremely relevant both in vivo (e.g., for oligomerization, immune response, induced folding, substrate binding, and molecular switches) and in vitro (e.g., for biosensing, catalysis, chromatography, and implantation). Here, we use a minimalist computational model system to investigate how polar/nonpolar patterns on a surface can induce the folding of an otherwise unstructured peptide. We show that a model peptide that exists in the bulk as a molten globular state consisting of many interconverting structures can fold into either a helix-coil-helix or an extended helix structure in the presence of a complementary designed patterned surface at low hydrophobicity (3.7%) or a uniform surface at high hydrophobicity (50%). However, we find that a carefully chosen surface pattern can bind to and catalyze the folding of a natively unfolded protein much more readily or effectively than a surface with a noncomplementary or uniform distribution of hydrophobic residues.

Reconciling theories of chaperonin accelerated folding with experimental evidence

For the last 20 years, a large volume of experimental and theoretical work has been undertaken to understand how chaperones like GroEL can assist protein folding in the cell. The most accepted explanation appears to be the simplest: GroEL, like most other chaperones, helps proteins fold by preventing aggregation. However, evidence suggests that, under some conditions, GroEL can play a more active role by accelerating protein folding. A large number of models have been proposed to explain how this could occur. Focused experiments have been designed and carried out using different protein substrates with conclusions that support many different mechanisms. In the current article, we attempt to see the forest through the trees. We review all suggested mechanisms for chaperonin-mediated folding and weigh the plausibility of each in light of what we now know about the most stringent, essential, GroEL-dependent protein substrates.

Stability of tethered proteins

The stability of tethered globular proteins under denaturing conditions was interrogated with a hydrophobic surface, since conventional structural methods like circular dichroism (CD) and fluorescence or infrared spectroscopy could not be used because of the presence of an opaque solid substrate and extremely low surface concentrations. For free protein in solution, CD spectra gave well-known unfolding denaturing curves for lysozyme (LYS) and ribonuclease A (RNase A). The unfolding process for covalently tethered LYS and RNase A was followed, with multimolecular force spectroscopy (using an atomic force microscope in force-mode), via the adhesion energy between a functionalized self-assembled monolayer (CH(3)-SAM) probe and the protein molecules covalently bound to a carboxylic SAM on a gold-coated glass coverslip. The adhesion energy passed through a maximum for the tethered proteins during excursions with temperature or chemical denaturants. The initial rise in adhesion energy on increasing the temperature or GuHCl concentration was due to increasing exposure of the unfolded hydrophobic core of the proteins to the CH(3)-SAM tip, while the decrease in adhesion energy at high temperature or large concentrations of denaturant is attributed to interprotein association with nearest neighbors. Attempts to recover their folded state upon cooling (or reducing GuHCl concentration) were unsuccessful. Also, dilution of surface-tethered LYS reduced the aggregation with nearest neighbors about 6-fold. These results are in qualitative agreement with Monte Carlo simulations on a simple two-letter lattice protein model, especially for low concentrations of grafted proteins.

Thermodynamics and kinetics of protein folding under confinement

Understanding the effects of confinement on protein stability and folding kinetics is important for describing protein folding in the cellular environment. We have investigated the effects of confinement on two structurally distinct proteins as a function of the dimension d(c) and characteristic size R of the confining boundary. We find that the stabilization of the folded state relative to bulk conditions is quantitatively described by R(-gamma(c)), where the exponent gamma(c) is approximately 5/3 independent of the dimension of confinement d(c) (cylindrical, planar, or spherical). Moreover, we find that the logarithm of the folding rates also scale as R(-gamma(c)), with deviations only being seen for very small confining geometries, where folding is downhill; for both stability and kinetics, the dominant effect is the change in the free energy of the unfolded state. A secondary effect on the kinetics is a slight destabilization of the transition state by confinement, although the contacts present in the confined transition state are essentially identical to the bulk case. We investigate the effect of confinement on the position-dependent diffusion coefficients D(Q) for dynamics along the reaction coordinate Q (fraction of native contacts). The diffusion coefficients only change in the unfolded state basin, where they are increased because of compaction.

Rattling the cage: computational models of chaperonin-mediated protein folding

Chaperonins are known to maintain the stability of the proteome by facilitating the productive folding of numerous misfolded or aggregation-prone proteins and are thus essential for cell viability. Despite their established importance, the mechanism by which chaperonins facilitate protein folding remains unknown. Computer simulation techniques are now being employed to complement experimental ones in order to shed light on this mystery. Here we review previous computational models of chaperonin-mediated protein folding in the context of the two main hypotheses for chaperonin function: iterative annealing and landscape modulation. We then discuss new results pointing to the importance of solvent (a previously neglected factor) in chaperonin activity. We conclude with our views on the future role of simulation in studying chaperonin activity as well as protein folding in other biologically relevant confined contexts.

Do chaperonins boost protein yields by accelerating folding or preventing aggregation?

The GroEL chaperonin has the ability to behave as an unfoldase, repeatedly denaturing proteins upon binding, which in turn can free them from kinetic traps and increase their folding rates. The complex formed by GroEL+GroES+ATP can also act as an infinite dilution cage, enclosing proteins within a protective container where they can fold without danger of aggregation. Controversy remains over which of these two properties is more critical to the GroEL/ES chaperonin's function. We probe the importance of the unfoldase nature of GroEL under conditions where aggregation is the predominant protein degradation pathway. We consider the effect of a hypothetical mutation to GroEL which increases the cycle frequency of GroEL/ES by increasing the rate of hydrolysis of GroEL-bound ATP. Using a simple kinetic model, we show that this modified chaperonin would be self-defeating: any potential reduction in folding time would be negated by an increase in time spent in the bulk, causing an increase in aggregation and a net decrease in protein folding yields.

Different mechanistic requirements for prokaryotic and eukaryotic chaperonins: a lattice study

MOTIVATION

The folding of many proteins in vivo and in vitro is assisted by molecular chaperones. A well-characterized molecular chaperone system is the chaperonin GroEL/GroES from Escherichia coli which has a homolog found in the eukaryotic cytosol called CCT. All chaperonins have a ring structure with a cavity in which the substrate protein folds. An interesting difference between prokaryotic and eukaryotic chaperonins is in the nature of the ATP-mediated conformational changes that their ring structures undergo during their reaction cycle. Prokaryotic chaperonins are known to exhibit a highly cooperative concerted change of their cavity surface while in eukaryotic chaperonins the change is sequential. Approximately 70% of proteins in eukaryotic cells are multi-domain whereas in prokaryotes single-domain proteins are more common. Thus, it was suggested that the different modes of action of prokaryotic and eukaryotic chaperonins can be explained by the need of eukaryotic chaperonins to facilitate folding of multi-domain proteins.

RESULTS

Using a 2D square lattice model, we generated two large populations of single-domain and double-domain substrate proteins. Chaperonins were modeled as static structures with a cavity wall with which the substrate protein interacts. We simulated both concerted and sequential changes of the cavity surfaces and demonstrated that folding of single-domain proteins benefits from concerted but not sequential changes whereas double-domain proteins benefit also from sequential changes. Thus, our results support the suggestion that the different modes of allosteric switching of prokaryotic and eukaryotic chaperonin rings have functional implications as it enables eukaryotic chaperonins to better assist multi-domain protein folding.

A tale of two tails: why are terminal residues of proteins exposed?

MOTIVATION

It is widely known that terminal residues of proteins (i.e. the N- and C-termini) are predominantly located on the surface of proteins and exposed to the solvent. However, there is no good explanation as to the forces driving this phenomenon. The common explanation that terminal residues are charged, and charged residues prefer to be on the surface, cannot explain the magnitude of the phenomenon. Here, we survey a large number of proteins from the PDB in order to explore, quantitatively, this phenomenon, and then we use a lattice model to study the mechanisms involved.

RESULTS

The location of the termini was examined for 425 small monomeric proteins (50-200 amino acids) and it was found that the average solvent accessibility of termini residues is 87.1% compared with 49.2% of charged residues and 35.9% of all residues. Using a cutoff of 50% of the maximal possible exposure, 80.3% of the N-terminal and 86.1% of the C-terminal residues are exposed compared to 32% for all residues. In addition, terminal residues are much more distant from the center of mass of their proteins than other residues. Using a 2D lattice, a large population of model proteins was studied on three levels: structural selection of compact structures, thermodynamic selection of conformations with a pronounced energy gap and kinetic selection of fast folding proteins using Monte-Carlo simulations. Progressively, each selection raises the proportion of proteins with termini on the surface, resulting in similar proportions to those observed for real proteins.

GroEL-mediated protein folding: making the impossible, possible

Protein folding is a spontaneous process that is essential for life, yet the concentrated and complex interior of a cell is an inherently hostile environment for the efficient folding of many proteins. Some proteins-constrained by sequence, topology, size, and function-simply cannot fold by themselves and are instead prone to misfolding and aggregation. This problem is so deeply entrenched that a specialized family of proteins, known as molecular chaperones, evolved to assist in protein folding. Here we examine one essential class of molecular chaperones, the large, oligomeric, and energy utilizing chaperonins or Hsp60s. The bacterial chaperonin GroEL, along with its co-chaperonin GroES, is probably the best-studied example of this family of protein-folding machine. In this review, we examine some of the general properties of proteins that do not fold well in the absence of GroEL and then consider how folding of these proteins is enhanced by GroEL and GroES. Recent experimental and theoretical studies suggest that chaperonins like GroEL and GroES employ a combination of protein isolation, unfolding, and conformational restriction to drive protein folding under conditions where it is otherwise not possible.

Folding on the chaperone: yield enhancement through loose binding

A variety of small cageless chaperones have been discovered that can assist protein folding without the consumption of ATP. These include mini-chaperones (catalytically active fragments of larger chaperones), as well as small proteins such as alpha-casein and detergents acting as "artificial chaperones." These chaperones all possess exposed hydrophobic patches on their surface that act as recognition sites for misfolded proteins. They lack the complexity of chaperonins (that encapsulate proteins in their inner rings) and their study can offer insight into the minimal requirements for chaperone function. We use molecular dynamics simulations to investigate how a cageless chaperone, modeled as a sphere of tunable hydrophobicity, can assist folding of a substrate protein. We find that under steady-state (non-stress) conditions, cageless chaperones that bind to a single substrate protein increase folding yields by reducing the time the substrate spends in an aggregation-prone state in a dual manner: (a) by competing for aggregation-prone hydrophobic sites on the surface of a protein, hence reducing the time the protein spends unprotected in the bulk and (b) by accelerating folding rates of the protein. In both cases, the chaperone must bind to and hold the protein loosely enough to allow the protein to change its conformation and fold while bound. Loose binding may enable small cageless chaperones to help proteins fold and avoid aggregation under steady-state conditions, even at low concentrations, without the consumption of ATP.

Hydrophobic collapse in (in silico) protein folding

A model of hydrophobic collapse, which is treated as the driving force for protein folding, is presented. This model is the superposition of three models commonly used in protein structure prediction: (1) 'oil-drop' model introduced by Kauzmann, (2) a lattice model introduced to decrease the number of degrees of freedom for structural changes and (3) a model of the formation of hydrophobic core as a key feature in driving the folding of proteins. These three models together helped to develop the idea of a fuzzy-oil-drop as a model for an external force field of hydrophobic character mimicking the hydrophobicity-differentiated environment for hydrophobic collapse. All amino acids in the polypeptide interact pair-wise during the folding process (energy minimization procedure) and interact with the external hydrophobic force field defined by a three-dimensional Gaussian function. The value of the Gaussian function usually interpreted as a probability distribution is treated as a normalized hydrophobicity distribution, with its maximum in the center of the ellipsoid and decreasing proportionally with the distance versus the center. The fuzzy-oil-drop is elastic and changes its shape and size during the simulated folding procedure.

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.

Mimicking the action of GroEL in molecular dynamics simulations: application to the refinement of protein structures

Bacterial chaperonin, GroEL, together with its co-chaperonin, GroES, facilitates the folding of a variety of polypeptides. Experiments suggest that GroEL stimulates protein folding by multiple cycles of binding and release. Misfolded proteins first bind to an exposed hydrophobic surface on GroEL. GroES then encapsulates the substrate and triggers its release into the central cavity of the GroEL/ES complex for folding. In this work, we investigate the possibility to facilitate protein folding in molecular dynamics simulations by mimicking the effects of GroEL/ES namely, repeated binding and release, together with spatial confinement. During the binding stage, the (metastable) partially folded proteins are allowed to attach spontaneously to a hydrophobic surface within the simulation box. This destabilizes the structures, which are then transferred into a spatially confined cavity for folding. The approach has been tested by attempting to refine protein structural models generated using the ROSETTA procedure for ab initio structure prediction. Dramatic improvements in regard to the deviation of protein models from the corresponding experimental structures were observed. The results suggest that the primary effects of the GroEL/ES system can be mimicked in a simple coarse-grained manner and be used to facilitate protein folding in molecular dynamics simulations. Furthermore, the results support the assumption that the spatial confinement in GroEL/ES assists the folding of encapsulated proteins.

Nanopore-protein interactions dramatically alter stability and yield of the native state in restricted spaces

We have studied the stability and the yield of the folded WW domains in a spherical nanopore to provide insights into the changes in the folding characteristics due to interactions of the polypeptide (SP) with the walls of the pore. Using different models for the interactions between the nanopore and the polypeptide chain we have obtained results that are relevant to a broad range of experiments. (a) In the temperature and the strength of the SP-pore interaction plane (lambda), there are four "phases," namely, the unfolded state, the native state, the molten globule phase (MG), and the surface interaction-stabilized (SIS) state. The MG and SIS states are populated at moderate and large values of lambda, respectively. For a fixed pore size, the folding rates vary non-monotonically as lambda is varied with a maximum at lambda approximately 1 at which the SP-nanopore interaction is comparable to the stability of the native state. At large lambda values, the WW domain is kinetically trapped in the SIS states. Using multiple sequence alignment, we conclude that similar folding mechanism should be observed in other WW domains as well. (b) To mimic the changes in the nature of the allosterically driven SP-GroEL interactions we consider two models for the dynamic Anfinsen cage (DAC). In DAC1, the SP-cavity interaction cycles between hydrophobic (lambda>0) and hydrophilic (lambda=0) with a period tau. The yield of the native state is a maximum for an optimum value of tau=tau(OPT). At tau=tau(OPT), the largest yield of the native state is obtained when tau(H) approximately tau(P) where tau(H)(tau(P)) is the duration for which the cavity is hydrophobic (hydrophilic). Thus, in order to enhance the native state yield, the cycling rate, for a given loading rate of the GroEL nanomachine, should be maximized. In DAC2, the volume of the cavity is doubled (as happens when ATP and GroES bind to GroEL) and the SP-pore interaction simultaneously changes from hydrophobic to hydrophilic. In this case, we find greater increase in yield of the native state compared to DAC1 at all values of tau.

The mechanism of GroEL/GroES folding/refolding of protein substrates revisited

The thermodynamics and kinetics of zinc-cytochrome c (ZnCyt c) interactions with Escherichia coli molecular chaperone GroEL (Chaperonin 60; Cpn60) are described. Zinc(II)-porphyrin represents a flexible fluorescent probe for thermodynamic complex formation between GroEL and ZnCyt c, as well as for stopped-flow fluorescence kinetic experiments. Data suggests that GroEL and GroEL/GroES-assisted refolding of unfolded ZnCyt c takes place by a mechanism that is quite close to the Anfinsen Cage hypothesis for molecular chaperone activity. However, even in the presence of ATP, GroEL/GroES-assisted refolding of ZnCyt c takes place at approximately half the rate of refolding of ZnCyt c alone. On the other hand, there is little evidence for refolding behaviour consistent with the Iterative Annealing hypothesis. This includes a complete lack of GroEL or GroEL/GroES-assisted enhancement of refolding rate constant k(2) associated with the unfolding of a putative misfolded state I (Zn) on the pathway to the native state. Reviewing our data in the light of data from other laboratories, we observe that all forward rate enhancements or reductions could be accounted for in terms of thermodynamic coupling (adjusting positions of refolding equilibria) due to binding interactions between GroEL and unfolded protein substrates, driven by thermodynamic considerations. Therefore, we propose that passive kinetic partitioning should be considered the core mechanism of the GroEL/GroES molecular chaperone machinery, wherein the core function is to bind unfolded protein substrates leading to a blockade of aggregation pathways and to increases in molecular flux through productive folding pathway(s).

Folding zones inside the ribosomal exit tunnel

Helicity of membrane proteins can be manifested inside the ribosome tunnel, but the determinants of compact structure formation inside the tunnel are largely unexplored. Using an extended nascent peptide as a molecular tape measure of the ribosomal tunnel, we have previously demonstrated helix formation inside the tunnel. Here, we introduce a series of consecutive polyalanines into different regions of the tape measure to monitor the formation of compact structure in the nascent peptide. We find that the formation of compact structure of the polyalanine sequence depends on its location. Calculation of free energies for the equilibria between folded and unfolded nascent peptides in different regions of the tunnel shows that there are zones of secondary structure formation inside the ribosomal exit tunnel. These zones may have an active role in nascent-chain compaction.

Water and molecular chaperones act as weak links of protein folding networks: energy landscape and punctuated equilibrium changes point towards a game theory of proteins

Water molecules and molecular chaperones efficiently help the protein folding process. Here we describe their action in the context of the energy and topological networks of proteins. In energy terms water and chaperones were suggested to decrease the activation energy between various local energy minima smoothing the energy landscape, rescuing misfolded proteins from conformational traps and stabilizing their native structure. In kinetic terms water and chaperones may make the punctuated equilibrium of conformational changes less punctuated and help protein relaxation. Finally, water and chaperones may help the convergence of multiple energy landscapes during protein-macromolecule interactions. We also discuss the possibility of the introduction of protein games to narrow the multitude of the energy landscapes when a protein binds to another macromolecule. Both water and chaperones provide a diffuse set of rapidly fluctuating weak links (low affinity and low probability interactions), which allow the generalization of all these statements to a multitude of networks.

Quality control in the bacterial periplasm

Studies of the mechanisms that Gram-negative bacteria use to sense and respond to stress have led to a greater understanding of protein folding in both cytoplasmic and extracytoplasmic locations. In response to stressful conditions, bacteria induce a variety of stress response systems, examples of which are the sigma(E) and Cpx systems in Escherichia coli. Induction of these stress response systems results in upregulation of several gene targets that have been shown to be important for protein folding under normal conditions. Here we review the identification of stress response systems and their corresponding gene targets in E. coli. In addition, we discuss the apparent redundancy of the folding factors in the periplasm, and we consider the potential importance of the functional overlap that exists.

Protein structural transitions and their functional role

Living cells are a collection of molecular machines which carry out many of the functions essential for the cell's existence, differentiation and reproduction. Most, though not all, of these machines are made up of proteins. Because of their complexity, an understanding of how they work requires a synergistic combination of experimental and theoretical studies. In this paper we outline our studies of two such protein machines. One is GroEL, the chaperone from Escherichia coli, which aids in protein folding; the other is F(1)-ATPase, a motor protein which synthesizes and hydrolyses ATP.

The unfolding action of GroEL on a protein substrate

A molecular dynamics simulation of the active unfolding of denatured rhodanese by the chaperone GroEL is presented. The compact denatured protein is bound initially to the cis cavity and forms stable contacts with several of the subunits. As the cis ring apical domains of GroEL undergo the transition from the closed to the more open (ATP-bound) state, they exert a force on rhodanese that leads to the increased unfolding of certain loops. The contacts between GroEL and rhodanese are analyzed and their variation during the GroEL transition is shown. The major contacts, which give rise to the stretching force, are found to be similar to those observed in crystal structures of peptides bound to the apical domains. The results of the simulation show that multidomain interactions play an essential role, in accord with experiments. Implications of the results for mutation experiments and for the action of GroEL are discussed.

Accelerated folding in the weak hydrophobic environment of a chaperonin cavity: creation of an alternate fast folding pathway

Recent experiments suggest that the folding of certain proteins can take place entirely within a chaperonin-like cavity. These substrate proteins experience folding rate enhancements without undergoing multiple rounds of ATP-induced binding and release from the chaperonin. Rather, they undergo only a single binding event, followed by sequestration into the chaperonin cage. The present work uses molecular dynamics simulations to investigate the folding of a highly frustrated protein within this chaperonin cavity. The chaperonin interior is modeled by a sphere with a lining of tunable degree of hydrophobicity. We demonstrate that a moderately hydrophobic environment, similar to the interior of the GroEL cavity upon complexion with ATP and GroES, is sufficient to accelerate the folding of a frustrated protein by more than an order of magnitude. Our simulations support a mechanism by which the moderately hydrophobic chaperonin environment provides an alternate pathway to the native state through a transiently bound intermediate state.

Mimicking the action of folding chaperones in molecular dynamics simulations: Application to the refinement of homology-based protein structures

A novel method for the refinement of misfolded protein structures is proposed in which the properties of the solvent environment are oscillated in order to mimic some aspects of the role of molecular chaperones play in protein folding in vivo. Specifically, the hydrophobicity of the solvent is cycled by repetitively altering the partial charges on solvent molecules (water) during a molecular dynamics simulation. During periods when the hydrophobicity of the solvent is increased, intramolecular hydrogen bonding and secondary structure formation are promoted. During periods of increased solvent polarity, poorly packed regions of secondary structures are destabilized, promoting structural rearrangement. By cycling between these two extremes, the aim is to minimize the formation of long-lived intermediates. The approach has been applied to the refinement of structural models of three proteins generated by using the ROSETTA procedure for ab initio structure prediction. A significant improvement in the deviation of the model structures from the corresponding experimental structures was observed. Although preliminary, the results indicate computationally mimicking some functions of molecular chaperones in molecular dynamics simulations can promote the correct formation of secondary structure and thus be of general use in protein folding simulations and in the refinement of structural models of small- to medium-size proteins.

Model study of protein unfolding by interfaces

We study interface-induced protein unfolding on hydrophobic and polar interfaces by means of a two-dimensional lattice model and an exhaustive enumeration ground-state structure search, for a set of model proteins of length 20 residues. We compare the effects of the two types of interfaces, and search for criteria that influence the retention of a protein's native-state structure upon adsorption. We find that the unfolding proceeds by a large, sudden loss of native contacts. The unfolding at polar interfaces exhibits similar behavior to that at hydrophobic interfaces but with a much weaker interface coupling strength. Further, we find that the resistance of proteins to unfolding in our model is positively correlated with the magnitude of the folding energy in the native-state structure, the thermal stability (or energy gap) for that structure, and the interface energy for native-state adsorption. We find these factors to be of roughly equal importance.

Effects of confinement in chaperonin assisted protein folding: rate enhancement by decreasing the roughness of the folding energy landscape

Chaperonins, such as the GroE complex of the bacteria Escherichia coli, assist the folding of proteins under non-permissive folding conditions by providing a cavity in which the newly translated or translocated protein can be encapsulated. Whether the chaperonin cage plays a passive role in protecting the protein from aggregation, or an active role in accelerating folding rates, remains a matter of debate. Here, we investigate the role of confinement in chaperonin mediated folding through molecular dynamics simulations. We designed a substrate protein with an alpha/beta sandwich fold, a common structural motif found in GroE substrate proteins and confined it to a spherical hydrophilic cage which mimicked the interior of the GroEL/ES cavity. The thermodynamics and kinetics of folding were studied over a wide range of temperature and cage radii. Confinement was seen to significantly raise the collapse temperature, T(c), as a result of the associated entropy loss of the unfolded state. The folding temperature, T(f), on the other hand, remained unaffected by encapsulation, a consequence of the folding mechanism of this protein that involves an initial collapse to a compact misfolded state prior to rearranging to the native state. Folding rates were observed to be either accelerated or retarded compared to bulk folding rates, depending on the temperature of the simulation. Rate enhancements due to confinement were observed only at temperatures above the temperature T(m), which corresponds to the temperature at which the protein folds fastest. For this protein, T(m) lies above the folding temperature, T(f), implying that encapsulation alone will not lead to a rate enhancement under conditions where the native state is stable (T<T(f)). For confinement to positively impact folding rates under physiological conditions, it is hence necessary for the protein to exhibit a folding transition above the temperature at which it exhibits its fastest folding rate (T(m)<T(f)). We designed a protein with this property by reducing the energetic frustration in the original alpha/beta sandwich substrate protein. The modified protein exhibited a twofold acceleration in folding rates upon encapsulation. This rate enhancement is due to a mechanistic change in folding involving the elimination, upon encapsulation, of accessible local energy minima corresponding to structures with large radii of gyration. For this protein, confinement hence plays more than the role of a passive cage, but rather adopts an active role, accelerating folding rates by decreasing the roughness of the energy landscape of the protein.

Folding with and without encapsulation by cis- and trans-only GroEL-GroES complexes

Although a cis mechanism of GroEL-mediated protein folding, occurring inside a hydrophilic chamber encapsulated by the co-chaperonin GroES, has been well documented, recently the GroEL-GroES-mediated folding of aconitase, a large protein (82 kDa) that could not be encapsulated, was described. This process required GroES binding to the ring opposite the polypeptide (trans) to drive release and productive folding. Here, we have evaluated this mechanism further using trans-only complexes in which GroES is closely tethered to one of the two GroEL rings, blocking polypeptide binding by that ring. In vitro, trans-only folded aconitase with kinetics identical to GroEL-GroES. Surprisingly, trans-only also folded smaller GroEL-GroES-dependent substrates, Rubisco and malate dehydrogenase, but at rates slower than the cis reaction. Remarkably, in vivo, a plasmid encoding a trans-only complex rescued a GroEL-deficient strain, but the colony size was approximately one-tenth that produced by wild-type GroEL-GroES. We conclude that a trans mechanism, involving rounds of binding to an open ring and direct release into the bulk solution, can be generally productive although, where size permits, cis encapsulation supports more efficient folding.

Directed evolution of substrate-optimized GroEL/S chaperonins

GroEL/S chaperonin ring complexes fold many unrelated proteins. To understand the basis and extent of the chaperonin substrate spectrum, we used rounds of selection and DNA shuffling to obtain GroEL/S variants that dramatically enhanced folding of a single substrate-green fluorescent protein (GFP). Changes in the substrate-optimized chaperonins increase the polarity of the folding cavity and alter the ATPase cycle. These findings reveal a surprising plasticity of GroEL/S, which can be exploited to aid folding of recombinant proteins. Our studies also reveal a conflict between specialization and generalization of chaperonins as increased GFP folding comes at the expense of the ability of GroEL/S to fold its natural substrates. This conflict and the nature of the ring structure may help explain the evolution of cellular chaperone systems.

Application of a chaperone-based refolding method to two- and three-dimensional off-lattice protein models

A model of protein-chaperone interaction as a two-phase (unfolding/refolding) iterative annealing mechanism able to promote structural segregation of hydrophobic and hydrophilic monomers and thereby facilitate access to nativelike states has recently been applied successfully to two 22-mers of the Honeycutt and Thirumalai BLN (hydrophobic, hydrophilic, neutral) heteropolymer model. This technique is here applied to a much wider data set: 94 8-mers of the off-lattice protein model originally presented in two dimensions by Stillinger and Head-Gordon, and later extended into three dimensions by Irbäck and Potthast; the model chaperone is shown to be equally successful, and by progressive elaboration of the chaperone model as in the earlier BLN model work, to be utilizing very similar underlying mechanisms. It is demonstrated that on average, contacts with the model chaperone give rise to a consistent movement in structure space in the direction of more nativelike structures; this method of global minimization does not therefore rely fundamentally on random search. Insofar as the responses to the chaperone of the two- and three-dimensional forms of the substrate model do differ, this can be interpreted as reflecting the different handling of hydrophilic monomers in the models-in particular, whether there is active repulsion between these and monomers of hydrophobic character. The chaperone-induced refolding method is also tested on a set of 220 9-mer chains of each version of the substrate model, where it is seen that the two-dimensional model, with its more clearly distinguished roles for the hydrophobic and hydrophilic monomers, shows a more favorable scaling behavior.

GroEL channels the folding of thioredoxin along one kinetic route

Many proteins display complex folding kinetics, which represent multiple parallel folding pathways emanating from multiple unfolded forms and converging to the unique native form. The small protein thioredoxin from Escherichia coli is one such protein. The effect of the chaperonin GroEL on modulating the complex energy landscape that separates the unfolded ensemble from the native state of thioredoxin has been studied. It is shown that while the fluorescence change accompanying folding occurs in five kinetic phases in the absence of GroEL, only the two slowest kinetic phases are discernible in the presence of saturating concentrations of GroEL. This result is shown to be consistent with only one out of several available folding routes being operational in the presence of GroEL. It is shown that native protein, which forms via fast as well as slow routes in the absence of GroEL, forms only via a slow route in its presence. The effect of GroEL on the folding of thioredoxin is shown to be the consequence of it binding differentially to the many folding-competent forms. While some of these forms can continue folding when bound to GroEL, others cannot. All molecules are then drawn into the operational folding route by the law of mass action. This observation indicates a new role for GroEL, which is to bias the energy landscape of a folding polypeptide towards fewer available pathways. It is suggested that such channeling might be a mechanism to avoid possible aggregation-prone routes available to a refolding polypeptide in vivo.

Global minimization of an off-lattice potential energy function using a chaperone-based refolding method

A global energy minimization method based on what is known about the mechanisms of the GroEL/GroES chaperonin system is applied to two 22-mers of an off-lattice protein model whose native states are beta-hairpins and which have structural similarity to short peptides known to interact strongly with the GroEL substrate binding domain. These model substrates have been used by other workers to test the effectiveness of a number of global minimization techniques, and are regarded as providing a significant challenge. The minimization method developed here is progressively elaborated from an initial simple form that targets exposed hydrophobic regions for unfolding to include a refolding phase that encourages the later recompactification of partly unfolded substrate; this refolding phase is seen to be crucial in the successful application of the method. The optimal handling of hydrophilic monomers within the model is also systematically explored, and it is seen that the best interpretation of their role is one that allows the chaperonin model to operate in "proofreading" mode whereby misfolded substrates are recognized by their surface exposure of a large proportion of hydrophobic monomers. The final version of the model allows native-like structures to be found quickly, on average for the two 22-mer substrates after 6 or 7 chaperone contacts. These results compare very favorably with those that have been obtained elsewhere using generic global minimization methods such as those based on thermal annealing. The paper concludes with a discussion of the place of the technique within the general category of hypersurface deformation methods for global minimization, and with suggestions as to how the chaperone-based method developed here could be elaborated so as to be effective on longer substrate chains that give rise to more complex tertiary structures in their native states.

Chaperonin-mediated protein folding

Molecular chaperones are required to assist folding of a subset of proteins in Escherichia coli. We describe a conceptual framework for understanding how the GroEL-GroES system assists misfolded proteins to reach their native states. The architecture of GroEL consists of double toroids stacked back-to-back. However, most of the fundamentals of the GroEL action can be described in terms of the single ring. A key idea in our framework is that, with coordinated ATP hydrolysis and GroES binding, GroEL participates actively by repeatedly unfolding the substrate protein (SP), provided that it is trapped in one of the misfolded states. We conjecture that the unfolding of SP becomes possible because a stretching force is transmitted to the SP when the GroEL particle undergoes allosteric transitions. Force-induced unfolding of the SP puts it on a higher free-energy point in the multidimensional energy landscape from which the SP can either reach the native conformation with some probability or be trapped in one of the competing basins of attraction (i.e., the SP undergoes kinetic partitioning). The model shows, in a natural way, that the time scales in the dynamics of the allosteric transitions are intimately coupled to folding rates of the SP. Several scenarios for chaperonin-assisted folding emerge depending on the interplay of the time scales governing the cycle. Further refinement of this framework may be necessary because single molecule experiments indicate that there is a great dispersion in the time scales governing the dynamics of the chaperonin cycle.

Multiple-sequence information provides protection against mis-specified potential energy functions in the lattice model of proteins

The neutral nets in the lattice models of proteins are composed of a group of similar sequences that encode for the same native structure. The maximal neutral nets in an HP and an AB lattice model are investigated by exhaustive enumeration in the conformation space. The result demonstrates that the performance of mis-specified potential functions can be improved significantly by averaging over sequences in the neutral nets.