Analysis of the factors that stabilize a designed two-stranded antiparallel beta-sheet

Autonomously folding beta-hairpins (two-strand antiparallel beta-sheets) have become increasingly valuable tools for probing the forces that control peptide and protein conformational preferences. We examine the effects of variations in sequence and solvent on the stability of a previously designed 12-residue peptide (1). This peptide adopts a beta-hairpin conformation containing a two-residue loop (D-Pro-Gly) and a four-residue interstrand sidechain cluster that is observed in the natural protein GB1. We show that the conformational propensity of the loop segment plays an important role in beta-hairpin stability by comparing 1 with (D)P--> N mutant 2. In addition, we show that the sidechain cluster contributes both to conformational stability and to folding cooperativity by comparing 1 with mutant 3, in which two of the four cluster residues have been changed to serine. Thermodynamic analysis suggests that the high loop-forming propensity of the (D)PG segment decreases the entropic cost of beta-hairpin formation relative to the more flexible NG segment, but that the conformational rigidity of (D)PG may prevent optimal contacts between the sidechains of the GB1-derived cluster. The enthalpic favorability of folding in these designed beta-hairpins suggests that they are excellent scaffolds for studying the fundamental mechanisms by which amino acid sidechains interact with one another in folded proteins.     

Buried hydrophobic side-chains essential for the folding of the parallel beta-helix domains of the P22 tailspike

The processive beta-strands and turns of a polypeptide parallel beta-helix represent one of the topologically simplest beta-sheet folds. The three subunits of the tailspike adhesin of phage P22 each contain 13 rungs of a parallel beta-helix followed by an interdigitated section of triple-stranded beta-helix. Long stacks of hydrophobic residues dominate the elongated buried core of these two beta-helix domains and extend into the core of the contiguous triple beta-prism domain. To test whether these side-chain stacks represent essential residues for driving the chain into the correct fold, each of three stacked phenylalanine residues within the buried core were substituted with less bulky amino acids. The mutant chains with alanine in place of phenylalanine were defective in intracellular folding. The chains accumulated exclusively in the aggregated inclusion body state regardless of temperature of folding. These severe folding defects indicate that the stacked phenylalanine residues are essential for correct parallel beta-helix folding. Replacement of the same phenylalanine residues with valine or leucine also impaired folding in vivo, but with less severity. Mutants were also constructed in a second buried stack that extends into the intertwined triple-stranded beta-helix and contiguous beta-prism regions of the protein. These mutants exhibited severe defects in later stages of chain folding or assembly, accumulating as misfolded but soluble multimeric species. The results indicate that the formation of the buried hydrophobic stacks is critical for the correct folding of the parallel beta-helix, triple-stranded beta-helix, and beta-prism domains in the tailspike protein.     

Factors involved in the stability of isolated beta-sheets: Turn sequence, beta-sheet twisting, and hydrophobic surface burial

We have recently reported on the design of a 20-residue peptide able to form a significant population of a three-stranded up-and-down antiparallel beta-sheet in aqueous solution. To improve our beta-sheet model in terms of the folded population, we have modified the sequences of the two 2-residue turns by introducing the segment DPro-Gly, a sequence shown to lead to more rigid type II' beta-turns. The analysis of several NMR parameters, NOE data, as well as Deltadelta(CalphaH), DeltadeltaC(beta), and Deltadelta(Cbeta) values, demonstrates that the new peptide forms a beta-sheet structure in aqueous solution more stable than the original one, whereas the substitution of the DPro residues by LPro leads to a random coil peptide. This agrees with previous results on beta-hairpin-forming peptides showing the essential role of the turn sequence for beta-hairpin folding. The well-defined beta-sheet motif calculated for the new designed peptide (pair-wise RMSD for backbone atoms is 0.5 +/- 0.1 A) displays a high degree of twist. This twist likely contributes to stability, as a more hydrophobic surface is buried in the twisted beta-sheet than in a flatter one. The twist observed in the up-and-down antiparallel beta-sheet motifs of most proteins is less pronounced than in our designed peptide, except for the WW domains. The additional hydrophobic surface burial provided by beta-sheet twisting relative to a "flat" beta-sheet is probably more important for structure stability in peptides and small proteins like the WW domains than in larger proteins for which there exists a significant contribution to stability arising from their extensive hydrophobic cores.     

