Folding under the influence

Fast closure of long loops at the initiation of the folding transition of globular proteins studied by time-resolved FRET-based methods

The protein folding problem would be considered “solved” when it will be possible to “read genes”, i.e., to predict the native fold of proteins, their dynamics, and the mechanism of fast folding based solely on sequence data. The long-term goal should be the creation of an algorithm that would simulate the stepwise mechanism of folding, which constrains the conformational space and in which random search for stable interactions is possible. Here, we focus attention on the initial phases of the folding transition starting with the compact disordered collapsed ensemble, in search of the initial sub-domain structural biases that direct the otherwise stochastic dynamics of the backbone. Our studies are designed to test the “loop hypothesis”, which suggests that fast closure of long loop structures by non-local interactions between clusters of mainly non-polar residues is an essential conformational step at the initiation of the folding transition of globular proteins. We developed and applied experimental methods based on time-resolved resonance excitation energy transfer (trFRET) measurements combined with fast mixing methods and studied the initial phases of the folding of

The loop hypothesis: contribution of early formed specific non-local interactions to the determination of protein folding pathways

The extremely fast and efficient folding transition (in seconds) of globular proteins led to the search for some unifying principles embedded in the physics of the folding polypeptides. Most of the proposed mechanisms highlight the role of local interactions that stabilize secondary structure elements or a folding nucleus as the starting point of the folding pathways, i.e., a "bottom-up" mechanism. Non-local interactions were assumed either to stabilize the nucleus or lead to the later steps of coalescence of the secondary structure elements. An alternative mechanism was proposed, an "up-down" mechanism in which it was assumed that folding starts with the formation of very few non-local interactions which form closed long loops at the initiation of folding. The possible biological advantage of this mechanism, the "loop hypothesis", is that the hydrophobic collapse is associated with ordered compactization which reduces the chance for degradation and misfolding. In the present review the experiments, simulations and theoretical consideration that either directly or indirectly support this mechanism are summarized. It is argued that experiments monitoring the time-dependent development of the formation of specifically targeted early-formed sub-domain structural elements, either long loops or secondary structure elements, are necessary. This can be achieved by the time-resolved FRET-based "double kinetics" method in combination with mutational studies. Yet, attempts to improve the time resolution of the folding initiation should be extended down to the sub-microsecond time regime in order to design experiments that would resolve the classes of proteins which first fold by local or non-local interactions.

A comparison of experimental and computational methods for mapping the interactions present in the transition state for folding of FKBP12

The folding pathway of FKBP12, a 107 residue α/β protein, has been characterised in detail using a combination of experimental and computational techniques. FKBP12 follows a two-state model of folding in which only the denatured and native states are significantly populated; no intermediate states are detected. The refolding rate constant in water is 4 s(-1) at 25 °C. Two different experimental strategies were employed for studying the transition state for folding. In the first case, a non-mutagenic approach was used and the unfolding and refolding of the wild-type protein measured as a function of experimental conditions such as temperature, denaturant, ligand and trifluoroethanol (TFE) concentration. These data suggest a compact transition state relative to the unfolded state with some 70% of the surface area buried. The ligand-binding site, whichis mainly formed by two long loops, is largely unstructured in the transition state. TFE experiments suggest that the α-helix may be formed in the transition state. The second experimental approach involved using protein engineering techniques with φ-value analysis. Residue-specific information on the structure and energetics of the transition state can be obtained by this method. 34 mutations were made at sites throughout the protein to probe the extent of secondary and tertiary structure in the transition state. In contrast to some other proteins of this size, no element of structure is fully formed in the transition state, instead, the transition state is similar to that found for smaller, single-domain proteins, such as chymotrypsin inhibitor 2 and the SH3 domainfrom α-spectrin. For FKBP12, the central three strands of the β-sheet (2, 4 and 5), comprise the most structured region of the transition state. In particular Val 101, which is one of the most highly buried residues and located in the middle of the central β-strand,makes approximately 60% of its native interactions. The outer β-strands, and the ends of the central β-strands are formed to a lesser degree. The short α-helix is largely unstructured in the transition state as are the loops. The data are consistent with a nucleation-condensation model of folding, the nucleus of which is formed by side chains within β-strands 2, 4 and 5 and the C-terminus of the α-helix. These residues are distant in the primary sequence, demonstrating the importance of tertiary interactions in the transition state. High-temperature molecular dynamic simulations on the unfoldingpathway of FKBP12 are in good agreement with the experimental results.

