Folding mechanisms of proteins with high sequence identity but different folds

The problem of how a protein folds from a linear chain of amino acids to the three-dimensional structure necessary for function is often investigated using proteins with a low degree of sequence identity that adopt different folds. The design of pairs of proteins with a high degree of sequence identity but different folds offers the opportunity for a complementary study; in two highly similar sequences, which residues are the most important in directing folding to a particular structure? Here we use molecular dynamics simulations to characterize the folding-unfolding pathways of a pair of proteins designed by Bryan and co-workers [Alexander, P. A., et al. (2005) Biochemistry 44, 14045-14054; He, Y. N., et al. (2005) Biochemistry 44, 14055-14061]. Despite being 59% identical, the two protein sequences fold to two different structures. The first sequence folds to the alpha+beta protein G structure and the second to the all-alpha-helical protein A structure. We show that the final protein structure is determined early along the folding pathway. In folding to the protein G structure, the single alpha-helix (alpha1) and the beta3-beta4 turn fold early. Formation of the hairpin turn essentially prevents folding to helical structure in this region of the protein. This early structure is then consolidated by formation of long-range hydrophobic interactions between alpha1 and the beta3-beta4 turn. The protein A sequence differs both in the residues that form the beta3-beta4 turn and also in many of the residues that form the early hydrophobic interactions in the protein G structure. Instead, in the protein A sequence, a more hierarchical mechanism is observed, with helices folding before many of the tertiary interactions are formed. We find that small, but critical, sequence differences determine the topology of the protein early along the folding pathway, which help to explain the process by which one fold can evolve into another.     

Conserved patterns and interactions in the unfolding transition state across SH3 domain structural homologues

Proteins with similar structures are generally assumed to arise from similar sequences. However, there are more cases than not where this is not true. The dogma is that sequence determines structure; how, then, can very different sequences fold to the same structure? Here, we employ high temperature unfolding simulations to probe the pathways and specific interactions that direct the folding and unfolding of the SH3 domain. The SH3 metafold in the Dynameomics Database consists of 753 proteins with the same structure, but varied sequences and functions. To investigate the relationship between sequence and structure, we selected 17 targets from the SH3 metafold with high sequence variability. Six unfolding simulations were performed for each target, transition states were identified, revealing two general folding/unfolding pathways at the transition state. Transition states were also expressed as mathematical graphs of connected chemical nodes, and it was found that three positions within the structure, independent of sequence, were consistently more connected within the graph than any other nearby positions in the sequence. These positions represent a hub connecting different portions of the structure. Multiple sequence alignment and covariation analyses also revealed certain positions that were more conserved due to packing constraints and stabilizing long‐range contacts. This study demonstrates that members of the SH3 domain with different sequences can unfold through two main pathways, but certain characteristics are conserved regardless of the sequence or unfolding pathway. While sequence determines structure, we show that disparate sequences can provide similar interactions that influence folding and lead to similar structures.     

Protein sequence and structure relationship ARMA spectral analysis: application to membrane proteins.

If it is assumed that the primary sequence determines the three-dimensional folded structure of a protein, then the regular folding patterns, such as alpha-helix, beta-sheet, and other ordered patterns in the three-dimensional structure must correspond to the periodic distribution of the physical properties of the amino acids along the primary sequence. An AutoRegressive Moving Average (ARMA) model method of spectral analysis is applied to analyze protein sequences represented by the hydrophobicity of their amino acids. The results for several membrane proteins of known structures indicate that the periodic distribution of hydrophobicity of the primary sequence is closely related to the regular folding patterns in a protein's three-dimensional structure. We also applied the method to the transmembrane regions of acetylcholine receptor alpha subunit and Shaker potassium channel for which no atomic resolution structure is available. This work is an extension of our analysis of globular proteins by a similar method.     

Benchmarking of TASSER in the ab initio limit