Force field influences in beta-hairpin folding simulations

All-atom force fields are now routinely used for more detailed understanding of protein folding mechanisms. However, it has been pointed out that use of all-atom force fields does not guarantee more accurate representations of proteins; in fact, sometimes it even leads to biased structural distributions. Indeed, several issues remain to be solved in force field developments, such as accurate treatment of implicit solvation for efficient conformational sampling and proper treatment of backbone interactions for secondary structure propensities. In this study, we first investigate the quality of several recently improved backbone interaction schemes in AMBER for folding simulations of a beta-hairpin peptide, and further study their influences on the peptide's folding mechanism. Due to the significant number of simulations needed for a thorough analysis of tested force fields, the implicit Poisson-Boltzmann solvent was used in all simulations. The chosen implicit solvent was found to be reasonable for studies of secondary structures based on a set of simulations of both alpha-helical and beta-hairpin peptides with the TIP3P explicit solvent as benchmark. Replica exchange molecular dynamics was also utilized for further efficient conformational sampling. Among the tested AMBER force fields, ff03 and a revised ff99 force field were found to produce structural and thermodynamic data in comparably good agreement with the experiment. However, detailed folding pathways, such as the order of backbone hydrogen bond zipping and the existence of intermediate states, are different between the two force fields, leading to force field-dependent folding mechanisms.     

Effects of turn residues in directing the formation of the beta-sheet and in the stability of the beta-sheet

The designed peptide (denoted 20-mer, sequence VFITS(D)PGKTYTEV(D)PGOKILQ) has been shown to form a three-strand antiparallel beta-sheet. It is generally believed that the (D)Pro-Gly segment has the propensity to adopt a type II' beta-turn, thereby promoting the formation of this beta-sheet. Here, we replaced (D)Pro-Gly with Asp-Gly, which should favor a type I' turn, to examine the influence of different type of turns on the stability of the beta-sheet. Contrary to our expectation, the mutant peptide, denoted P6D, forms a five-residue type I turn plus a beta-bulge between the first two strands due to a one amino-acid frameshift in the hydrogen bonding network and side-chain inversion of the first beta-strand. In contrast, the same kind of substitution at (D)Pro-14 in the double mutant, denoted P6DP14D, does not yield the same effect. These observations suggest that the SDGK sequence disfavors the type I' conformation while the VDGO sequence favors a type I' turn, and that the frameshift in the first strand provides a way for the peptide to accommodate a disfavored turn sequence by protruding a bulge in the formation of the beta-hairpin. Thus, different types of turns can affect the stability of a beta-structure.     

A partially folded intermediate species of the beta-sheet protein apo-pseudoazurin is trapped during proline-limited folding

The folding of apo-pseudoazurin, a 123-residue, predominantly beta-sheet protein with a complex Greek key topology, has been investigated using several biophysical techniques. Kinetic analysis of refolding using far- and near-ultraviolet circular dichroism (UV CD) shows that the protein folds slowly to the native state with rate constants of 0.04 and 0.03 min(-1), respectively, at pH 7.0 and at 15 degrees C. This process has an activation enthalpy of approximately 90 kJ/mole and is catalyzed by cyclophilin A, indicating that folding is limited by trans-cis proline isomerization, presumably around the Xaa-Pro 20 bond that is in the cis isomer in the native state. Before proline isomerization, an intermediate accumulates during folding. This species has a substantial signal in the far-UV CD, a nonnative signal in the near-UV CD, exposed hydrophobic surfaces (judged by 1-anilino naphthalenesulphonate binding), a noncooperative denaturation transition, and a dynamic structure (revealed by line broadening on the nuclear magnetic resonance time scale). We compare the properties of this intermediate with partially folded states of other proteins and discuss its role in folding of this complex Greek key protein.     

