Redesigning the folding energetics of a model three-helix bundle protein by site-directed mutagenesis

Because of their limited size and complexity, de novo designed proteins are particularly useful for the detailed investigation of folding thermodynamics and mechanisms. Here, we describe how subtle changes in the hydrophobic core of a model three-helix bundle protein (GM-0) alter its folding energetics. To explore the folding tolerance of GM-0 toward amino acid sequence variability, two mutant proteins (GM-1 and GM-2) were generated. In the mutants, cavities were created in the hydrophobic core of the protein by either singly (GM-1; L35A variant) or doubly (GM-2; L35A/I39A variant) replacing large hydrophobic side chains by smaller Ala residues. The folding of GM-0 is characterized by two partially folded intermediate states exhibiting characteristics of molten globules, as evidenced by pressure-unfolding and pressure-assisted cold denaturation experiments. In contrast, the folding energetics of both mutants, GM-1 and GM-2, exhibit only one folding intermediate. Our results support the view that simple but biologically important folding motifs such as the three-helix bundle can reveal complex folding plasticity, and they point to the role of hydrophobic packing as a determinant of the overall stability and folding thermodynamic of the helix bundle.     

Folding kinetics of designer proteins. Application of the diffusion-collision model to a de novo designed four-helix bundle.

A folding algorithm is described, based on the diffusion-collision model, combining static and dynamic calculational methods. The algorithm is applied to predict the basic structure and schematic folding pathways of an artificial four-helix bundle.     

Characterization of the folding kinetics of a three-helix bundle protein via a minimalist Langevin model.

We use a simple off-lattice Langevin model of protein folding to characterize the folding and unfolding of a fast-folding, 46 residue three-helix bundle. Under conditions at which the C-terminal helix is 30 % stable, we observe a clear three-state folding mechanism. In the on-pathway intermediate state, the middle and C-terminal helices are folded and in contact with each other, while the N-terminal region remains disordered. Nevertheless, under these conditions this intermediate is thermodynamically unstable relative to its unfolded state. The first and highest folding barrier corresponds to the organization of the hinge between the middle and C-terminal helices. A subsequent major barrier corresponds to the organization of the hinge between the middle and N-terminal helices. Hyperstabilizing the hinge regions leads to twice the folding rate that is obtained from hyperstabilizing the helices, even though much fewer contacts are involved in hinge hyperstabilization than in helix hyperstabilization. Unfolding follows single-exponential kinetics, even at temperatures only slightly above the folding transition temperature.     

Local secondary structure content predicts folding rates for simple,two-state proteins

Many single-domain proteins exhibit two-state folding kinetics, with folding rates that span more than six orders of magnitude. A quantity of much recent interest for such proteins is their contact order, the average separation in sequence between contacting residue pairs. Numerous studies have reached the surprising conclusion that contact order is well-correlated with the logarithm of the folding rate for these small, well-characterized molecules. Here, we investigate the physico-chemical basis for this finding by asking whether contact order is actually a composite number that measures the fraction of local secondary structure in the protein; viz. turns, helices, and hairpins. To pursue this question, we calculated the secondary structure content for 24 two-state proteins and obtained coefficients that predict their folding rates. The predicted rates correlate strongly with experimentally determined rates, comparable to the correlation with contact order. Further, these predicted folding rates are correlated strongly with contact order. Our results suggest that the folding rate of two-state proteins is a function of their local secondary structure content, consistent with the hierarchic model of protein folding. Accordingly, it should be possible to utilize secondary structure prediction methods to predict folding rates from sequence alone.     

Sub-microsecond protein folding

We have investigated the structure, equilibria, and folding kinetics of an engineered 35-residue subdomain of the chicken villin headpiece, an ultrafast-folding protein. Substitution of two buried lysine residues by norleucine residues stabilizes the protein by 1 kcal/mol and increases the folding rate sixfold, as measured by nanosecond laser T-jump. The folding rate at 300 K is (0.7 micros)(-1) with little or no temperature dependence, making this protein the first sub-microsecond folder, with a rate only twofold slower than the theoretically predicted speed limit. Using the 70 ns process to obtain the effective diffusion coefficient, the free energy barrier height is estimated from Kramers theory to be less than approximately 1 kcal/mol. X-ray crystallographic determination at 1A resolution shows no significant change in structure compared to the single-norleucine-substituted molecule and suggests that the increased stability is electrostatic in origin. The ultrafast folding rate, very accurate X-ray structure, and small size make this engineered villin subdomain an ideal system for simulation by atomistic molecular dynamics with explicit solvent.     

