Thermodynamics and kinetics of folding of common-type acylphosphatase: comparison to the highly homologous muscle isoenzyme

A short structural extension dictates the early stages of folding of a PDZ domain

PDZ domains are highly abundant protein-protein interaction modules in human. One of the most extensively characterized PDZ domain, the third PDZ domain from PSD-95 (PDZ3), contains an α-helical C-terminal extension that has a key role in the function of the domain. Here we compared the folding of PDZ3 with a truncated variant (PDZ3Δα3), lacking the additional helix, by means of the so-called Φ-value analysis, an experimental technique that allows inferring the structure of folding transition states. Experiments reveal subtle but detectable differences in the folding of PDZ3Δα3 versus PDZ3, as probed by structural characterization of the folding transition states. These differences appear more remarkable in the early stages of folding, with a detectable shift of the folding nucleus. The presented results allow demonstrating that the native state exerts a weak bias at the early stages of folding, which appear to be characterized by alternative pathways.

Free-energy landscape of two-state protein acylphosphatase with large contact order revealed by force-dependent folding and unfolding dynamics

Acylphosphatase (AcP) is a small protein with 98 amino acid residues that catalyzes the hydrolysis of carboxyl-phosphate bonds. AcP is a typical two-state protein with slow folding rate due to its relatively large contact order in the native structure. The mechanical properties and unfolding behavior of AcP has been studied by atomic force microscope. Here using stable magnetic tweezers, we measured the force-dependent folding rates within a force range 1-3 pN, and unfolding rates 15-40 pN. The obtained unfolding rates show different force sensitivities at forces below and above ∼27 pN, which determines a free-energy landscape with two energy barriers. Our results indicate that the free-energy landscape of small globule proteins have general Bactrian camel shape, and large contact order of the native state produces a high barrier dominate at low forces.

Protein folding in vitro and in the cell: From a solitary journey to a team effort

Correct protein folding is essential for the health and function of living organisms. Yet, it is not well understood how unfolded proteins reach their native state and avoid aggregation, especially within the cellular milieu. Some proteins, especially small, single-domain and apparent two-state folders, successfully attain their native state upon dilution from denaturant. Yet, many more proteins undergo misfolding and aggregation during this process, in a concentration-dependent fashion. Once formed, native and aggregated states are often kinetically trapped relative to each other. Hence, the early stages of protein life are absolutely critical for proper kinetic channeling to the folded state and for long-term solubility and function. This review summarizes current knowledge on protein folding/aggregation mechanisms in buffered solution and within the bacterial cell, highlighting early stages. Remarkably, teamwork between nascent chain, ribosome, trigger factor and Hsp70 molecular chaperones enables all proteins to overcome aggregation propensities and reach a long-lived bioactive state.

PFDB: A standardized protein folding database with temperature correction

We constructed a standardized protein folding kinetics database (PFDB) in which the logarithmic rate constants of all listed proteins are calculated at the standard temperature (25 °C). A temperature correction based on the Eyring–Kramers equation was introduced for proteins whose folding kinetics were originally measured at temperatures other than 25 °C. We verified the temperature correction by comparing the logarithmic rate constants predicted and experimentally observed at 25 °C for 14 different proteins, and the results demonstrated improvement of the quality of the database. PFDB consists of 141 (89 two-state and 52 non-two-state) single-domain globular proteins, which has the largest number among the currently available databases of protein folding kinetics. PFDB is thus intended to be used as a standard for developing and testing future predictive and theoretical studies of protein folding. PFDB can be accessed from the following link: http://lee.kias.re.kr/~bala/PFDB.

Exploring the Denatured State Ensemble by Single-Molecule Chemo-Mechanical Unfolding: The Effect of Force, Temperature, and Urea

While it is widely appreciated that the denatured state of a protein is a heterogeneous conformational ensemble, there is still debate over how this ensemble changes with environmental conditions. Here, we use single-molecule chemo-mechanical unfolding, which combines force and urea using the optical tweezers, together with traditional protein unfolding studies to explore how perturbants commonly used to unfold proteins (urea, force, and temperature) affect the denatured-state ensemble. We compare the urea m-values, which report on the change in solvent accessible surface area for unfolding, to probe the denatured state as a function of force, temperature, and urea. We find that while the urea- and force-induced denatured states expose similar amounts of surface area, the denatured state at high temperature and low urea concentration is more compact. To disentangle these two effects, we use destabilizing mutations that shift the T_m and C_m. We find that the compaction of the denatured state is related to changing temperature as the different variants of acyl-coenzyme A binding protein have similar m-values when they are at the same temperature but different urea concentration. These results have important implications for protein folding and stability under different environmental conditions.