Autonomous folding of a peptide corresponding to the N-terminal beta-hairpin from ubiquitin.

The N-terminal 17 residues of ubiquitin have been shown by 1H NMR to fold autonomously into a beta-hairpin structure in aqueous solution. This structure has a specific, native-like register, though side-chain contacts differ in detail from those observed in the intact protein. An autonomously folding hairpin has previously been identified in the case of streptococcal protein G, which is structurally homologous with ubiquitin, but remarkably, the two are not in topologically equivalent positions in the fold. This suggests that the organization of folding may be quite different for proteins sharing similar tertiary structures. Two smaller peptides have also been studied, corresponding to the isolated arms of the N-terminal hairpin of ubiquitin, and significant differences from simple random coil predictions observed in the spectra of these subfragments, suggestive of significant limitation of the backbone conformational space sampled, presumably as a consequence of the strongly beta-structure favoring composition of the sequences. This illustrates the ability of local sequence elements to express a propensity for beta-structure even in the absence of actual sheet formation. Attempts were made to estimate the population of the folded state of the hairpin, in terms of a simple two-state folding model. Using published "random coil" values to model the unfolded state, and values derived from native ubiquitin for the putative unique, folded state, it was found that the apparent population varied widely for different residues and with different NMR parameters. Use of the spectra of the subfragment peptides to provide a more realistic model of the unfolded state led to better agreement in the estimates that could be obtained from chemical shift and coupling constant measurements, while making it clear that some other approaches to population estimation could not give meaningful results, because of the tendency to populate the beta-region of conformational space even in the absence of the hairpin structure.     

The roles of turn formation and cross-strand interactions in fibrillization of peptides derived from the OspA single-layer beta-sheet

We previously demonstrated that a beta-hairpin peptide, termed BH(9-10), derived from a single-layer beta-sheet of Borrelia OspA protein, formed a native-like beta-turn in trifluoroethanol (TFE) solution, and it assembled into amyloid-like fibrils at higher TFE concentrations. This peptide is highly charged, and fibrillization of such a hydrophilic peptide is quite unusual. In this study, we designed a circularly permutated peptide of BH(9-10), termed BH(10-9). When folded into their respective beta-hairpin structures found in OspA, these peptides would have identical cross-strand interactions but different turns connecting the strands. NMR study revealed that BH(10-9) had little propensity to form a turn structure both in aqueous and TFE solutions. At higher TFE concentration, BH(10-9) precipitated with a concomitant alpha-to-beta conformational conversion, in a similar manner to the BH(9-10) fibrillization. However, the BH(10-9) precipitates were nonfibrillar aggregation. The precipitation kinetics of BH(10-9) was exponential, consistent with a first-order molecular assembly reaction, while the fibrillization of BH(9-10) showed sigmoidal kinetics, indicative of a two-step reaction consisting of nucleation and molecular assembly. The correlation between native-like turn formation and fibrillization of our peptide system strongly suggests that BH(9-10) adopts a native-like beta-hairpin conformation in the fibrils. Remarkably, seeding with the preformed BH(10-9) precipitates changed the two-step BH(9-10) fibrillization to a one-step molecular assembly reaction, and disrupted the BH(9-10) fibril structure, indicating interactions between the BH(10-9) aggregates and the BH(9-10) peptide. Our results suggest that, in these peptides, cross-strand interactions are the driving force for molecular assembly, and turn formation limits modes of peptide assembly.     

Energy landscape and dynamics of the beta-hairpin G peptide and its isomers: Topology and sequences

