Position effect of cross-strand side-chain interactions on beta-hairpin formation

Roles of Amphipathicity and Hydrophobicity in the Micelle-Driven Structural Switch of a 14-mer Peptide Core from a Choline-Binding Repeat

Choline-binding repeats (CBRs) are ubiquitous sequences with a β-hairpin core that are found in the surface proteins of several microorganisms such as S. pneumoniae (pneumococcus). Previous studies on a 14-mer CBR sequence derived from the pneumoccal LytA autolysin (LytA_239-252 peptide) have demonstrated a switch behaviour for this peptide, so that it acquires a stable, native-like β-hairpin conformation in aqueous solution but is reversibly transformed into an amphipathic α-helix in the presence of detergent micelles. With the aim of understanding the factors responsible for this unusual β-hairpin to α-helix transition, and to specifically assess the role of peptide hydrophobicity and helical amphipathicity in the process, we designed a series of LytA_239-252 variants affecting these two parameters and studied their interaction with dodecylphosphocholine (DPC) micelles by solution NMR, circular dichroism and fluorescence spectroscopies. Our results indicate that stabilising cross-strand interactions become essential for β-hairpin stability in the absence of optimal turn sequences. Moreover, both amphipathicity and hydrophobicity display comparable importance for helix stabilisation of CBR-derived peptides in micelles, indicating that these sequences represent a novel class of micelle/membrane-interacting peptides.

Nature of aryl-tyrosine interactions contribute to β-hairpin scaffold stability: NMR evidence for alternate ring geometry

The specific contribution of the acidic-aromatic β-sheet favouring amino acid tyrosine to the stability of short octapeptide β-hairpin structures is presented here. Solution NMR analysis in near-apolar environments suggests the energetically favourable mode of interaction to be T-shaped face-to-edge (FtE) and that a Trp-Tyr interacting pair is the most stabilizing. Alternate aryl geometries also exist in solution, which readily equilibrate between a preferred π···π conformation to an aromatic-amide conformation, without any change in the backbone structure. While the phenolic ring is readily accommodated at the "edge" of FtE aryl interactions, it exhibits an overall lowered contribution to scaffold stability in the "face" orientation. Such differential tyrosine interactions are key to its dual nature in proteins.

Micelle-Triggered β-Hairpin to α-Helix Transition in a 14-Residue Peptide from a Choline-Binding Repeat of the Pneumococcal Autolysin LytA

Choline-binding modules (CBMs) have a ββ-solenoid structure composed of choline-binding repeats (CBR), which consist of a β-hairpin followed by a short linker. To find minimal peptides that are able to maintain the CBR native structure and to evaluate their remaining choline-binding ability, we have analysed the third β-hairpin of the CBM from the pneumococcal LytA autolysin. Circular dichroism and NMR data reveal that this peptide forms a highly stable native-like β-hairpin both in aqueous solution and in the presence of trifluoroethanol, but, strikingly, the peptide structure is a stable amphipathic α-helix in both zwitterionic (dodecylphosphocholine) and anionic (sodium dodecylsulfate) detergent micelles, as well as in small unilamellar vesicles. This β-hairpin to α-helix conversion is reversible. Given that the β-hairpin and α-helix differ greatly in the distribution of hydrophobic and hydrophilic side chains, we propose that the amphipathicity is a requirement for a peptide structure to interact and to be stable in micelles or lipid vesicles. To our knowledge, this "chameleonic" behaviour is the only described case of a micelle-induced structural transition between two ordered peptide structures.

Design of monomeric water-soluble β-hairpin and β-sheet peptides

Since the first report in 1993 (JACS 115, 5887-5888) of a peptide able to form a monomeric β-hairpin structure in aqueous solution, the design of peptides forming either β-hairpins (two-stranded antiparallel β-sheets) or three-stranded antiparallel β-sheets has become a field of growing interest and activity. These studies have yielded great insights into the principles governing the stability and folding of β-hairpins and antiparallel β-sheets. This chapter provides an overview of the reported β-hairpin/β-sheet peptides focussed on the applied design criteria, reviews briefly the factors contributing to β-hairpin/β-sheet stability, and describes a protocol for the de novo design of β-sheet-forming peptides based on them. Guidelines to select appropriate turn and strand residues and to avoid self-association are provided. The methods employed to check the success of new designed peptides are also summarized. Since NMR is the best technique to that end, NOEs and chemical shifts characteristic of β-hairpins and three-stranded antiparallel β-sheets are given.

