Searching for multiple folding pathways of a nearly symmetrical protein: temperature dependent phi-value analysis of the B domain of protein A

The B domain of protein A (BdpA) is a popular paradigm for simulating protein folding pathways. The discrepancies between so many simulations and subsequent experimental testing may be attributable to the protein being highly symmetrical: changing experimental conditions could perturb the subtle interplay between the effects of symmetry in the native structure and the effects of asymmetry from specific interactions in a given sequence. If the protein folds via multiple pathways, perturbations, such as temperature, denaturant concentration, and mutation, should change the flux of micro pathways, leading to changes in the bulk properties of the transition state. We tested this hypothesis by conducting a Phi-analysis of BdpA as a function of temperature from 25.0 degrees C to 60.0 degrees C. The Phi-values had no significant dependence on temperature and the values at 55.0 degrees C (denaturing conditions) are very similar to those at 25.0 degrees C (folding conditions), indicating the structure of the transition state does not significantly change although the experimental conditions are considerably altered. The results suggest that BdpA folds via a single dominant folding pathway.     

Fast and faster: a designed variant of the B-domain of protein A folds in 3 microsec

We have introduced the mutation glycine 29 to alanine, designed to increase the rate of protein folding, into the B-domain of protein A (BdpA). From NMR lineshape analysis, we find the G29A mutation increases the folding rate constant by threefold; the folding time is 3 microsec. Although wild-type BdpA folds extremely fast, simple-point mutations can still speed up the folding; thus, the folding rate is not evolutionarily maximized. The short folding time of G29A BdpA (the shortest time yet reported) makes it an attractive candidate for an all-atom molecular dynamics simulation that could potentially show a complete folding reaction starting from an extended chain. We also constructed a fluorescent variant of BdpA by mutating phenylalanine 13 to tryptophan, allowing fluorescence-based time-resolved temperature-jump measurements. Temperature jumps and NMR complement each other, and give a very complete picture of the folding kinetics.     

Phage P22 tailspike protein: removal of head-binding domain unmasks effects of folding mutations on native-state thermal stability.

A shortened, recombinant protein comprising residues 109-666 of the tailspike endorhamnosidase of Salmonella phage P22 was purified from Escherichia coli and crystallized. Like the full-length tailspike, the protein lacking the amino-terminal head-binding domain is an SDS-resistant, thermostable trimer. Its fluorescence and circular dichroism spectra indicate native structure. Oligosaccharide binding and endoglycosidase activities of both proteins are identical. A number of tailspike folding mutants have been obtained previously in a genetic approach to protein folding. Two temperature-sensitive-folding (tsf) mutations and the four known global second-site suppressor (su) mutations were introduced into the shortened protein and found to reduce or increase folding yields at high temperature. The mutational effects on folding yields and subunit folding kinetics parallel those observed with the full-length protein. They mirror the in vivo phenotypes and are consistent with the substitutions altering the stability of thermolabile folding intermediates. Because full-length and shortened tailspikes aggregate upon thermal denaturation, and their denaturant-induced unfolding displays hysteresis, kinetics of thermal unfolding were measured to assess the stability of the native proteins. Unfolding of the shortened wild-type protein in the presence of 2% SDS at 71 degrees C occurs at a rate of 9.2 x 10(-4) s(-1). It reflects the second kinetic phase of unfolding of the full-length protein. All six mutations were found to affect the thermal stability of the native protein. Both tsf mutations accelerate thermal unfolding about 10-fold. Two of the su mutations retard thermal unfolding up to 5-fold, while the remaining two mutations accelerate unfolding up to 5-fold. The mutational effects can be rationalized on the background of the recently determined crystal structure of the protein.     

Evidences for the unfolding mechanism of three‐dimensional domain swapping

The full or partial unfolding of proteins is widely believed to play an essential role in three‐dimensional domain swapping. However, there is little research that has rigorously evaluated the association between domain swapping and protein folding/unfolding. Here, we examined a kinetic model in which domain swapping occurred via the denatured state produced by the complete unfolding of proteins. The relationships between swapping kinetics and folding/unfolding thermodynamics were established, which were further adopted as criteria to show that the proposed mechanism dominates in three representative proteins: Cyanovirin‐N (CV‐N), the C‐terminal domain of SARS‐CoV main protease (M^pro‐C), and a single mutant of oxidized thioredoxin (Trx_W28A^ox).     