We have investigated free energy landscape [MM/PBSA + normal modes entropy] of permutations in the G peptide (41-56) from the protein G B1 domain by studying six isomers corresponding to moving the hydrophobic cluster along the beta-strands (toward the turn: T1, AGEWTYDDKTFTVTET; T2, GEDTWDYATFTVTKTE; T3, GEDDWTYATFTVTKTE; toward the end: E1, WTYDDAGETKTFTVT; E2, WEYTGDDATKTETFTV; E3, WTYEGDDATKTETFTV). The free energy terms include molecular mechanics energy, Poisson-Boltzmann electrostatic solvation energy, surface area solvation energy, and conformational entropy estimated by using normal mode analysis. From the wild type to T1, then T3, and finally T2, we see a progressively changing energy landscape, toward a less stable beta-hairpin structure. Moving the hydrophobic cluster outside toward the end region causes a greater change in the energy landscape. alpha-Helical instead of a beta-hairpin structure was the most stable form for the E2 isomer. However, no matter how much the sequence changes, for all variants studied, ideal "native" beta-hairpin topologies remain as minima (regardless of whether global or local) in the energy landscape. In general, we find that the energy landscape is dependent on the hydrophobic cluster topology and on the sequence. Our present study indicates that the key is the relative conformational energies of the different conformations. Changes in the sequence strongly modulate the relative stabilities of topologically similar regions in the energy landscape, rather than redefine the topology space. This finding is consistent with a population redistribution in the process of protein folding. The limited variation of topological space, compared with the number of possible sequence changes, may relate to the observation that the number of known protein folds are far less than the sequential allowance.     

Importance of secondary structural specificity determinants in protein folding: insertion of a native beta-sheet sequence into an alpha-helical coiled-coil

To examine how a short secondary structural element derived from a native protein folds when in a different protein environment, we inserted an 11-residue beta-sheet segment (cassette) from human immunoglobulin fold, Fab new, into an alpha-helical coiled-coil host protein (cassette holder). This de novo design protein model, the structural cassette mutagenesis (SCM) model, allows us to study protein folding principles involving both short- and long-range interactions that affect secondary structure stability and conformation. In this study, we address whether the insertion of this beta-sheet cassette into the alpha-helical coiled-coil protein would result in conformational change nucleated by the long-range tertiary stabilization of the coiled-coil, therefore overriding the local propensity of the cassette to form beta-sheet, observed in its native immunoglobulin fold. The results showed that not only did the nucleating helices of the coiled-coil on either end of the cassette fail to nucleate the beta-sheet cassette to fold with an alpha-helical conformation, but also the entire chimeric protein became a random coil. We identified two determinants in this cassette that prevented coiled-coil formation: (1) a tandem dipeptide NN motif at the N-terminal of the beta-sheet cassette, and (2) the hydrophilic Ser residue, which would be buried in the hydrophobic core if the coiled-coil structure were to fold. By amino acid substitution of these helix disruptive residues, that is, either the replacement of the NN motif with high helical propensity Ala residues or the substitution of Ser with Leu to enhance hydrophobicity, we were able to convert the random coil chimeric protein into a fully folded alpha-helical coiled-coil. We hypothesized that this NN motif is a "secondary structural specificity determinant" which is very selective for one type of secondary structure and may prevent neighboring residues from adopting an alternate protein fold. These sequences with secondary structural specificity determinants have very strong local propensity to fold into a specific secondary structure and may affect overall protein folding by acting as a folding initiation site.     

Characterization of two potentially universal turn motifs that shape the repeated five-residues fold--crystal structure of a lumenal pentapeptide repeat protein from Cyanothece 51142