Folding 19 proteins to their native state and stability of large proteins from a coarse-grained model

We develop an intermediate resolution model, where the backbone is modeled with atomic resolution but the side chain with a single bead, by extending our previous model (Proteins (2013) DOI: 10.1002/prot.24269) to properly include proline, preproline residues and backbone rigidity. Starting from random configurations, the model properly folds 19 proteins (including a mutant 2A3D sequence) into native states containing β sheet, α helix, and mixed α/β. As a further test, the stability of H-RAS (a 169 residue protein, critical in many signaling pathways) is investigated: The protein is stable, with excellent agreement with experimental B-factors. Despite that proteins containing only α helices fold to their native state at lower backbone rigidity, and other limitations, which we discuss thoroughly, the model provides a reliable description of the dynamics as compared with all atom simulations, but does not constrain secondary structures as it is typically the case in more coarse-grained models. Further implications are described.

Tryptophan residues: scarce in proteins but strong stabilizers of β-hairpin peptides

Tryptophan plays important roles in protein stability and recognition despite its scarcity in proteins. Except as fluorescent groups, they have been used rarely in peptide design. Nevertheless, Trp residues were crucial for the stability of some designed minimal proteins. In 2000, Trp-Trp pairs were shown to contribute more than any other hydrophobic interaction to the stability of β-hairpin peptides. Since then, Trp-Trp pairs have emerged as a paradigm for the design of stable β-hairpins, such as the Trpzip peptides. Here, we analyze the nature of the stabilizing capacity of Trp-Trp pairs by reviewing the β-hairpin peptides containing Trp-Trp pairs described up to now, the spectroscopic features and geometry of the Trp-Trp pairs, and their use as binding sites in β-hairpin peptides. To complete the overview, we briefly go through the other relevant β-hairpin stabilizing Trp-non-Trp interactions and illustrate the use of Trp in the design of short peptides adopting α-helical and mixed α/β motifs. This review is of interest in the field of rational design of proteins, peptides, peptidomimetics, and biomaterials.

LpxA: a natural nanotube

UDP-N-acetylglucosamine 3-O-acyltransferase is a protein with a left-handed parallel beta-helix, which is a natural nanotube. They are associated with unusual high stability. To identify the reason behind the structural stability of beta-helical nanotubular structure, we have performed a total of 4 mus molecular dynamics simulations of the protein in implicit solvent at four different temperatures and monitored the unfolding pathway. The correlation in movement between different regions of the nanotubular structure has been identified from the dynamical cross-correlation map and contribution of some specific residues towards unfolding transition has been identified by principal component analysis. Difference in stability of the three loop regions has also been characterized. Construction of the unfolding conformational energy landscape identifies the probable intermediates that can appear in the unfolding pathway of the protein.

Sequence inversion and phenylalanine surrogates at the beta-turn enhance the antibiotic activity of gramicidin S

A series of gramicidin S (GS) analogues have been synthesized where the Phe (i + 1) and Pro (i + 2) residues of the beta-turn have been swapped while the respective chiralities (D-, L-) at each position are preserved, and Phe is replaced by surrogates with aromatic side chains of diverse size, orientation, and flexibility. Although most analogues preserve the beta-sheet structure, as assessed by NMR, their antibiotic activities turn out to be highly dependent on the bulkiness and spatial arrangement of the aromatic side chain. Significant increases in microbicidal potency against both Gram-positive and Gram-negative pathogens are observed for several analogues, resulting in improved therapeutic profiles. Data indicate that seemingly minor replacements at the GS beta-turn can have significant impact on antibiotic activity, highlighting this region as a hot spot for modulating GS plasticity and activity.

The inherent flexibility of peptides and protein fragments quantitized by CD in conjunction with CCA+

