Acceleration of the folding of acylphosphatase by stabilization of local secondary structure

The addition of trifluoroethanol or hexafluoroisopropanol converts the apparent two-state folding of acylphosphatase, a small alpha/beta protein, into a multistate mechanism where secondary structure accumulates significantly in the denatured state before folding to the native state. This results in a marked acceleration of folding as revealed by following the intrinsic fluorescence and circular dichroism changes upon folding. The folding rate is at a maximum when the secondary-structure content of the denatured state corresponds to that of the native state, while further stabilization of secondary structure decreases the folding rate. These findings indicate that stabilization of intermediate structure can either enhance or retard folding depending on its nature and content of native-like interactions.     

Stabilisation of alpha-helices by site-directed mutagenesis reveals the importance of secondary structure in the transition state for acylphosphatase folding

The effects of stabilising mutations on the folding process of common-type acylphosphatase have been investigated. The mutations were designed to increase the helical propensity of the regions of the polypeptide chain corresponding to the two alpha-helices of the native protein. Various synthetic peptides incorporating the designed mutations were produced and their helical content estimated by circular dichroism. The most substantial increase in helical content is found for the peptide carrying five mutations in the second alpha-helix. Acylphosphatase variants containing the corresponding mutations display, to different extents, enhanced conformational stabilities as indicated by equilibrium urea denaturation experiments monitored by changes of intrinsic fluorescence. All the protein variants studied here refold with apparent two-state kinetics. Mutations in the first alpha-helix are responsible for a small increase in the refolding rate, accompanied by a marked decrease in the unfolding rate. On the other hand, multiple mutations in the second helix result in a considerable increase in the refolding rate without any significant effect on the unfolding rate. Addition of trifluoroethanol was found to accelerate the folding of the acylphosphatase variants, the extent of the acceleration being inversely proportional to the intrinsic rate of folding of the corresponding mutant. The trifluoroethanol-induced acceleration is far less marked for those variants whose alpha-helical structure is efficiently stabilised by amino acid replacements. This observation suggests that trifluoroethanol acts in a similar manner to the stabilising mutations in promoting native-like secondary structure. Analysis of the kinetic data indicates that the second helix is fully consolidated in the transition state for folding of acylphosphatase, whereas the first helix is only partially formed. These data suggest that the second helix is an important element in the folding process of the protein.     

Contact order revisited: influence of protein size on the folding rate

Guided by the recent success of empirical model predicting the folding rates of small two-state folding proteins from the relative contact order (CO) of their native structures, by a theoretical model of protein folding that predicts that logarithm of the folding rate decreases with the protein chain length L as L(2/3), and by the finding that the folding rates of multistate folding proteins strongly correlate with their sizes and have very bad correlation with CO, we reexamined the dependence of folding rate on CO and L in attempt to find a structural parameter that determines folding rates for the totality of proteins. We show that the Abs_CO = CO x L, is able to predict rather accurately folding rates for both two-state and multistate folding proteins, as well as short peptides, and that this Abs_CO scales with the protein chain length as L(0.70 +/- 0.07) for the totality of studied single-domain proteins and peptides.     

A critical assessment of the topomer search model of protein folding using a continuum explicit-chain model with extensive conformational sampling

Recently, a series of closely related theoretical constructs termed the "topomer search model" (TSM) has been proposed for the folding mechanism of small, single-domain proteins. A basic assumption of the proposed scenarios is that the rate-limiting step in folding is an essentially unbiased, diffusive search for a conformational state called the native topomer defined by an overall native-like topological pattern. Successes in correlating TSM-predicted folding rates with that of real proteins have been interpreted as experimental support for the model. To better delineate the physics entailed, key TSM concepts are examined here using extensive Langevin dynamics simulations of continuum C(alpha) chain models. The theoretical native topomers of four experimentally well-studied two-state proteins are characterized. Consistent with the TSM perspective, we found that the sizes of the native topomers increase with experimental folding rate. However, a careful determination of the corresponding probabilities that the native topomers are populated during a random search fails to reproduce the previously predicted folding rates. Instead, our results indicate that an unbiased TSM search for the native topomer amounts to a Levinthal-like process that would take an impossibly long average time to complete. Furthermore, intraprotein contacts in all four native topomers considered exhibit no apparent correlation with the experimental phi-values determined from the folding kinetics of these proteins. Thus, the present findings suggest that certain basic, generic yet essential energetic features in protein folding are not accounted for by TSM scenarios to date.     

Modulation of the multistate folding of designed TPR proteins through intrinsic and extrinsic factors