Spatial ranges of driving forces are a key determinant of protein folding cooperativity and rate diversity

The physical basis of two-state-like folding transitions and the tremendous diversity in folding rates is elucidated by directly simulating the folding kinetics of 52 representative proteins. Relative to the results from a common modeling approach, the diversity of the simulated folding rates can be increased from ~10(2.1) to the experimental ~10(6.0) by a modest decrease in the spatial range of the attractive potential. The required theoretical range is consistent with desolvation physics and is notably much more permissive than that needed for two-state-like homopolymer collapse.

Probing the protein-folding mechanism using denaturant and temperature effects on rate constants

Protein folding has been extensively studied, but many questions remain regarding the mechanism. Characterizing early unstable intermediates and the high-free-energy transition state (TS) will help answer some of these. Here, we use effects of denaturants (urea, guanidinium chloride) and temperature on folding and unfolding rate constants and the overall equilibrium constant as probes of surface area changes in protein folding. We interpret denaturant kinetic m-values and activation heat capacity changes for 13 proteins to determine amounts of hydrocarbon and amide surface buried in folding to and from TS, and for complete folding. Predicted accessible surface area changes for complete folding agree in most cases with structurally determined values. We find that TS is advanced (50-90% of overall surface burial) and that the surface buried is disproportionately amide, demonstrating extensive formation of secondary structure in early intermediates. Models of possible pre-TS intermediates with all elements of the native secondary structure, created for several of these proteins, bury less amide and hydrocarbon surface than predicted for TS. Therefore, we propose that TS generally has both the native secondary structure and sufficient organization of other regions of the backbone to nucleate subsequent (post-TS) formation of tertiary interactions. The approach developed here provides proof of concept for the use of denaturants and other solutes as probes of amount and composition of the surface buried in coupled folding and other large conformational changes in TS and intermediates in protein processes.

Structural and thermodynamic investigations on the aggregation and folding of acylphosphatase by molecular dynamics simulations and solvation free energy analysis

Protein engineering method to study the mutation effects on muscle acylphosphatase (AcP) has been actively applied to describe kinetics and thermodynamics associated with AcP aggregation as well as folding processes. Despite the extensive mutation experiments, the molecular origin and the structural motifs for aggregation and folding kinetics as well as thermodynamics of AcP have not been rationalized at the atomic resolution. To this end, we have investigated the mutation effects on the structures and thermodynamics for the aggregation and folding of AcP by using the combination of fully atomistic, explicit-water molecular dynamics simulations, and three-dimensional reference interaction site model theory. The results indicate that the A30G mutant with the fastest experimental aggregation rate displays considerably decreased α1-helical contents as well as disrupted hydrophobic core compared to the wild-type AcP. Increased solvation free energy as well as hydrophobicity upon A30G mutation is achieved due to the dehydration of hydrophilic side chains in the disrupted α1-helix region of A30G. In contrast, the Y91Q mutant with the slowest aggregation rate shows a non-native H-bonding network spanning the mutation site to hydrophobic core and α1-helix region, which rigidifies the native state protein conformation with the enhanced α1-helicity. Furthermore, Y91Q exhibits decreased solvation free energy and hydrophobicity compared to wild type due to more exposed and solvated hydrophilic side chains in the α1-region. On the other hand, the experimentally observed slower folding rates in both mutants are accompanied by decreased helicity in α2-helix upon mutation. We here provide the atomic-level structures and thermodynamic quantities of AcP mutants and rationalize the structural origin for the changes that occur upon introduction of those mutations along the AcP aggregation and folding processes.

Mechanical unfolding of acylphosphatase studied by single-molecule force spectroscopy and MD simulations

