The folding pathway of the cell-cycle regulatory protein p13suc1: clues for the mechanism of domain swapping

Dissecting the structural and functional consequences of the evolutionary proline-glycine deletion in the wing 1 region of the forkhead domain of human FoxP1

The human FoxP transcription factors dimerize via three-dimensional domain swapping, a unique feature among the human Fox family, as result of evolutionary sequence adaptations in the forkhead domain. This is the case for the conserved glycine and proline residues in the wing 1 region, which are absent in FoxP proteins but present in most of the Fox family. In this work, we engineered both glycine (G) and proline-glycine (PG) insertion mutants to evaluate the deletion events in FoxP proteins in their dimerization, stability, flexibility, and DNA-binding ability. We show that the PG insertion only increases protein stability, whereas the single glycine insertion decreases the association rate and protein stability and promotes affinity to the DNA ligand.

From covalent transition states in chemistry to noncovalent in biology: from β- to Φ-value analysis of protein folding

Solving the mechanism of a chemical reaction requires determining the structures of all the ground states on the pathway and the elusive transition states linking them. 2024 is the centenary of Brønsted's landmark paper that introduced the β-value and structure-activity studies as the only experimental means to infer the structures of transition states. It involves making systematic small changes in the covalent structure of the reactants and analysing changes in activation and equilibrium-free energies. Protein engineering was introduced for an analogous procedure, Φ-value analysis, to analyse the noncovalent interactions in proteins central to biological chemistry. The methodology was developed first by analysing noncovalent interactions in transition states in enzyme catalysis. The mature procedure was then applied to study transition states in the pathway of protein folding - 'part (b) of the protein folding problem'. This review describes the development of Φ-value analysis of transition states and compares and contrasts the interpretation of β- and Φ-values and their limitations. Φ-analysis afforded the first description of transition states in protein folding at the level of individual residues. It revealed the nucleation-condensation folding mechanism of protein domains with the transition state as an expanded, distorted native structure, containing little fully formed secondary structure but many weak tertiary interactions. A spectrum of transition states with various degrees of structural polarisation was then uncovered that spanned from nucleation-condensation to the framework mechanism of fully formed secondary structure. Φ-analysis revealed how movement of the expanded transition state on an energy landscape accommodates the transition from framework to nucleation-condensation mechanisms with a malleability of structure as a unifying feature of folding mechanisms. Such movement follows the rubric of analysis of classical covalent chemical mechanisms that began with Brønsted. Φ-values are used to benchmark computer simulation, and Φ and simulation combine to describe folding pathways at atomic resolution.

Single-molecule optical tweezers reveals folding steps of the domain swapping mechanism of a protein

Domain swapping is a mechanism of protein oligomerization by which two or more subunits exchange structural elements to generate an intertwined complex. Numerous studies support a diversity of swapping mechanisms in which structural elements can be exchanged at different stages of the folding pathway of a subunit. Here, we used single-molecule optical tweezers technique to analyze the swapping mechanism of the forkhead DNA-binding domain of human transcription factor FoxP1. FoxP1 populates folded monomers in equilibrium with a swapped dimer. We generated a fusion protein linking two FoxP1 domains in tandem to obtain repetitive mechanical folding and unfolding trajectories. Thus, by stretching the same molecule several times, we detected either the independent folding of each domain or the elusive swapping step between domains. We found that a swapped dimer can be formed directly from fully or mostly folded monomer. In this situation, the interaction between the monomers in route to the domain-swapped dimer is the rate-limiting step. This approach is a useful strategy to test the different proposed domain swapping mechanisms for proteins with relevant physiological functions.

The protonation state of an evolutionarily conserved histidine modulates domainswapping stability of FoxP1

Forkhead box P (FoxP) proteins are members of the versatile Fox transcription factors, which control the timing and expression of multiple genes for eukaryotic cell homeostasis. Compared to other Fox proteins, they can form domain-swapped dimers through their DNA-binding –forkhead– domains, enabling spatial reorganization of distant chromosome elements by tethering two DNA molecules together. Yet, domain swapping stability and DNA binding affinity varies between different FoxP proteins. Experimental evidence suggests that the protonation state of a histidine residue conserved in all Fox proteins is responsible for pH-dependent modulation of these interactions. Here, we explore the consequences of the protonation state of another histidine (H59), only conserved within FoxM/O/P subfamilies, on folding and dimerization of the forkhead domain of human FoxP1. Dimer dissociation kinetics and equilibrium unfolding experiments demonstrate that protonation of H59 leads to destabilization of the domain-swapped dimer due to an increase in free energy difference between the monomeric and transition states. This pH–dependence is abolished when H59 is mutated to alanine. Furthermore, anisotropy measurements and molecular dynamics evidence that H59 has a direct impact in the local stability of helix H3. Altogether, our results highlight the relevance of H59 in domain swapping and folding stability of FoxP1.

