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相似文献   

Autonomously Folding Protein Fragments Reveal Differences in the Energy Landscapes of Homologous RNases H

An important approach to understanding how a protein sequence encodes its energy landscape is to compare proteins with different sequences that fold to the same general native structure. In this work, we compare E. coli and T. thermophilus homologs of the protein RNase H. Using protein fragments, we create equilibrium mimics of two different potential partially-folded intermediates (I_core and I_core+1) hypothesized to be present on the energy landscapes of these two proteins. We observe that both T. thermophilus RNase H (ttRNH) fragments are folded and have distinct stabilities, indicating that both regions are capable of autonomous folding and that both intermediates are present as local minima on the ttRNH energy landscape. In contrast, the two E. coli RNase H (ecRNH) fragments have very similar stabilities, suggesting that the presence of additional residues in the I_core+1 fragment does not affect the folding or structure as compared to I_core. NMR experiments provide additional evidence that only the I_core intermediate is populated by ecRNH. This is one of the biggest differences that has been observed between the energy landscapes of these two proteins. Additionally, we used a FRET experiment in the background of full-length ttRNH to specifically monitor the formation of the I_core+1 intermediate. We determine that the ttRNH I_core+1 intermediate is likely the intermediate populated prior to the rate-limiting barrier to global folding, in contrast to E. coli RNase H for which I_core is the folding intermediate. This result provides new insight into the nature of the rate-limiting barrier for the folding of RNase H.     

Slow and Bimolecular Folding of a De Novo Designed Monomeric Protein DS119

De novo protein design offers a unique means to test and advance our understanding of how proteins fold. However, most current design methods are native structure eccentric and folding kinetics has rarely been considered in the design process. Here, we show that a de novo designed mini-protein DS119, which folds into a βαβ structure, exhibits unusually slow and concentration-dependent folding kinetics. For example, the folding time for 50 μM of DS119 was estimated to be ∼2 s. Stopped-flow fluorescence resonance energy transfer experiments further suggested that its folding was likely facilitated by a transient dimerization process. Taken together, these results highlight the need for consideration of the entire folding energy landscape in de novo protein design and provide evidence suggesting nonnative interactions can play a key role in protein folding.     

Förster distances for fluorescence resonant energy transfer between mCherry and other visible fluorescent proteins

We present, for the red fluorescent protein mCherry acting as both fluorescence resonant energy transfer (FRET) donor and acceptor, Förster critical distance (r₀) values with five important visible fluorescent protein (VFP) variants as well as with itself. The pair EYFP-mCherry exhibits an r₀ of 5.66 nm, equaling or exceeding any combination of VFPs reported previously. Moreover, mCherry should be an excellent chromophore for homo-FRET with an r₀ of 5.10 nm for energy transfer between two mCherry moieties. Finally, mCherry exhibits higher r₀ values than does DsRed. These characteristics, combined with mCherry’s rapid folding and excellent spectral properties, suggest that mCherry constitutes a valuable long-wavelength hetero-FRET acceptor and probe for homo-FRET experiments.     

Effects of ionic strength on the folding and stability of SAMP1, a ubiquitin-like halophilic protein

Our knowledge of the folding behavior of proteins from extremophiles is limited at this time. These proteins may more closely resemble the primordial proteins selected in early evolution under extreme conditions. The small archaeal modifier protein 1 (SAMP1) studied in this report is an 87-residue protein with a β-grasp fold found in the halophile Haloferax volcanii from the Dead Sea. To gain insight into the effects of salt on the stability and folding mechanism of SAMP1, we conducted equilibrium and kinetic folding experiments as a function of sodium chloride concentration. The results revealed that increasing ionic strength accelerates refolding and slows down unfolding of SAMP1, giving rise to a pronounced salt-induced stabilization. With increasing NaCl concentration, the rate of folding observed via a combination of continuous-flow (0.1–2 ms time range) and stopped-flow measurements (>2 ms) exhibited a >100-fold increase between 0.1 and 1.5 M NaCl and leveled off at higher concentrations. Using the Linderström-Lang smeared charge formalism to model electrostatic interactions in ground and transition states encountered during folding, we showed that the observed salt dependence is dominated by Debye-Hückel screening of electrostatic repulsion among numerous negatively charged residues. Comparisons are also drawn with three well-studied mesophilic members of the β-grasp superfamily: protein G, protein L, and ubiquitin. Interestingly, the folding rate of SAMP1 in 3 M sodium chloride is comparable to that of protein G, ubiquitin, and protein L at lower ionic strength. The results indicate the important role of electrostatic interactions in protein folding and imply that proteins have evolved to minimize unfavorable charge-charge interactions under their specific native conditions.     