Single-molecule manipulation methods provide a powerful means to study protein transitions. Here we combined single-molecule force spectroscopy and steered molecular-dynamics simulations to study the mechanical properties and unfolding behavior of the small enzyme acylphosphatase (AcP). We find that mechanical unfolding of AcP occurs at relatively low forces in an all-or-none fashion and is decelerated in the presence of a ligand, as observed in solution measurements. The prominent energy barrier for the transition is separated from the native state by a distance that is unusually long for alpha/beta proteins. Unfolding is initiated at the C-terminal strand (beta(T)) that lies at one edge of the beta-sheet of AcP, followed by unraveling of the strand located at the other. The central strand of the sheet and the two helices in the protein unfold last. Ligand binding counteracts unfolding by stabilizing contacts between an arginine residue (Arg-23) and the catalytic loop, as well as with beta(T) of AcP, which renders the force-bearing units of the protein resistant to force. This stabilizing effect may also account for the decelerated unfolding of ligand-bound AcP in the absence of force.

What lessons can be learned from studying the folding of homologous proteins?

The studies of the folding of structurally related proteins have proved to be a very important tool for investigating protein folding. Here we review some of the insights that have been gained from such studies. Our highlighted studies show just how such an investigation should be designed and emphasise the importance of the synergy between experiment and theory. We also stress the importance of choosing the right system carefully, exploiting the excellent structural and sequence databases at our disposal.

Urea denatured state ensembles contain extensive secondary structure that is increased in hydrophobic proteins

The goal of this article is to gain a better understanding of the denatured state ensemble (DSE) of proteins through an experimental and computational study of their denaturation by urea. Proteins unfold to different extents in urea and the most hydrophobic proteins have the most compact DSE and contain almost as much secondary structure as folded proteins. Proteins that unfold to the greatest extent near pH 7 still contain substantial amounts of secondary structure. At low pH, the DSE expands due to charge-charge interactions and when the net charge per residue is high, most of the secondary structure is disrupted. The proteins in the DSE appear to contain substantial amounts of polyproline II conformation at high urea concentrations. In all cases considered, including staph nuclease, the extent of unfolding by urea can be accounted for using the data and approach developed in the laboratory of Wayne Bolen (Auton et al., Proc Natl Acad Sci 2007; 104:15317-15323).

Agitation and high ionic strength induce amyloidogenesis of a folded PDZ domain in native conditions

Amyloid fibril formation is a distinctive hallmark of a number of degenerative diseases. In this process, protein monomers self-assemble to form insoluble structures that are generally referred to as amyloid fibrils. We have induced in vitro amyloid fibril formation of a PDZ domain by combining mechanical agitation and high ionic strength under conditions otherwise close to physiological (pH 7.0, 37 degrees C, no added denaturants). The resulting aggregates enhance the fluorescence of the thioflavin T dye via a sigmoidal kinetic profile. Both infrared spectroscopy and circular dichroism spectroscopy detect the formation of a largely intermolecular beta-sheet structure. Atomic force microscopy shows straight, rod-like fibrils that are similar in appearance and height to mature amyloid-like fibrils. Under these conditions, before aggregation, the protein domain adopts an essentially native-like structure and an even higher conformational stability (DeltaG(U-F)(H2O)). These results show a new method for converting initially folded proteins into amyloid-like aggregates. The methodological approach used here does not require denaturing conditions; rather, it couples agitation with a high ionic strength. Such an approach offers new opportunities to investigate protein aggregation under conditions in which a globular protein is initially folded, and to elucidate the physical forces that promote amyloid fibril formation.

The folding process of acylphosphatase from Escherichia coli is remarkably accelerated by the presence of a disulfide bond

The acylphosphatase from Escherichia coli (EcoAcP) is the first AcP so far studied with a disulfide bond. A mutational variant of the enzyme lacking the disulfide bond has been produced by substituting the two cysteine residues with alanine (EcoAcP mutational variant C5A/C49A, mutEcoAcP). The native states of the two protein variants are similar, as shown by far-UV and near-UV circular dichroism and dynamic light-scattering measurements. From unfolding experiments at equilibrium using intrinsic fluorescence and far-UV circular dichroism as probes, EcoAcP shows an increased conformational stability as compared with mutEcoAcP. The wild-type protein folds according to a two-state model with a very fast rate constant (k(F)(H2O)=72,600 s(-1)), while mutEcoAcP folds ca 1500-fold slower, via the accumulation of a partially folded species. The correlation between the hydrophobicity of the polypeptide chain and the folding rate, found previously in the AcP-like structural family, is maintained only when considering the mutant but not the wild-type protein, which folds much faster than expected from this correlation. Similarly, the correlation between the relative contact order and the folding rate holds only for mutEcoAcP. The correlation also holds for EcoAcP, provided the relative contact order value is recalculated by considering the disulfide bridge as an alternate path for the backbone to determine the shortest sequence separation between contacting residues. These results indicate that the presence of a disulfide bond in a protein is an important determinant of the folding rate and allows its contribution to be determined in quantitative terms.