The genome of the diurnal cyanobacterium Cyanothece sp. PCC 51142 has recently been sequenced and observed to contain 35 pentapeptide repeat proteins (PRPs). These proteins, while present throughout the prokaryotic and eukaryotic kingdoms, are most abundant in cyanobacteria. The sheer number of PRPs in cyanobacteria coupled with their predicted location in every cellular compartment argues for important, yet unknown, physiological and biochemical functions. To gain biochemical insights, the crystal structure for Rfr32, a 167-residue PRP with an N-terminal 29-residue signal peptide, was determined at 2.1 A resolution. The structure is dominated by 21 tandem pentapeptide repeats that fold into a right-handed quadrilateral beta-helix, or Rfr-fold, as observed for the tandem pentapeptide repeats in the only other PRP structure, the mycobacterial fluoroquinoline resistance protein MfpA from Mycobacterium tuberculosis. Sitting on top of the Rfr-fold are two short, antiparallel alpha-helices, bridged with a disulfide bond, that perhaps prevent edge-to-edge aggregation at the C terminus. Analysis of the main-chain (Phi,Psi) dihedral orientations for the pentapeptide repeats in Rfr32 and MfpA makes it possible to recognize the structural details for the two distinct types of four-residue turns adopted by the pentapeptide repeats in the Rfr-fold. These turns, labeled type II and type IV beta-turns, may be universal motifs that shape the Rfr-fold in all PRPs.     

Crystal structures and increased stabilization of the protein G variants with switched folding pathways NuG1 and NuG2

We recently described two protein G variants (NuG1 and NuG2) with redesigned first hairpins that were almost twice as stable, folded 100-fold faster, and had a switched folding mechanism relative to the wild-type protein. To test the structural accuracy of our design algorithm and to provide insights to the dramatic changes in the kinetics and thermodynamics of folding, we have now determined the crystal structures of NuG1 and NuG2 to 1.8 A and 1.85 A, respectively. We find that they adopt hairpin structures that are closer to the computational models than to wild-type protein G; the RMSD of the NuG1 hairpin to the design model and the wild-type structure are 1.7 A and 5.1 A, respectively. The crystallographic B factor in the redesigned first hairpin of NuG1 is systematically higher than the second hairpin, suggesting that the redesigned region is somewhat less rigid. A second round of structure-based design yielded new variants of NuG1 and NuG2, which are further stabilized by 0.5 kcal/mole and 0.9 kcal/mole.     

Fast protein folding on downhill energy landscape

Proteins fold in a time range of microseconds to minutes despite the large amount of possible conformers. Molecular dynamics simulations of a three-stranded antiparallel beta-sheet peptide (for a total of 12.6 microsec and 72 folding events) show that at the melting temperature the unfolded state ensemble contains many more conformers than those sampled during a folding event.     

The role of strand 1 of the C β-sheet in the structure and function of α1-antitrypsin

Serpins inhibit cognate serine proteases involved in a number of important processes including blood coagulation and inflammation. Consequently, loss of serpin function or stability results in a number of disease states. Many of the naturally occurring mutations leading to disease are located within strand 1 of the C beta-sheet of the serpin. To ascertain the structural and functional importance of each residue in this strand, which constitutes the so-called distal hinge of the reactive center loop of the serpin, an alanine scanning study was carried out on recombinant alpha(1)-antitrypsin Pittsburgh mutant (P1 = Arg). Mutation of the P10' position had no effect on its inhibitory properties towards thrombin. Mutations to residues P7' and P9' caused these serpins to have an increased tendency to act as substrates rather than inhibitors, while mutations at P6' and P8' positions caused the serpin to behave almost entirely as a substrate. Mutations at the P6' and P8' residues of the C beta-sheet, which are buried in the hydrophobic core in the native structure, caused the serpin to become highly unstable and polymerize much more readily. Thus, P6' and P8' mutants of alpha(1)-antitrypsin had melting temperatures 14 degrees lower than wild-type alpha(1)-antitrypsin. These results indicate the importance of maintaining the anchoring of the distal hinge to both the inhibitory mechanism and stability of serpins, the inhibitory mechanism being particularly sensitive to any perturbations in this region. The results of this study allow more informed analysis of the effects of mutations found at these positions in disease-associated serpin variants.     

Elongation of the BH8 β-hairpin peptide: Electrostatic interactions in β-hairpin formation and stability