Biophysical analysis of the MHR motif in folding and domain swapping of the HIV capsid protein C-terminal domain

Infection by human immunodeficiency virus (HIV) depends on the function, in virion morphogenesis and other stages of the viral cycle, of a highly conserved structural element, the major homology region (MHR), within the carboxyterminal domain (CTD) of the capsid protein. In a modified CTD dimer, MHR is swapped between monomers. While no evidence for MHR swapping has been provided by structural models of retroviral capsids, it is unknown whether it may occur transiently along the virus assembly pathway. Whatever the case, the MHR-swapped dimer does provide a novel target for the development of anti-HIV drugs based on the concept of trapping a nonnative capsid protein conformation. We have carried out a thermodynamic and kinetic characterization of the domain-swapped CTD dimer in solution. The analysis includes a dissection of the role of conserved MHR residues and other amino acids at the dimerization interface in CTD folding, stability, and dimerization by domain swapping. The results revealed some energetic hotspots at the domain-swapped interface. In addition, many MHR residues that are not in the protein hydrophobic core were nevertheless found to be critical for folding and stability of the CTD monomer, which may dramatically slow down the swapping reaction. Conservation of MHR residues in retroviruses did not correlate with their contribution to domain swapping, but it did correlate with their importance for stable CTD folding. Because folding is required for capsid protein function, this remarkable MHR-mediated conformational stabilization of CTD may help to explain the functional roles of MHR not only during immature capsid assembly but in other processes associated with retrovirus infection. This energetic dissection of the dimerization interface in MHR-swapped CTD may also facilitate the design of anti-HIV compounds that inhibit capsid assembly by conformational trapping of swapped CTD dimers.

The Landscape of Intertwined Associations in Homooligomeric Proteins

We present an overview of the full repertoire of intertwined associations in homooligomeric proteins. This overview summarizes recent findings on the different categories of intertwined associations in known protein structures, their assembly modes, the properties of their interfaces, and their structural plasticity. Furthermore, the current body of knowledge on the so-called three-dimensional domain-swapped systems is reexamined in the context of the wider landscape of intertwined homooligomers, with a particular focus on the mechanistic aspects that underpin intertwined self-association processes in proteins. Insights gained from this integrated overview into the physical and biological roles of intertwining are highlighted.

Regions which are Responsible for Swapping are also Responsible for Folding and Misfolding

Domain swapping is a term used to describe a process when two or more protein chains exchange identical structural elements. Some cases of amyloid formation can be explained through a domain swapping mechanism therefore this deserves theoretical consideration and studying. It has been demonstrated that diverse proteins in sequence and structure are able to oligomerize via domain swapping. This allows us to suggest that the exchangeable regions are important in folding and misfolding processes of proteins, i.e. the residues from the swapping regions are typically incorporated into the native structure early during its formation. The modeling of folding of the proteins with swapped domains demonstrates that the regions exchanged in the oligomeric form in most cases are also responsible for folding and misfolding. For 11 out of 17 proteins, swapping regions intersect with the predicted amyloidogenic regions. Moreover, for 10 out of 17 proteins, high Φ-values (>0.5) belong to residues from the swapping regions. Our data confirm that the exchangeable regions are important in folding, misfolding, and domain swapping processes of the proteins, therefore the suggestion that domain swapping can serve as a mechanism for functional interconversion between monomers and oligomers is likely to be correct.

Folding intermediate and folding nucleus for I-->N and U-->I-->N transitions in apomyoglobin: contributions by conserved and nonconserved residues