Structure, conformational stability, and enzymatic properties of acylphosphatase from the hyperthermophile Sulfolobus solfataricus

The structure of AcP from the hyperthermophilic archaeon Sulfolobus solfataricus has been determined by (1)H-NMR spectroscopy and X-ray crystallography. Solution and crystal structures (1.27 A resolution, R-factor 13.7%) were obtained on the full-length protein and on an N-truncated form lacking the first 12 residues, respectively. The overall Sso AcP fold, starting at residue 13, displays the same betaalphabetabetaalphabeta topology previously described for other members of the AcP family from mesophilic sources. The unstructured N-terminal tail may be crucial for the unusual aggregation mechanism of Sso AcP previously reported. Sso AcP catalytic activity is reduced at room temperature but rises at its working temperature to values comparable to those displayed by its mesophilic counterparts at 25-37 degrees C. Such a reduced activity can result from protein rigidity and from the active site stiffening due the presence of a salt bridge between the C-terminal carboxylate and the active site arginine. Sso AcP is characterized by a melting temperature, Tm, of 100.8 degrees C and an unfolding free energy, DeltaG(U-F)H2O, at 28 degrees C and 81 degrees C of 48.7 and 20.6 kJ mol(-1), respectively. The kinetic and structural data indicate that mesophilic and hyperthermophilic AcP's display similar enzymatic activities and conformational stabilities at their working conditions. Structural analysis of the factor responsible for Sso AcP thermostability with respect to mesophilic AcP's revealed the importance of a ion pair network stabilizing particularly the beta-sheet and the loop connecting the fourth and fifth strands, together with increased density packing, loop shortening and a higher alpha-helical propensity.

Crystal structure and structural stability of acylphosphatase from hyperthermophilic archaeon Pyrococcus horikoshii OT3

To elucidate the structural basis for the high stability of acylphosphatase (AcP) from Pyrococcus horikoshii OT3, we determined its crystal structure at 1.72 A resolution. P. horikoshii AcP possesses high stability despite its approximately 30% sequence identity with eukaryotic enzymes that have moderate thermostability. The overall fold of P. horikoshii AcP was very similar to the structures of eukaryotic counterparts. The crystal structure of P. horikoshii AcP shows the same fold betaalphabetabetaalphabeta topology and the conserved putative catalytic residues as observed in eukaryotic enzymes. Comparison with the crystal structure of bovine common-type AcP and that of D. melanogaster AcP (AcPDro2) as representative of eukaryotic AcP revealed some significant characteristics in P. horikoshii AcP that likely play important roles in structural stability: (1) shortening of the flexible N-terminal region and long loop; (2) an increased number of ion pairs on the protein surface; (3) stabilization of the loop structure by hydrogen bonds. In P. horikoshii AcP, two ion pair networks were observed one located in the loop structure positioned near the C-terminus, and other on the beta-sheet. The importance of ion pairs for structural stability was confirmed by site-directed mutation and denaturation induced by guanidium chloride.

Desolvation is a likely origin of robust enthalpic barriers to protein folding