The role of sequence and structure in protein folding kinetics; the diffusion-collision model applied to proteins L and G

The diffusion-collision model (DCM) is applied to the folding kinetics of protein L and protein G. In the DCM, the two proteins are treated as consisting of two beta-hairpins and one alpha-helix, so that they are isomorphous with the three-helix bundle DCM model. In the absence of sequence dependent factors, both proteins would fold in the same way in the DCM, with the coalescence of the N-terminal hairpin and the helix slightly favored over the C-terminal hairpin and the helix because the former are closer together than the latter. However, sequence dependent factors make the N-terminal hairpin of protein L and the C-terminal hairpin of protein G more stable in the ensemble of unfolded conformations. This difference in the stabilities gives rise to the difference in the calculated folding behavior, in agreement with experiment.     

Retinal proteins as model systems for membrane protein folding

Experimental folding studies of membrane proteins are more challenging than water-soluble proteins because of the higher hydrophobicity content of membrane embedded sequences and the need to provide a hydrophobic milieu for the transmembrane regions. The first challenge is their denaturation: due to the thermodynamic instability of polar groups in the membrane, secondary structures in membrane proteins are more difficult to disrupt than in soluble proteins. The second challenge is to refold from the denatured states. Successful refolding of membrane proteins has almost always been from very subtly denatured states. Therefore, it can be useful to analyze membrane protein folding using computational methods, and we will provide results obtained with simulated unfolding of membrane protein structures using the Floppy Inclusions and Rigid Substructure Topography (FIRST) method. Computational methods have the advantage that they allow a direct comparison between diverse membrane proteins. We will review here both, experimental and FIRST studies of the retinal binding proteins bacteriorhodopsin and mammalian rhodopsin, and discuss the extension of the findings to deriving hypotheses on the mechanisms of folding of membrane proteins in general. This article is part of a Special Issue entitled: Retinal Proteins—You can teach an old dog new tricks.     

Brownian dynamics simulation of protein folding: a study of the diffusion-collision model

The dynamic aspects of protein folding are described by a series of diffusion-collision steps involving structural units (microdomains) of various sizes that combine to form the protein in its native state. A method is introduced for obtaining the rate constants for the basic diffusion-collision step by use of Brownian dynamics. The method is applied to an investigation of the folding dynamics of two α-helices connected by a flexible (random-coil) polypeptide chain. The results of this full three-dimensional treatment are compared with simplified model calculations for the diffusion-collision step. Of particular interest are the nature of the collision dynamics and the role of the intervening peptide chain.     

The folding transition state of the cold shock protein is strongly polarized

It has been shown recently that an 11-residue peptide fragment of transthyretin, TTR(105-115), can form amyloid fibrils in vitro by adopting an extended beta-strand conformation. We used molecular dynamics simulations on systems of TTR(105-115) peptides, for a total length of about 5 micros, to explore the process of self-assembly and the structures of the resulting aggregates. Our results suggest that an antiparallel association of the beta-strands is more probable than a parallel one and that the central residues (T106-L111) in a beta-strand have a high propensity to form inter-peptide hydrogen bonds. The study of the dynamics of self-association indicated that, for this peptide, trajectories leading to conformations with high alpha-helical content are off-pathway from those leading to aggregates with high beta-structure content. We also show that the diverse oligomeric structures that form spontaneously in the molecular dynamics simulations are, to a large extent, compatible with solid-state NMR experimental measurements, including chemical shifts, on fully formed fibrils. The strategy that we present may therefore be used in the design of new experiments to determine the structure of amyloid fibrils, such as those involving site-specific isotope labelling schemes to measure key inter-atomic distances.     

Rapid collapse precedes the fast two-state folding of the cold shock protein

The cold shock protein Bc-Csp folds very rapidly in a reaction that is well described by a kinetic two-state mechanism without intermediates. We measured the shortening of six intra-protein distances during folding by F?rster resonance energy transfer (FRET) in combination with stopped-flow experiments. Single tryptophan residues were engineered into the protein as the donors, and single 5-(((acetylamino)ethyl)amino)naphthalene-1-sulfonate (AEDANS) residues were placed as the acceptors at solvent-exposed sites of Bc-Csp. Their R0 value of about 22 A was well suited for following distance changes during the folding of this protein with a high sensitivity. The mutagenesis and the labeling did not alter the refolding kinetics. The changes in energy transfer during folding were monitored by both donor and acceptor emission and reciprocal effects were found. In two cases the donor-acceptor distances were similar in the unfolded and the folded state and, as a consequence, the kinetic changes in energy transfer upon folding were very small. For four donor/acceptor pairs we found that > or =50% of the increase in energy transfer upon folding occurred prior to the rate-limiting step of folding. This reveals that about half of the shortening of the intra-molecular distances upon folding has occurred already before the rate-limiting step and suggests that the fast two-state folding reaction of Bc-Csp is preceded by a very rapid collapse.     

