Native structural propensity in cellular retinoic acid-binding protein I 64-88: the role of locally encoded structure in the folding of a beta-barrel protein

A central question in protein folding is the relative importance of locally encoded structure and cooperative interactions among residues distant in sequence. We have been exploring this question in a predominantly β-sheet protein, since β-structure formation clearly relies on both local and global sequence information. We present evidence that a 24-residue peptide corresponding to two linked hairpins of cellular retinoic acid-binding protein I (CRABP I) adopts significant native structure in aqueous solution. Prior work from our laboratory showed that the two turns contained in this fragment (turns III and IV) had the highest tendency of any of the eight turns in this anti-parallel β-barrel to fold into native turns. In addition, the primary sequence of these two turns is well conserved throughout the structural family to which CRABP I belongs, and residues in the turns and their associated hairpins participate in a network of conserved long-range interactions. We propose that the strong local-sequence biases within the chain segment comprising turns III and IV favor longer-range interactions that are crucial to the folding and native-state stability of CRABP I, and may play a similar role in related intracellular lipid-binding proteins (iLBPs).     

Nanosecond time scale folding dynamics of a pentapeptide in water

Reverse turns, four-residue sections of polypeptides where the chain changes direction by about 180 degrees, are thought to be important protein folding initiation structures. However, the time scale and mechanism for their formation have yet to be determined experimentally. To develop a microscopic picture of the formation of protein folding initiation structures, we have carried out a pair of 2.2-ns molecular dynamics simulations of Tyr-Pro-Gly-Asp-Val, a peptide which is known to form a high population of reverse turns in water. In the first simulation, which was started with the peptide in an ideal type II reverse turn involving the first four residues, the turn unfolded after about 1.4 ns. After about 0.6 ns in the second simulation, which was started with the peptide in a fully extended conformation, the peptide folded into a type II turn which had a transient existence before unfolding. The peptide remained unfolded for another 0.9 ns before folding into a type I turn involving the last four residues. The type I turn lasted for about 0.2 ns before unfolding. Thus, these simulations showed that protein folding initiation structures can form and dissolve on the nanosecond time scale. Furthermore, the atomic-level detail of the simulations allowed us to identify some of the interactions which can stabilize the folded structures. The type II turns were stabilized by either a salt bridge between the terminal groups or a backbone-C-terminal group hydrogen bond, and the type I turns were stabilized by a hydrophobic interaction between the proline and valine-side chains.     

前导肽介导的蛋白折叠机制及其对脂肪酶影响的研究进展

大多数蛋白质的形成过程主要由合成前体蛋白和合成功能蛋白两个步骤组成.在这个过程中,前导肽能够辅助蛋白质折叠或抑制它的活性.前导肽作为脂肪酶结构中重要的一段多肽链,通常作为分子内分子伴侣来辅助脂肪酶的折叠,同时该序列上包括糖基化位点在内的一些特殊位点,对酶的活性、极端环境稳定性、甲醇耐受性和底物特异性等性质具有重要影响....     

Reduced tendency to form a beta turn in peptides from the P22 tailspike protein correlates with a temperature-sensitive folding defect

A family of mutants of the P22 bacteriophage tailspike protein has been characterized as temperature sensitive for folding (tsf) by King and co-workers [King, J. (1986) Bio/Technology 4, 297-303]. There is substantial evidence that the tsf mutations alter the folding pathway but not the stability of the final folded protein. Several point mutations are known to cause the tsf phenotype; most of these occur in regions of the tailspike sequence likely to take up reverse turns. Hence, it has been hypothesized that the correct folding of the P22 tailspike protein requires formation of turns and that the mutations causing tsf phenotypes interfere at this critical stage. We have tested this hypothesis by study of isolated peptides corresponding to a region of the P22 tailspike harboring a tsf mutation. Comparison of the tendencies of wild-type and tsf sequences to adopt turn conformations was achieved by the synthesis of peptides with flanking cysteine residues and the use of a thiol-disulfide exchange assay. We find that the wild-type sequence, either as a decapeptide (Ac-CVKFPGIETC-CONH2) or as a dodecapeptide (Ac-CYVKFPGIETLC-CONH2), has a 3-5-fold greater tendency for its termini to approach closely enough to form the intramolecular disulfide than do the peptide sequences corresponding to the tsf mutant sequences, which have a Gly----Arg substitution (Ac-CVKFPRIETC-CONH2 or Ac-CYVKFPRIETLC-CONH2). A peptide with a D-Arg substituted for the Gly has a slightly higher turn propensity than does the wild type. Together with data from nuclear magnetic resonance analysis of the oxidized peptides, this suggests that a type II beta turn is favored by the wild-type sequence. Our results on isolated peptides from the P22 tailspike protein support the model for its folding that includes reverse turn formation as a critical step.     