Finding the lowest free energy conformation of a protein is an NP-hard problem: Proof and implications

The protein folding problem and the notion of NP-completeness and NP-hardness are discussed. A lattice model is suggested to capture the essece of protein folding. For this model we present a proof that finding the lowest free energy conformation belongs to the class of NP-hard problems. The implications of the proof are discussed and we suggest that the natural folding process cannot be considered as a search for the global free energy minimum. However, we suggest an explanation as to why, for many proteins, the native functional conformation maycoincide with the lowest free energy conformation.     

On Simplified Global Nonlinear Function for Fitness Landscape: A Case Study of Inverse Protein Folding

The construction of fitness landscape has broad implication in understanding molecular evolution, cellular epigenetic state, and protein structures. We studied the problem of constructing fitness landscape of inverse protein folding or protein design, with the aim to generate amino acid sequences that would fold into an a priori determined structural fold which would enable engineering novel or enhanced biochemistry. For this task, an effective fitness function should allow identification of correct sequences that would fold into the desired structure. In this study, we showed that nonlinear fitness function for protein design can be constructed using a rectangular kernel with a basis set of proteins and decoys chosen a priori. The full landscape for a large number of protein folds can be captured using only 480 native proteins and 3,200 non-protein decoys via a finite Newton method. A blind test of a simplified version of fitness function for sequence design was carried out to discriminate simultaneously 428 native sequences not homologous to any training proteins from 11 million challenging protein-like decoys. This simplified function correctly classified 408 native sequences (20 misclassifications, 95% correct rate), which outperforms several other statistical linear scoring function and optimized linear function. Our results further suggested that for the task of global sequence design of 428 selected proteins, the search space of protein shape and sequence can be effectively parametrized with just about 3,680 carefully chosen basis set of proteins and decoys, and we showed in addition that the overall landscape is not overly sensitive to the specific choice of this set. Our results can be generalized to construct other types of fitness landscape.     

Direct Observation of Parallel Folding Pathways Revealed Using a Symmetric Repeat Protein System

Although progress has been made to determine the native fold of a polypeptide from its primary structure, the diversity of pathways that connect the unfolded and folded states has not been adequately explored. Theoretical and computational studies predict that proteins fold through parallel pathways on funneled energy landscapes, although experimental detection of pathway diversity has been challenging. Here, we exploit the high translational symmetry and the direct length variation afforded by linear repeat proteins to directly detect folding through parallel pathways. By comparing folding rates of consensus ankyrin repeat proteins (CARPs), we find a clear increase in folding rates with increasing size and repeat number, although the size of the transition states (estimated from denaturant sensitivity) remains unchanged. The increase in folding rate with chain length, as opposed to a decrease expected from typical models for globular proteins, is a clear demonstration of parallel pathways. This conclusion is not dependent on extensive curve-fitting or structural perturbation of protein structure. By globally fitting a simple parallel-Ising pathway model, we have directly measured nucleation and propagation rates in protein folding, and have quantified the fluxes along each path, providing a detailed energy landscape for folding. This finding of parallel pathways differs from results from kinetic studies of repeat-proteins composed of sequence-variable repeats, where modest repeat-to-repeat energy variation coalesces folding into a single, dominant channel. Thus, for globular proteins, which have much higher variation in local structure and topology, parallel pathways are expected to be the exception rather than the rule.     

Folding Steps in the Fibrillation of Functional Amyloid: Denaturant Sensitivity Reveals Common Features in Nucleation and Elongation