Experimental data from global analyses of temperature (T) and denaturant dependence of the folding rates of small proteins led to an intrinsic enthalpic folding barrier hypothesis: to a good approximation, the T-dependence of folding rate under constant native stability conditions is Arrhenius. Furthermore, for a given protein, the slope of isostability folding rate versus 1/T is essentially independent of native stability. This hypothesis implies a simple relationship between chevron and Eyring plots of folding that is easily discernible when both sets of rates are expressed as functions of native stability. Using experimental data in the literature, we verify the predicted chevron-Eyring relationship for 14 proteins and determine their intrinsic enthalpic folding barriers, which vary approximately from 15 kcal/mol to 40 kcal/mol for different proteins. These enthalpic barriers do not appear to correlate with folding rates, but they exhibit correlation with equilibrium unfolding enthalpy at room temperature. Intrinsic enthalpic barriers with similarly high magnitudes apply as well to at least two cases of peptide-peptide and peptide-protein association, suggesting that these barriers are a hallmark of certain general and fundamental kinetic processes during folding and binding. Using a class of explicit-chain C(alpha) protein models with constant elementary enthalpic desolvation barriers between C(alpha) positions, we show that small microscopic pairwise desolvation barriers, which are a direct consequence of the particulate nature of water, can act cooperatively to give rise to a significant overall enthalpic barrier to folding. This theoretical finding provides a physical rationalization for the high intrinsic enthalpic barriers in protein folding energetics. Ramifications of entropy-enthalpy compensation in hydrophobic association for the height of enthalpic desolvation barrier are discussed.

Prediction of protein folding rates from the amino acid sequence-predicted secondary structure

We present a method for predicting folding rates of proteins from their amino acid sequences only, or rather, from their chain lengths and their helicity predicted from their sequences. The method achieves 82% correlation with experiment over all 64 "two-state" and "multistate" proteins (including two artificial peptides) studied up to now.

Early collapse is not an obligate step in protein folding

The dimensions and secondary structure content of two proteins which fold in a two-state manner are measured within milliseconds of denaturant dilution using synchrotron-based, stopped-flow small-angle X-ray scattering and far-UV circular dichroism spectroscopy. Even upon a jump to strongly native conditions, neither ubiquitin nor common-type acylphosphatase contract prior to the major folding event. Circular dichroism and fluorescence indicate that negligible amounts of secondary and tertiary structures form in the burst phase. Thus, for these two denatured states, collapse and secondary structure formation are not energetically downhill processes even under aqueous, low-denaturant conditions. In addition, water appears to be as good a solvent as that with high concentrations of denaturant, when considering the over-all dimensions of the denatured state. However, the removal of denaturant does subtly alter the distribution of backbone dihedral phi,psi angles, most likely resulting in a shift from the polyproline II region to the helical region of the Ramachandran map. We consider the thermodynamic origins of these behaviors along with implications for folding mechanisms and computer simulations thereof.

Selection of antibody fragments specific for an alpha-helix region of acylphosphatase

The native state of common-type acylphosphatase (AcP) elicits two alpha-helices spanning residues 22-32 and 55-67 in the protein sequence. A peptide corresponding to the second alpha-helix (helix-2) of the protein was used to select phage antibodies consisting of a single chain fragment variable. The selection was performed in the presence of trifluoroethanol, a cosolvent known to induce the formation of helical structure in peptides and proteins. Phage scFv antibodies capable of binding the peptide specifically in a trifluoroethanol-induced alpha-helical conformation were isolated by affinity selection (biopanning). Some of these scFvs were also able to bind the native protein but not the peptide in a non-helical unstructured state. This indicates that the structural determinant recognized by the selected antibodies is the alpha-helical conformation of this specific region, rather than simply its amino acid sequence. This study shows that phage display libraries can be used to raise antibodies one can use as reagents able to target regions of a protein with a specific native-like secondary structure.

An Alternative Route for the Folding of Large RNAs: Apparent Two-state Folding by a Group II Intron Ribozyme

Despite a growing literature on the folding of RNA, our understanding of tertiary folding in large RNAs derives from studies on a small set of molecular examples, with primary focus on group I introns and RNase P RNA. To broaden the scope of RNA folding models and to better understand group II intron function, we have examined the tertiary folding of a ribozyme (D135) that is derived from the self-splicing ai5gamma intron from yeast mitochondria. The D135 ribozyme folds homogeneously and cooperatively into a compact, well-defined tertiary structure that includes all regions critical for active-site organization and substrate recognition. When D135 was treated with increasing concentrations of Mg(2+) and then subjected to hydroxyl radical footprinting, similar Mg(2+) dependencies were seen for internalization of all regions of the molecule, suggesting a highly cooperative folding behavior. In this work, we show that global folding and compaction of the molecule have the same magnesium dependence as the local folding previously observed. Furthermore, urea denaturation studies indicate highly cooperative unfolding of the ribozyme that is governed by thermodynamic parameters similar to those for forward folding. In fact, D135 folds homogeneously and cooperatively from the unfolded state to its native, active structure, thereby demonstrating functional reversibility in RNA folding. Taken together, the data are consistent with two-state folding of the D135 ribozyme, which is surprising given the size and multi-domain structure of the RNA. The findings establish that the accumulation of stable intermediates prior to formation of the native state is not a universal feature of RNA folding and that there is an alternative paradigm in which the folding landscape is relatively smooth, lacking rugged features that obstruct folding to the native state.