Systematic circular permutation of an entire protein reveals essential folding elements

The importance of chain connectivity in determining protein function and stability can be examined by breaking the peptide backbone using a technique such as circular permutation. Cleavage at certain positions results in a complete loss of the ability of the protein to fold. When such cleavage sites occur sequentially in the primary structure, we call the region a 'folding element', a new concept that could assist in our understanding of the protein folding problem. The folding elements of dihydrofolate reductase have been assigned by conducting a systematic circular permutation analysis in which the peptide backbone was sequentially broken between every pair of residues in the protein. The positions of folding elements do not appear to correspond to secondary structure motifs, substrate or coenzyme binding sites, or accessible surface area. However, almost all of the amino acid residues known to be involved in early folding events are located within the folding elements.     

Turn scanning by site-directed mutagenesis: application to the protein folding problem using the intestinal fatty acid binding protein

We have systematically mutated residues located in turns between beta-strands of the intestinal fatty acid binding protein (IFABP), and a glycine in a half turn, to valine and have examined the stability, refolding rate constants and ligand dissociation constants for each mutant protein. IFABP is an almost all beta-sheet protein exhibiting a topology comprised of two five-stranded sheets surrounding a large cavity into which the fatty acid ligand binds. A glycine residue is located in seven of the eight turns between the antiparallel beta-strands and another in a half turn of a strand connecting the front and back sheets. Mutations in any of the three turns connecting the last four C-terminal strands slow the folding and decrease stability with the mutation between the last two strands slowing folding dramatically. These data suggest that interactions between the last four C-terminal strands are highly cooperative, perhaps triggered by an initial hydrophobic collapse. We suggest that this trigger is collapse of the highly hydrophobic cluster of amino acids in the D and E strands, a region previously shown to also affect the last stage of the folding process (Kim et al., 1997). Changing the glycine in the strand between the front and back sheets also results in a unstable, slow folding protein perhaps disrupting the D-E strand interactions. For most of the other turn mutations there was no apparent correlation between stability and refolding rate constants. In some turns, the interaction between strands, rather than the turn type, appears to be critical for folding while in others, turn formation itself appears to be a rate limiting step. Although there is no simple correlation between turn formation and folding kinetics, we propose that turn scanning by mutagenesis will be a useful tool for issues related to protein folding.     

A model study of protein nascent chain and cotranslational folding using hydrophobic-polar residues

To study protein nascent chain folding during biosynthesis, we investigate the folding behavior of models of hydrophobic and polar (HP) chains at growing length using both two-dimensional square lattice model and an optimized three-dimensional 4-state discrete off-lattice model. After enumerating all possible sequences and conformations of HP heteropolymers up to length N = 18 and N = 15 in two and three-dimensional space, respectively, we examine changes in adopted structure, stability, and tolerance to single point mutation as the nascent chain grows. In both models, we find that stable model proteins have fewer folded nascent chains during growth, and often will only fold after reaching full length. For the few occasions where partial chains of stable proteins fold, these partial conformations on average are very similar to the corresponding parts of the final conformations at full length. Conversely, we find that sequences with fewer stable nascent chains and sequences with native-like folded nascent chains are more stable. In addition, these stable sequences in general can have many more point mutations and still fold into the same conformation as the wild type sequence. Our results suggest that stable proteins are less likely to be trapped in metastable conformations during biosynthesis, and are more resistant to point-mutations. Our results also imply that less stable proteins will require the assistance of chaperone and other factors during nascent chain folding. Taken together with other reported studies, it seems that cotranslational folding may not be a general mechanism of in vivo protein folding for small proteins, and in vitro folding studies are still relevant for understanding how proteins fold biologically.     

