How protein thermodynamics and folding mechanisms are altered by the chaperonin cage: molecular simulations

Confinement effects on the kinetics and thermodynamics of protein dimerization

In the cell, protein complexes form by relying on specific interactions between their monomers. Excluded volume effects due to molecular crowding would lead to correlations between molecules even without specific interactions. What is the interplay of these effects in the crowded cellular environment? We study dimerization of a model homodimer when the mondimers are free and when they are tethered to each other. We consider a structured environment: Two monomers first diffuse into a cavity of size L and then fold and bind within the cavity. The folding and binding are simulated by using molecular dynamics based on a simplified topology based model. The confinement in the cell is described by an effective molecular concentration C approximately L(-3). A two-state coupled folding and binding behavior is found. We show the maximal rate of dimerization occurred at an effective molecular concentration C(op) approximately = 1 mM, which is a relevant cellular concentration. In contrast, for tethered chains the rate keeps at a plateau when C < C(op) but then decreases sharply when C > C(op). For both the free and tethered cases, the simulated variation of the rate of dimerization and thermodynamic stability with effective molecular concentration agrees well with experimental observations. In addition, a theoretical argument for the effects of confinement on dimerization is also made.

Effect of the ordered water on protein folding: an off-lattice Gō-like model study

Recent experiments and numerical simulations have shown that the water molecules confined on the surfaces of some substrates, including the surfaces of cellular components in tissues and cells, form icelike ordered structures. If a protein folds in an environment with those icelike ordered water molecules, its behavior may be different from that in bulk water. Here, the effect of this ordered water environment on protein folding is studied by using an off-lattice Gō-like model. It is found that the ordered water environment significantly improves the native state stability and greatly speeds up the folding rate of the proteins.

Monte Carlo simulations of proteins in cages: influence of confinement on the stability of intermediate states

We study the folding of small proteins inside confining potentials. Proteins are described using an effective potential model that contains the Ramachandran angles as degrees of freedom and does not need any a priori information about the native state. Hydrogen bonds, dipole-dipole-, and hydrophobic interactions are taken explicitly into account. An interesting feature displayed by this potential is the presence of metastable intermediates between the unfolded and native states. We consider different types of confining potentials to describe proteins folding inside cages with repulsive or attractive walls. Using the Wang-Landau algorithm, we determine the density of states and analyze in detail the thermodynamical properties of the confined proteins for different sizes of the cages. We show that confinement dramatically reduces the phase space available to the protein and that the presence of intermediate states can be controlled by varying the properties of the confining potential. Cages with strongly attractive walls destabilize the intermediate states and lead to a two-state folding into a configuration that is less stable than the native structure. However, cages with slightly attractive walls enhance the stability of native structure and induce a folding process, which occurs through intermediate configurations.

Development of free-energy-based models for chaperonin containing TCP-1 mediated folding of actin

A free-energy-based approach is used to describe the mechanism through which chaperonin-containing TCP-1 (CCT) folds the filament-forming cytoskeletal protein actin, which is one of its primary substrates. The experimental observations on the actin folding and unfolding pathways are collated and then re-examined from this perspective, allowing us to determine the position of the CCT intervention on the actin free-energy folding landscape. The essential role for CCT in actin folding is to provide a free-energy contribution from its ATP cycle, which drives actin to fold from a stable, trapped intermediate I3, to a less stable but now productive folding intermediate I2. We develop two hypothetical mechanisms for actin folding founded upon concepts established for the bacterial type I chaperonin GroEL and extend them to the much more complex CCT system of eukaryotes. A new model is presented in which CCT facilitates free-energy transfer through direct coupling of the nucleotide hydrolysis cycle to the phases of actin substrate maturation.

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.

Dynamics of a self-avoiding polymer chain in slit, tube, and cube confinements

Monte Carlo simulations are presented for the observation of the dynamics of a self-avoiding polymer in three types of confinement: slit, tube, and cube. We pay special attention to the parameter regime where the characteristic confinement dimension is smaller than the radius of gyration of the unconfined polymer. On the basis of the bond-fluctuation model, we measured the rotation time of the end-to-end vector of the polymer, the diffusion time for the center of the polymer to move a distance comparable to the root mean square end-to-end distance, and the looping time for the ends of the polymer to approach each other from an open position. As functions of the confinement width and polymer length, these three time scales are discussed in light of scaling theories.

Characterization of the unfolded state of repeat proteins

The unfolded state ensemble of proteins has been described as a structurally featureless state. While this approach is supported by the fact that many unfolded proteins follow the scaling law behavior of a random coil, there is evidence that the unfolded states of various proteins are stabilized by native or non-native interactions. Recently, the existence of extensive non-native structure was reported for a repeat protein, which resulted in a scaling law exponent that is significantly smaller than that of a random polymer [Cortajarena et al., J. Mol. Biol. 382(1), 203-212 (2008)]. It was concluded that the high compactness of this protein stems from a significant fraction of interacting PP(II) helical segments in the unfolded state. In this study, we aim at providing possible molecular understanding of this anomalous compactness of the unfolded state and to investigate its origin. Using a hierarchy of computational models, we ask whether in general the unfolded state of a repeat protein is likely to be intrinsically more compact than the unfolded state of globular proteins, or whether this phenomenon depends mostly on the occurrence of a specific sequence that promotes PP(II) conformations. Our results suggest that the formation of the PP(II) conformation is indeed essential, yet the recurring sequence of repeat proteins promotes the interactions between these PP(II) segments and the formation of non-native interactions in the unfolded state.

Folding of elongated proteins: conventional or anomalous?