Functional bacterial amyloids (FuBA) are intrinsically disordered proteins (IDPs) which rapidly and efficiently aggregate, forming extremely stable fibrils. The conversion from IDP to amyloid is evolutionarily optimized and likely couples folding to association. Many FuBA contain several imperfect repeat sequences which contribute to the stability of mature FuBA fibrils. Aggregation can be considered an intermolecular extension of the process of intramolecular protein folding which has traditionally been studied using chemical denaturants. Here we employ denaturants to investigate folding steps during fibrillation of CsgA and FapC. We quantify protein compactification (i.e. the extent of burial of otherwise exposed surface area upon association of proteins) during different stages of fibrillation based on the dependence of fibrillation rate constants on the denaturant concentration (m-values) determined from fibrillation curves. For both proteins, urea mainly affects nucleation and elongation (not fragmentation), consistent with the fact that these steps involve both intra- and intermolecular association. The two steps have similar m-values, indicating that activation steps in nucleation and elongation involve the same level of folding. Surprisingly, deletion of two or three repeats from FapC leads to larger m-values (i.e. higher compactification) during the activation step of fibril growth. This observation is extended by SAXS analysis of the fibrils which indicates that weakening of the amyloidogenic core caused by repeat deletions causes a larger portion of normally unstructured regions of the protein to be included into the amyloid backbone. We conclude that the sensitivity of fibrillation to denaturants can provide useful insight into molecular mechanisms of aggregation.     

Slow and Bimolecular Folding of a De Novo Designed Monomeric Protein DS119

De novo protein design offers a unique means to test and advance our understanding of how proteins fold. However, most current design methods are native structure eccentric and folding kinetics has rarely been considered in the design process. Here, we show that a de novo designed mini-protein DS119, which folds into a βαβ structure, exhibits unusually slow and concentration-dependent folding kinetics. For example, the folding time for 50 μM of DS119 was estimated to be ∼2 s. Stopped-flow fluorescence resonance energy transfer experiments further suggested that its folding was likely facilitated by a transient dimerization process. Taken together, these results highlight the need for consideration of the entire folding energy landscape in de novo protein design and provide evidence suggesting nonnative interactions can play a key role in protein folding.     

Twin-Arginine Translocation of Active Human Tissue Plasminogen Activator in Escherichia coli

When eukaryotic proteins with multiple disulfide bonds are expressed at high levels in Escherichia coli, the efficiency of thiol oxidation and isomerization is typically not sufficient to yield soluble products with native structures. Even when such proteins are secreted into the oxidizing periplasm or expressed in the cytoplasm of cells carrying mutations in the major intracellular disulfide bond reduction systems (e.g., trxB gor mutants), correct folding can be problematic unless a folding modulator is simultaneously coexpressed. In the present study we explored whether the bacterial twin-arginine translocation (Tat) pathway could serve as an alternative expression system for obtaining appreciable levels of recombinant proteins which exhibit complex patterns of disulfide bond formation, such as full-length human tissue plasminogen activator (tPA) (17 disulfides) and a truncated but enzymatically active version of tPA containing nine disulfides (vtPA). Remarkably, targeting of both tPA and vtPA to the Tat pathway resulted in active protein in the periplasmic space. We show here that export by the Tat translocator is dependent upon oxidative protein folding in the cytoplasm of trxB gor cells prior to transport. Whereas previous efforts to produce high levels of active tPA or vtPA in E. coli required coexpression of the disulfide bond isomerase DsbC, we observed that Tat-targeted vtPA and tPA reach a native conformation without thiol-disulfide oxidoreductase coexpression. These results demonstrate that the Tat system may have inherent and unexpected benefits compared with existing expression strategies, making it a viable alternative for biotechnology applications that hinge on protein expression and secretion.     

Measuring “Unmeasurable” Folding Kinetics of Proteins by Single-Molecule Force Spectroscopy