Generation of a consensus protein domain dictionary

MOTIVATION: The discovery of new protein folds is a relatively rare occurrence even as the rate of protein structure determination increases. This rarity reinforces the concept of folds as reusable units of structure and function shared by diverse proteins. If the folding mechanism of proteins is largely determined by their topology, then the folding pathways of members of existing folds could encompass the full set used by globular protein domains. RESULTS: We have used recent versions of three common protein domain dictionaries (SCOP, CATH and Dali) to generate a consensus domain dictionary (CDD). Surprisingly, 40% of the metafolds in the CDD are not composed of autonomous structural domains, i.e. they are not plausible independent folding units. This finding has serious ramifications for bioinformatics studies mining these domain dictionaries for globular protein properties. However, our main purpose in deriving this CDD was to generate an updated CDD to choose targets for MD simulation as part of our dynameomics effort, which aims to simulate the native and unfolding pathways of representatives of all globular protein consensus folds (metafolds). Consequently, we also compiled a list of representative protein targets of each metafold in the CDD. Availability and implementation: This domain dictionary is available at www.dynameomics.org.     

Unusual cold denaturation of a small protein domain

A thermal unfolding study of the 45-residue α-helical domain UBA(2) using circular dichroism is presented. The protein is highly thermostable and exhibits a clear cold unfolding transition with the onset near 290 K without denaturant. Cold denaturation in proteins is rarely observed in general and is quite unique among small helical protein domains. The cold unfolding was further investigated in urea solutions, and a simple thermodynamic model was used to fit all thermal and urea unfolding data. The resulting thermodynamic parameters are compared to those of other small protein domains. Possible origins of the unusual cold unfolding of UBA(2) are discussed.     

Exhaustive enumeration of protein domain families

Domains are considered as the basic units of protein folding, evolution, and function. Decomposing each protein into modular domains is thus a basic prerequisite for accurate functional classification of biological molecules. Here, we present ADDA, an automatic algorithm for domain decomposition and clustering of all protein domain families. We use alignments derived from an all-on-all sequence comparison to define domains within protein sequences based on a global maximum likelihood model. In all, 90% of domain boundaries are predicted within 10% of domain size when compared with the manual domain definitions given in the SCOP database. A representative database of 249,264 protein sequences were decomposed into 450,462 domains. These domains were clustered on the basis of sequence similarities into 33,879 domain families containing at least two members with less than 40% sequence identity. Validation against family definitions in the manually curated databases SCOP and PFAM indicates almost perfect unification of various large domain families while contamination by unrelated sequences remains at a low level. The global survey of protein-domain space by ADDA confirms that most large and universal domain families are already described in PFAM and/or SMART. However, a survey of the complete set of mobile modules leads to the identification of 1479 new interesting domain families which shuffle around in multi-domain proteins. The data are publicly available at ftp://ftp.ebi.ac.uk/pub/contrib/heger/adda.     

Rapid cooperative two-state folding of a miniature alpha-beta protein and design of a thermostable variant

There is a great deal of interest in developing small stably folded miniature proteins. A limited number of these molecules have been described, however they typically have not been characterized in depth. In particular, almost no detailed studies of the thermodynamics and folding kinetics of these proteins have been reported. Here we describe detailed studies of the thermodynamics and kinetics of folding of a 39 residue mixed alpha-beta protein (NTL9(1-39)) derived from the N-terminal domain of the ribosomal protein L9. The protein folds cooperatively and rapidly in a two-state fashion to a native state typical of those found for normal globular proteins. At pH 5.4 in 20mM sodium acetate, 100mM NaCl the temperature of maximum stability is 6 degrees C, the t(m) is 65.3 degrees C, deltaH degrees (t(m)) is between 24.6 kcalmol(-1) and 26.3 kcalmol(-1), and deltaC(p) degrees is 0.38 kcalmol(-1)deg(-1). The thermodynamic parameters are in the range expected on the basis of per residue values determined from databases of globular proteins. H/2H exchange measurements reveal a set of amides that exchange via global unfolding, exactly as expected for a normal cooperatively folded globular protein. Kinetic measurements show that folding is two-state folding. The folding rate is 640 s(-1) and the value of deltaG degrees calculated from the folding and unfolding rates is in excellent agreement with the equilibrium value. A designed thermostable variant, generated by mutating K12 to M, was characterized and found to have a t(m) of 82 degrees C. Equilibrium and kinetic measurements demonstrate that its folding is cooperative and two-state.     