An elongated version of the de novo designed beta-hairpin peptide, BH8, has allowed us to gain insight into the role of electrostatic interactions in beta-hairpin stability. A Lys-Glu electrostatic pair has been introduced by adding a residue at the beginning and at the end of the N-terminal and C-terminal strands, respectively, of the beta-hairpin structure, in both orientations. The two resulting peptides and controls having Ala residues at these positions and different combinations of Ala with Lys, or Glu residues, have been analyzed by nuclear magnetic resonance (NMR), under different pH and ionic strength conditions. All of the NMR parameters, in particular the conformational shift analysis of Calpha protons and the coupling constants, (3)J(HNalpha), correlate well and the population estimates are in reasonable agreement among the different methods used. In the most structured peptides, we find an extension of the beta-hairpin structure comprising the two extra residues. Analysis of the pH and salt dependence shows that ionic pairs contribute to beta-hairpin stability. The interaction is electrostatic in nature and can be screened by salt. There is also an important salt-independent contribution of negatively charged groups to the stability of this family of beta-hairpin peptides.     

Influence of hPin1 WW N-terminal domain boundaries on function, protein stability, and folding

An N-terminally truncated and cooperatively folded version (residues 6-39) of the human Pin1 WW domain (hPin1 WW hereafter) has served as an excellent model system for understanding triple-stranded beta-sheet folding energetics. Here we report that the negatively charged N-terminal sequence (Met1-Ala-Asp-Glu-Glu5) previously deleted, and which is not conserved in highly homologous WW domain family members from yeast or certain fungi, significantly increases the stability of hPin1 WW (approximately 4 kJ mol(-1) at 65 degrees C), in the context of the 1-39 sequence based on equilibrium measurements. N-terminal truncations and mutations in conjunction with a double mutant cycle analysis and a recently published high-resolution X-ray structure of the hPin1 cis/trans-isomerase suggest that the increase in stability is due to an energetically favorable ionic interaction between the negatively charged side chains in the N terminus of full-length hPin1 WW and the positively charged epsilon-ammonium group of residue Lys13 in beta-strand 1. Our data therefore suggest that the ionic interaction between Lys13 and the charged N terminus is the optimal solution for enhanced stability without compromising function, as ascertained by ligand binding studies. Kinetic laser temperature-jump relaxation studies reveal that this stabilizing interaction has not formed to a significant extent in the folding transition state at near physiological temperature, suggesting a differential contribution of the negatively charged N-terminal sequence to protein stability and folding rate. As neither the N-terminal sequence nor Lys13 are highly conserved among WW domains, our data further suggest that caution must be exercised when selecting domain boundaries for WW domains for structural, functional, or thermodynamic studies.     

Structural and dynamic characterization of an unfolded state of poplar apo-plastocyanin formed under nondenaturing conditions

Plastocyanin is a predominantly beta-sheet protein containing a type I copper center. The conformational ensemble of a denatured state of apo-plastocyanin formed in solution under conditions of low salt and neutral pH has been investigated by multidimensional heteronuclear NMR spectroscopy. Chemical shift assignments were obtained by using three-dimensional triple-resonance NMR experiments to trace through-bond heteronuclear connectivities along the backbone and side chains. The (3)J(HN,Halpha) coupling constants, (15)N-edited proton-proton nuclear Overhauser effects (NOEs), and (15)N relaxation parameters were also measured for the purpose of structural and dynamic characterization. Most of the residues corresponding to beta-strands in the folded protein exhibit small upfield shifts of the (13)C(alpha) and (13)CO resonances relative to random coil values, suggesting a slight preference for backbone dihedral angles in the beta region of (phi,psi) space. This is further supported by the presence of strong sequential d(alphaN)(i, i + 1) NOEs throughout the sequence. The few d(NN)(i, i + 1) proton NOEs that are observed are mostly in regions that form loops in the native plastocyanin structure. No medium or long-range NOEs were observed. A short sequence, between residues 59 and 63, was found to populate a nonnative helical conformation in the unfolded state, as indicated by the shift of the (13)C(alpha), (13)CO, and (1)H(alpha) resonances relative to random coil values and by the decreased values of the (3)J(HN,Halpha) coupling constants. The (15)N relaxation parameters indicate restriction of motions on a nanosecond timescale in this region. Intriguingly, this helical conformation is present in a sequence that is close to but not in the same location as the single short helix in the native folded protein. The results are consistent with earlier NMR studies of peptide fragments of plastocyanin and confirm that the regions of the sequence that form beta-strands in the native protein spontaneously populate the beta-region of (phi,psi) space under folding conditions, even in the absence of stabilizing tertiary interactions. We conclude that the state of apo-plastocyanin present under nondenaturing conditions is a noncompact unfolded state with some evidence of nativelike and nonnative local structuring that may be initiation sites for folding of the protein.     