Syncytial phenotype of C-terminally truncated herpes simplex virus type 1 gB is associated with diminished membrane interactions

The cytoplasmic domain of glycoprotein B (gB) from herpes simplex virus type 1 (HSV-1) is an important regulator of membrane fusion. C-terminal truncations of the cytoplasmic domain lead to either hyperfusion or fusion-null phenotypes. Currently, neither the structure of the cytoplasmic domain nor its mechanism of fusion regulation is known. Here we show, for the first time, that the full-length cytoplasmic domain of HSV-1 gB associates stably with lipid membranes, preferentially binding to membranes containing anionic head groups. This interaction involves a large increase in helical content. However, the truncated cytoplasmic domains associated with the hyperfusion phenotype show a small increase in helical structure and a diminished association with lipid membranes, whereas the one associated with the fusion-null phenotype shows no increase in helical structure and only a minimal association with lipid membranes. We hypothesize that stable binding to lipid membranes is an important part of the mechanism by which the cytoplasmic domain negatively regulates membrane fusion. Moreover, our experiments with truncated cytoplasmic domains point to two specific regions that are critical for membrane interactions. Taken together, our work provides several important new insights into the architecture of the cytoplasmic domain of HSV-1 gB and its interaction with lipid membranes.

Ubiquitin: a small protein folding paradigm

For the past twenty years, the small, 76-residue protein ubiquitin has been used as a model system to study protein structure, stability, folding and dynamics. In this time, ubiquitin has become a paradigm for both the experimental and computational folding communities. The folding energy landscape is now uniquely characterised with a plethora of information available on not only the native and denatured states, but partially structured states, alternatively folded states and locally unfolded states, in addition to the transition state ensemble. This Perspective focuses on the experimental characterisation of ubiquitin using a comprehensive range of biophysical techniques.

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.

Trifluoroethanol increases the stability of Delta(5)-3-ketosteroid isomerase. 15N NMR relaxation studies

In the equilibrium unfolding process of Delta(5)-3-ketosteroid isomerase from Pseudomonas testosteroni by urea, it was observed that the enzyme stability increases by 2.5 kcal/mol in the presence of 5% trifluoroethanol (TFE). To elucidate the increased enzyme stability by TFE, the backbone dynamics of Delta(5)-3-ketosteroid isomerase were studied in the presence and absence of 5% TFE by (15)N NMR relaxation measurements, and the motional parameters (S(2), tau(e), and R(ex)) were extracted from the relaxation data using the model-free formalism. The presence of 5% TFE causes little change or a slight increase in the order parameters (S(2)) for a number of residues, which are located mainly in the dimer interface region. However, the majority of the residues exhibit reduced order parameters in the presence of 5% TFE, indicating that high frequency (pico- to nanosecond) motions are generally enhanced by TFE. The results suggest that the entropy can be an important factor for the enzyme stability, and the increase in entropy by TFE is partially responsible for the increased stability of Delta(5)-3-ketosteroid isomerase.

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 building block folding model and the kinetics of protein folding

Here we show that qualitatively, the building blocks folding model accounts for three-state versus the two-state protein folding. Additionally, it is consistent with the faster versus slower folding rates of the two-state proteins. Specifically, we illustrate that the building blocks size, their mode of associations in the native structure, the number of ways they can combinatorially assemble, their population times and the way they are split in the iterative, step-by-step structural dissection which yields the anatomy trees, explain a broad range of folding rates. We further show that proteins with similar general topologies may have different folding pathways, and hence different folding rates. On the other hand, the effect of mutations resembles that of changes in conditions, shifting the population times and hence the energy landscapes. Hence, together with the secondary structure type and the extent of local versus non-local interactions, a coherent, consistent rationale for folding kinetics can be outlined, in agreement with experimental results. Given the native structure of a protein, these guidelines enable a qualitative prediction of the folding kinetics. We further describe these in the context of the protein folding energy landscape. Quantitatively, in principle, the diffusion-collision model for the building block association can be used. However, the folding rates of the building blocks and traps in their formation and association, need to be considered.