Effects of alanine substitutions in alpha-helices of sperm whale myoglobin on protein stability.

The peptide backbones in folded native proteins contain distinctive secondary structures, alpha-helices, beta-sheets, and turns, with significant frequency. One question that arises in folding is how the stability of this secondary structure relates to that of the protein as a whole. To address this question, we substituted the alpha-helix-stabilizing alanine side chain at 16 selected sites in the sequence of sperm whale myoglobin, 12 at helical sites on the surface of the protein, and 4 at obviously internal sites. Substitution of alanine for bulky side chains at internal sites destabilizes the protein, as expected if packing interactions are disrupted. Alanine substitutions do not uniformly stabilize the protein, either in capping positions near the ends of helices or at mid-helical sites near the surface of myoglobin. When corrected for the extent of exposure of each side chain replaced by alanine at a mid-helix position, alanine replacement still has no clear effect in stabilizing the native structure. Thus linkage between the stabilization of secondary structure and tertiary structure in myoglobin cannot be demonstrated, probably because of the relatively small free energy differences between side chains in stabilizing isolated helix. By contrast, about 80% of the variance in free energy observed can be accounted for by the loss in buried surface area of the native residue substituted by alanine. The differential free energy of helix stabilization does not account for any additional variation.     

Role of local sequence in the folding of cellular retinoic abinding protein I: structural propensities of reverse turns

The experiments described here explore the role of local sequence in the folding of cellular retinoic acid binding protein I (CRABP I). This is a 136-residue, 10-stranded, antiparallel beta-barrel protein with seven beta-hairpins and is a member of the intracellular lipid binding protein (iLBP) family. The relative roles of local and global sequence information in governing the folding of this class of proteins are not well-understood. In question is whether the beta-turns are locally defined by short-range interactions within their sequences, and are thus able to play an active role in reducing the conformational space available to the folding chain, or whether the turns are passive, relying upon global forces to form. Short (six- and seven-residue) peptides corresponding to the seven CRABP I turns were analyzed by circular dichroism and NMR for their tendencies to take up the conformations they adopt in the context of the native protein. The results indicate that two of the peptides, encompassing turns III and IV in CRABP I, have a strong intrinsic bias to form native turns. Intriguingly, these turns are on linked hairpins in CRABP I and represent the best-conserved turns in the iLBP family. These results suggest that local sequence may play an important role in narrowing the conformational ensemble of CRABP I during folding.     

Polymer models of protein stability, folding, and interactions

The unfolded state and flexible linkers in the folded structure play essential roles in protein stability and folding and protein-protein interactions. Intrinsic to these roles is the fact that unfolded proteins and flexible linkers sample many different conformations. Polymer models may capture this and complement experiments in elucidating the contributions of the unfolded state and flexible linkers. Here I review what can be predicted from these models and how well these predictions match experiments. For example, Gaussian chain models give quantitatively reasonable predictions of the effects of residual charge-charge interactions in the unfolded state and qualitatively reasonable results for the effects of spatial confinement and macromolecular crowding on protein stability. A wormlike chain model has met with success in quantifying the effects of flexible linkers in binding affinity enhancement and in regulatory switches. In future developments, more realistic models may emerge from molecular dynamics simulations, and these models will guide experiments to advance our understanding of the unfolded state and flexible linkers.     

Transient structure formation in unfolded acyl-coenzyme A-binding protein observed by site-directed spin labelling