Complex folding kinetics of a multidomain protein

Spectrin domains are three-helix bundles, commonly found in large tandem arrays. Equilibrium studies have shown that spectrin domains are significantly stabilized by their neighbors. In this work we show that domain:domain interactions can also have profound effects on their kinetic behavior. We have studied the folding of a tandem pair of spectrin domains (R1617) using a combination of single- and double-jump stopped flow experiments (monitoring folding by both circular dichroism and fluorescence). Mutant proteins were also used to investigate the complex folding kinetics. We find that, although the domains fold and unfold individually, there is a single rate-determining step for both folding and unfolding of the protein. This is consistent with the equilibrium observation of cooperative folding of the entire two-domain protein. The results may have important biological implications. Not only will the protein fold more efficiently during cotranslational folding, but the ability of the multidomain protein to withstand thermal unfolding in the cell will be dramatically increased. This study suggests that caution has to be exercised when extrapolating from single domains to larger proteins with a number of independently folding modules arranged in tandem. The multidomain protein spectrin is certainly more than "the sum of its parts".     

Understanding the folding and stability of a designed WW domain protein with replica exchange molecular dynamics simulations

WW domain proteins are usually regarded as simple models for understanding the folding mechanism of β-sheet. CC45 is an artificial protein that is capable of folding into the same structure as WW domain. In this article, the replica exchange molecular dynamics simulations are performed to investigate the folding mechanism of CC45. The analysis of thermal stability shows that β-hairpin 1 is more stable than β-hairpin 2 during the unfolding process. Free energy analysis shows that the unfolding of this protein substantially proceeds through solvating the smaller β-hairpin 2, followed by the unfolding of β-hairpin 1. We further propose the unfolding process of CC45 and the folding mechanism of two β-hairpins. These results are similar to the previous folding studies of formin binding protein 28 (FBP28). Compared with FBP28, it is found that CC45 has more aromatic residues in N-terminal loop, and these residues contact with C-terminal loop to form the outer hydrophobic core, which increases the stability of CC45. Knowledge about the stability and folding behaviour of CC45 may help in understanding the folding mechanisms of the β-sheet and in designing new WW domains.     

Asymmetric effect of domain interactions on the kinetics of folding in yeast phosphoglycerate kinase

The aim of this work is to shed more light on the effect of domain-domain interactions on the kinetics and the pathway of protein folding. A model protein system consisting of several single-tryptophan variants of the two-domain yeast phosphoglycerate kinase (PGK) and its individual domains was studied. Refolding was initiated from the guanidine-unfolded state by stopped-flow or manual mixing and monitored by tryptophan fluorescence from 1 msec to 1000 sec. Denaturant titrations of both individual domains showed apparent two-state unfolding transitions. Refolding kinetics of the individual domains from different denaturant concentrations, however, revealed the presence of intermediate structures during titration for both domains. Refolding of the same domains within the complete protein showed that domain-domain interactions direct the folding of both domains, but in an asymmetric way. Folding of the N domain was already altered within 1 msec, while detectable changes in the folding of the C domain occurred only 60-100 msec after initiating refolding. All mutants showed a hyperfluorescent kinetic intermediate. Both the disappearance of this intermediate and the completion of the folding were significantly faster in the individual N domain than in the complete protein. On the contrary, folding of the individual C domain was slower than in the complete protein. The presence of the C domain directs the refolding of the N domain along a completely different pathway than that of the individual N domain, while folding of the individual C domain follows the same path as within the complete protein.     

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

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

The native state conformational ensemble of the SH3 domain from alpha-spectrin.