Replacement of proline with valine does not remove an apparent proline isomerization-dependent folding event in CRABP I

Site-directed mutagenesis has frequently been used to replace proline with other amino acids in order to determine if proline isomerization is responsible for a slow phase during refolding. Replacement of Pro 85 with alanine in cellular retinoic acid binding protein I (CRABP-I) abolished the slowest refolding phase, suggesting that this phase is due to proline isomerization in the unfolded state. To further test this assumption, we mutated Pro 85 to valine, which is the conservative replacement in the two most closely related proteins in the family (cellular retinoic acid binding protein II and cellular retinol binding protein I). The mutant protein was about 1 kcal/mole more stable than wild type. Retinoic acid bound equally well to wild type and P85V-CRABP I, confirming the functional integrity of this mutation. The refolding and unfolding kinetics of the wild-type and mutant proteins were characterized by stopped flow fluorescence and circular dichroism. The mutant P85V protein refolded with three kinetic transitions, the same number as wild-type protein. This result conflicts with the P85A mutant, which lost the slowest refolding rate. The P85V mutation also lacked a kinetic unfolding intermediate found for wild-type protein. These data suggest that proline isomerization may not be responsible for the slowest folding phase of CRABP I. As such, the loss of a slow refolding phase upon mutation of a proline residue may not be diagnostic for proline isomerization effects on protein folding.     

Energy landscape of a peptide consisting of alpha-helix, 3(10)-helix, beta-turn, beta-hairpin, and other disordered conformations

The energy landscape of a peptide [Ace-Lys-Gln-Cys-Arg-Glu-Arg-Ala-Nme] in explicit water was studied with a multicanonical molecular dynamics simulation, and the AMBER parm96 force field was used for the energy calculation. The peptide was taken from the recognition helix of the DNA-binding protein, c-MYB: A rugged energy landscape was obtained, in which the random-coil conformations were dominant at room temperature. The CD spectra of the synthesized peptide revealed that it is in the random state at room temperature. However, the 300 K canonical ensemble, Q(300K), contained alpha-helix, 3(10)-helix, beta-turn, and beta-hairpin structures with small but notable probabilities of existence. The complete alpha-helix, imperfect alpha-helix, and random-coil conformations were separated from one another in the conformational space. This means that the peptide must overcome energy barriers to form the alpha-helix. The overcoming process may correspond to the hydrogen-bond rearrangements from peptide-water to peptide-peptide interactions. The beta-turn, imperfect 3(10)-helix, and beta-hairpin structures, among which there are no energy barriers at 300 K, were embedded in the ensemble of the random-coil conformations. Two types of beta-hairpin with different beta-turn regions were observed in Q(300K). The two beta-hairpin structures may have different mechanisms for the beta-hairpin formation. The current study proposes a scheme that the random state of this peptide consists of both ordered and disordered conformations. In contrast, the energy landscape obtained from the parm94 force field was funnel like, in which the peptide formed the helical conformation at room temperature and random coil at high temperature.     

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.