Kinetic investigation on the wild-type apomyoglobin and its 12 mutants with substitutions of hydrophobic residues by Ala was performed using stopped-flow fluorescence. Characteristics of the kinetic intermediate I and the folding nucleus were derived solely from kinetic data, namely, the slow-phase folding rate constants and the burst-phase amplitudes of Trp fluorescence intensity. This allowed us to pioneer the phi-analysis for apomyoglobin. As shown, these mutations drastically destabilized the native state N and produced minor (for conserved residues of G, H helices) or even negligible (for nonconserved residues of B, C, D, E helices) destabilizing effect on the state I. On the other hand, conserved residues of A, G, H helices made a smaller contribution to stability of the folding nucleus at the rate-limiting I-->N transition than nonconserved residues of B, D, E helices. Thus, conserved side chains of the A-, G-, H-residues become involved in the folding nucleus before crossing the main barrier, whereas nonconserved side chains of the B-, D-, E-residues join the nucleus in the course of the I-->N transition.

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.

Domain swapping in materials design

Peptide self-assembly can be used as a bottom-up approach to material fabrication. Although many different types of materials can be prepared from peptides, hydrogels are perhaps one of the most common. Gels typically result from the self-assembly of peptides into fibrillar networks. Controlling the structural morphology of these fibrils and the networks they form allows direct control over a given material's bulk properties. However, exerting this control is extremely difficult as the mechanistic rules that govern peptide self-assembly are far from being established. Conversely, several amyloidogenic proteins have been shown to self-assemble into fibrils using a mechanism known as domain swapping. Here, discrete units of secondary structure or even whole domains are exchanged (swapped) among discrete proteins during self-assembly to form extended networks with precise structural control. This review discusses several common mechanistic variations of domain swapping using naturally occurring proteins as examples. The possibility of using these principles to design peptides capable of controlled assembly and fibril formation leading to materials with targeted properties is explored.

Phi-analysis of the folding of the B domain of protein A using multiple optical probes

We examined the co-operativity of ultra-fast folding of a protein and whether the Phi-value analysis of its transition state depended on the location of the optical probe. We incorporated in turn a tryptophan residue into each of the three helices of the B domain of Protein A. Each Trp mutant of the three-helix bundle protein was used as a pseudo-wild-type parent for Phi-analysis in which the intrinsic Trp fluorescence probed the formation of each helix during the transition state. Apart from local effects in the immediate vicinity of the probe, the three separate sets of Phi-values were in excellent agreement, demonstrating the overall co-operativity of folding and the robustness of the Phi-analysis. The transition state of folding of Protein A contains the second helix being well formed with many stabilizing tertiary hydrophobic interactions. In contrast, the first and the third helices are more poorly structured in the transition state. The mechanism of folding thus involves the concurrent formation of secondary and tertiary interactions, and is towards the nucleation-condensation extreme in the nucleation-condensation-framework continuum of mechanism, with helix 2 being the nucleus. We provide an error analysis of Phi-values derived purely from the kinetics of two-state chevron plots.

Folding rate prediction using n-order contact distance for proteins with two- and three-state folding kinetics

It is a challenging task to understand the relationship between sequences and folding rates of proteins. Previous studies are found that one of contact order (CO), long-range order (LRO), total contact distance (TCD), chain topology parameter (CTP), and effective length (Leff) has a significant correlation with folding rate of proteins. In this paper, we introduce a new parameter called n-order contact distance (nOCD) and use it to predict folding rate of proteins with two- and three-state folding kinetics. A good linear correlation between the folding rate logarithm lnkf and nOCD with n=1.2, alpha=0.6 is found for two-state folders (correlation coefficient is -0.809, P-value<0.0001) and n=2.8, alpha=1.5 for three-state folders (correlation coefficient is -0.816, P-value<0.0001). However, this correlation is completely absent for three-state folders with n=1.2, alpha=0.6 (correlation coefficient is 0.0943, P-value=0.661) and for two-state folders with n=2.8, alpha=1.5 (correlation coefficient is -0.235, P-value=0.2116). We also find that the average number of contacts per residue Pm in the interval of m for two-state folders is smaller than that for three-state folders. The probability distribution P(gamma) of residue having gamma pairs of contacts fits a Gaussian distribution for both two- and three-state folders. We observe that the correlations between square radius of gyration S2 and number of residues for two- and three-state folders are both good, and the correlation coefficient is 0.908 and 0.901, and the slope of the fitting line is 1.202 and 0.795, respectively. Maybe three-state folders are more compact than two-state folders. Comparisons with nTCD and nCTP are also made, and it is found that nOCD is the best one in folding rate prediction.