Paramagnetic relaxation has been used to monitor the formation of structure in the folding peptide chain of guanidinium chloride-denatured acyl-coenzyme A-binding protein. The spin label (1-oxyl-2,2,5,5-tetramethyl-3-pyrroline-3-methyl)methanesulfonate (MTSL) was covalently bound to a single cysteine residue introduced into five different positions in the amino acid sequence. It was shown that the formation of structure in the folding peptide chain at conditions where 95% of the sample is unfolded brings the relaxation probe close to a wide range of residues in the peptide chain, which are not affected in the native folded structure. It is suggested that the experiment is recording the formation of many discrete and transient structures in the polypeptide chain in the preface of protein folding. Analysis of secondary chemical shifts shows a high propensity for alpha-helix formation in the C-terminal part of the polypeptide chain, which forms an alpha-helix in the native structure and a high propensity for turn formation in two regions of the polypeptide that form turns in the native structure. The results contribute to the idea that native-like structural elements form transiently in the unfolded state, and that these may be of importance to the initiation of protein folding.     

Accurate computer-based design of a new backbone conformation in the second turn of protein L.

The rational design of loops and turns is a key step towards creating proteins with new functions. We used a computational design procedure to create new backbone conformations in the second turn of protein L. The Protein Data Bank was searched for alternative turn conformations, and sequences optimal for these turns in the context of protein L were identified using a Monte Carlo search procedure and an energy function that favors close packing. Two variants containing 12 and 14 mutations were found to be as stable as wild-type protein L. The crystal structure of one of the variants has been solved at a resolution of 1.9 A, and the backbone conformation in the second turn is remarkably close to that of the in silico model (1.1 A RMSD) while it differs significantly from that of wild-type protein L (the turn residues are displaced by an average of 7.2 A). The folding rates of the redesigned proteins are greater than that of the wild-type protein and in contrast to wild-type protein L the second beta-turn appears to be formed at the rate limiting step in folding.     

Sequence swapping does not result in conformation swapping for the beta4/beta5 and beta8/beta9 beta-hairpin turns in human acidic fibroblast growth factor

The beta-turn is the most common type of nonrepetitive structure in globular proteins, comprising ~25% of all residues; however, a detailed understanding of effects of specific residues upon beta-turn stability and conformation is lacking. Human acidic fibroblast growth factor (FGF-1) is a member of the beta-trefoil superfold and contains a total of five beta-hairpin structures (antiparallel beta-sheets connected by a reverse turn). beta-Turns related by the characteristic threefold structural symmetry of this superfold exhibit different primary structures, and in some cases, different secondary structures. As such, they represent a useful system with which to study the role that turn sequences play in determining structure, stability, and folding of the protein. Two turns related by the threefold structural symmetry, the beta4/beta5 and beta8/beta9 turns, were subjected to both sequence-swapping and poly-glycine substitution mutations, and the effects upon stability, folding, and structure were investigated. In the wild-type protein these turns are of identical length, but exhibit different conformations. These conformations were observed to be retained during sequence-swapping and glycine substitution mutagenesis. The results indicate that the beta-turn structure at these positions is not determined by the turn sequence. Structural analysis suggests that residues flanking the turn are a primary structural determinant of the conformation within the turn.     

Structural and sequence features of two residue turns in beta‐hairpins

Beta‐turns in beta‐hairpins have been implicated as important sites in protein folding. In particular, two residue β‐turns, the most abundant connecting elements in beta‐hairpins, have been a major target for engineering protein stability and folding. In this study, we attempted to investigate and update the structural and sequence properties of two residue turns in beta‐hairpins with a large data set. For this, 3977 beta‐turns were extracted from 2394 nonhomologous protein chains and analyzed. First, the distribution, dihedral angles and twists of two residue turn types were determined, and compared with previous data. The trend of turn type occurrence and most structural features of the turn types were similar to previous results, but for the first time Type II turns in beta‐hairpins were identified. Second, sequence motifs for the turn types were devised based on amino acid positional potentials of two‐residue turns, and their distributions were examined. From this study, we could identify code‐like sequence motifs for the two residue beta‐turn types. Finally, structural and sequence properties of beta‐strands in the beta‐hairpins were analyzed, which revealed that the beta‐strands showed no specific sequence and structural patterns for turn types. The analytical results in this study are expected to be a reference in the engineering or design of beta‐hairpin turn structures and sequences. Proteins 2014; 82:1721–1733. © 2014 Wiley Periodicals, Inc.     