A large scale test of computational protein design: folding and stability of nine completely redesigned globular proteins

A previously developed computer program for protein design, RosettaDesign, was used to predict low free energy sequences for nine naturally occurring protein backbones. RosettaDesign had no knowledge of the naturally occurring sequences and on average 65% of the residues in the designed sequences differ from wild-type. Synthetic genes for ten completely redesigned proteins were generated, and the proteins were expressed, purified, and then characterized using circular dichroism, chemical and temperature denaturation and NMR experiments. Although high-resolution structures have not yet been determined, eight of these proteins appear to be folded and their circular dichroism spectra are similar to those of their wild-type counterparts. Six of the proteins have stabilities equal to or up to 7kcal/mol greater than their wild-type counterparts, and four of the proteins have NMR spectra consistent with a well-packed, rigid structure. These encouraging results indicate that the computational protein design methods can, with significant reliability, identify amino acid sequences compatible with a target protein backbone.

Comparison of the folding processes of distantly related proteins. Importance of hydrophobic content in folding

The N-terminal domain of HypF from Escherichia coli (HypF-N) is a 91 residue protein module sharing the same folding topology and a significant sequence identity with two extensively studied human proteins, muscle and common-type acylphosphatases (mAcP and ctAcP). With the aim of learning fundamental aspects of protein folding from the close comparison of so similar proteins, the folding process of HypF-N has been studied using stopped-flow fluorescence. While mAcP and ctAcP fold in a two-state fashion, HypF-N was found to collapse into a partially folded intermediate before reaching the fully folded conformation. Formation of a burst-phase intermediate is indicated by the roll over in the Chevron plot at low urea concentrations and by the large jump of intrinsic and 8-anilino-1-naphtalenesulphonic acid-derived fluorescence immediately after removal of denaturant. Furthermore, HypF-N was found to fold rapidly with a rate constant that is approximately two and three orders of magnitudes faster than ctAcP and mAcP, respectively. Differences between the bacterial protein and the two human counterparts were also found as to the involvement of proline isomerism in their respective folding processes. The results clearly indicate that features that are often thought to be relevant in protein folding are not highly conserved in the evolution of the acylphosphatase superfamily. The large difference in folding rate between mAcP and HypF-N cannot be entirely accounted for by the difference in relative contact order or related topological metrics. The analysis shows that the higher folding rate of HypF-N is in part due to the relatively high hydrophobic content of this protein. This conclusion, which is also supported by the highly significant correlation found between folding rate and hydrophobic content within a group of proteins displaying the topology of HypF-N and AcPs, suggests that the average hydrophobicity of a protein sequence is an important determinant of its folding rate.

Prediction of folding rates and transition-state placement from native-state geometry

A variety of experimental and theoretical studies have established that the folding process of monomeric proteins is strongly influenced by the topology of the native state. In particular, folding times have been shown to correlate well with the contact order, a measure of contact locality. Our investigation focuses on identifying additional topologic properties that correlate with experimentally measurable quantities, such as folding rates and transition-state placement, for both two- and three-state folders. The validation against data from 40 experiments shows that a particular topological property that measures the interdependence of contacts, termed cliquishness or clustering coefficient, can account with statistically significant accuracy both for the transition state placement and especially for folding rates. The observed correlations can be further improved by optimally combining the distinct topological information captured by cliquishness and contact order.

Evidence for sequential barriers and obligatory intermediates in apparent two-state protein folding