The complexity of the mechanisms by which proteins fold has been shown by many studies to be governed by their native-state topologies. This was manifested in the ability of the native topology-based model to capture folding mechanisms and the success of folding rate predictions based on various topological measures, such as the contact order. However, while the finer details of topological complexity have been thoroughly examined and related to folding kinetics, simpler characteristics of the protein, such as its overall shape, have been largely disregarded. In this study, we investigated the folding of proteins with an unusual elongated geometry that differs substantially from the common globular structure. To study the effect of the elongation degree on the folding kinetics, we used repeat proteins, which become more elongated as they include more repeating units. Some of these have apparently anomalous experimental folding kinetics, with rates that are often less than expected on the basis of rates for globular proteins possessing similar topological complexity. Using experimental folding rates and a larger set of rates obtained from simulations, we have shown that as the protein becomes increasingly elongated, its folding kinetics becomes slower and deviates more from the rate expected on the basis of topology measures fitted for globular proteins. The observed slow kinetics is a result of a more complex pathway in which stable intermediates composed of several consecutive repeats can appear. We thus propose a novel measure, an elongation-sensitive contact order, that takes into account both the extent of elongation and the topological complexity of the protein. This new measure resolves the apparent discrimination between the folding of globular and elongated repeat proteins. Our study extends the current capabilities of folding-rate predictions by unifying the kinetics of repeat and globular proteins.

Lattice model of equilibrium polymerization. VII. Understanding the role of "cooperativity" in self-assembly

Cooperativity is an emergent many-body phenomenon related to the degree to which elementary entities (particles, molecules, organisms) collectively interact to form larger scale structures. From the standpoint of a formal mean field description of chemical reactions, the cooperativity index m, describing the number of elements involved in this structural self-organization, is the order of the reaction. Thus, m for molecular self-assembly is the number of molecules in the final organized structure, e.g., spherical micelles. Although cooperativity is crucial for regulating the thermodynamics and dynamics of self-assembly, there is a limited understanding of this aspect of self-assembly. We analyze the cooperativity by calculating essential thermodynamic properties of the classical mth order reaction model of self-assembly (FAm model), including universal scaling functions describing the temperature and concentration dependence of the order parameter and average cluster size. The competition between self-assembly and phase separation is also described. We demonstrate that a sequential model of thermally activated equilibrium polymerization can quantitatively be related to the FAm model. Our analysis indicates that the essential requirement for "cooperative" self-assembly is the introduction of constraints (often nonlocal) acting on the individual assembly events to regulate the thermodynamic free energy landscape and, thus, the thermodynamic sharpness of the assembly transition. An effective value of m is defined for general self-assembly transitions, and we find a general tendency for self-assembly to become a true phase transition as m-->infinity. Finally, various quantitative measures of self-assembly cooperativity are discussed in order to identify experimental signatures of cooperativity in self-assembling systems and to provide a reliable metric for the degree of transition cooperativity.

Simulations of proteins with inhomogeneous degrees of freedom: The effect of thermostats

Proteins with a long flexible polymeric tail attached at their N- or C-terminus are studied using molecular dynamics (MD) simulations of a coarse-grained model for protein folding where the temperature is regulated by either the Berendsen or the Langevin thermostat. These thermostats show different abilities to regulate the temperature of these systems that include flexible and more rigid regions. In the simulations with the Berendsen thermostat, the flexible tail is significantly hotter than the protein, both in its folded and unfolded states. Upon weakening the strength of the Berendsen thermostat, the temperature gradient between the fast and the slow degrees of freedom is significantly decreased, yet linkage between the temperatures of the flexible tail and the protein remains. The Langevin thermostat is proven to regulate the temperature of these inhomogenous systems reliably, without discriminating between the slow and fast degrees of freedom. The Langevin thermostat is less sensitive than is the Berendsen thermostat to the strength of the coupling between the protein system and the thermal bath. Our study calls for special care in choosing the thermostat for MD simulations of systems with inhomogenous degrees of freedom. Using the Berendsen thermostat with strong coupling would result in mistaken thermodynamic descriptions of such systems.

Effects of confinement and crowding on folding of model proteins

We perform molecular dynamics simulations for a simple coarse-grained model of crambin placed inside of a softly repulsive sphere of radius R. The confinement makes folding at the optimal temperature slower and affects the folding scenarios, but both effects are not dramatic. The influence of crowding on folding are studied by placing several identical proteins within the sphere, denaturing them, and then by monitoring refolding. If the interactions between the proteins are dominated by the excluded volume effects, the net folding times are essentially like for a single protein. An introduction of inter-proteinic attractive contacts hinders folding when the strength of the attraction exceeds about a half of the value of the strength of the single protein contacts. The bigger the strength of the attraction, the more likely is the occurrence of aggregation and misfolding.

Potential for modulation of the hydrophobic effect inside chaperonins

Despite the spontaneity of some in vitro protein-folding reactions, native folding in vivo often requires the participation of barrel-shaped multimeric complexes known as chaperonins. Although it has long been known that chaperonin substrates fold upon sequestration inside the chaperonin barrel, the precise mechanism by which confinement within this space facilitates folding remains unknown. We examine the possibility that the chaperonin mediates a favorable reorganization of the solvent for the folding reaction. We discuss the effect of electrostatic charge on solvent-mediated hydrophobic forces in an aqueous environment. Based on these physical arguments, we construct a simple, phenomenological theory for the thermodynamics of density and hydrogen-bond order fluctuations in liquid water. Within the framework of this model, we investigate the effect of confinement inside a chaperonin-like cavity on the configurational free energy of water by calculating solvent free energies for cavities corresponding to the different conformational states in the ATP-driven catalytic cycle of the prokaryotic chaperonin GroEL. Our findings suggest that one function of chaperonins may involve trapping unfolded proteins and subsequently exposing them to a microenvironment in which the hydrophobic effect, a crucial thermodynamic driving force for folding, is enhanced.

Effect of glycosylation on protein folding: a close look at thermodynamic stabilization