Surface amino acids as sites of temperature-sensitive folding mutations in the P22 tailspike protein

Temperature-sensitive folding (tsf) mutations in the gene for the thermostable P22 tailspike interfere with the polypeptide chain folding and association pathway at restrictive temperature without altering the thermostability of the protein once correctly folded and assembled at permissive temperature. Though the native proteins matured at permissive temperature are biologically active, many of them display alterations in electrophoretic mobility. The native forms of 15 of these tsf mutant proteins have been purified and characterized. The purified proteins differed in electrophoretic mobility and isoelectric point from wild type but did not show evidence of major conformational alterations. The results suggest that the electrophoretic variations conferred by the 15 tsf amino acid substitutions are due to changes in the net charge at solvent-accessible sites in the native form of the mutant protein. During the maturation of the chains at restrictive temperature, these sites influence the conformation of intermediates in chain folding and association. The amino acid sequences at these sites resemble those found at turns in polypeptide chains. The isolation of tsf mutations requires that the mature structure of the tailspike accommodates the mutant amino acid substitution without loss of function. The solvent-accessible sites are probably at the surface of this structural protein. This would explain how bulky mutant substitutions, such as arginines for glycines, are accommodated in the native tailspike structure. Such sites, stabilizing intermediates in the folding pathway and located on the surface of the mature protein, probably represent a general class of conformational substrates for tsf mutations.     

Local sequence information in cellular retinoic acid-binding protein I: specific residue roles in beta-turns

We have recently shown that two of the beta-turns (III and IV) in the ten-stranded, beta-clam protein, cellular retinoic acid-binding protein I (CRABP I), are favored in short peptide fragments, arguing that they are encoded by local interactions (K. S. Rotondi and L. M. Gierasch, Biochemistry, 2003, Vol. 42, pp. 7976-7985). In this paper we examine these turns in greater detail to dissect the specific local interactions responsible for their observed native conformational biases. Conformations of peptides corresponding to the turn III and IV fragments were examined under conditions designed to selectively disrupt stabilizing interactions, using pH variation, chaotrope addition, or mutagenesis to probe specific side-chain influences. We find that steric constraints imposed by excluded volume effects between near neighbor residues (i,i+2), favorable polar (i,i+2) interactions, and steric permissiveness of glycines are the principal factors accounting for the observed native bias in these turns. Longer-range stabilizing interactions across the beta-turns do not appear to play a significant role in turn stability in these short peptides, in contrast to their importance in hairpins. Additionally, our data add to a growing number of examples of the 3:5 type I turn with a beta-bulge as a class of turns with high propensity to form locally defined structure. Current work is directed at the interplay between the local sequence information in the turns and more long-range influences in the mechanism of folding of this predominantly beta-sheet protein.     

Importance of context in protein folding: secondary structural propensities versus tertiary contact-assisted secondary structure formation

Molecular dynamics simulations can be used to reveal the detailed conformational behaviors of peptides and proteins. By comparing fragment and full-length protein simulations, we can investigate the role of each peptide segment in the folding process. Here, we take advantage of information regarding the helix formation process from our previous simulations of barnase and protein A as well as new simulations of four helical fragments from these proteins at three different temperatures, starting with both helical and extended structures. Segments with high helical propensity began the folding process by tethering the chain through side chain interactions involving either polar interactions, such as salt bridges, or hydrophobic staples. These tethers were frequently nonnative (i.e., not i --> i + 4 spacing) and provided a scaffold for other residues, thereby limiting the conformational search. The helical structure then propagated on both sides of the tether. Segments with low stability and propensity formed later in the folding process and utilized contacts with other portions of the protein when folding. These helices formed via a tertiary contact-assisted mechanism, primarily via hydrophobic contacts between residues distant in sequence. Thus, segments with different helical propensities appear to play different roles during protein folding. Furthermore, the active role of nonlocal side chains in helix formation highlights why we must move beyond simple hierarchical models of protein folding.     