The folding/unfolding equilibrium of the alpha-spectrin SH3 domain has been measured by NMR-detected hydrogen/deuterium exchange and by differential scanning calorimetry. Protection factors against exchange have been obtained under native conditions for more than half of the residues in the domain. Most protected residues are located at the beta-strands, the short 3(10) helix, and part of the long RT loop, whereas the loops connecting secondary structure elements show no measurable protection. Apparent stability constants per residue and their corresponding Gibbs energies have been calculated from the exchange experiments. The most stable region of the SH3 domain is defined by the central portions of the beta-strands. The peptide binding region, on the other hand, is composed of a highly stable region (residues 53-57) and a highly unstable region, the loop between residues 34-41 (n-Src loop). All residues in the domain have apparent Gibbs energies lower than the global unfolding Gibbs energy measured by differential scanning calorimetry, indicating that under our experimental conditions the amide exchange of all residues in the SH3 domain occurs primarily via local unfolding reactions. A structure-based thermodynamic analysis has allowed us to predict correctly the thermodynamics of the global unfolding of the domain and to define the ensemble of conformational states that quantitatively accounts for the observed pattern of hydrogen exchange protection. These results demonstrate that under native conditions the SH3 domain needs to be considered as an ensemble of conformations and that the hydrogen exchange data obtained under those conditions cannot be interpreted by a two-state equilibrium. The observation that specific regions of a protein are able to undergo independent local folding/unfolding reactions indicates that under native conditions the scale of cooperative interactions is regional rather than global.     

Analysis of the thermodynamics of binding of an SH3 domain to proline-rich peptides using a chimeric fusion protein

A complete understanding of the thermodynamic determinants of binding between SH3 domains and proline-rich peptides is crucial to the development of rational strategies for designing ligands for these important domains. Recently we engineered a single-chain chimeric protein by fusing the α-spectrin Src homology region 3 (SH3) domain to the decapeptide APSYSPPPPP (p41). This chimera mimics the structural and energetic features of the interaction between SH3 domains and proline-rich peptides. Here we show that analysing the unfolding thermodynamics of single-point mutants of this chimeric fusion protein constitutes a very useful approach to deciphering the thermodynamics of SH3-ligand interactions. To this end, we investigated the contribution of each proline residue of the ligand sequence to the SH3-peptide interaction by producing six single Pro-Ala mutants of the chimeric protein and analysing their unfolding thermodynamics by differential scanning calorimetry (DSC). Structural analyses of the mutant chimeras by circular dichroism, fluorescence and NMR together with NMR-relaxation measurements indicate conformational flexibility at the binding interface, which is strongly affected by the different Pro-Ala mutations. An analysis of the DSC thermograms on the basis of a three-state unfolding model has allowed us to distinguish and separate the thermodynamic magnitudes of the interaction at the binding interface. The model assumes equilibrium between the “unbound” and “bound” states at the SH3-peptide binding interface. The resulting thermodynamic magnitudes classify the different proline residues according to their importance in the interaction as P2∼P7∼P10 > P9∼P6 > P8, which agrees well with Lim's model for the interaction between SH3 domains and proline-rich peptides. In addition, the thermodynamic signature of the interaction is the same as that usually found for this type of binding, with a strong enthalpy-entropy compensation for all the mutants. This compensation appears to derive from an increase in conformational flexibility concomitant to the weakening of the interactions at the binding interface. We conclude that our approach, based on DSC and site-directed mutagenesis analysis of chimeric fusion proteins, may serve as a suitable tool to analyse the energetics of weak biomolecular interactions such as those involving SH3 domains.     

p23 and HSP20/alpha-crystallin proteins define a conserved sequence domain present in other eukaryotic protein families

We identified families of proteins characterized by the presence of a domain similar to human p23 protein, which include proteins such as Sgt1, involved in the yeast kinetochore assembly; melusin, involved in specific interactions with the cytoplasmic integrin beta1 domain; Rar1, related to pathogenic resistance in plants, and to development in animals; B5+B5R flavo-hemo cytochrome NAD(P)H oxidoreductase type B in humans and mice; and NudC, involved in nucleus migration during mitosis. We also found that p23 and the HSP20/alpha-crystallin family of heat shock proteins, which share the same three-dimensional folding, show a pattern of conserved residues that points to a common origin in the evolution of both protein domains. The p23 and HSP20/alpha-crystallin phylogenetic relationship and their similar role in chaperone activity suggest a common function, probably involving protein-protein interaction, for those proteins containing p23-like domains.     