ECD spectroscopy is traditionally used for rapid, non-atomic level structure analysis of natural products such as peptides and proteins. Unlike globular proteins, peptides less frequently adopt a single 3D-fold in a time average manner. Moreover, they exhibit an ensemble of conformers composed of a multitude of substantially different structures. In principle, both ECD- and vibrational circular dichroism (VCD)-spectroscopy are sensitive enough to pick up structural information on these dynamic ensembles. However, the interpretation of the raw spectral data of these highly dynamic molecular systems can be cumbersome. The herein presented Convex Constraint Analysis Plus method, or CCA+ for short (http://www.chem.elte.hu/departments/protnmr/cca/), provides a unique opportunity for spectral ensemble analysis of peptides, glycopeptides, peptidomimetics, and other foldamers. The precision and accuracy of the approach is presented here through different peptide model systems. An interesting temperature and pH dependent folding and unfolding of a miniprotein (e.g. Tc5b variant) is also described. Analysis of CD spectra sets strongly affected by solvent and ion type is also introduced to account for severe environmental-induced structure influencing effect(s). The deconvolution makes always possible the quantitative data analysis even when the interpretation of the deconvolution resulted in pure CD curves is complex.

Therapeutic index of gramicidin S is strongly modulated by D-phenylalanine analogues at the beta-turn

Analogues of the cationic antimicrobial peptide gramicidin S (GS), cyclo(Val-Orn-Leu-D-Phe-Pro)2, with d-Phe residues replaced by different (restricted mobility, mostly) surrogates have been synthesized and used in SAR studies against several pathogenic bacteria. While all D-Phe substitutions are shown by NMR to preserve the overall beta-sheet conformation, they entail subtle structural alterations that lead to significant modifications in biological activity. In particular, the analogue incorporating D-Tic (1,2,3,4-tetrahydroisoquinoline-3-carboxylic acid) shows a modest but significant increase in therapeutic index, mostly due to a sharp decrease in hemolytic effect. The fact that NMR data show a shortened distance between the D-Tic aromatic ring and the Orn delta-amino group may help explain the improved antibiotic profile of this analogue.

Dependence of folding dynamics and structural stability on the location of a hydrophobic pair in beta-hairpins

We study the dependence of folding time, nucleation site, and stability of a model beta-hairpin on the location of a cross-strand hydrophobic pair, using a coarse-grained off-lattice model with the aid of Monte Carlo simulations. Our simulations have produced 6500 independent folding trajectories dynamically, forming the basis for extensive statistical analysis. Four folding pathways, zipping-out, middle-out, zipping-in, and reptation, have been closely monitored and discussed in all seven sequences studied. A hydrophobic pair placed near the beta-turn or in the middle section effectively speed up folding; a hydrophobic pair placed close to the terminal ends or next to the beta-turn encourages stability of the entire chain.

NMR analysis of aromatic interactions in designed peptide beta-hairpins

Designed octapeptide beta-hairpins containing a central (D)Pro-Gly segment have been used as a scaffold to place the aromatic residues Tyr and Trp at various positions on the antiparallel beta-strands. Using a set of five peptide hairpins, aromatic interactions have been probed across antiparallel beta-sheets, in the non-hydrogen bonding position (Ac-L-Y-V-(D)P-G-L-Y/W-V-OMe: peptides 1 and 2), diagonally across the strands (Boc-Y/W-L-V-(D)P-G-W-L-V-OMe: peptides 3 and 6), and along the strands at positions i and i + 2 (Boc-L-L-V-(D)P-G-Y-L-W-OMe: peptide 4). Two peptides served as controls (Boc-L-L-V-(D)P-G-Y-W-V-OMe: peptide 5; Boc-L-Y-V-(D)P-G-L-L-V-OMe: peptide 7) for aromatic interactions. All studies have been carried out using solution NMR methods in CDCl(3) + 10% DMSO-d(6) and have been additionally examined in CD(3)OH for peptides 1 and 2. Inter-ring proton-proton nuclear Overhauser effects (NOEs) and upfield shifted aromatic proton resonances have provided firm evidence for specific aromatic interactions. Calculated NMR structures for peptides 1 and 2, containing aromatic pairs at facing non-hydrogen bonded positions, revealed that T-shaped arrangements of the interacting pairs of rings are favored, with ring current effects leading to extremely upfield chemical shifts and temperature dependences for specific aromatic protons. Anomalous far-UV CD spectra appeared to be a characteristic feature in peptides where the two aromatic residues are spatially proximal. The observation of the close approach of aromatic rings in organic solvents suggests that interactions of an electrostatic nature may be favored. This situation may be compared to the case of aqueous solutions, where clustering of aromatic residues is driven by solvophobic (hydrophobic) forces.

Calorimetric evidence for a two-state unfolding of the beta-hairpin peptide trpzip4

beta-Sheets are a common secondary structural element found in proteins. The difficulty in studying beta-sheet folding and stability is that their formation is often dependent on the tertiary structural environment within the protein. However, the discovery of water-soluble beta-hairpin peptides has allowed them to be used as model systems because they represent the smallest units of beta-sheet structure independent of tertiary structural context. Trpzip4 has been used as a model beta-hairpin peptide to study beta-hairpin folding and stability because it is highly soluble in aqueous solutions, maintains its monomeric state, and shows reversible cooperative thermal unfolding. The previously determined thermodynamic parameters for trpzip4 thermal unfolding vary depending on the spectroscopic probe used, which questions the assumption that trpzip4 unfolds in a two-state manner. Here we provide direct calorimetric evidence that the unfolding of trpzip4 follows a two-state unfolding mode. Furthermore, the thermal unfolding of trpzip4 monitored using near- and far-UV-CD yielded thermodynamic parameters similar to those determined calorimetrically, providing additional evidence for a two-state unfolding mode.

β-Hairpin formation in aqueous solution and in the presence of trifluoroethanol: A1H and13C nuclear magnetic resonance conformational study of designed peptides

In order to check our current knowledge on the principles involved in beta-hairpin formation, we have modified the sequence of a 3:5 beta-hairpin forming peptide with two different purposes, first to increase the stability of the formed 3:5 beta-hairpin, and second to convert the 3:5 beta-hairpin into a 2:2 beta-hairpin. The conformational behavior of the designed peptides was investigated in aqueous solution and in 30% trifluoroethanol (TFE) by analysis of the following nuclear magnetic resonance (NMR) parameters: nuclear Overhauser effect (NOE) data, and C(alpha)H, (13)C(alpha), and (13)C(beta) conformational shifts. From the differences in the ability to adopt beta-hairpin structures in these peptides, we have arrived to the following conclusions: (i) beta-Hairpin population increases with the statistical propensity of residues to occupy each turn position. (ii) The loop length, and in turn, the beta-hairpin type, can be modified as a function of the type of turn favored by the loop sequence. These two conclusions reinforce previous results about the importance of beta-turn sequence in beta-hairpin folding. (iii) Side-chain packing on each face of the beta-sheet may play a major role in beta-hairpin stability; hence simplified analysis in terms of isolated pair interactions and intrinsic beta-sheet propensities is insufficient. (iv) Contributions to beta-hairpin stability of turn and strand sequences are not completely independent. (v) The burial of hydrophobic surface upon beta-hairpin formation that, in turn, depends on side-chain packing also contributes to beta-hairpin stability. (vi) As previously observed, TFE stabilizes beta-hairpin structures, but the extent of the contribution of different factors to beta-hairpin formation is sometimes different in aqueous solution and in 30% TFE.

Quantifying amino acid conformational preferences and side-chain-side-chain interactions in beta-hairpins

The intrinsic conformational biases of individual amino acids and their interstrand side-chain-side-chain (SC-SC) interactions both contribute to the stability of beta-sheets. The relative magnitudes of these effects have been difficult to assess in the context of folded proteins, where tertiary contacts complicate the quantitative analysis of local effects. We now report the results of such an analysis in a much simpler system, a short, stabilized beta-hairpin structure where intrastrand (conformational) and interstrand (SC-SC) influences can be distinguished in the absence of competing protein tertiary interactions. A comprehensive comparison of all pairwise combinations of 11 N-terminal and 7 C-terminal amino acids within an 8-residue, @-tide-stabilized [in which @ denotes the 1,2-dihydro-3(6H)-pyridinyl unit] beta-hairpin reveals distinct differences between the various pairings and shows that the intrastrand and interstrand effects are of comparable magnitude in contributing to the stability of the folded forms over the unfolded forms.

Characterization of a possible amyloidogenic precursor in glutamine-repeat neurodegenerative diseases

Several neurodegenerative diseases are linked to expanded repeats of glutamine residues, which lead to the formation of amyloid fibrils and neuronal death. The length of the repeats correlates with the onset of Huntington's disease, such that healthy individuals have <38 residues and individuals with >38 repeats exhibit symptoms. Because it is difficult to obtain atomic-resolution structural information for poly(l-glutamine) (polyQ) in aqueous solution experimentally, we performed molecular dynamics simulations to investigate the conformational behavior of this homopolymer. In simulations of 20-, 40-, and 80-mer polyQ, we observed the formation of the "alpha-extended chain" conformation, which is characterized by alternating residues in the alpha(L) and alpha(R) conformations to yield a sheet. The structural transition from disordered random-coil conformations to the alpha-extended chain conformation exhibits modest length and temperature dependence, in agreement with the experimental observation that aggregation depends on length and temperature. We propose that fibril formation in polyQ may occur through an alpha-sheet structure, which was proposed by Pauling and Corey. Also, we propose an atomic-resolution model of how the inhibitory peptide QBP1 (polyQ-binding peptide 1) may bind to polyQ in an alpha-extended chain conformation to inhibit fibril formation.

IR Study of Cross-Strand Coupling in a β-Hairpin Peptide Using Isotopic Labels

Model beta-hairpin peptides can be used to develop understanding of fundamental elements of beta-sheet secondary structure formation and stability. We have studied two 13C-labeled variants of a beta-hairpin peptide modified from a design originally proposed by Gellman: Arg-Tyr-Val-Glu-Val-Aib-Gly-Lys-Lys-Ile-Leu-Gln. (In this peptide, the two italicized residues form a beta-turn, while 13C-labels are on the amide C=O of Val3, Lys8 in HBG-L and Val3, Ile10 in HBG-S.) Both these peptides are labeled on opposite strands of the hairpin, but differ in the labeling pattern. One (HBG-L) forms a large (14-atom) H-bonded ring of labeled C=Os, while the other (HBG-S) forms a small (10-atom) H-bonded ring. These impact the amide I infrared spectra, with HBG-L having a 13C frequency and intensity higher than that of HBG-S, in good agreement with our spectral simulations based on quantum mechanically derived force fields. The thermal behavior of both peptides yields a broad thermal transition and lacks an isosbestic point. The 13C band for HBG-L has the largest intensity change with temperature, distinct from the 12C change and the HBG-S 13C change.

Applications of model β‐hairpin peptides

In recent years, beta-hairpin peptides have been studied in great detail. Much of the focus has been on the thermodynamic stability of beta-hairpin structure. Structural measurements have been conducted with nuclear magnetic resonance, with additional information obtained from circular dichroism, Fourier transform infrared, and molecular dynamic simulation studies. Point mutations, both in the beta-strands and in the turn region, have systematically explored the role of turn sequence, side-chain-side-chain interactions, intramolecular hydrogen bonding, and beta-strand length on beta-hairpin peptide conformational stability. In addition to studying the elements of structural stability independently, the cooperative nature of the individual components to combine to form the overall structure has also been investigated. Because the beta-hairpin peptides often spontaneously form their conformation, they have begun to serve as models for studying peptide binding and therapeutic agents.

The role of the unfolded state in hairpin stability

The effects of a T3S mutation on the stability of a 3:5 beta-hairpin forming peptide (YITNSNGTWT) are investigated. Molecular dynamics simulations in explicit water indicate that the wild-type peptide forms a stable hairpin whereas the T3S mutant does not, in agreement with the experimental data. Thermodynamic integration calculations for the mutation of Thr to Ser suggest that the free-energy changes in the folded state are small, but the corresponding changes in the unfolded state are large and favorable. One of the main reasons for the difference appears to be the formation of a stable cluster involving the Tyr1 and Ser3 hydroxyl groups and their interaction with the C-terminal carboxylate group, which was observed after unfolding of the T3S mutant. Further analysis of the side-chain preferences of Thr and Ser indicate that the corresponding cluster in the wild-type peptide is unstable due to the high preference of the Thr chi1 dihedral for g+ states, which appeared to be incompatible with formation of a stable cluster. The results suggest that one should consider the nature of the unfolded state when attempting to fully explain the effects of mutations on hairpin stability.

Sequence dependence of beta-hairpin structure: comparison of a salt bridge and an aromatic interaction

A comparison of the contributions and position dependence of cross-strand electrostatic and aromatic side-chain interactions to beta-sheet stability has been performed by using nuclear magnetic resonance in a well-folded beta-hairpin peptide of the general sequence XRTVXVdPGOXITQX. Phe-Phe and Glu-Lys pairs were varied at the internal and terminal non-hydrogen-bonded position, and the resulting stability was measured by the effects on alpha-hydrogen and aromatic hydrogen chemical shifts. It was determined that the introduction of a Phe-Phe pair resulted in a more folded peptide, regardless of position, and a more tightly folded core. Substitution of the Glu-Lys pair at the internal position results in a less folded peptide and increased fraying at the terminal residues. Upfield shifting of the aromatic hydrogens provided evidence for an edge-face aromatic interaction, regardless of position of the Phe-Phe pair. In peptides with two Phe-Phe pairs, substitution with Glu-Lys at either position resulted in a weakening of the aromatic interaction and a subsequent decrease in peptide stability. Thermal denaturation of the peptides containing Phe-Phe indicates that the aromatic interaction is enthalpically favored, whereas the folding of hairpins with cross-strand Glu-Lys pairs was less enthalpically favorable but entropically more favorable.

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.

Expected and unexpected results from combined β-hairpin design elements

A model beta-hairpin dodecapeptide [EFGWVpGKWTIK] was designed by including a favorable D-ProGly Type II' beta-turn sequence and a Trp-zip interaction, while also incorporating a beta-strand unfavorable glycine residue in the N-terminal strand. This peptide is highly folded and monomeric in aqueous solution as determined by combined analysis with circular dichroism and 1H NMR spectroscopy. A peptide representing the folded conformation of the model beta-hairpin [cyclic(EFGWVpGKWTIKpG)] and a linear peptide representing the unfolded conformation [EFGWVPGKWTIK] yield unexpected relative deviations between the CD and 1H NMR spectroscopic results that are attributed to variations in the packing interactions of the aromatic side chains. Mutational analysis of the model beta-hairpin indicates that the Trp-zip interaction favors folding and stability relative to an alternate hydrophobic cluster between Trp and Tyr residues [EFGYVpGKWTIK]. The significance of select diagonal interactions in the model beta-hairpin was tested by rearranging the cross-strand hydrophobic interactions to provide a folded peptide [EWFGIpGKTYWK] displaying evidence of an unusual backbone conformation at the hydrophobic cluster. This unusual conformation does not appear to be a result of the glycine residue in the beta-strand, as replacement with a serine results in a peptide [EWFSIpGKTYWK] with a similar and seemingly characteristic CD spectrum. However, an alternate arrangement of hydrophobic residues with a Trp-zip interaction in a similar position to the parent beta-hairpin [EGFWVpGKWITK] results in a folded beta-hairpin conformation. The differences between side chain packing of these peptides precludes meaningful thermodynamic analysis and illustrates the caution necessary when interpreting beta-hairpin folding thermodynamics that are driven, at least in part, by aromatic cross strand interactions.

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.

The energy landscape of unsolvated peptides: the role of context in the stability of alanine/glycine helices

Ion mobility measurements have been used to examine the conformations present for unsolvated Ac-(AG)(7)A+H(+) and (AG)(7)A+H(+) peptides (Ac = acetyl, A = alanine, and G = glycine) over a broad temperature range (100-410 K). The results are compared to those recently reported for Ac-A(4)G(7)A(4)+H(+) and A(4)G(7)A(4)+H(+), which have the same compositions but different sequences. Ac-(AG)(7)A+H(+) shows less conformational diversity than Ac-A(4)G(7)A(4)+H(+); it is much less helical than Ac-A(4)G(7)A(4)+H(+) at the upper end of the temperature range studied, and at low temperatures, one of the two Ac-A(4)G(7)A(4)+H(+) features assigned to helical conformations is missing for Ac-(AG)(7)A+H(+). Molecular dynamics simulations suggest that the different conformational preferences are not due to differences in the stabilities of the helical states, but differences in the nonhelical states: it appears that Ac-(AG)(7)A+H(+) is more flexible and able to adopt lower energy globular conformations (compact random looking three-dimensional structures) than Ac-A(4)G(7)A(4)+H(+). The helix to globule transition that occurs for Ac-(AG)(7)A+H(+) at around 250-350 K is not a direct (two-state) process, but a creeping transition that takes place through at least one and probably several intermediates.

Insights into the determinants of beta-sheet stability: 1H and 13C NMR conformational investigation of three-stranded antiparallel beta-sheet-forming peptides

In a previous study we designed a 20-residue peptide able to adopt a significant population of a three-stranded antiparallel beta-sheet in aqueous solution (de Alba et al. [1999]Protein Sci.8, 854-865). In order to better understand the factors contributing to beta-sheet folding and stability we designed and prepared nine variants of the parent peptide by substituting residues at selected positions in its strands. The ability of these peptides to form the target motif was assessed on the basis of NMR parameters, in particular NOE data and 13Calpha conformational shifts. The populations of the target beta-sheet motif were lower in the variants than in the parent peptide. Comparative analysis of the conformational behavior of the peptides showed that, as expected, strand residues with low intrinsic beta-sheet propensities greatly disfavor beta-sheet folding and that, as already found in other beta-sheet models, specific cross-strand side chain-side chain interactions contribute to beta-sheet stability. More interestingly, the performed analysis indicated that the destabilization effect of the unfavorable strand residues depends on their location at inner or edge strands, being larger at the latter. Moreover, in all the cases examined, favorable cross-strand side chain-side chain interactions were not strong enough to counterbalance the disfavoring effect of a poor beta-sheet-forming residue, such as Gly.

Thermodynamic analysis of beta-hairpin-forming peptides from the thermal dependence of (1)H NMR chemical shifts

The temperature dependence of the (1)H chemical shifts of six designed peptides previously shown to adopt beta-hairpin structures in aqueous solution has been analyzed in terms of two-state (beta-hairpin left arrow over right arrow coil) equilibrium. The stability of the beta-hairpins formed by these peptides, as derived from their T(m) (midpoint transition temperature) values, parallels in general their ability to adopt those structures as deduced from independent NMR parameters: NOEs, Deltadelta(C)(alpha)(H), Deltadelta(C)(alpha), and Deltadelta(C)(beta) values. The observed T(m) values are dependent on the particular position within the beta-hairpin that is probed, indicating that their folding to a beta-hairpin conformation deviates from a "true" two-state transition. To obtain individual T(m) values for each hairpin region in each peptide, a simplified model of a successive uncoupled two-state equilibrium covering the entire process has been applied. The distribution of T(m) values obtained for the different beta-hairpin regions (turn, strands, backbone, side chains) in the six analyzed peptides reveals a similar pattern. A model for beta-hairpin folding is proposed on the basis of this pattern and the reasonable assumption that regions showing higher T(m) values are the last ones to unfold and, presumably, the first to form. With this assumption, the analysis suggests that turn formation is the first event in beta-hairpin folding. This is consistent with previous results on the essential role of the turn sequence in beta-hairpin folding.

Pescador: the PEptides in Solution ConformAtion Database: Online Resource

In recent years a large body of data has been obtained from Nuclear Magnetic Resonance and Circular Dichroism experiments on the influence of the amino acid sequence and various other parameters on the conformational state of peptides in solution. Interpreting the experimental data in terms of the conformational populations of the peptides remains a key problem, for which current solutions leave appreciable room for improvement. Considering that making this body of data available for surveys and analysis should be instrumental in tackling the problem, we undertook the development of Pescador: The 'PEptides in Solution ConformAtion Database: Online Resource'. Pescador contains data from NMR and CD spectroscopy on peptides in solution as well as information on the structural parameters derived from these data. It also features specialized Web-based tools for data deposition, and means for readily accessing the stored information for analysis purposes. To illustrate the use of the database in deriving information for the conformational analysis of peptides, we show how the alpha proton delta-values stored in Pescador and measured by NMR for different peptides in different laboratories can be used to derive a new set of 'random coil' chemical shift values. Firstly, we show these values to be very similar to those obtained experimentally for model peptides in water, and their variation with increasing Tri-Fluoro-Ethanol (TFE) concentration is similar to that reported for model peptides. We show, furthermore, that the chemical shift data in Pescador can be used to derive correction factors that take into account effects of neighboring residues. These correction factors compare favorably with those recently derived from a series of model GGXGG peptides (Schwarzinger et al., 2001). These encouraging results suggest that, as the quantity of NMR data on peptide deposited in Pescador increases, surveys of these data should be a valuable means of deriving key parameters for the analysis of peptide conformation.

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.

Length-dependent stability and strand length limits in antiparallel beta -sheet secondary structure

Designed peptides that fold autonomously to specific conformations in aqueous solution are useful for elucidating protein secondary structural preferences. For example, autonomously folding model systems have been essential for establishing the relationship between alpha-helix length and alpha-helix stability, which would be impossible to probe with alpha-helices embedded in folded proteins. Here, we use designed peptides to examine the effect of strand length on antiparallel beta-sheet stability. alpha-Helices become more stable as they grow longer. Our data show that a two-stranded beta-sheet ("beta-hairpin") becomes more stable when the strands are lengthened from five to seven residues, but that further strand lengthening to nine residues does not lead to further beta-hairpin stabilization for several extension sequences examined. (In one case, all-threonine extension, there may be an additional stabilization on strand lengthening from seven to nine residues.) These results suggest that there may be an intrinsic limit to strand length for most sequences in antiparallel beta-sheet secondary structure.

Computer-aided design of beta-sheet peptides

The design of beta-sheet proteins is still a challenge in the field of de novo protein design. Here, we have tested the validity of automatic design methods to create and/or improve beta-sheet peptides and proteins. We chose Betanova, a three-stranded beta-sheet peptide, as target system, and, as an automatic design tool, a protein design algorithm called PERLA (protein engineering rotamer library algorithm). PERLA was used to define both stabilising and destabilising single- and multiple-residue mutations of Betanova. Conformational analysis by NMR spectroscopy and far-UV circular dichroism (CD) allowed us to evaluate population differences among the set of designed peptides. Some of the new mutants are approximately 1 kcal/mol more stable than the wild-type peptide. Comparison of the scale of predicted and observed stabilities demonstrates that they are in good agreement for most peptides studied. Our results show that automatic design algorithms can be successfully applied to the design of beta-sheet peptides.

13C(alpha) and 13C(beta) chemical shifts as a tool to delineate beta-hairpin structures in peptides

Unravelling the factors that contribute to the formation and the stability of beta-sheet structure in peptides is a subject of great current interest. A beta-hairpin, the smallest beta-sheet motif, consists of two antiparallel hydrogen-bonded beta-strands linked by a loop region. We have performed a statistical analysis on protein beta-hairpins showing that the most abundant types of beta-hairpins, 2:2, 3:5 and 4:4, have characteristic patterns of 13C(alpha) and 13C(beta) conformational shifts, as expected on the basis of their phi and psi angles. This fact strongly supports the potential value of 13C(alpha) and 13C(beta) conformational shifts as a means to identify beta-hairpin motifs in peptides. Their usefulness was confirmed by analysing the patterns of 13C(alpha) and 13C(beta) conformational shifts in 13 short peptides, 10-15 residues long, that adopt beta-hairpin structures in aqueous solution. Furthermore, we have investigated their potential as a method to quantify beta-hairpin populations in peptides.