Tetratricopeptide repeats (TPRs) are a class of all alpha-helical repeat proteins that are comprised of 34-aa helix-turn-helix motifs. These stack together to form nonglobular structures that are stabilized by short-range interactions from residues close in primary sequence. Unlike globular proteins, they have few, if any, long-range nonlocal stabilizing interactions. Several studies on designed TPR proteins have shown that this modular structure is reflected in their folding, that is, modular multistate folding is observed as opposed to two-state folding. Here we show that TPR multistate folding can be suppressed to approximate two-state folding through modulation of intrinsic stability or extrinsic environmental variables. This modulation was investigated by comparing the thermodynamic unfolding under differing buffer regimes of two distinct series of consensus-designed TPR proteins, which possess different intrinsic stabilities. A total of nine proteins of differing sizes and differing consensus TPR motifs were each thermally and chemically denatured and their unfolding monitored using differential scanning calorimetry (DSC) and CD/fluorescence, respectively. Analyses of both the DSC and chemical denaturation data show that reducing the total stability of each protein and repeat units leads to observable two-state unfolding. These data highlight the intimate link between global and intrinsic repeat stability that governs whether folding proceeds by an observably two-state mechanism, or whether partial unfolding yields stable intermediate structures which retain sufficient stability to be populated at equilibrium.     

Folding barrier in horse cytochrome c: support for a classical folding pathway

Native-state structures and conformations of ferrocytochrome c, nitrosylcytochrome c, and carbonmonoxycytochrome c are very similar. They are, however, immensely different from each other in terms of thermodynamic stability. The dramatic destabilization of ferrocytochrome c to the extent of 12 kcal mol(-1) produces no effect on the folding rate, and this is so in spite of the fact that all three test-tube variants fold in an apparent two-state manner. For all three proteins the folding barrier is early in time, sizable in energy, and is of the same magnitude (approximately 6.5 kcal mol(-1)). These results raise some challenges to the "new view" of protein folding. An early transition state, the search for which consumes most of the observed folding time, is suggested.     

Structure is lost incrementally during the unfolding of barstar

Coincidental equilibrium unfolding transitions observed by multiple structural probes are taken to justify the modeling of protein unfolding as a two-state, N <==> U, cooperative process. However, for many of the large number of proteins that undergo apparently two-state equilibrium unfolding reactions, folding intermediates are detected in kinetic experiments. The small protein barstar is one such protein. Here the two-state model for equilibrium unfolding has been critically evaluated in barstar by estimating the intramolecular distance distribution by time-resolved fluorescence resonance energy transfer (TR-FRET) methods, in which fluorescence decay kinetics are analyzed by the maximum entropy method (MEM). Using a mutant form of barstar containing only Trp 53 as the fluorescence donor and a thionitrobenzoic acid moiety attached to Cys 82 as the fluorescence acceptor, the distance between the donor and acceptor has been shown to increase incrementally with increasing denaturant concentration. Although other probes, such as circular dichroism and fluorescence intensity, suggest that the labeled protein undergoes two-state equilibrium unfolding, the TR-FRET probe clearly indicates multistate equilibrium unfolding. Native protein expands progressively through a continuum of native-like forms that achieve the dimensions of a molten globule, whose heterogeneity increases with increasing denaturant concentration and which appears to be separated from the unfolded ensemble by a free energy barrier.     

Probing the self-association, intermolecular contacts, and folding propensity of amelogenin

Amelogenins are an intrinsically disordered protein family that plays a major role in the development of tooth enamel, one of the most highly mineralized materials in nature. Monomeric porcine amelogenin possesses random coil and residual secondary structures, but it is not known which sequence regions would be conformationally attractive to potential enamel matrix targets such as other amelogenins (self-assembly), other matrix proteins, cell surfaces, or biominerals. To address this further, we investigated recombinant porcine amelogenin (rP172) using "solvent engineering" techniques to simultaneously promote native-like structure and induce amelogenin oligomerization in a manner that allows identification of intermolecular contacts between amelogenin molecules. We discovered that in the presence of 2,2,2-trifluoroethanol (TFE) significant folding transitions and stabilization occurred primarily within the N- and C-termini, while the polyproline Type II central domain was largely resistant to conformational transitions. Seven Pro residues (P2, P127, P130, P139, P154, P157, P162) exhibited conformational response to TFE, and this indicates these Pro residues act as folding enhancers in rP172. The remaining Pro residues resisted TFE perturbations and thus act as conformational stabilizers. We also noted that TFE induced rP172 self-association via the formation of intermolecular contacts involving P4-H6, V19-P33, and E40-T58 regions of the N-terminus. Collectively, these results confirm that the N- and C-termini of amelogenin are conformationally responsive and represent potential interactive sites for amelogenin-target interactions during enamel matrix mineralization. Conversely, the Pro, Gln central domain is resistant to folding and this may have important functional significance for amelogenin.     

Transition state contact orders correlate with protein folding rates