Correlation of folding kinetics with the number and isomerization states of prolines in three homologous proteins of the RNase family

Several studies attribute the slower phases in protein folding to prolyl isomerizations, and several others do not. A correlation exists between the number of prolines in a protein and the complexity of the mechanism with which it folds. In this study, we have demonstrated a direct correlation between the number of cis-prolyl bonds in a native protein and the complexity with which it folds via slower phases by studying the folding of three structurally homologous proteins of the ribonuclease family, namely RNase A, onconase and angiogenin, which differ in the number and isomerization states of their proline residues.     

Probing the folding and unfolding dynamics of secondary and tertiary structures in a three-helix bundle protein

Fast relaxation kinetics studies of the B-domain of staphylococcal protein A were performed to characterize the folding and unfolding of this small three-helix bundle protein. The relaxation kinetics were initiated using a laser-induced temperature jump and probed using time-resolved infrared spectroscopy. The kinetics monitored within the amide I' absorbance of the polypeptide backbone exhibit two distinct kinetics phases with nanosecond and microsecond relaxation times. The fast kinetics relaxation time is close to the diffusion limits placed on protein folding reactions. The fast kinetics phase is dominated by the relaxation of the solvated helix (nu = 1632 cm(-1)), which reports on the fast relaxation of the individual helices. The slow kinetics phase follows the cooperative relaxation of the native helical bundle core that is monitored by both solvated (nu = 1632 cm(-1)) and buried helical IR bands (nu = 1652 cm(-1)). The folding rates of the slow kinetics phase calculated over an extended temperature range indicate that the core formation of this protein follows a pathway that is energetically downhill. The unfolding rates are much more strongly temperature-dependent indicating an activated process with a large energy barrier. These results provide significant insight into the primary process of protein folding and suggest that fast formation of helices can drive the folding of helical proteins.     

Outlining folding nuclei in globular proteins

Our theoretical approach for prediction of folding/unfolding nuclei in three-dimensional protein structures is based on a search for free energy saddle points on networks of protein unfolding pathways. Under some approximations, this search is performed rapidly by dynamic programming and results in prediction of Phi values, which can be compared with those found experimentally. In this study, we optimize some details of the model (specifically, hydrogen atoms are taken into account in addition to heavy atoms), and compare the theoretically obtained and experimental Phi values (which characterize involvement of residues in folding nuclei) for all 17 proteins, where Phi values are now known for many residues. We show that the model provides good Phi value predictions for proteins whose structures have been determined by X-ray analysis (the average correlation coefficient is 0.65), with a more limited success for proteins whose structures have been determined by NMR techniques only (the average correlation coefficient is 0.34), and that the transition state free energies computed from the same model are in a good anticorrelation with logarithms of experimentally measured folding rates at mid-transition (the correlation coefficient is -0.73).     

The folding mechanism of a two-domain protein: folding kinetics and domain docking of the gene-3 protein of phage fd

The gene-3 protein (G3P) of filamentous phages is essential for the infection of Escherichia coli. The carboxy-terminal domain anchors this protein in the phage coat, whereas the two amino-terminal domains N1 and N2 protrude from the phage surface. We analyzed the folding mechanism of the two-domain fragment N1-N2 of G3P (G3P(*)) and the interplay between folding and domain assembly. For this analysis, a variant of G3P(*) was used that contained four stabilizing mutations (IIHY-G3P(*)). The observed refolding kinetics extend from 10 ms to several hours. Domain N1 refolds very rapidly (with a time constant of 9.4 ms at 0.5 M guanidinium chloride, 25 degrees C) both as a part of IIHY-G3P(*) and as an isolated protein fragment. The refolding of domain N2 is slower and involves two reactions with time constants of seven seconds and 42 seconds. These folding reactions of the individual domains are followed by a very slow, spectroscopically silent docking process, which shows a time constant of 6200 seconds. This reaction was detected by a kinetic unfolding assay for native molecules. Before docking, N1 and N2 unfold fast and independently, after docking they unfold slowly in a correlated fashion. A high energy barrier is thus created by domain docking, which protects G3P kinetically against unfolding. The slow domain docking is possibly important for the infection of E.coli by the phage. Upon binding to the F pilus, the N2 domain separates from N1 and the binding site for TolA on domain N1 is exposed. Since domain reassembly is so slow, this binding site remains accessible until pilus retraction has brought N1 close to TolA on the bacterial surface.     