A framework for describing topological frustration in models of protein folding

In a natively folded protein of moderate or larger size, the protein backbone may weave through itself in complex ways, raising questions about what sequence of events might have to occur in order for the protein to reach its native configuration from the unfolded state. A mathematical framework is presented here for describing the notion of a topological folding barrier, which occurs when a protein chain must pass through a hole or opening, formed by other regions of the protein structure. Different folding pathways encounter different numbers of such barriers and therefore different degrees of frustration. A dynamic programming algorithm finds the optimal theoretical folding path and minimal degree of frustration for a protein based on its natively folded configuration. Calculations over a database of protein structures provide insights into questions such as whether the path of minimal frustration might tend to favor folding from one or from many sites of folding nucleation, or whether proteins favor folding around the N terminus, thereby providing support for the hypothesis that proteins fold co-translationally. The computational methods are applied to a multi-disulfide bonded protein, with computational findings that are consistent with the experimentally observed folding pathway. Attention is drawn to certain complex protein folds for which the computational method suggests there may be a preferred site of nucleation or where folding is likely to proceed through a relatively well-defined pathway or intermediate. The computational analyses lead to testable models for protein folding.     

Principles of protein folding--a perspective from simple exact models.

General principles of protein structure, stability, and folding kinetics have recently been explored in computer simulations of simple exact lattice models. These models represent protein chains at a rudimentary level, but they involve few parameters, approximations, or implicit biases, and they allow complete explorations of conformational and sequence spaces. Such simulations have resulted in testable predictions that are sometimes unanticipated: The folding code is mainly binary and delocalized throughout the amino acid sequence. The secondary and tertiary structures of a protein are specified mainly by the sequence of polar and nonpolar monomers. More specific interactions may refine the structure, rather than dominate the folding code. Simple exact models can account for the properties that characterize protein folding: two-state cooperativity, secondary and tertiary structures, and multistage folding kinetics--fast hydrophobic collapse followed by slower annealing. These studies suggest the possibility of creating "foldable" chain molecules other than proteins. The encoding of a unique compact chain conformation may not require amino acids; it may require only the ability to synthesize specific monomer sequences in which at least one monomer type is solvent-averse.     

Structure and dynamics of an acid-denatured protein G mutant

NMR studies of protein denatured states provide insights into potential initiation sites for folding that may be too transient to be observed kinetically. We have characterized the structure and dynamics of the acid-denatured state of protein G by using a F30H mutant of G(B1) which is on the margin of stability. At 5 degrees C, F30H-G(B1) is greater than 95% folded at pH 7.0 and is greater than 95% unfolded at pH 4.0. This range of stability is useful because the denatured state can be examined under relatively mild conditions which are optimal for folding G(B1). We have assigned almost all backbone (15)N, H(N), and H(alpha) resonances in the acid-denatured state. Chemical shift, coupling constant, and NOE data indicate that the denatured state has considerably more residual structure when studied under these mild conditions than in the presence of chemical denaturants. The acid-denatured state populates nativelike conformations with both alpha-helical and beta-hairpin characteristics. To our knowledge, this is the first example of a denatured state with NOE and coupling constant evidence for beta-hairpin character. A number of non-native turn structures are also detected, particularly in the region corresponding to the beta1-beta2 hairpin of the folded state. Steady-state ?(1)H-(15)N? NOE results demonstrate restricted backbone flexibility in more structured regions of the denatured protein. Overall, our studies suggest that regions of the helix, the beta3-beta4 hairpin, and the beta1-beta2 turn may serve as potential initiation sites for folding of G(B). Furthermore, residual structure in acid-denatured F30H-G(B1) is more extensive than in peptide fragments corresponding to the beta1-beta2, alpha-helix, and beta3-beta4 regions, suggesting additional medium-to-long-range interactions in the full-length polypeptide chain.