Folding and unfolding are fundamental biological processes in cell and are important for the biological functions of proteins. Characterizing the folding and unfolding kinetics of proteins is important for understanding the energetic landscape leading to the active native conformations of these molecules. However, the thermal or chemical-induced unfolding of many proteins is irreversible in vitro, precluding characterization of the folding kinetics of such proteins, just as it is impossible to “un-boil” an egg. Irreversible unfolding often manifests as irreversible aggregation of unfolded polypeptide chains, which typically occurs between denatured protein molecules in response to the exposure of hydrophobic residues to solvent. An example of such a protein where thermal denaturation results in irreversible aggregation is the β-1,4 endoxylanase from Bacillus circulans (BCX). Here, we report the use of single-molecule atomic force microscopy to directly measure the folding kinetics of BCX in vitro. By mechanically unfolding BCX, we essentially allowed only one unfolded molecule to exist in solution at a given time, effectively eliminating the possibility for aggregation. We found that BCX can readily refold back to the native state, allowing us to measure its folding kinetics for the first time. Our results demonstrate that single-molecule force-spectroscopy-based methods can adequately tackle the challenge of “un-boiling eggs”, providing a general methodology to characterize the folding kinetics of many proteins that suffer from irreversible denaturation and thus cannot be characterized using traditional equilibrium methodologies.     

Toward resolution of ambiguity for the unfolded state

The unfolded states in proteins and nucleic acids remain weakly understood despite their importance in folding processes; misfolding diseases (Parkinson's and Alzheimer's); natively unfolded proteins (as many as 30% of eukaryotic proteins, according to Fink); and the study of ribozymes. Research has been hindered by the inability to quantify the residual (native) structure present in an unfolded protein or nucleic acid. Here, a scaling model is proposed to quantify the molar degree of folding and the unfolded state. The model takes a global view of protein structure and can be applied to a number of analytic methods and to simulations. Three examples are given of application to small-angle scattering from pressure-induced unfolding of SNase, from acid-unfolded cytochrome c, and from folding of Azoarcus ribozyme. These examples quantitatively show three characteristic unfolded states for proteins, the statistical nature of a protein folding pathway, and the relationship between extent of folding and chain size during folding for charge-driven folding in RNA.     

Role of Protein Stabilizers on the Conformation of the Unfolded State of Cytochrome c and Its Early Folding Kinetics: INVESTIGATION AT SINGLE MOLECULAR RESOLUTION*

An insight into the conformation and dynamics of unfolded and early intermediate states of a protein is essential to understand the mechanism of its aggregation and to design potent inhibitor molecules. Fluorescence correlation spectroscopy has been used to study the effects of several model protein stabilizers on the conformation of the unfolded state and early folding dynamics of tetramethyl rhodamine-labeled cytochrome c from Saccharomyces cerevisiae at single molecular resolution. Special attention has been given to arginine, which is a widely used stabilizer for improving refolding yield of different proteins. The value of the hydrodynamic radius (r_H) obtained by analyzing the intensity fluctuations of the diffusing molecules has been found to increase in a two-state manner as the protein is unfolded by urea. The results further show that the presence of arginine and other protein stabilizers favors a relatively structured conformation of the unfolded states (r_H of 29 Å) over an extended one (r_H of 40 Å), which forms in their absence. Also, the time constant of a kinetic component (τ_R) of about 30 μs has been observed by analyzing the correlation functions, which represents formation of a collapsed state. This time constant varies with urea concentration representing an inverted Chevron plot that shows a roll-over and behavior in the absence of arginine. To the best of our knowledge, this is one of the first applications of fluorescence correlation spectroscopy to study direct folding kinetics of a protein.     

Transition states in protein folding kinetics: modeling phi-values of small beta-sheet proteins

Small single-domain proteins often exhibit only a single free-energy barrier, or transition state, between the denatured and the native state. The folding kinetics of these proteins is usually explored via mutational analysis. A central question is which structural information on the transition state can be derived from the mutational data. In this article, we model and structurally interpret mutational Φ-values for two small β-sheet proteins, the PIN and the FBP WW domains. The native structure of these WW domains comprises two β-hairpins that form a three-stranded β-sheet. In our model, we assume that the transition state consists of two conformations in which either one of the hairpins is formed. Such a transition state has been recently observed in molecular dynamics folding-unfolding simulations of a small designed three-stranded β-sheet protein. We obtain good agreement with the experimental data 1), by splitting up the mutation-induced free-energy changes into terms for the two hairpins and for the small hydrophobic core of the proteins; and 2), by fitting a single parameter, the relative degree to which hairpins 1 and 2 are formed in the transition state. The model helps us to understand how mutations affect the folding kinetics of WW domains, and captures also negative Φ-values that have been difficult to interpret.     

Chaperones GroEL/GroES Accelerate the Refolding of a Multidomain Protein through Modulating On-pathway Intermediates