Alternative binding proteins: affibody binding proteins developed from a small three-helix bundle scaffold

In recent years, classical antibody-based affinity reagents have been challenged by novel types of binding proteins developed by combinatorial protein engineering principles. One of these classes of binding proteins of non-Ig origin are the so-called affibody binding proteins, functionally selected from libraries of a small (6 kDa), non-cysteine three-helix bundle domain used as a scaffold. During the first 10 years since they were first described, high-affinity affibody binding proteins have been selected towards a large number of targets for use in a variety of applications, such as bioseparation, diagnostics, functional inhibition, viral targeting and in vivo tumor imaging/therapy. The small size offers the possibility to produce functional affibody binding proteins also by chemical synthesis production routes, which has been found to be advantageous for the site-specific introduction of various labels and radionuclide chelators.     

Conformational stability and folding mechanisms of dimeric proteins

The folding of multisubunit proteins is of tremendous biological significance since the large majority of proteins exist as protein-protein complexes. Extensive experimental and computational studies have provided fundamental insights into the principles of folding of small monomeric proteins. Recently, important advances have been made in extending folding studies to multisubunit proteins, in particular homodimeric proteins. This review summarizes the equilibrium and kinetic theory and models underlying the quantitative analysis of dimeric protein folding using chemical denaturation, as well as the experimental results that have been obtained. Although various principles identified for monomer folding also apply to the folding of dimeric proteins, the effects of subunit association can manifest in complex ways, and are frequently overlooked. Changes in molecularity typically give rise to very different overall folding behaviour than is observed for monomeric proteins. The results obtained for dimers have provided key insights pertinent to understanding biological assembly and regulation of multisubunit proteins. These advances have set the stage for future advances in folding involving protein-protein interactions for natural multisubunit proteins and unnatural assemblies involved in disease.     

Folding of RNase T1 is decelerated by a specific tertiary contact in a folding intermediate.

The replacement of tryptophan 59 of ribonuclease T1 by a tyrosine residue does not change the stability of the protein. However, it leads to a strong acceleration of a major, proline-limited reaction that is unusually slow in the refolding of the wild-type protein. The distribution of fast- and slow-folding species and the kinetic mechanism of slow folding are not changed by the mutation. Trp-59 is in close contact to Pro-39 in native RNase T1 and probably also in an intermediate that forms rapidly during folding. We suggest that this specific interaction interferes with the trans----cis reisomerization of the Tyr-38-Pro-39 bond at the stage of a native-like folding intermediate. The steric hindrance is abolished either by changing Trp-59 to a less bulky residue, such as tyrosine, or, by a destabilization of folding intermediates at increased concentrations of denaturant. Under such conditions folding of the wild-type protein and of the W59Y variant no longer differ. These results provide strong support for the proposal that trans----cis isomerization of Pro-39 is responsible for the major, very slow refolding reaction of RNase T1. They also indicate that specific tertiary interactions in folding intermediates do exist, but do not necessarily facilitate folding. They can have adverse effects and decelerate rate-limiting steps by trapping partially folded structures.     

Extending the folding nucleus of ubiquitin with an independently folding beta-hairpin finger: hurdles to rapid folding arising from the stabilisation of local interactions

The N-terminal beta-hairpin sequence of ubiquitin has been implicated as a folding nucleation site. To extend and stabilise the ubiquitin folding nucleus, we have inserted an autonomously folding 14-residue peptide sequence beta4 which in isolation forms a highly populated beta-hairpin (>70%) stabilised by local interactions. NMR structural analysis of the ubiquitin mutant (Ubeta4) shows that the hairpin finger is fully structured and stabilises ubiquitin by approximately 8kJmol(-1). Protein engineering and kinetic (phi(F)-value) analysis of a series of Ubeta4 mutants shows that the hairpin extension of Ubeta4 is also significantly populated in the transition state (phi(F)-values >0.7) and has the effect of templating the formation of native contacts in the folding nucleus of ubiquitin. However, at low denaturant concentrations the chevron plot of Ubeta4 shows a small deviation from linearity (roll-over effect), indicative of the population of a compact collapsed state, which appears to arise from over-stabilisation of local interactions. Destabilising mutations within the native hairpin sequence and within the engineered hairpin extension, but not elsewhere, eliminate this non-linearity and restore apparent two-state behaviour. The pitfall to stabilising local interactions is to present hurdles to the rapid and efficient folding of small proteins down a smooth folding funnel by trapping partially folded or misfolded states that must unfold or rearrange before refolding.