Glycosylation is one of the most common posttranslational modifications to occur in protein biosynthesis, yet its effect on the thermodynamics and kinetics of proteins is poorly understood. A minimalist model based on the native protein topology, in which each amino acid and sugar ring was represented by a single bead, was used to study the effect of glycosylation on protein folding. We studied in silico the folding of 63 engineered SH3 domain variants that had been glycosylated with different numbers of conjugated polysaccharide chains at different sites on the protein's surface. Thermal stabilization of the protein by the polysaccharide chains was observed in proportion to the number of attached chains. Consistent with recent experimental data, the degree of thermal stabilization depended on the position of the glycosylation sites, but only very weakly on the size of the glycans. A thermodynamic analysis showed that the origin of the enhanced protein stabilization by glycosylation is destabilization of the unfolded state rather than stabilization of the folded state. The higher free energy of the unfolded state is enthalpic in origin because the bulky polysaccharide chains force the unfolded ensemble to adopt more extended conformations by prohibiting formation of a residual structure. The thermodynamic stabilization induced by glycosylation is coupled with kinetic stabilization. The effects introduced by the glycans on the biophysical properties of proteins are likely to be relevant to other protein polymeric conjugate systems that regularly occur in the cell as posttranslational modifications or for biotechnological purposes.

Effect of glycosylation on protein folding: a close look at thermodynamic stabilization

Glycosylation is one of the most common posttranslational modifications to occur in protein biosynthesis, yet its effect on the thermodynamics and kinetics of proteins is poorly understood. A minimalist model based on the native protein topology, in which each amino acid and sugar ring was represented by a single bead, was used to study the effect of glycosylation on protein folding. We studied in silico the folding of 63 engineered SH3 domain variants that had been glycosylated with different numbers of conjugated polysaccharide chains at different sites on the protein's surface. Thermal stabilization of the protein by the polysaccharide chains was observed in proportion to the number of attached chains. Consistent with recent experimental data, the degree of thermal stabilization depended on the position of the glycosylation sites, but only very weakly on the size of the glycans. A thermodynamic analysis showed that the origin of the enhanced protein stabilization by glycosylation is destabilization of the unfolded state rather than stabilization of the folded state. The higher free energy of the unfolded state is enthalpic in origin because the bulky polysaccharide chains force the unfolded ensemble to adopt more extended conformations by prohibiting formation of a residual structure. The thermodynamic stabilization induced by glycosylation is coupled with kinetic stabilization. The effects introduced by the glycans on the biophysical properties of proteins are likely to be relevant to other protein polymeric conjugate systems that regularly occur in the cell as posttranslational modifications or for biotechnological purposes.

An entropic perspective of protein stability on surfaces

The interaction of proteins with surfaces regulates numerous processes in nature, science, and technology. In many applications, it is desirable to place proteins on surfaces in an active state, and tethering represents one manner in which to accomplish this. However, a clear understanding of how tether placement and design affects protein activity is lacking. Available theoretical models predict that proteins will be stabilized when tethered to substrates. Such models suggest that the surface reduces the number of states accessible to the unfolded state of the protein, thereby reducing the entropic cost of folding on the surface compared to the bulk case. Recent studies, however, have shown that this stabilization is not always seen. The purpose of this article is to determine the validity of the theory with a thorough thermodynamic analysis of the folding of peptides attached to surfaces. Configuration-temperature-density-of-states Monte Carlo simulations are used to examine the behavior of four different peptides of different secondary and tertiary structure. It is found that the surface does reduce the entropic cost of folding for tethered peptides, as the theory suggests. This effect, however, does not always translate into improved stability because the surface may also have a destabilizing enthalpic effect. The theory neglects this effect and assumes that the enthalpy of folding is the same on and off the surface. Both the enthalpic and entropic contributions to the stability are found to be topology- and tether-placement-specific; we show that stability cannot be predicted a priori. A detailed analysis of the folding of protein A shows how the same protein can be both stabilized and destabilized on a surface depending upon how the tethering enhances or hinders the ability of the peptide to form correct tertiary structures.

Molecular simulation of protein dynamics in nanopores. I. Stability and folding

Discontinuous molecular dynamics simulations, together with the protein intermediate resolution model, an intermediate-resolution model of proteins, are used to carry out several microsecond-long simulations and study folding transition and stability of alpha-de novo-designed proteins in slit nanopores. Both attractive and repulsive interaction potentials between the proteins and the pore walls are considered. Near the folding temperature T(f) and in the presence of the attractive potential, the proteins undergo a repeating sequence of folding/partially folding/unfolding transitions, with T(f) decreasing with decreasing pore sizes. The unfolded states may even be completely adsorbed on the pore's walls with a negative potential energy. In such pores the energetic effects dominate the entropic effects. As a result, the unfolded state is stabilized, with a folding temperature T(f) which is lower than its value in the bulk and that, compared with the bulk, the folding rate decreases. The opposite is true in the presence of a repulsive interaction potential between the proteins and the walls. Moreover, for short proteins in very tight pores with attractive walls, there exists an unfolded state with only one alpha-helical hydrogen bond and an energy nearly equal to that of the folded state. The proteins have, however, high entropies, implying that they cannot fold onto their native structure, whereas in the presence of repulsive walls the proteins do attain their native structure. There is a pronounced asymmetry between the two termini of the protein with respect to their interaction with the pore walls. The effect of a variety of factors, including the pore size and the proteins' length, as well as the temperature, is studied in detail.

Essential role of the chaperonin folding compartment in vivo

The GroEL/GroES chaperonin system of Escherichia coli forms a nano-cage allowing single protein molecules to fold in isolation. However, as the chaperonin can also mediate folding independently of substrate encapsulation, it remained unclear whether the folding cage is essential in vivo. To address this question, we replaced wild-type GroEL with mutants of GroEL having either a reduced cage volume or altered charge properties of the cage wall. A stepwise reduction in cage size resulted in a gradual loss of cell viability, although the mutants bound non-native protein efficiently. Strikingly, a mild reduction in cage size increased the yield and the apparent rate of green fluorescent protein folding, consistent with the view that an effect of steric confinement can accelerate folding. As shown in vitro, the observed acceleration of folding was dependent on protein encapsulation by GroES but independent of GroES cycling regulated by the GroEL ATPase. Altering the net-negative charge of the GroEL cage wall also strongly affected chaperonin function. Based on these findings, the GroEL/GroES compartment is essential for protein folding in vivo.