Protein oligomerization through domain swapping: role of inter-molecular interactions and protein concentration

Domain swapping has been shown to be an important mechanism controlling multiprotein assembly and has been suggested recently as a possible mechanism underlying protein aggregation. Understanding oligomerization via domain swapping is therefore of theoretical and practical importance. By using a symmetrized structure-based (Gō) model, we demonstrate that in the free-energy landscape of domain swapping, a large free-energy barrier separates monomeric and domain-swapped dimeric configurations. We investigate the effect of finite monomer concentration, by implementing a new semi-analytical method, which involves computing the second virial coefficient, a thermodynamic indicator of inter-molecular interactions. This method, together with the symmetrized structure-based (Gō) model, minimizes the need for expensive many-protein simulations, providing a convenient framework to investigate concentration effect. Finally, we perform direct simulations of domain-swapped trimer formation, showing that this modeling approach can be used for higher-order oligomers.

Role of Conformational Heterogeneity in Domain Swapping and Adapter Function of the Cks Proteins

Cks proteins are adapter molecules that coordinate the assembly of multiprotein complexes. They share the ability to domain swap by exchanging a beta-strand, beta4. Here we use NMR spectroscopy and molecular dynamics simulations to investigate the dynamic properties of human Cks1 and its response on assembly with components of the SCF(Skp2) ubiquitin ligation machinery. In the NMR experiment with the free form of Cks1, a subset of residues displayed elevated R2 values and the cross-peaks of neighboring residues were missing from the spectrum, indicating a substantial conformational exchange contribution on the microsecond to millisecond time scale. Strikingly the region of greatest conformational variability was the beta4-strand that domain swaps to form the dimer. Binding of the ligand common to all Cks proteins, Cdk2, suppressed the conformational heterogeneity. This response was specific to Cdk2 binding; in contrast, binding of Skp2, a ligand unique to human Cks1, did not alter the dynamic behavior. Short time (<5 ns) molecular dynamics simulations indicate that residues of Cks1 that form the binding site for phosphorylated ligands are considerably more flexible in the free form of Cks1 than they are in the Cdk2-Cks1 complex. A cooperative interaction between Cdk2 and Cks1 is suggested, which reduces the configurational entropy of Cks1 and therefore facilitates phosphoprotein binding. Indications of an unusual dynamic behavior of strand beta4 in the free form of Cks1 were obtained from longer time scale (50 ns) dynamics simulations. A spontaneous reversible unzipping of hydrogen bonds between beta4 and beta2 was observed, suggesting an early intermediate structure for unfolding and/or domain swapping. We propose that the dynamic properties of the beta-sheet and its modification upon ligand binding underlie the domain swapping ability and the adapter function of Cks proteins.

The Folding Mechanism of a Dimeric β-Barrel Domain

The dimeric beta-barrel domain is an unusual topology, shared only by two viral origin binding proteins, where secondary, tertiary and quaternary structure are coupled, and where the dimerization interface is composed of two four-stranded half-beta-barrels. The folding of the DNA binding domain of the E2 transcriptional regulator from human papillomavirus, strain-16, takes place through a stable and compact monomeric intermediate, with 31% the stability of the folded dimeric domain. Double jump multiple wavelength experiments allowed the reconstruction of the fluorescence spectrum of the monomeric intermediate at 100 milliseconds, indicating that tryptophan residues, otherwise buried in the folded state, are accessible to the solvent. Burial of surface area as well as differential behavior to ionic strength and pH with respect to the native ground state, plus the impossibility of having over 2500 A2 of surface area of the half-barrel exposed to the solvent, indicates that the formation of a non-native compact tertiary structure precedes the assembly of native quaternary structure. The monomeric intermediate can dimerize, albeit with a weaker affinity (approximately 1 microM), to yield a non-native dimeric intermediate, which rearranges to the native dimer through a parallel folding channel, with a unimolecular rate-limiting step. Folding pathways from either acid or urea unfolded states are identical, making the folding model robust. Unfolding takes place through a major phase accounting for apparently all the secondary structure change, with identical rate constant to that of the fluorescence unfolding experiment. In contrast to the folding direction, no unfolding intermediate was found.