Validation of helical tilt angles in the solution NMR structure of the Z domain of Staphylococcal protein A by combined analysis of residual dipolar coupling and NOE data

Staphylococcal protein A (SpA) is a virulence factor from Staphylococcus aureus that is able to bind to immunoglobulins. The 3D structures of its immunoglobulin (Ig) binding domains have been extensively studied by NMR and X-ray crystallography, and are often used as model structures in developing de novo or ab initio strategies for predicting protein structure. These small three-helix-bundle structures, reported in free proteins or Ig-bound complexes, have been determined previously using medium- to high-resolution data. Although the location and relative orientation of the three helices in most of these published 3D domain structures are consistent, there are significant differences among the reported structures regarding the tilt angle of the first helix (helix 1). We have applied residual dipolar coupling data, together with nuclear Overhauser enhancement and scalar coupling data, in refining the NMR solution structure of an engineered IgG-binding domain (Z domain) of SpA. Our results demonstrate that the three helices are almost perfectly antiparallel in orientation, with the first helix tilting slightly away from the other two helices. We propose that this high-accuracy structure of the Z domain of SpA is a more suitable target for theoretical predictions of the free domain structure than previously published lower-accuracy structures of protein A domains.     

Thermodynamic analysis of the folding of the streptococcal protein G IgG-binding domains B1 and B2: why small proteins tend to have high denaturation temperatures.

We have cloned, expressed, and characterized two naturally occurring variations of the IgG-binding domain of streptococcal protein G. The domain is a stable cooperative folding unit of 56 amino acids, which maintains a unique folded structure without disulfide cross-links or tight ligand binding. We have studied the thermodynamics of the unfolding reaction for the two versions of this domain, designated B1 and B2, which differ by six amino acids. They have denaturation temperatures of 87.5 degrees C and 79.4 degrees C, respectively at pH 5.4, as determined by differential scanning calorimetry. Thermodynamic state functions for the unfolding reaction (delta G, delta H, delta S, and delta Cp) have been determined and reveal several interesting insights into the behavior of very small proteins. First, though the B1 domain has a heat denaturation point close to 90 degrees C, it is not unusually stable at physiologically relevant temperatures (delta G = 25 kJ/mol at 37 degrees C). This behavior occurs because the stability profile (delta G vs temperature) is flat and shallow due to the small delta S and delta Cp for unfolding. Related to this point is the second observation that small changes in the free energy of unfolding of the B-domain due to mutation or change in solvent conditions lead to large shifts in the heat denaturation temperature. Third, the magnitude and relative contributions of hydrophobic vs nonhydrophobic forces (per amino acid residue) to the total free energy of folding of the B-domain are remarkably typical of other globular proteins of much larger size.     

An Overlapping Region between the Two Terminal Folding Units of the Outer Surface Protein A (OspA) Controls Its Folding Behavior

Although many naturally occurring proteins consist of multiple domains, most studies on protein folding to date deal with single-domain proteins or isolated domains of multi-domain proteins. Studies of multi-domain protein folding are required for further advancing our understanding of protein folding mechanisms. Borrelia outer surface protein A (OspA) is a β-rich two-domain protein, in which two globular domains are connected by a rigid and stable single-layer β-sheet. Thus, OspA is particularly suited as a model system for studying the interplays of domains in protein folding. Here, we studied the equilibria and kinetics of the urea-induced folding–unfolding reactions of OspA probed with tryptophan fluorescence and ultraviolet circular dichroism. Global analysis of the experimental data revealed compelling lines of evidence for accumulation of an on-pathway intermediate during kinetic refolding and for the identity between the kinetic intermediate and a previously described equilibrium unfolding intermediate. The results suggest that the intermediate has the fully native structure in the N-terminal domain and the single layer β-sheet, with the C-terminal domain still unfolded. The observation of the productive on-pathway folding intermediate clearly indicates substantial interactions between the two domains mediated by the single-layer β-sheet. We propose that a rigid and stable intervening region between two domains creates an overlap between two folding units and can energetically couple their folding reactions.     

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.