Trifluoroethanol-assisted protein folding: fragment 53--103 of bovine alpha-lactalbumin

Fragment 53--103 of bovine alpha-lactalbumin, prepared by limited peptic digestion of the protein at low pH, is a 51-residue polypeptide chain crosslinked by two disulfide bonds encompassing helix C (residues 86--98) of the native protein. Refolding of the fully reduced fragment (four--SH groups) is expected to lead to three fully oxidized isomers, the native (61--77, 73--91) and the two misfolded species named ribbon (61--91, 73--77) and beads (61--73, 77--91) isomers. The fragment with correct disulfide bonds was formed in approx. 30% yield when refolding was conducted in aqueous solution at neutral pH in the presence of the redox system constituted by reduced and oxidized glutathione. On the other hand, when the reaction was conducted in 30% (v/v) trifluoroethanol (TFE), the oxidative refolding to the native isomer was almost quantitative. To provide an explanation of the beneficial effect of TFE in promoting the correct oxidative folding, the conformational features of the various fragment species were analyzed by far-UV circular dichroism measurements. The fully reduced fragment is largely unfolded in water, but it becomes helical in aqueous TFE. Correctly refolded fragment is produced most when the helical contents of the reduced and oxidized fragment in aqueous TFE are roughly equal. It is proposed that 30% TFE promotes a native-like format of the fragment and thus an efficient and correct pairing of disulfides. Higher concentrations of TFE, instead, promote some non-native helical secondary structure in the fragment species, thus hampering correct folding.

Point mutations and sequence variability in proteins: redistributions of preexisting populations

Here we study the effect of point mutations in proteins on the redistributions of the conformational substates. We show that regardless of the location of a mutation in the protein structure and of its type, the observed movements of the backbone recur largely at the same positions in the structures. Despite the different interactions that are disrupted and formed by the residue substitution, not only are the conformations very similar, but the regions that move are also the same, regardless of their sequential or spatial distance from the mutation. This observation leads us to conclude that, apart from some extreme cases, the details of the interactions are not critically important in determining the protein conformation or in specifying which parts of the protein would be more prone to take on different local conformations in response to changes in the sequence. This finding further illustrates why proteins manifest a robustness toward many mutational events. This nonuniform distribution of the conformer population is consistently observed in a variety of protein structural types. Topology is critically important in determining folding pathways, kinetics, building block cutting, and anatomy trees. Here we show that topology is also very important in determining which regions of the protein structure will respond to sequence changes, regardless of the sequential or spatial location of the mutation.

The duplication of an eight-residue helical stretch in Staphylococcal nuclease is not helical: a model for evolutionary change

A common method of evolutionary change is gene duplication, followed by other events that lead to new function, decoration of folds, oligomerization, or other changes. As part of a study on the potential for evolutionary change created by duplicated sequences, we have carried out a crystallographic study on a mutant of Staphylococcal nuclease in which residues 55-62 have been duplicated in a wild-type variant termed PHS. In the parental protein (PHS) these residues form the first two turns of a helix running from residue 54 to 68 (hereafter designated as helix I). The crystal structure of the mutant is very similar to that of the parental, with helix I being unaltered. The duplicated residues are accommodated by expanding an existing loop N-terminal to helix I. In addition, circular dichroism (CD) studies have been carried out on a parental peptide containing helix I with six flanking residues at each terminus (residues 48-74) and on the same peptide expanded by the duplication, as a function of 2,2,2-trifluoroethanol (TFE) concentration. Each peptide possesses only modest helical propensity in solution. Our data, which is different from what was observed in T4 lysozyme, show that the conformation of the duplicated sequence is determined by a balance of sequential and longer-range effects. Thus duplicating sequence need not mean duplicating structure. Proteins 2000;40:465-472.

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.