Computational studies of the reversible domain swapping of p13suc1

A minimalist representation of protein structures using a Go-like potential for interactions is implemented to investigate the mechanisms of the domain swapping of p13suc1, a protein that exists in two native conformations: a monomer and a domain-swapped dimer formed by the exchange of a beta-strand. Inspired by experimental studies which showed a similarity of the transition states for folding of the monomer and the dimer, in this study we justify this similarity in molecular descriptions. When intermediates are populated in the simulations, formation of a domain-swapped dimer initiates from the ensemble of unfolded monomers, given by the fact that the dimer formation occurs at the folding/unfolding temperature of the monomer (T(f)). It is also shown that transitions, leading to a dimer, involve the presence of two intermediates, one of them has a dimeric form and the other is monomeric; the latter is much more populated than the former. However, at temperatures lower than T(f), the population of intermediates decreases. It is argued that the two folded forms may coexist in absence of intermediates at a temperature much lower than T(f). Computational simulations enable us to find a mechanism, "lock-and-dock", for domain swapping of p13suc1. To explore the route toward dimer formation, the folding of unstructured monomers must be retarded by first locking one of the free ends of each chain. Then, the other free termini could follow and dock at particular regions, where most intrachain contacts are formed, and thus define the transition states of the dimer. The simulations also showed that a decrease in the maximum distance between monomers increased their stability, which is explained based on confinement arguments. Although the simulations are based on models extracted from the native structure of the monomer and the dimer of p13suc1, the mechanism of the domain-swapping process could be general, not only for p13suc1.

Non-linear effects of temperature and urea on the thermodynamics and kinetics of folding and unfolding of hisactophilin

Extensive measurements and analysis of thermodynamic stability and kinetics of urea-induced unfolding and folding of hisactophilin are reported for 5-50 degrees C, at pH 6.7. Under these conditions hisactophilin has moderate thermodynamic stability, and equilibrium and kinetic data are well fit by a two-state transition between the native and the denatured states. Equilibrium and kinetic m values decrease with increasing temperature, and decrease with increasing denaturant concentration. The betaF values at different temperatures and urea concentrations are quite constant, however, at about 0.7. This suggests that the transition state for hisactophilin unfolding is native-like and changes little with changing solution conditions, consistent with a narrow free energy profile for the transition state. The activation enthalpy and entropy of unfolding are unusually low for hisactophilin, as is also the case for the corresponding equilibrium parameters. Conventional Arrhenius and Eyring plots for both folding and unfolding are markedly non-linear, but these plots become linear for constant DeltaG/T contours. The Gibbs free energy changes for structural changes in hisactophilin have a non-linear denaturant dependence that is comparable to non-linearities observed for many other proteins. These non-linearities can be fit for many proteins using a variation of the Tanford model, incorporating empirical quadratic denaturant dependencies for Gibbs free energies of transfer of amino acid constituents from water to urea, and changes in fractional solvent accessible surface area of protein constituents based on the known protein structures. Noteworthy exceptions that are not well fit include amyloidogenic proteins and large proteins, which may form intermediates. The model is easily implemented and should be widely applicable to analysis of urea-induced structural transitions in proteins.

Domain swapping is a consequence of minimal frustration

The same energy landscape principles associated with the folding of proteins into their monomeric conformations should also describe how these proteins oligomerize into domain-swapped conformations. We tested this hypothesis by using a simplified model for the epidermal growth factor receptor pathway substrate 8 src homology 3 domain protein, both of whose monomeric and domain-swapped structures have been solved. The model, which we call the symmetrized Gō-type model, incorporates only information regarding the monomeric conformation in an energy function for the dimer to predict the domain-swapped conformation. A striking preference for the correct domain-swapped structure was observed, indicating that overall monomer topology is a main determinant of the structure of domain-swapped dimers. Furthermore, we explore the free energy surface for domain swapping by using our model to characterize the mechanism of oligomerization.

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).

Intermediates Control Domain Swapping during Folding of p13