Despite a vast amount information on the interplay of GroEL, GroES, and ATP in chaperone-assisted folding, the molecular details on the conformational dynamics of folding polypeptide during its GroEL/GroES-assisted folding cycle is quite limited. Practically no such studies have been reported to date on large proteins, which often have difficulty folding in vitro. The effect of the GroEL/GroES chaperonin system on the folding pathway of an 82-kDa slow folding protein, malate synthase G (MSG), was investigated. GroEL bound to the burst phase intermediate of MSG and accelerated the slowest kinetic phase associated with the formation of native topology in the spontaneous folding pathway. GroEL slowly induced conformational changes on the bound burst phase intermediate, which was then transformed into a more folding-compatible form. Subsequent addition of ATP or GroES/ATP to the GroEL-MSG complex led to the formation of the native state via a compact intermediate with the rate several times faster than that of spontaneous refolding. The presence of GroES doubled the ATP-dependent reactivation rate of bound MSG by preventing multiple cycles of its GroEL binding and release. Because GroES bound to the trans side of GroEL-MSG complex, it may be anticipated that confinement of the substrate underneath the co-chaperone is not required for accelerating the rate in the assisted folding pathway. The potential role of GroEL/GroES in assisted folding is most likely to modulate the conformation of MSG intermediates that can fold faster and thereby eliminate the possibility of partial aggregation caused by the slow folding intermediates during its spontaneous refolding pathway.     

Observation of Continuous Contraction and a Metastable Misfolded State during the Collapse and Folding of a Small Protein

To obtain proper insight into how structure develops during a protein folding reaction, it is necessary to understand the nature and mechanism of the polypeptide chain collapse reaction, which marks the initiation of folding. Here, the time-resolved fluorescence resonance energy transfer technique, in which the decay of the fluorescence light intensity with time is used to determine the time evolution of the distribution of intra-molecular distances, has been utilized to study the folding of the small protein, monellin. It is seen that when folding begins, about one-third of the protein molecules collapse into a molten globule state (I_MG), from which they relax by continuous further contraction to transit to the native state. The larger fraction gets trapped into a metastable misfolded state. Exit from this metastable state occurs via collapse to the lower free energy I_MG state. This exit is slow, on a time scale of seconds, because of activation energy barriers. The trapped misfolded molecules as well as the I_MG molecules contract continuously and slowly as structure develops. A phenomenological model of Markovian evolution of the polymer chain undergoing folding, incorporating these features, has been developed, which fits well the experimentally observed time evolution of distance distributions. The observation that the “wrong turn” to the misfolded state has not been eliminated by evolution belies the common belief that the folding of functional protein sequences is very different from that of a random heteropolymer, and the former have been selected by evolution to fold quickly.     

How Difficult Is It to Fold a Knotted Protein? In Silico Insights from Surface-Tethered Folding Experiments

We explore the effect of surface tethering on the folding process of a lattice protein that contains a trefoil knot in its native structure via Monte Carlo simulations. We show that the outcome of the tethering experiment depends critically on which terminus is used to link the protein to a chemically inert plane. In particular, if surface tethering occurs at the bead that is closer to the knotted core the folding rate becomes exceedingly slow and the protein is not able to find the native structure in all the attempted folding trajectories. Such low folding efficiency is also apparent from the analysis of the probability of knot formation, p_knot, as a function of nativeness. Indeed, p_knot increases abruptly from ∼0 to ∼1 only when the protein has more than 80% of its native contacts formed, showing that a highly compact conformation must undergo substantial structural re-arrangement in order to get effectively knotted. When the protein is surface tethered by the bead that is placed more far away from the knotted core p_knot is higher than in the other folding setups (including folding in the bulk), especially if conformations are highly native-like. These results show that the mobility of the terminus closest to the knotted core is critical for successful folding of trefoil proteins, which, in turn, highlights the importance of a knotting mechanism that is based on a threading movement of this terminus through a knotting loop. The results reported here predict that if this movement is blocked, knotting occurs via an alternative mechanism, the so-called spindle mechanism, which is prone to misfolding. Our simulations show that in the three considered folding setups the formation of the knot is typically a late event in the folding process. We discuss the implications of our findings for co-translational folding of knotted trefoils.