A significant number of protein sequences in a given proteome have no obvious evolutionarily related protein in the database of solved protein structures, the PDB. Under these conditions, ab initio or template-free modeling methods are the sole means of predicting protein structure. To assess its expected performance on proteomes, the TASSER structure prediction algorithm is benchmarked in the ab initio limit on a representative set of 1129 nonhomologous sequences ranging from 40 to 200 residues that cover the PDB at 30% sequence identity and which adopt alpha, alpha + beta, and beta secondary structures. For sequences in the 40-100 (100-200) residue range, as assessed by their root mean square deviation from native, RMSD, the best of the top five ranked models of TASSER has a global fold that is significantly close to the native structure for 25% (16%) of the sequences, and with a correct identification of the structure of the protein core for 59% (36%). In the absence of a native structure, the structural similarity among the top five ranked models is a moderately reliable predictor of folding accuracy. If we classify the sequences according to their secondary structure content, then 64% (36%) of alpha, 43% (24%) of alpha + beta, and 20% (12%) of beta sequences in the 40-100 (100-200) residue range have a significant TM-score (TM-score > or = 0.4). TASSER performs best on helical proteins because there are less secondary structural elements to arrange in a helical protein than in a beta protein of equal length, since the average length of a helix is longer than that of a strand. In addition, helical proteins have shorter loops and dangling tails. If we exclude these flexible fragments, then TASSER has similar accuracy for sequences containing the same number of secondary structural elements, irrespective of whether they are helices and/or strands. Thus, it is the effective configurational entropy of the protein that dictates the average likelihood of correctly arranging the secondary structure elements.     

Identification and characterization of key substructures involved in the early folding events of a (beta/alpha)8-barrel protein as studied by experimental and computational methods

