Andrew I. Jewett

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.

Cooperativity, Smooth Energy Landscapes and the Origins of Topology-dependent Protein Folding Rates

The relative folding rates of simple, single-domain proteins, proteins whose folding energy landscapes are smooth, are highly dispersed and strongly correlated with native-state topology. In contrast, the relative folding rates of small, Gō-potential lattice polymers, which also exhibit smooth energy landscapes, are poorly dispersed and insignificantly correlated with native-state topology. Here, we investigate this discrepancy in light of a recent, quantitative theory of two-state folding kinetics, the topomer search model. This model stipulates that the topology-dependence of two-state folding rates is a direct consequence of the extraordinarily cooperative equilibrium folding of simple proteins. We demonstrate that traditional Gō polymers lack the extreme cooperativity that characterizes the folding of naturally occurring, two-state proteins and confirm that the folding rates of a diverse set of Gō 27-mers are poorly dispersed and effectively uncorrelated with native state topology. Upon modestly increasing the cooperativity of the Gō-potential, however, significantly increased dispersion and strongly topology-dependent kinetics are observed. These results support previous arguments that the cooperative folding of simple, single-domain proteins gives rise to their topology-dependent folding rates. We speculate that this cooperativity, and thus, indirectly, the topology-rate relationship, may have arisen in order to generate the smooth energetic landscapes upon which rapid folding can occur.

Sequence periodicity and secondary structure propensity in model proteins

We explore the question of whether local effects (originating from the amino acids intrinsic secondary structure propensities) or nonlocal effects (reflecting the sequence of amino acids as a whole) play a larger role in determining the fold of globular proteins. Earlier circular dichroism studies have shown that the pattern of polar, non polar amino acids (nonlocal effect) dominates over the amino acid intrinsic propensity (local effect) in determining the secondary structure of oligomeric peptides. In this article, we present a coarse grained computational model that allows us to quantitatively estimate the role of local and nonlocal factors in determining both the secondary and tertiary structure of small, globular proteins. The amino acid intrinsic secondary structure propensity is modeled by a dihedral potential term. This dihedral potential is parametrized to match with experimental measurements of secondary structure propensity. Similarly, the magnitude of the attraction between hydrophobic residues is parametrized to match the experimental transfer free energies of hydrophobic amino acids. Under these parametrization conditions, we systematically explore the degree of frustration a given polar, non polar pattern can tolerate when the secondary structure intrinsic propensities are in opposition to it. When the parameters are in the biophysically relevant range, we observe that the fold of small, globular proteins is determined by the pattern of polar, non polar amino acids regardless of their instrinsic secondary structure propensities. Our simulations shed new light on previous observations that tertiary interactions are more influential in determining protein structure than secondary structure propensity. The fact that this can be inferred using a simple polymer model that lacks most of the biochemical details points to the fundamental importance of binary patterning in governing folding.

The effect of surface tethering on the folding of the src-SH3 protein domain

The effect of surface tethering on the folding mechanism of the src-SH3 protein domain was investigated using a coarse-grained Gō-type protein model. The protein was tethered at various locations along the protein chain and the thermodynamics and kinetics of folding were studied using replica exchange and constant temperature Langevin dynamics. Our simulations reveal that tethering in a structured part of the transition state can dramatically alter the folding mechanism, while tethering in an unstructured part leaves the folding mechanism unaltered as compared to bulk folding. Interestingly, there is only modest correlation between the tethering effect on the folding mechanism and its effect on thermodynamic stability and folding rates. We suggest locations on the protein at which tethering could be performed in single-molecule experiments so as to leave the folding mechanism unaltered from the bulk.

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 chaperonins 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.

Relative stability of de novo four–helix bundle proteins: Insights from coarse grained molecular simulations

We use a recently developed coarse‐grained computational model to investigate the relative stability of two different sets of de novo designed four–helix bundle proteins. Our simulations suggest a possible explanation for the experimentally observed increase in stability of the four–helix bundles with increasing sequence length. In details, we show that both short subsequences composed only by polar residues and additional nonpolar residues inserted, via different point mutations in ad hoc positions, seem to play a significant role in stabilizing the four–helix bundle conformation in the longer sequences. Finally, we propose an additional mutation that rescues a short amino acid sequence that would otherwise adopt a compact misfolded state. Our work suggests that simple computational models can be used as a complementary tool in the design process of de novo proteins.