We have used molecular dynamics simulations restrained by experimental phi values derived from protein engineering experiments to determine the structures of the transition state ensembles of ten proteins that fold with two-state kinetics. For each of these proteins we then calculated the average contact order in the transition state ensemble and compared it with the corresponding experimental folding rate. The resulting correlation coefficient is similar to that computed for the contact orders of the native structures, supporting the use of native state contact orders for predicting folding rates. The native contacts in the transition state also correlate with those of the native state but are found to be about 30% lower. These results show that, despite the high levels of heterogeneity in the transition state ensemble, the large majority of contributing structures have native-like topologies and that the native state contact order captures this phenomenon.     

Non-redox-active small-molecules can accelerate oxidative protein folding by novel mechanisms

Multi-disulfide-bond-containing proteins acquire their native structures through an oxidative folding reaction which involves formation of native disulfide bonds through thiol-disulfide exchange reactions between cysteines and disulfides coupled to a conformational folding event. Oxidative folding rates of the four-disulfide-bond-containing protein bovine pancreatic ribonuclease A (RNase A) in the presence of the synthetic redox-active molecule, (+/-)-trans-1,2-bis(2-mercaptoacetamido)cyclohexane (BMC), and in combination with non-redox-active trimethylamine-N-oxide (TMAO), and trifluorethanol were determined by HPLC analysis. The data indicate that regeneration of RNase A is enhanced 2-fold by BMC (50 microM) and 3-fold upon addition of TMAO (0.2 M) and TFE (3% v/v) relative to control experiments performed in the absence of small-molecules. Examination of the native tendency of the fully-reduced polypeptide and the stability of key folding intermediates suggests that the increased oxidative folding rate can be attributed to native-like elements induced within the fully-reduced polypeptide and the stabilization of native-like species by added non-redox-active molecules.     

Critical nucleation size in the folding of small apparently two-state proteins

For apparently two-state proteins, we found that the size (number of folded residues) of a transition state is mostly encoded by the topology, defined by total contact distance (TCD) of the native state, and correlates with its folding rate. This is demonstrated by using a simple procedure to reduce the native structures of the 41 two-state proteins with native TCD as a constraint, and is further supported by analyzing the results of eight proteins from protein engineering studies. These results support the hypothesis that the major rate-limiting process in the folding of small apparently two-state proteins is the search for a critical number of residues with the topology close to that of the native state.     

Trifluoroethanol-induced conformational transitions of proteins: insights gained from the differences between alpha-lactalbumin and ribonuclease A

The trifluoroethanol (TFE)-induced structural changes of two proteins widely used in folding experiments, bovine alpha-lactalbumin, and bovine pancreatic ribonuclease A, have been investigated. The experiments were performed using circular dichroism spectroscopy in the far- and near-UV region to monitor changes in the secondary and tertiary structures, respectively, and dynamic light scattering to measure the hydrodynamic dimensions and the intermolecular interactions of the proteins in different conformational states. Both proteins behave rather differently under the influence of TFE: alpha-lactalbumin exhibits a molten globule state at low TFE concentrations before it reaches the so-called TFE state, whereas ribonuclease A is directly transformed into the TFE state at TFE concentrations above 40% (v/v). The properties of the TFE-induced states are compared with those of equilibrium and kinetic intermediate states known from previous work to rationalize the use of TFE in yielding information about the folding of proteins. Additionally, we report on the properties of TFE/water and TFE/buffer mixtures derived from dynamic light scattering investigations under conditions used in our experiments.     

Contact order dependent protein folding rates: kinetic consequences of a cooperative interplay between favorable nonlocal interactions and local conformational preferences

Physical mechanisms underlying the empirical correlation between relative contact order (CO) and folding rate among naturally occurring small single-domain proteins are investigated by evaluating postulated interaction schemes for a set of three-dimensional 27mer lattice protein models with 97 different CO values. Many-body interactions are constructed such that contact energies become more favorable when short chain segments sequentially adjacent to the contacting residues adopt native-like conformations. At a given interaction strength, this scheme leads to folding rates that are logarithmically well correlated with CO (correlation coefficient r = 0.914) and span more than 2.5 orders of magnitude, whereas folding rates of the corresponding Gō models with additive contact energies have much less logarithmic correlation with CO and span only approximately one order of magnitude. The present protein chain models also exhibit calorimetric cooperativity and linear chevron plots similar to that observed experimentally for proteins with apparent simple two-state folding/unfolding kinetics. Thus, our findings suggest that CO-dependent folding rates of real proteins may arise partly from a significant positive coupling between nonlocal contact favorabilities and local conformational preferences.     

First principles prediction of protein folding rates