The 13-kDa protein p13(suc1) has two folded states, a monomer and a structurally similar domain-swapped dimer formed by exchange of a beta-strand. The refolding reaction of p13(suc1) is multiphasic, and in this paper we analyze the kinetics as a function of denaturant and protein concentration and compare the behavior of wild type and a set of mutants previously designed with dimerization propensities that span 9 orders of magnitude. We show that the folding reactions of wild type and all mutants produce the monomer predominantly despite their very different equilibrium behavior. However, the addition of low concentrations of denaturant in the refolding buffer leads to thermodynamic control of the folding reaction with products that correspond to the wild type and mutant equilibrium dimerization propensities. We present evidence that the kinetic control in the absence of urea arises because of the population of the folding intermediates. Intermediates are usually considered to be detrimental to folding because they slow down the reaction; however, our work shows that intermediates buffer the monomeric folding pathway against the effect of mutations that favor the nonfunctional, dimeric state at equilibrium.

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.

Chain length is the main determinant of the folding rate for proteins with three-state folding kinetics

We demonstrate that chain length is the main determinant of the folding rate for proteins with the three-state folding kinetics. The logarithm of their folding rate in water (k(f)) strongly anticorrelates with their chain length L (the correlation coefficient being -0.80). At the same time, the chain length has no correlation with the folding rate for two-state folding proteins (the correlation coefficient is -0.07). Another significant difference of these two groups of proteins is a strong anticorrelation between the folding rate and Baker's "relative contact order" for the two-state folders and the complete absence of such correlation for the three-state folders.

The unfolding story of three-dimensional domain swapping

Three-dimensional domain swapping is the event by which a monomer exchanges part of its structure with identical monomers to form an oligomer where each subunit has a similar structure to the monomer. The accumulating number of observations of this phenomenon in crystal structures has prompted speculation as to its biological relevance. Domain swapping was originally proposed to be a mechanism for the emergence of oligomeric proteins and as a means for functional regulation, but also to be a potentially harmful process leading to misfolding and aggregation. We highlight experimental studies carried out within the last few years that have led to a much greater understanding of the mechanism of domain swapping and of the residue- and structure-specific features that facilitate the process. We discuss the potential biological implications of domain swapping in light of these findings.

Weak cooperativity in the core causes a switch in folding mechanism between two proteins of the cks family

The human protein ckshs1 (cks1) is a 79 residue alpha/beta protein with low thermodynamic and kinetic stability. Its folding mechanism was probed by mutation at sites throughout the structure. Many of the mutations caused changes in the slope of the unfolding arm of the chevron plot. The effects can be rationalised in terms of either transition-state movement or native-state "breathing", and in either case, the magnitude of the effect enables the sequence of events in the folding reaction to be determined. Those sites that fold early exhibit a small perturbation, whilst those sites that fold late exhibit a large perturbation. The results show that cks1 folds sequential pairs of beta-strands first; beta1/beta2 and beta3/beta4. Subsequently, these pairs pack against each other and onto the alpha-helical region to form the core. The folding process of cks1 contrasts with that of the homologue, suc1. The 113 residue suc1 has the same beta-sheet core structure but, additionally, two large insertions that confer much greater thermodynamic and kinetic stability. The more extensive network of tertiary interactions in suc1 provides sufficient enthalpic gain to overcome the entropic cost of forming the core and thus tips the balance in favour of non-local interactions: the non-local, central beta-strand pair, beta2/beta4, forms first and the periphery strands pack on later. Moreover, the greater cooperativity of the core of suc1 protects its folding from perturbation and consequently the slope of the unfolding arm of the chevron plot is much less sensitive to mutation.

Sequential unfolding of ankyrin repeats in tumor suppressor p16

The ANK repeat is a ubiquitous 33-residue motif that adopts a beta hairpin helix-loop-helix fold. Multiple tandem repeats stack in a linear manner to produce an elongated structure that is stabilized predominantly by short-range interactions between residues close in sequence. The tumor suppressor p16(INK4) consists of four repeats and represents the minimal ANK folding unit. We found from Phi value analysis that p16 unfolded sequentially. The two N-terminal ANK repeats, which are distorted from the canonical ANK structure in all INK4 proteins and which are important for functional specificity, were mainly unstructured in the rate-limiting transition state for folding/unfolding, while the two C-terminal repeats were fully formed. A sequential unfolding mechanism could have implications for the cellular fate of wild-type and cancer-associated mutant p16 proteins.

The cross-reactive calcium-binding pollen allergen, Phl p 7, reveals a novel dimer assembly