A number of studies have examined the structural properties of late folding intermediates of (beta/alpha)8-barrel proteins involved in tryptophan biosynthesis, whereas there is little information available about the early folding events of these proteins. To identify the contiguous polypeptide segments important to the folding of the (beta/alpha)8-barrel protein Escherichia coli N-(5'-phosphoribosyl)anthranilate isomerase, we structurally characterized fragments and circularly permuted forms of the protein. We also simulated thermal unfolding of the protein using molecular dynamics. Our fragmentation experiments demonstrate that the isolated (beta/alpha)(1-4)beta5 fragment is almost as stable as the full-length protein. The far and near-UV CD spectra of this fragment are indicative of native-like secondary and tertiary structures. Structural analysis of the circularly permutated proteins shows that if the protein is cleaved within the two N-terminal betaalpha modules, the amount of secondary structure is unaffected, whereas, when cleaved within the central (beta/alpha)(3-4)beta5 segment, the protein simply cannot fold. An ensemble of the denatured structures produced by thermal unfolding simulations contains a persistent local structure comprised of beta3, beta4 and beta5. The presence of this three-stranded beta-barrel suggests that it may be an important early-stage folding intermediate. Interactions found in (beta/alpha)(3-4)beta5 may be essential for the early events of ePRAI folding if they provide a nucleation site that directs folding.     

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.     

Solution NMR structures of IgG binding domains with artificially evolved high levels of sequence identity but different folds

We describe here the solution NMR structures of two IgG binding domains with highly homologous sequences but different three-dimensional structures. The proteins, G311 and A219, are derived from the IgG binding domains of their wild-type counterparts, protein G and protein A, respectively. Through a series of site-directed mutations and phage display selections, the sequences of G311 and A219 were designed to converge to a point of high-level sequence identity while keeping their respective wild-type tertiary folds. Structures of both artificially evolved sequences were determined by NMR spectroscopy. The main chain fold of G311 can be superimposed on the wild-type alpha/beta protein G structure with a backbone rmsd of 1.4 A, and the A219 structure can be overlaid on the wild-type three-alpha-helix protein A fold also with a backbone rmsd of 1.4 A. The structure of G311, in particular, accommodates a large number of mutational changes without undergoing a change in the overall fold of the main chain. The structural differences are maintained despite a high level (59%) of sequence identity. These proteins serve as starting points for further experiments that will probe basic concepts of protein folding and conformational switching.     

Effect of single amino acid replacements on the folding and stability of dihydrofolate reductase from Escherichia coli

The role of the secondary structure in the folding mechanism of dihydrofolate reductase from Escherichia coli was probed by studying the effects of amino acid replacements in two alpha helices and two strands of the central beta sheet on the folding and stability. The effects on stability could be qualitatively understood in terms of the X-ray structure for the wild-type protein by invoking electrostatic, hydrophobic, or hydrogen-bonding interactions. Kinetic studies focused on the two slow reactions that are thought to reflect the unfolding/refolding of two stable native conformers to/from their respective folding intermediates [Touchette, N. A., Perry, K. M., & Matthews, C. R. (1986) Biochemistry 25, 5445-5452]. Replacements at three different positions in helix alpha B selectively alter the relaxation time for unfolding while a single replacement in helix alpha C selectively alters the relaxation time for refolding. This behavior is characteristic of mutations that change the stability of the protein but do not affect the rate-limiting step. In striking contrast, replacements in strands beta F and beta G can affect both unfolding and refolding relaxation times. This behavior shows that these mutations alter the rate-limiting step in these native-to-intermediate folding reactions. It is proposed that the intermediates have an incorrectly formed beta sheet whose maturation to the structure found in the native conformation is one of the slow steps in folding.     

Structural homology of lens crystallins. II. Homology expressed by correlation coefficients and hydropathy profiles

Bovine lens alpha A- and alpha B-crystallin polypeptides show extensive sequence homology with each other, but apparently none with beta Bp- and gamma 2-crystallin. Despite only 30% sequence homology, the latter two proteins are assumed to have a strong correspondence in tertiary structure, consisting of four structurally similar folding units of antiparallel beta-sheet. We have tested for internal structural repeats in all crystallins, and structural homology between crystallins, by comparing various physical properties of the amino acid residues, such as bulkiness and propensity to form beta-sheet and beta-turn structure. Two procedures used a combination of five physical parameters to calculate correlation coefficients. The 4-fold structural repeat in gamma 2-crystallin and the internal duplication in beta Bp-crystallin were readily detectable, as was also the strong structural homology between corresponding folding units in beta Bp- and gamma 2-crystallin. However, for alpha-crystallin polypeptides, no conclusive support was obtained for either a four-unit or a six-unit folding, the two models previously considered by us. The third procedure compared smoothened hydropathy plots, representing hydrophilic and hydrophobic regions along the polypeptide sequences. Hydropathy profiles were found to show strong correspondence, particularly between alpha B-crystallin and beta Bp-crystallin. These observations support a similar 4-fold folding pattern for all bovine crystallins. A possible role in subunit interactions of the N-terminal folding unit, which has hydrophobic surface characteristics in both alpha- and beta-crystallin polypeptides, is proposed.     

A detailed unfolding pathway of a (beta/alpha)8-barrel protein as studied by molecular dynamics simulations

The (beta/alpha)(8)-barrel is the most common protein fold. Similar structural properties for folding intermediates of (beta/alpha)(8)-barrel proteins involved in tryptophan biosynthesis have been reported in a number of experimental studies; these intermediates have the last two beta-strands and three alpha-helices partially unfolded, with other regions of the polypeptide chain native-like in conformation. To investigate the detailed folding/unfolding pathways of these (beta/alpha)(8)-barrel proteins, temperature-induced unfolding simulations of N-(5'-phosphoribosyl)anthranilate isomerase from Escherichia coli were carried out using a special-purpose parallel computer system. Unfolding simulations at five different temperatures showed a sequential unfolding pathway comprised of several events. Early events in unfolding involved disruption of the last two strands and three helices, producing an intermediate ensemble similar to those detected in experimental studies. Then, denaturation of the first two betaalpha units and separation of the sixth strand from the fifth took place independently. The remaining central betaalphabetaalphabeta module persisted the longest during all simulations, suggesting an important role for this module as the incipient folding scaffold. Our simulations also predicted the presence of a nucleation site, onto which several hydrophobic residues condensed forming the foundation for the central betaalphabetaalphabeta module.     

Beta-sheet proteins with nearly identical structures have different folding intermediates

The folding mechanisms of two proteins in the family of intracellular lipid binding proteins, ileal lipid binding protein (ILBP) and intestinal fatty acid binding protein (IFABP), were examined. The structures of these all-beta-proteins are very similar, with 123 of the 127 amino acids of ILBP having backbone and C(beta) conformations nearly identical to those of 123 of the 131 residues of IFABP. Despite this structural similarity, the sequences of these proteins have diverged, with 23% sequence identity and an additional 16% sequence similarity. The folding process was completely reversible, and no significant concentrations of intermediates were observed by circular dichroism or fluorescence at equilibrium for either protein. ILBP was less stable than IFABP with a midpoint of 2. 9 M urea compared to 4.0 M urea for IFABP. Stopped-flow kinetic studies showed that both the folding and unfolding of these proteins were not monophasic, suggesting that either multiple paths or intermediate states were present during these processes. Proline isomerization is unlikely to be the cause of the multiphasic kinetics. ILBP had an intermediate state with molten globule-like spectral properties, whereas IFABP had an intermediate state with little if any secondary structure during folding and unfolding. Double-jump experiments showed that these intermediates appear to be on the folding path for each protein. The folding mechanisms of these proteins were markedly different, suggesting that the different sequences of these two proteins dictate different paths through the folding landscape to the same final structure.     

Differential stability of beta-sheets and alpha-helices in beta-lactamase: a high temperature molecular dynamics study of unfolding intermediates.

beta-Lactamase, which catalyzes beta-lactam antibiotics, is prototypical of large alpha/beta proteins with a scaffolding formed by strong noncovalent interactions. Experimentally, the enzyme is well characterized, and intermediates that are slightly less compact and having nearly the same content of secondary structure have been identified in the folding pathway. In the present study, high temperature molecular dynamics simulations have been carried out on the native enzyme in solution. Analysis of these results in terms of root mean square fluctuations in cartesian and [phi, psi] space, backbone dihedral angles and secondary structural hydrogen bonds forms the basis for an investigation of the topology of partially unfolded states of beta-lactamase. A differential stability has been observed for alpha-helices and beta-sheets upon thermal denaturation to putative unfolding intermediates. These observations contribute to an understanding of the folding/unfolding processes of beta-lactamases in particular, and other alpha/beta proteins in general.     

Topology, stability, sequence, and length: defining the determinants of two-state protein folding kinetics

The fastest simple, single domain proteins fold a million times more rapidly than the slowest. Ultimately this broad kinetic spectrum is determined by the amino acid sequences that define these proteins, suggesting that the mechanisms that underlie folding may be almost as complex as the sequences that encode them. Here, however, we summarize recent experimental results which suggest that (1) despite a vast diversity of structures and functions, there are fundamental similarities in the folding mechanisms of single domain proteins and (2) rather than being highly sensitive to the finest details of sequence, their folding kinetics are determined primarily by the large-scale, redundant features of sequence that determine a protein's gross structural properties. That folding kinetics can be predicted using simple, empirical, structure-based rules suggests that the fundamental physics underlying folding may be quite straightforward and that a general and quantitative theory of protein folding rates and mechanisms (as opposed to unfolding rates and thus protein stability) may be near on the horizon.     

Crystal structure of the disulfide bond-deficient azurin mutant C3A/C26A: how important is the S-S bond for folding and stability?

Azurin has a beta-barrel fold comprising eight beta-strands and one alpha helix. A disulfide bond between residues 3 and 26 connects the N-termini of beta strands beta1 and beta3. Three mutant proteins lacking the disulfide bond were constructed, C3A/C26A, C3A/C26I and a putative salt bridge (SB) in the C3A/S25R/C26A/K27R mutant. All three mutants exhibit spectroscopic properties similar to the wild-type protein. Furthermore, the crystal structure of the C3A/C26A mutant was determined at 2.0 A resolution and, in comparison to the wild-type protein, the only differences are found in the immediate proximity of the mutation. The mutants lose the 628 nm charge-transfer band at a temperature 10-22 degrees C lower than the wild-type protein. The folding of the zinc loaded C3A/C26A mutant was studied by guanidine hydrochloride (GdnHCl) induced denaturation monitored both by fluorescence and CD spectroscopy. The midpoint in the folding equilibrium, at 1.3 M GdnHCl, was observed using both CD and fluorescence spectroscopy. The free energy of folding determined from CD is -24.9 kJ.mol-1, a destabilization of approximately 20 kJ.mol-1 compared to the wild-type Zn2+-protein carrying an intact disulfide bond, indicating that the disulfide bond is important for giving azurin its stable structure. The C3A/C26I mutant is more stable and the SB mutant is less stable than C3A/C26A, both in terms of folding energy and thermal denaturation. The folding intermediate of the wild-type Zn2+-azurin is not observed for the disulfide-deficient C3A/C26A mutant. The rate of unfolding for the C3A/C26A mutant is similar to that of the wild-type protein, suggesting that the site of the mutation is not involved in an early unfolding reaction.     

The importance of surface loops for stabilizing an eightfold beta alpha barrel protein.

An important step in understanding how a protein folds is to determine those regions of the sequence that are critical to both its stability and its folding pathway. We chose phosphoribosyl anthranilate isomerase from Escherichia coli, which is a monomeric representative of the (beta alpha)8 barrel family of proteins, to construct a variant that carries an internal tandem duplication of the fifth beta alpha module. This (beta alpha)9 variant was enzymically active and therefore must have a wild-type (beta alpha)8 core. It had a choice a priori to fold to three different folding frames, which are distinguished by carrying the duplicated segment as an insert into one out of three different loops. Steady-state kinetic constants, the fluorescence properties of a crucial tryptophan residue, and limited proteolysis showed that the stable (beta alpha)9 variant carries the insertion between beta-strand 5 and alpha-helix 5. This preference can be explained by the important role of loops between alpha helices and beta strands in stabilizing the structure of the enzyme.     

The equilibrium unfolding pathway of a (beta/alpha)8 barrel

The (beta/alpha)(8) barrel is the most commonly occurring fold among enzymes. A key step towards rationally engineering (beta/alpha)(8) barrel proteins is to understand their underlying structural organization and folding energetics. Using misincorporation proton-alkyl exchange (MPAX), a new tool for solution structural studies of large proteins, we have performed a native-state exchange analysis of the prototypical (beta/alpha)(8) barrel triosephosphate isomerase. Three cooperatively unfolding subdomains within the structure are identified, as well as two partially unfolded forms of the protein. The C-terminal domain coincides with domains reported to exist in four other (beta/alpha)(8) barrels, but the two N-terminal domains have not been observed previously. These partially unfolded forms may represent sequential intermediates on the folding pathway of triosephosphate isomerase. The methods reported here should be applicable to a variety of other biological problems involving protein conformational changes.     

The tetrameric α-helical membrane protein GlpF unfolds via a dimeric folding intermediate

Many membrane proteins appear to be present and functional in higher-order oligomeric states. While few studies have analyzed the thermodynamic stability of α-helical transmembrane (TM) proteins under equilibrium conditions in the past, oligomerization of larger polytopic monomers has essentially not yet been studied. However, it is vital to study the folding of oligomeric membrane proteins to improve our understanding of the general mechanisms and pathways of TM protein folding. To investigate the folding and stability of the aquaglyceroporin GlpF from Escherichia coli, unfolding of the protein in mixed micelles was monitored by steady-state fluorescence and circular dichroism spectroscopy as well as by seminative sodium dodecyl sulfate-polyacrylamide gel electrophoresis analyses. On the basis of our results, it appears most likely that GlpF unfolds in a two-step process, involving the equilibrium of tetrameric, dimeric, and monomeric GlpF species. A kinetic analysis also indicates an intermediate along the kinetic GlpF unfolding pathway, and thus, two phases are involved in GlpF unfolding. While three-state unfolding pathways and a dimeric folding intermediate are not uncommon for water-soluble proteins, a stable (un)folding intermediate with a decreased oligomeric structure has not been detected or reported for any α-helical membrane protein.     

Redesigning the hydrophobic core of a four-helix-bundle protein.

Rationally redesigned variants of the 4-helix-bundle protein Rop are described. The novel proteins have simplified, repacked, hydrophobic cores and yet reproduce the structure and native-like physical properties of the wild-type protein. The repacked proteins have been characterized thermodynamically and their equilibrium and kinetic thermal and chemical unfolding properties are compared with those of wild-type Rop. The equilibrium stability of the repacked proteins to thermal denaturation is enhanced relative to that of the wild-type protein. The rate of chemically induced folding and unfolding of wild-type Rop is extremely slow when compared with other small proteins. Interestingly, although the repacked proteins are more thermally stable than the wild type, their rates of chemically induced folding and unfolding are greatly increased in comparison to wild type. Perhaps as a consequence of this, their equilibrium stabilities to chemical denaturants are slightly reduced in comparison to the wild type.     

Secondary structure-based profiles: use of structure-conserving scoring tables in searching protein sequence databases for structural similarities

The profile method, for detecting distantly related proteins by sequence comparison, has been extended to incorporate secondary structure information from known X-ray structures. The sequence of a known structure is aligned to sequences of other members of a given folding class. From the known structure, the secondary structure (alpha-helix, beta-strand or "other") is assigned to each position of the aligned sequences. As in the standard profile method, a position-dependent scoring table, termed a profile, is calculated from the aligned sequences. However, rather than using the standard Dayhoff mutation table in calculating the profile, we use distinct amino acid mutation tables for residues in alpha-helices, beta-strands or other secondary structures to calculate the profile. In addition, we also distinguish between internal and external residues. With this new secondary structure-based profile method, we created a profile for eight-stranded, antiparallel beta barrels of the insecticyanin folding class. It is based on the sequences of retinol-binding protein, insecticyanin and beta-lactoglobulin. Scanning the sequence database with this profile, it was possible to detect the sequence of avidin. The structure of streptavidin is known, and it appears to be distantly related to the antiparallel beta barrels. Also detected is the sequence of complement component C8, which we therefore predict to be a member of this folding class.     

Computer simulations of protein folding by targeted molecular dynamics

We have performed 128 folding and 45 unfolding molecular dynamics runs of chymotrypsin inhibitor 2 (CI2) with an implicit solvation model for a total simulation time of 0.4 microseconds. Folding requires that the three-dimensional structure of the native state is known. It was simulated at 300 K by supplementing the force field with a harmonic restraint which acts on the root-mean-square deviation and allows to decrease the distance to the target conformation. High temperature and/or the harmonic restraint were used to induce unfolding. Of the 62 folding simulations started from random conformations, 31 reached the native structure, while the success rate was 83% for the 66 trajectories which began from conformations unfolded by high-temperature dynamics. A funnel-like energy landscape is observed for unfolding at 475 K, while the unfolding runs at 300 K and 375 K as well as most of the folding trajectories have an almost flat energy landscape for conformations with less than about 50% of native contacts formed. The sequence of events, i.e., secondary and tertiary structure formation, is similar in all folding and unfolding simulations, despite the diversity of the pathways. Previous unfolding simulations of CI2 performed with different force fields showed a similar sequence of events. These results suggest that the topology of the native state plays an important role in the folding process.