Experimental studies have demonstrated that many small, single-domain proteins fold via simple two-state kinetics. We present a first principles approach for predicting these experimentally determined folding rates. Our approach is based on a nucleation-condensation folding mechanism, where the rate-limiting step is a random, diffusive search for the native tertiary topology. To estimate the rates of folding for various proteins via this mechanism, we first determine the probability of randomly sampling a conformation with the native fold topology. Next, we convert these probabilities into folding rates by estimating the rate that a protein samples different topologies during diffusive folding. This topology-sampling rate is calculated using the Einstein diffusion equation in conjunction with an experimentally determined intra-protein diffusion constant. We have applied our prediction method to the 21 topologically distinct small proteins for which two-state rate data is available. For the 18 beta-sheet and mixed alpha-beta native proteins, we predict folding rates within an average factor of 4, even though the experimental rates vary by a factor of approximately 4 x 10(4). Interestingly, the experimental folding rates for the three four-helix bundle proteins are significantly underestimated by this approach, suggesting that proteins with significant helical content may fold by a faster, alternative mechanism. This method can be applied to any protein for which the structure is known and hence can be used to predict the folding rates of many proteins prior to experiment.     

Comparison of the effects of 2,2,2-trifluoroethanol on peptide and protein structure and function

The co-solvent 2,2,2-trifluoroethanol (TFE) has been often used to aid formation of secondary structure in solution peptides or alternately as a denaturant within protein folding studies. Hen egg white lysozyme (HEWL) and a synthetic model peptide defining HEWL helix-4 were used as comparative model systems to systematically investigate the effect of increasing TFE concentrations on the structure of proteins and peptides. HEWL was analyzed using NMR, far-UV CD and fluorescence spectroscopy; with correlation of these results towards changes in enzymatic activity and the helix-4 peptide was analysed using NMR. Data illustrates two conflicting modes of interaction: Low TFE concentrations stabilize tertiary structure, observed from an increase in the number of NMR NOE contacts. Higher TFE concentrations denatured HEWL with the loss of lysozyme tertiary structure. The effects of TFE upon secondary structural elements within HEWL are distinct from those observed for the helix-4 peptide. This illustrates a dissimilar interaction of TFE towards both protein and peptide at equivalent TFE concentrations. The concentration that TFE promotes stabilization over denaturation is likely to be protein dependent although the structural action can be extrapolated to other protein systems with implications for the use of TFE in structural stability studies.     

Does trifluoroethanol affect folding pathways and can it be used as a probe of structure in transition states?

Nonaqueous co-solvents, particularly 2,2,2-trifluoroethanol (TFE), have been used as tools to study protein folding. By analyzing FKBP12, an alpha/beta-protein that folds with two-state kinetics, we have been able to address three key questions concerning the use of TFE. First, does TFE perturb the folding pathway? Second, can the observed changes in the rate of folding and unfolding in TFE be attributed to a change in free energy of a single state? Finally, can TFE be used to infer information on secondary structure formation in the transition state? Protein engineering experiments on FKBP12, coupled with folding and unfolding experiments in 0% and 9.6% TFE, conclusively show that TFE does not perturb the folding pathway of this protein. Our results also suggest that the changes in folding and unfolding rates observed in 9.6% TFE are due to a global effect of TFE on the protein, rather than the stabilization of any elements of secondary structure in the transition state. Thus, studies with TFE and other co-solvents can be accurately interpreted only when combined with other techniques.     

Folding thermodynamics and kinetics of the leucine-rich repeat domain of the virulence factor Internalin B

Although the folding of alpha-helical repeat proteins has been well characterized, much less is known about the folding of repeat proteins containing beta-sheets. Here we investigate the folding thermodynamics and kinetics of the leucine-rich repeat (LRR) domain of Internalin B (InlB), an extracellular virulence factor from the bacterium Lysteria monocytogenes. This domain contains seven tandem leucine-rich repeats, of which each contribute a single beta-strand that forms a continuous beta-sheet with neighboring repeats, and an N-terminal alpha-helical capping motif. Despite its modular structure, InlB folds in an equilibrium two-state manner, as reflected by the identical thermodynamic parameters obtained by monitoring its sigmoidal urea-induced unfolding transition by different spectroscopic probes. Although equilibrium two-state folding is common in alpha-helical repeat proteins, to date, InlB is the only beta-sheet-containing repeat protein for which this behavior is observed. Surprisingly, unlike other repeat proteins exhibiting equilibrium two-state folding, InlB also folds by a simple two-state kinetic mechanism lacking intermediates, aside from the effects of prolyl isomerization on the denatured state. However, like other repeat proteins, InlB also folds significantly more slowly than expected from contact order. When plotted against urea, the rate constants for the fast refolding and single unfolding phases constitute a linear chevron that, when fitted with a kinetic two-state model, yields thermodynamic parameters matching those observed for equilibrium folding. Based on these kinetic parameters, the transition state is estimated to comprise 40% of the total surface area buried upon folding, indicating that a large fraction of the native contacts are formed in the rate-limiting step to folding.