Kinetic model for the coupling between allosteric transitions in GroEL and substrate protein folding and aggregation

The bacterial chaperonin GroEL and the co-chaperonin GroES assist in the folding of a number of structurally unrelated substrate proteins (SPs). In the absence of chaperonins, SP folds by the kinetic partitioning mechanism (KPM), according to which a fraction of unfolded molecules reaches the native state directly, while the remaining fraction gets trapped in a potentially aggregation-prone misfolded state. During the catalytic reaction cycle, GroEL undergoes a series of allosteric transitions (T<-->R-->R"-->T) triggered by SP capture, ATP binding and hydrolysis, and GroES binding. We developed a general kinetic model that takes into account the coupling between the rates of the allosteric transitions and the folding and aggregation of the SP. Our model, in which the GroEL allosteric rates and SP-dependent folding and aggregation rates are independently varied without prior assumption, quantitatively fits the GroEL concentration-dependent data on the yield of native ribulose bisphosphate carboxylase/oxygenase (Rubisco) as a function of time. The extracted kinetic parameters for the GroEL reaction cycle are consistent with the available values from independent experiments. In addition, we also obtained physically reasonable parameters for the kinetic steps in the reaction cycle that are difficult to measure. If experimental values for GroEL allosteric rates are used, the time-dependent changes in native-state yield at eight GroEL concentrations can be quantitatively fit using only three SP-dependent parameters. The model predicts that the differences in the efficiencies (as measured by yields of the native state) of GroEL, single-ring mutant (SR1), and variants of SR1, in the rescue of mitochondrial malate dehydrogenase, citrate synthase, and Rubisco, are related to the large variations in the allosteric transition rates. We also show that GroEL/S mutants that efficiently fold one SP at the expense of all others are due to a decrease in the rate of a key step in the reaction cycle, which implies that wild-type GroEL has evolved as a compromise between generality and specificity. We predict that, under maximum loading conditions and saturating ATP concentration, the efficiency of GroEL (using parameters for Rubisco) depends predominantly on the rate of R-->R" transition, while the equilibrium constant of the T<-->R has a small effect only. Both under sub- and superstoichiometric GroEL concentrations, enhanced efficiency is achieved by rapid turnover of the reaction cycle, which is in accord with the predictions of the iterative annealing mechanism. The effects are most dramatic at substoichiometric conditions (most relevant for in vivo situations) when SP aggregation can outcompete capture of SP by chaperonins.

In silico chaperonin-like cycle helps folding of proteins for structure prediction

Currently, one of the most serious problems in protein-folding simulations for de novo structure prediction is conformational sampling of medium-to-large proteins. In vivo, folding of these proteins is mediated by molecular chaperones. Inspired by the functions of chaperonins, we designed a simple chaperonin-like simulation protocol within the framework of the standard fragment assembly method: in our protocol, the strength of the hydrophobic interaction is periodically modulated to help the protein escape from misfolded structures. We tested this protocol for 38 proteins and found that, using a certain defined criterion of success, our method could successfully predict the native structures of 14 targets, whereas only those of 10 targets were successfully predicted using the standard protocol. In particular, for non-alpha-helical proteins, our method yielded significantly better predictions than the standard approach. This chaperonin-inspired protocol that enhanced de novo structure prediction using folding simulations may, in turn, provide new insights into the working principles underlying the chaperonin system.

Helix formation inside a nanotube: possible influence of backbone-water hydrogen bonding by the confining surface through modulation of water activity

Recent molecular dynamics simulations of Sorin and Pande [J. Am. Chem. Soc. 128, 6316 (2006)] in explicit solvent found that helix formation of an alanine peptide is disfavored inside a nanotube relative to that in bulk solution. Here, we present a theory to quantitatively rationalize their simulation results. The basic idea is that the nonpolar inner surface of the nanotube creates a depletion layer and raises the activity of the confined water. The raised water activity, in turn, stabilizes the coil state through hydrogen bonding with the backbone amides and carbonyls. We account for the influence of water activity on helix formation within the Lifson-Roig theory. With physically reasonable parameters, the dependence of the helical content on the diameter of the nanotube obtained in the simulations is well reproduced.

Chaperonins: The hunt for the Group II mechanism

Chaperonins are multi-subunit complexes that enhance the efficiency of protein-folding reactions by capturing protein substrates in their central cavities. They occur in all prokaryotic and eukaryotic cell types and, alone amongst molecular chaperones, chaperonin knockouts are always lethal. Chaperonins come in two forms; the Group I are found in bacteria, mitochondria and plastids [W.A. Fenton, A.L. Horwich, Q. Rev. Biophys. 36 (2003) 229-256, [1]] and the Group II in the eukaryotic cytoplasm and in archaea [N.J. Cowan, S.A. Lewis, Adv. Protein Chem. 59 (2001) 73-104, [2]]. Both use energy derived from ATP binding and hydrolysis to drive a series of structural rearrangements that enable them to capture, engulf and then release polypeptide chains that have either not yet acquired the native, biologically active state or have been denatured in the cell.

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.

GroEL stimulates protein folding through forced unfolding

Many proteins cannot fold without the assistance of chaperonin machines like GroEL and GroES. The nature of this assistance, however, remains poorly understood. Here we demonstrate that unfolding of a substrate protein by GroEL enhances protein folding. We first show that capture of a protein on the open ring of a GroEL-ADP-GroES complex, GroEL's physiological acceptor state for non-native proteins in vivo, leaves the substrate protein in an unexpectedly compact state. Subsequent binding of ATP to the same GroEL ring causes rapid, forced unfolding of the substrate protein. Notably, the fraction of the substrate protein that commits to the native state following GroES binding and protein release into the GroEL-GroES cavity is proportional to the extent of substrate-protein unfolding. Forced protein unfolding is thus a central component of the multilayered stimulatory mechanism used by GroEL to drive protein folding.

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.

Comparative modeling of proteins

Three-dimensional analysis of protein structures is proving to be one of the most fruitful modes of biological and medical discovery in the early 21st century, providing fundamental insight into many (perhaps most) biochemical functions of relevance to the cause and treatment of diseases. Fully realizing such insight, however, would require analysis of too many distinct proteins for thorough laboratory analysis of all proteins to be feasible, thus, any method capable of accurate, efficient in silico structure prediction should prove highly expeditious. The technique generally acknowledged to provide the most accurate protein structure predictions, called comparative modeling, has, thus, attracted substantial attention and is the focus of this chapter. Although other reviews have reported on the method development and research history of comparative modeling, our discussion herein focuses on the general philosophy of the method and specific strategies for successfully achieving reliable and accurate models. The chapter, thus, relates aspects of template selection, sequence alignment, spatial alignment, loop and gap modeling, side chain modeling, structural refinement, and validation.

Macromolecular crowding and confinement: biochemical, biophysical, and potential physiological consequences

Expected and observed effects of volume exclusion on the free energy of rigid and flexible macromolecules in crowded and confined systems, and consequent effects of crowding and confinement on macromolecular reaction rates and equilibria are summarized. Findings from relevant theoretical/simulation and experimental literature published from 2004 onward are reviewed. Additional complexity arising from the heterogeneity of local environments in biological media, and the presence of nonspecific interactions between macromolecules over and above steric repulsion, are discussed. Theoretical and experimental approaches to the characterization of crowding- and confinement-induced effects in systems approaching the complexity of living organisms are suggested.

Large multiprotein structures modeling and simulation: the need for mesoscopic models

Recent observational techniques based upon confocal microscopy make it possible to observe cells at a scale that has never been probed before: the mesoscopic scale. In the eukaryotic cell nucleus, many objects demonstrating phenomena occurring at this scale, such as nuclear bodies, are current subjects of investigations. But from a modeling perspective, this scale has not been widely explored, and hence there is a lack of suitable models for such studies. By reviewing higher and lower scale modeling techniques, we analyze their relevance in the context of mesoscale phenomena. We emphasize important characteristics that should be included in a mesoscopic model: an explicit continuous three-dimensional space with discrete simplified molecules that still have the characteristics of steric volume exclusion and realistic distant interaction forces. Then we present 3DSPI, a model dedicated to studies of nuclear bodies based on a simple formalism inspired from molecular dynamics and coarse-grained models: particles interacting through a potential energy function and driven by an overdamped Langevin equation. Finally, we present the features expected to be included in the model, pointing out the difficulties that might arise.

Chaperone function in organic co-solvents: experimental characterization and modeling of a hyperthermophilic chaperone subunit from Methanocaldococcus jannaschii

Molecular chaperones play a central role in maintaining protein structure within a cell. Previously, we determined that the gene encoding a molecular chaperone, a thermosome, from the hyperthermophilic archaeon Methanocaldococcus jannaschii is upregulated upon lethal heat shock. We have recombinantly expressed this thermosome (rTHS) and show here that it is both stable and fully functional in aqueous solutions containing water-miscible organic co-solvents. Based on circular dichroism the secondary structure of rTHS was not affected by one-hour exposures to a variety of co-solvents including 30% v/v acetonitrile (ACN) and 50% methanol (MeOH). By contrast, the secondary structure of a mesophilic homologue, GroEL/GroES (GroE), was substantially disrupted. rTHS reduced the aggregation of ovalbumin and citrate synthase in 30% ACN, assisted refolding of citrate synthase upon solvent-inactivation, and stabilized citrate synthase and glutamate dehydrogenase in the direct presence of co-solvents. Apparent total turnover numbers of these enzymes in denaturing solutions increased by up to 2.5-fold in the presence of rTHS. Mechanistic models are proposed to help ascertain specific conditions that could enhance or limit organic solvent-induced chaperone activity. These models suggest that thermodynamic stability and the reversibility of enzyme unfolding play key roles in the effectiveness of enzyme recovery by rTHS.

Manipulating biopolymer dynamics by anisotropic nanoconfinement

How the geometry of nanosized confinement affects dynamics of biomaterials is interesting yet poorly understood. An elucidation of structural details upon nanosized confinement may benefit manufacturing pharmaceuticals in biomaterial sciences and medicine. The behavior of biopolymers in nanosized confinement is investigated using coarse-grained models and molecular simulations. Particularly, we address the effects of shapes of a confinement on protein-folding dynamics by measuring folding rates and dissecting structural properties of the transition states in nanosized spheres and ellipsoids. We find that when the form of a confinement resembles the geometrical properties of the transition states, the rates of folding kinetics are most enhanced. This knowledge of shape selectivity in identifying optimal conditions for reactions will have a broad impact in nanotechnology and pharmaceutical sciences.

Conformational sampling of peptides in cellular environments

Biological systems provide a complex environment that can be understood in terms of its dielectric properties. High concentrations of macromolecules and cosolvents effectively reduce the dielectric constant of cellular environments, thereby affecting the conformational sampling of biomolecules. To examine this effect in more detail, the conformational preference of alanine dipeptide, poly-alanine, and melittin in different dielectric environments is studied with computer simulations based on recently developed generalized Born methodology. Results from these simulations suggest that extended conformations are favored over alpha-helical conformations at the dipeptide level at and below dielectric constants of 5-10. Furthermore, lower-dielectric environments begin to significantly stabilize helical structures in poly-alanine at epsilon = 20. In the more complex peptide melittin, different dielectric environments shift the equilibrium between two main conformations: a nearly fully extended helix that is most stable in low dielectrics and a compact, V-shaped conformation consisting of two helices that is preferred in higher dielectric environments. An additional conformation is only found to be significantly populated at intermediate dielectric constants. Good agreement with previous studies of different peptides in specific, less-polar solvent environments, suggest that helix stabilization and shifts in conformational preferences in such environments are primarily due to a reduced dielectric environment rather than specific molecular details. The findings presented here make predictions of how peptide sampling may be altered in dense cellular environments with reduced dielectric response.

Folding kinetics of a lattice protein via a forward flux sampling approach

We implement a forward flux sampling approach [R. J. Allen et al., J. Chem. Phys. 124, 194111 (2006)] for calculating transition rate constants and for sampling paths of protein folding events. The algorithm generates trajectories for the transition between the unfolded and folded states as chains of partially connected paths, which can be used to obtain the transition-state ensemble and the properties that characterize these intermediates. We apply this approach to Monte Carlo simulations of a model lattice protein in open space and in confined spaces of varying dimensions. We study the effect of confinement on both protein thermodynamic stability and folding kinetics; the former by mapping free-energy landscapes and the latter by the determination of rate constants and mechanistic details of the folding pathway. Our results show that, for the range of temperatures where the native state is stable, confinement of a protein destabilizes the unfolded state by reducing its entropy, resulting in increased thermodynamic stability of the folded state. Relative to the folding in open space, we find that the kinetics can be accelerated at temperatures above the temperature at which the unconfined protein folds fastest and that the rate constant increases with the number of constrained dimensions. By examining the statistical properties of the transition-state ensemble, we detect signs of a classical nucleation folding mechanism for a core of native contacts formed at an early stage of the process. This nucleus acts as folding foci and is composed of those residues that have higher probability to form native contacts in the transition-state intermediates, which can vary depending on the confinement conditions of the system.

Effects of Crowding and Confinement on the Structures of the Transition State Ensemble in Proteins

Characterization of the structures of the transition state ensemble is a key step in describing the folding reaction. Using two variants of a coarse-grained model of the three-stranded beta-sheet WW domain and a fully automated progress variable clustering (PVC) algorithm, we have dissected the effect of macromolecular crowding and confinement on the changes in the transition state structures in comparison to bulk. Each amino acid is represented using a Calpha atom and a side chain. The distance between the Calpha atom and center of mass of the side chain is taken to be its effective van der Waals radius. For the bulk case, we predict using the PVC algorithm, which does not assume knowledge of the underlying folding reaction coordinate, that there are two classes of structures in the transition state ensemble (TSE). The structures in both of the classes are compact. The dominant cluster is more structured than the structures in the less populated class. In accord with bulk experiments, the residues in strands beta2 and beta3 and the interactions at the beta2-beta3 interface are structured. When only excluded volume interactions between the crowding particles and the WW domain are taken into account or when the protein is confined to an inert spherical pore, the overall structure of the TSE does not change dramatically. However, in this entropy dominated regime, the width of the TSE decreases and the structures become more oblate and less spherical as the volume fraction of crowding particle increases or when the pore radius decreases. It suggests that the shape changes, which are computed using the moment of inertia tensor, may represent the slow degrees of freedom during the folding process. When non-native interactions between side chains and interactions with the cavity of the pores are taken into account, the TSE becomes considerably broader. Although the topology in the transition has a fold similar to the native state, the structures are far more plastic than in the bulk. The TSE is sensitive to the size of the pore as well as interactions between the pore and the protein. The differences between the two cases (confinement in an inert pore and when pore-protein interactions are considered) arise due to the increased importance of enthalpic interactions in the transition state as the strength of the protein-pore interaction increases.

Protein folding under confinement: a role for solvent

Although most experimental and theoretical studies of protein folding involve proteins in vitro, the effects of spatial confinement may complicate protein folding in vivo. In this study, we examine the folding dynamics of villin (a small fast folding protein) with explicit solvent confined to an inert nanopore. We have calculated the probability of folding before unfolding (P(fold)) under various confinement regimes. Using P(fold) correlation techniques, we observed two competing effects. Confining protein alone promotes folding by destabilizing the unfolded state. In contrast, confining both protein and solvent gives rise to a solvent-mediated effect that destabilizes the native state. When both protein and solvent are confined we see unfolding to a compact unfolded state different from the unfolded state seen in bulk. Thus, we demonstrate that the confinement of solvent has a significant impact on protein kinetics and thermodynamics. We conclude with a discussion of the implications of these results for folding in confined environments such as the chaperonin cavity in vivo.

Nondestructive mechanical release of ordered polymer microfiber arrays from porous templates

The fabrication of one-dimensional (1D) nanostructures and microstructures inside the pores of porous templates is intensively investigated. The release of these structures is commonly accomplished by etching and destroying the templates. The 1D nanostructures and microstructures tend to condense because of the occurrence of capillary forces during drying of the specimens. It is shown that highly ordered arrays of polymer microfibers can be easily detached from silanized porous templates by mechanical lift-off. This procedure leaves the templates intact, thus allowing their recycling, and does not involve the use of solutions or solvents, thus circumventing condensation. Therefore, mechanical lift-off may enable the up-scaling of template-based approaches to the fabrication of highly ordered assemblies of 1D nanostructures and microstructures.

Coupling between allosteric transitions in GroEL and assisted folding of a substrate protein

Escherichia coli chaperonin, GroEL, helps proteins fold under nonpermissive conditions. During the reaction cycle, GroEL undergoes allosteric transitions in response to binding of a substrate protein (SP), ATP, and the cochaperonin GroES. Using coarse-grained representations of the GroEL and GroES structures, we explore the link between allosteric transitions and the folding of a model SP, a de novo-designed four-helix bundle protein, with low spontaneous yield. The ensemble of GroEL-bound SP is less structured than the bulk misfolded structures. Upon binding, which kinetically occurs in two stages, the SP loses not only native tertiary contacts but also experiences a decrease in helical content. During multivalent binding and the subsequent ATP-driven transition of GroEL the SP undergoes force-induced stretching. Upon encapsulation, which occurs upon GroES binding, the SP finds itself in a "hydrophilic" cavity in which it can reach the folded conformation. Surprisingly, we find that the yield of the native state in the expanded GroEL cavity is relatively small even after it remains in it for twice the spontaneous folding time. Thus, in accord with the iterative annealing mechanism, multiple rounds of binding, partial unfolding, and release of the SP are required to enhance the yield of the folded SP.

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.

Hydrophobic and Ionic Interactions in Nanosized Water Droplets

A number of situations such as protein folding in confined spaces, lubrication in tight spaces, and chemical reactions in confined spaces require an understanding of water-mediated interactions. As an illustration of the profound effects of confinement on hydrophobic and ionic interactions, we investigate the solvation of methane and methane decorated with charges in spherically confined water droplets. Free energy profiles for a single methane molecule in droplets, ranging in diameter (D) from 1 to 4 nm, show that the droplet surfaces are strongly favorable as compared to the interior. From the temperature dependence of the free energy in D = 3 nm, we show that this effect is entropically driven. The potentials of mean force (PMFs) between two methane molecules show that the solvent separated minimum in the bulk is completely absent in confined water, independent of the droplet size since the solute particles are primarily associated with the droplet surface. The tendency of methanes with charges (M(q+) and M(q-) with q(+) = |q(-)| = 0.4e, where e is the electronic charge) to be pinned at the surface depends dramatically on the size of the water droplet. When D = 4 nm, the ions prefer the interior whereas for D < 4 nm the ions are localized at the surface, but with much less tendency than for methanes. Increasing the ion charge to e makes the surface strongly unfavorable. Reflecting the charge asymmetry of the water molecule, negative ions have a stronger preference for the surface compared to positive ions of the same charge magnitude. With increasing droplet size, the PMFs between M(q+) and M(q-) show decreasing influence of the boundary owing to the reduced tendency for surface solvation. We also show that as the solute charge density decreases the surface becomes less unfavorable. The implications of our results for the folding of proteins in confined spaces are outlined.

Computer simulation of polypeptides in a confinement

A coarse-grained model of polypeptide chains confined in a slit formed by two parallel impenetrable surfaces was studied. The chains were flexible heteropolymers (polypeptides) built of two kinds of united atoms-hydrophobic and hydrophilic. The positions of the united atoms were restricted to the vertices of a [310] lattice. The force field consisted of a rigorous excluded volume, a long-distance potential between a pair of amino-acid residues and a local preference for forming secondary structure (helices). The properties of the chains were studied at a wide range of temperatures from good to bad solvent conditions. Monte-Carlo simulations were carried out using the algorithm based on the chain's local changes of conformation and employing the Replica Exchange technique. The influence of the chain length, the distances between the confining surfaces, the temperature and the force field on the dimension and the structure of chains were studied. It was shown that the presence of the confinement chain complicates the process of the chain collapse to low-temperature structures. For some conditions, one can find a rapid decrease of chain size and a second transition indicated by the rapid decrease of the total energy of the system.

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.

Molecular simulations of cotranslational protein folding: fragment stabilities, folding cooperativity, and trapping in the ribosome

Although molecular simulation methods have yielded valuable insights into mechanistic aspects of protein refolding in vitro, they have up to now not been used to model the folding of proteins as they are actually synthesized by the ribosome. To address this issue, we report here simulation studies of three model proteins: chymotrypsin inhibitor 2 (CI2), barnase, and Semliki forest virus protein (SFVP), and directly compare their folding during ribosome-mediated synthesis with their refolding from random, denatured conformations. To calibrate the methodology, simulations are first compared with in vitro data on the folding stabilities of N-terminal fragments of CI2 and barnase; the simulations reproduce the fact that both the stability and thermal folding cooperativity increase as fragments increase in length. Coupled simulations of synthesis and folding for the same two proteins are then described, showing that both fold essentially post-translationally, with mechanisms effectively identical to those for refolding. In both cases, confinement of the nascent polypeptide chain within the ribosome tunnel does not appear to promote significant formation of native structure during synthesis; there are however clear indications that the formation of structure within the nascent chain is sensitive to location within the ribosome tunnel, being subject to both gain and loss as the chain lengthens. Interestingly, simulations in which CI2 is artificially stabilized show a pronounced tendency to become trapped within the tunnel in partially folded conformations: non-cooperative folding, therefore, appears in the simulations to exert a detrimental effect on the rate at which fully folded conformations are formed. Finally, simulations of the two-domain protease module of SFVP, which experimentally folds cotranslationally, indicate that for multi-domain proteins, ribosome-mediated folding may follow different pathways from those taken during refolding. Taken together, these studies provide a first step toward developing more realistic methods for simulating protein folding as it occurs in vivo.

Structural features of the GroEL-GroES nano-cage required for rapid folding of encapsulated protein

GroEL and GroES form a chaperonin nano-cage for proteins up to approximately 60 kDa to fold in isolation. Here we explored the structural features of the chaperonin cage critical for rapid folding of encapsulated substrates. Modulating the volume of the GroEL central cavity affected folding speed in accordance with confinement theory. Small proteins (approximately 30 kDa) folded more rapidly as the size of the cage was gradually reduced to a point where restriction in space slowed folding dramatically. For larger proteins (approximately 40-50 kDa), either expanding or reducing cage volume decelerated folding. Additionally, interactions with the C-terminal, mildly hydrophobic Gly-Gly-Met repeat sequences of GroEL protruding into the cavity, and repulsion effects from the negatively charged cavity wall were required for rapid folding of some proteins. We suggest that by combining these features, the chaperonin cage provides a physical environment optimized to catalyze the structural annealing of proteins with kinetically complex folding pathways.

Protein translocation through a tunnel induces changes in folding kinetics: a lattice model study

Compaction of a nascent polypeptide chain inside the ribosomal exit tunnel, before it leaves the ribosome, has been proposed to accelerate the folding of newly synthesized proteins following their release from the ribosome. Thus, we used Kinetic Monte Carlo simulations of a minimalist on-lattice model to explore the effect that polypeptide translocation through a variety of channels has on protein folding kinetics. Our results demonstrate that tunnel confinement promotes faster folding of a well-designed protein relative to its folding in free space by displacing the unfolded state towards more compact structures that are closer to the transition state. Since the tunnel only forbids rarely visited, extended configurations, it has little effect on a "poorly designed" protein whose unfolded state is largely composed of low-energy, compact, misfolded configurations. The beneficial effect of the tunnel depends on its width; for example, a too-narrow tunnel enforces unfolded states with limited or no access to the transition state, while a too-wide tunnel has no effect on the unfolded state entropy. We speculate that such effects are likely to play an important role in the folding of some proteins or protein domains in the cellular environment and might dictate whether a protein folds co-translationally or post-translationally.

Molecular simulations of cotranslational protein folding: fragment stabilities, folding cooperativity, and trapping in the ribosome

Although molecular simulation methods have yielded valuable insights into mechanistic aspects of protein refolding in vitro, they have up to now not been used to model the folding of proteins as they are actually synthesized by the ribosome. To address this issue, we report here simulation studies of three model proteins: chymotrypsin inhibitor 2 (CI2), barnase, and Semliki forest virus protein (SFVP), and directly compare their folding during ribosome-mediated synthesis with their refolding from random, denatured conformations. To calibrate the methodology, simulations are first compared with in vitro data on the folding stabilities of N-terminal fragments of CI2 and barnase; the simulations reproduce the fact that both the stability and thermal folding cooperativity increase as fragments increase in length. Coupled simulations of synthesis and folding for the same two proteins are then described, showing that both fold essentially post-translationally, with mechanisms effectively identical to those for refolding. In both cases, confinement of the nascent polypeptide chain within the ribosome tunnel does not appear to promote significant formation of native structure during synthesis; there are however clear indications that the formation of structure within the nascent chain is sensitive to location within the ribosome tunnel, being subject to both gain and loss as the chain lengthens. Interestingly, simulations in which CI2 is artificially stabilized show a pronounced tendency to become trapped within the tunnel in partially folded conformations: non-cooperative folding, therefore, appears in the simulations to exert a detrimental effect on the rate at which fully folded conformations are formed. Finally, simulations of the two-domain protease module of SFVP, which experimentally folds cotranslationally, indicate that for multi-domain proteins, ribosome-mediated folding may follow different pathways from those taken during refolding. Taken together, these studies provide a first step toward developing more realistic methods for simulating protein folding as it occurs in vivo.

Nanotube Confinement Denatures Protein Helices

In striking contrast to simple polymer physics theory, which does not account for solvent effects, we find that physical confinement of solvated biopolymers decreases solvent entropy, which in turn leads to a reduction in the organized structural content of the polymer. Since our theory is based on a fundamental property of water-protein statistical mechanics, we expect it to have broad implications in many biological and material science contexts.

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.

Translocation boost protein-folding efficiency of double-barreled chaperonins

Incorrect folding of proteins in living cells may lead to malfunctioning of the cell machinery. To prevent such cellular disasters from happening, all cells contain molecular chaperones that assist nonnative proteins in folding into the correct native structure. One of the most studied chaperone complexes is the GroEL-GroES complex. The GroEL part has a "double-barrel" structure, which consists of two cylindrical chambers joined at the bottom in a symmetrical fashion. The hydrophobic rim of one of the GroEL chambers captures nonnative proteins. The GroES part acts as a lid that temporarily closes the filled chamber during the folding process. Several capture-folding-release cycles are required before the nonnative protein reaches its native state. Here we report molecular simulations that suggest that translocation of the nonnative protein through the equatorial plane of the complex boosts the efficiency of the chaperonin action. If the target protein is correctly folded after translocation, it is released. However, if it is still nonnative, it is likely to remain trapped in the second chamber, which then closes to start a reverse translocation process. This shuttling back and forth continues until the protein is correctly folded. Our model provides a natural explanation for the prevalence of double-barreled chaperonins. Moreover, we argue that internal folding is both more efficient and safer than a scenario where partially refolded proteins escape from the complex before being recaptured.

Structural transitions of confined model proteins: molecular dynamics simulation and experimental validation

Proteins fold in a confined space not only in vivo, i.e., folding assisted by molecular chaperons and chaperonins in a crowded cellular medium, but also in vitro as in production of recombinant proteins. Despite extensive work on protein folding in bulk, little is known about how and to what extent the thermodynamics and kinetics of protein folding are altered by confinement. In this work, we use a Gō-like off-lattice model to investigate the folding and stability of an all beta-sheet protein in spherical cages of different sizes and surface hydrophobicity. We find whereas extreme confinement inhibits correct folding, a hydrophilic cage stabilizes the protein due to restriction of the unfolded configurations. In a hydrophobic cage, however, strong attraction from the cage surface destabilizes the confined protein because of competition between self-aggregation and adsorption of hydrophobic residues. We show that the kinetics of protein collapse and folding is strongly correlated with both the cage size and the surface hydrophobicity. It is demonstrated that a cage of moderate size and hydrophobicity optimizes both the folding yield and kinetics of structural transitions. To support the simulation results, we have also investigated the refolding of hen-egg lysozyme in the presence of cetyltrimethylammoniumbromide (CTAB) surfactants that provide an effective confinement of the proteins by micellization. The influence of the surfactant hydrophobicity on the structural and biological activity of the protein is determined with circular dichroism spectrum, fluorescence emission spectrum, and biological activity assay. It is shown that, as predicted by coarse-grained simulations, CTAB micelles facilitate the collapse of denatured lysozyme, whereas the addition of beta-cyclodextrin-grafted-PNIPAAm, a weakly hydrophobic stripper, dissociates CTAB micelles and promotes the conformational rearrangement and thereby gives an improved recovery of lysozyme activity.