Many small proteins fold fast and without detectable intermediates. This is frequently taken as evidence against the importance of partially folded states, which often transiently accumulate during folding of larger proteins. To get insight into the properties of free energy barriers in protein folding we analyzed experimental data from 23 proteins that were reported to show non-linear activation free-energy relationships. These non-linearities are generally interpreted in terms of broad transition barrier regions with a large number of energetically similar states. Our results argue against the presence of a single broad barrier region. They rather indicate that the non-linearities are caused by sequential folding pathways with consecutive distinct barriers and a few obligatory high-energy intermediates. In contrast to a broad barrier model the sequential model gives a consistent picture of the folding barriers for different variants of the same protein and when folding of a single protein is analyzed under different solvent conditions. The sequential model is also able to explain changes from linear to non-linear free energy relationships and from apparent two-state folding to folding through populated intermediates upon single point mutations or changes in the experimental conditions. These results suggest that the apparent discrepancy between two-state and multi-state folding originates in the relative stability of the intermediates, which argues for the importance of partially folded states in protein folding.

Folding and aggregation are selectively influenced by the conformational preferences of the alpha-helices of muscle acylphosphatase

The native state of human muscle acylphosphatase (AcP) presents two alpha-helices. In this study we have investigated folding and aggregation of a number of protein variants having mutations aimed at changing the propensity of these helical regions. Equilibrium and kinetic measurements of folding indicate that only helix-2, spanning residues 55-67, is largely stabilized in the transition state for folding therefore playing a relevant role in this process. On the contrary, the aggregation rate appears to vary only for the variants in which the propensity of the region corresponding to helix-1, spanning residues 22-32, is changed. Mutations that stabilize the first helix slow down the aggregation process while those that destabilize it increase the aggregation rate. AcP variants with the first helix destabilized aggregate with rates increased to different extents depending on whether the introduced mutations also alter the propensity to form beta-sheet structure. The fact that the first alpha-helix is important for aggregation and the second helix is important for folding indicates that these processes are highly specific. This partitioning does not reflect the difference in intrinsic alpha-helical propensities of the two helices, because helix-1 is the one presenting the highest propensity. Both processes of folding and aggregation do not therefore initiate from regions that have simply secondary structure propensities favorable for such processes. The identification of the regions involved in aggregation and the understanding of the factors that promote such a process are of fundamental importance to elucidate the principles by which proteins have evolved and for successful protein design.

The Thermal Stability of the Fusarium solani pisi Cutinase as a Function of pH

We have investigated the thermal stability of the Fusarium solani pisi cutinase as a function of pH, in the range from pH 2-12. Its highest enzymatic activity coincides with the pH-range at which it displays its highest thermal stability. The unfolding of the enzyme as a function of pH was investigated by microcalorimetry. The ratio between the calorimetric enthalpy (DeltaH(cal)) and the van't Hoff enthalpy (DeltaH(v)) obtained, is far from unity, indicating that cutinase does not exhibit a simple two state unfolding behaviour. The role of pH on the electrostatic contribution to the thermal stability was assessed using TITRA. We propose a molecular interpretation for the pH-variation in enzymatic activity.

Initial denaturing conditions influence the slow folding phase of acylphosphatase associated with proline isomerization

The folding kinetics of human common-type acylphosphatase (cAcP) from its urea- and TFE-denatured states have been determined by stopped-flow fluorescence techniques. The refolding reaction from the highly unfolded state formed in urea is characterized by double exponential behavior that includes a slow phase associated with isomerism of the Gly53-Pro54 peptide bond. However, this slow phase is absent when refolding is initiated by dilution of the highly a-helical denatured state formed in the presence of 40% trifluoroethanol (TFE). NMR studies of a peptide fragment corresponding to residues Gly53-Gly69 of cAcP indicate that only the native-like trans isomer of the Gly-Pro peptide bond is significantly populated in the presence of TFE, whereas both the cis and trans isomers are found in an approximately 1:9 ratio for the peptide bond in aqueous solution. Molecular modeling studies in conjunction with NMR experiments suggest that the trans isomer of the Gly53-Pro54 peptide bond is stabilized in TFE by the formation of a nonnative-like hydrogen bond between the CO group of Gly53 and the NH group of Lys57. These results therefore reveal that a specific nonnative interaction in the denatured state can increase significantly the overall efficiency of refolding.

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.

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.