The timothy grass pollen allergen Phl p 7 assembles most of the IgE epitopes of a novel family of 2 EF-hand calcium-binding proteins and therefore represents a diagnostic marker allergen and vaccine candidate for immunotherapy. Here we report the first three-dimensional structure of a representative of the 2 EF-hand allergen family, Phl p 7, in the calcium-bound form. The protein occurs as a novel dimer assembly with unique features: in contrast to well known EF-hand proteins such as calmodulin, parvalbumin or the S100 proteins, Phl p 7 adopts an extended conformation. Two protein monomers assemble in a head-to-tail arrangement with domain-swapped EF-hand pairing. The intertwined dimer adopts a barrel-like structure with an extended hydrophobic cavity providing a ligand-binding site. Calcium binding acts as a conformational switch between an open and a closed dimeric form of Phl p 7. These findings are interesting in the context of lipid- and calcium-dependent pollen tube growth. Furthermore, the structure of Phl p 7 allows for the rational development of vaccine strategies for treatment of sensitized allergic patients.

Simple physical models connect theory and experiment in protein folding kinetics

Our understanding of the principles underlying the protein-folding problem can be tested by developing and characterizing simple models that make predictions which can be compared to experimental data. Here we extend our earlier model of folding free energy landscapes, in which each residue is considered to be either folded as in the native state or completely disordered, by investigating the role of additional factors representing hydrogen bonding and backbone torsion strain, and by using a hybrid between the master equation approach and the simple transition state theory to evaluate kinetics near the free energy barrier in greater detail. Model calculations of folding phi-values are compared to experimental data for 19 proteins, and for more than half of these, experimental data are reproduced with correlation coefficients between r=0.41 and 0.88; calculations of transition state free energy barriers correlate with rates measured for 37 single domain proteins (r=0.69). The model provides insight into the contribution of alternative-folding pathways, the validity of quasi-equilibrium treatments of the folding landscape, and the magnitude of the Arrhenius prefactor for protein folding. Finally, we discuss the limitations of simple native-state-based models, and as a more general test of such models, provide predictions of folding rates and mechanisms for a comprehensive set of over 400 small protein domains of known structure.

The structure of the transition state for folding of domain-swapped dimeric p13suc1

suc1 has two native states, a monomer and a domain-swapped dimer, in which one molecule exchanges a beta strand with an identical partner. Thus, monomer and dimer have the same structures but are topologically distinct. Importantly, residues that exchange are part of the folding nucleus of the monomer and therefore forming these interactions in the dimer would be expected to incur a large entropic cost. Here we present the transition state for folding/unfolding of domain-swapped dimeric suc1 and compare it with its monomeric counterpart. The same overall structure is observed in the two transition states but the phi values are consistently higher for the domain-swapped dimer. Thus, a greater entropic penalty for bringing together the key interactions in the dimer is overcome by mobilizing more contacts in the transition state, thereby achieving a greater enthalpic gain.

Solution NMR study of the monomeric form of p13suc1 protein sheds light on the hinge region determining the affinity for a phosphorylated substrate

Cyclin-dependent kinase subunit (CKS) proteins bind to cyclin-dependent kinases and target various proteins to phosphorylation and proteolysis during cell division. Crystal structures showed that CKS can exist both in a closed monomeric conformation when bound to the kinase and in an inactive C-terminal beta-strand-exchanged conformation. With the exception of the hinge loop, however, both crystal structures are identical, and no new protein interface is formed in the dimer. Protein engineering studies have pinpointed the crucial role of the proline 90 residue of the p13(suc1) CKS protein from Schizosaccharomyces pombe in the monomer-dimer equilibrium and have led to the concept of a loaded molecular spring of the beta-hinge motif. Mutation of this hinge proline into an alanine stabilizes the protein and prevents the occurrence of swapping. However, other mutations further away from the hinge as well as ligand binding can equally shift the equilibrium between monomer and dimer. To address the question of differential affinity through relief of the strain, here we compare the ligand binding of the monomeric form of wild-type S. pombe p13(suc1) and its hinge mutant P90A in solution by NMR spectroscopy. We indeed observed a 5-fold difference in affinity with the wild-type protein being the most strongly binding. Our structural study further indicates that both wild-type and the P90A mutant proteins adopt in solution the closed conformation but display different dynamic properties in the C-terminal beta-sheet involved in domain swapping and protein interactions.

Protein folding and three-dimensional domain swapping: a strained relationship?

Many proteins function as multimeric assemblies into which the folded individual promoters organize as higher order structures. An oligomerization mechanism that appears to impose the coordination of events during folding and oligomer assembly is three-dimensional domain swapping. Recent studies have focused on revealing the structural basis of domain swapping and a possible role for domain swapping in the regulation of protein aggregation and activity.

Three-dimensional domain swapping in p13suc1 occurs in the unfolded state and is controlled by conserved proline residues

p13suc1 has two native states, a monomer and a domain-swapped dimer. We show that their folding pathways are connected by the denatured state, which introduces a kinetic barrier between monomer and dimer under native conditions. The barrier is lowered under conditions that speed up unfolding, thereby allowing, to our knowledge for the first time, a quantitative dissection of the energetics of domain swapping. The monomer-dimer equilibrium is controlled by two conserved prolines in the hinge loop that connects the exchanging domains. These two residues exploit backbone strain to specifically direct dimer formation while preventing higher-order oligomerization. Thus, the loop acts as a loaded molecular spring that releases tension in the monomer by adopting its alternative conformation in the dimer. There is an excellent correlation between domain swapping and aggregation, suggesting they share a common mechanism. These insights have allowed us to redesign the domain-swapping propensity of suc1 from a fully monomeric to a fully dimeric protein.

Mechanisms of protein folding

The strong correlation between protein folding rates and the contact order suggests that folding rates are largely determined by the topology of the native structure. However, for a given topology, there may be several possible low free energy paths to the native state and the path that is chosen (the lowest free energy path) may depend on differences in interaction energies and local free energies of ordering in different parts of the structure. For larger proteins whose folding is assisted by chaperones, such as the Escherichia coli chaperonin GroEL, advances have been made in understanding both the aspects of an unfolded protein that GroEL recognizes and the mode of binding to the chaperonin. The possibility that GroEL can remove non-native proteins from kinetic traps by unfolding them either during polypeptide binding to the chaperonin or during the subsequent ATP-dependent formation of folding-active complexes with the co-chaperonin GroES has also been explored.

Conversion of two-state to multi-state folding kinetics on fusion of two protein foldons

Chymotrypsin inhibitor 2 (CI2) is the archetypal single-foldon protein that folds in simple two-state kinetics without the accumulation of a folding intermediate. To model the effects of fusion of single foldons to give a multi-foldon protein, we engineered a "double-CI2" protein, in which another CI2 polypeptide was inserted into the loop region of the parent CI2. CD and HSQC spectra demonstrated that while the double-CI2 protein adopted two kinds of native conformations, CI2-like structure was almost preserved in both the domains of double-CI2. In the folding kinetic studies, double-CI2 exhibited a remarkable rollover of the observed folding rates at low denaturant concentrations, indicating that double-CI2 accumulated a kinetic folding intermediate. The different folding mechanisms between WT-CI2 and double-CI2 support the present view that protein size or number of domains is an important determinant for formation of folding intermediates.

Sequence conservation provides the best prediction of the role of proline residues in p13suc1

The unique nature of the proline side-chain imposes severe constraints on the polypeptide backbone, and thus it seems likely that it plays a special structural or functional role in the architecture of proteins. We have investigated the role of proline residues in suc1, a member of the cyclin-dependent kinase (cks) family of proteins, whose known function is to bind to and regulate the activity of the major mitotic cdk. The effect on stability of mutation to alanine of all but two of the eight proline residues is correlated with their conservation within the family. The remaining two proline residues are located in the hinge loop between two beta-strands that mediates a domain-swapping process involving exchange of a beta-strand between two monomers to form a dimer pair. Mutation of these proline residues to alanine stabilises the protein. cdk binding is unaffected by these mutations, but dimerisation is altered. We propose, therefore, that the double-proline motif is conserved for the purpose of domain swapping, which suggests that this phenomenon plays a role in the function of cks proteins. Thus, the conservation of the proline residues is a good indicator of their roles in suc1, either in the stabilisation of the native state or in performing functions that are as yet unknown. In addition, the strain resulting from two of the proline residues was relieved successfully by mutation of the preceeding residue to glycine, suggesting a general method for designing more stable proteins.