Molecular dynamics simulations of certain mutant peptide models from staphylococcal nuclease reveal that initial hydrophobic collapse associated with turn propensity drives β‐hairpin folding

An important nucleation event during the folding of staphylococcal nuclease involves the formation of a β‐hairpin by the sequence ²¹DTVKLMYKGQPMTFR³⁵. Earlier studies show that the turn sequence ‘YKGQP’ has an important role in the folding of this β‐hairpin. To understand the active or passive nature of the turn sequence ‘YKGQP’ in the folding of the aforementioned β‐hairpin sequence, we studied glycine mutant peptides Ac‐²DTVKLMYGGQPMTFR¹⁶‐NMe (K9G:15), Ac‐²DTVKLMYKGGPMTFR¹⁶‐NMe (Q11G:15), Ac‐²DTVKLMYGGGPMTFR¹⁶‐NMe (K9G/Q11G:15), and Ac‐²DTVKLMGGGGGMTFR¹⁶‐NMe (penta‐G:15) by using molecular dynamics simulations, starting with two different unfolded states, polyproline II and extended conformational forms. Further, 5mer mutant turn peptides Ac‐²YGGQP⁶‐NMe (K3G:5), Ac‐²YKGGP⁶‐NMe (Q5G:5), Ac‐²YGGGP⁶‐NMe (K3G/Q5G:5), and Ac‐²GGGGG⁶‐NMe (penta‐G:5) were also studied individually. Our results show that an initial hydrophobic collapse and loop closure occurs in all 15mer mutants, but only K9G:15 mutant forms a stable native‐like β‐hairpin. In the other 15mer mutants, the hydrophobic collapsed state would not proceed to β‐hairpin formation. Of the different simulations performed for the penta‐G:15 mutant, in only one simulation a nonnative β‐hairpin conformation is sampled with highly flexible loop region (⁸GGGGG¹²), which has no specific conformational preference as a 5mer. While the sequence ‘YGGQP’ in the K3G:5 simulation shows relatively higher β‐turn propensity, the presence of this sequence in K9G:15 peptide seems to be driving the β‐hairpin formation. Thus, these results seem to suggest that for the formation of a stable β‐hairpin, the initial hydrophobic collapse is to be assisted by a turn propensity. Initial hydrophobic collapse alone is not sufficient to guide β‐hairpin formation. Copyright © 2013 European Peptide Society and John Wiley & Sons, Ltd.     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. III. Dynamics of long‐range hydrophobic interactions

A 20‐residue peptide, IG(42–61), derived from the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptoccocus was studied using circular dichroism, nuclear magnetic resonance (NMR) spectroscopy at various temperatures and by differential scanning calorimetry (DSC). Unlike other related peptides studied so far, this peptide displays two heat capacity peaks in DSC measurements (at a scanning rate of 1.5 deg/min at a peptide concentration of 0.07 mM), which suggests a three‐state folding/unfolding process. The results from DSC and NMR measurements suggest the formation of a dynamic network of hydrophobic interactions stabilizing the structure, which resembles a β‐hairpin shape over a wide range of temperatures (283–313 K). Our results show that IG (42–61) possesses a well‐organized three‐dimensional structure stabilized by long‐range hydrophobic interactions (Tyr50 ··· Phe57 and Trp48 ··· Val59) at T = 283 K and (Trp48 ··· Val59) at 305 and 313 K. The mechanism of β‐hairpin folding and unfolding, as well as the influence of peptide length on its conformational properties, are also discussed. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

β‐Hairpin folding and stability: molecular dynamics simulations of designed peptides in aqueous solution

The structural properties of a 10‐residue and a 15‐residue peptide in aqueous solution were investigated by molecular dynamics simulation. The two designed peptides, SYINSDGTWT and SESYINSDGTWTVTE, had been studied previously by NMR at 278 K and the resulting model structures were classified as 3:5 β‐hairpins with a type I + G1 β‐bulge turn. In simulations at 278 K, starting from the NMR model structure, the 3:5 β‐hairpin conformers proved to be stable over the time period evaluated (30 ns). Starting from an extended conformation, simulations of the decapeptide at 278 K, 323 K and 353 K were also performed to study folding. Over the relatively short time scales explored (30 ns at 278 K and 323 K, 56 ns at 353 K), folding to the 3:5 β‐hairpin could only be observed at 353 K. At this temperature, the collapse to β‐hairpin‐like conformations is very fast. The conformational space accessible to the peptide is entirely dominated by loop structures with different degrees of β‐hairpin character. The transitions between different types of ordered loops and β‐hairpins occur through two unstructured loop conformations stabilized by a single side‐chain interaction between Tyr2 and Trp9, which facilitates the changes of the hydrogen‐bond register. In agreement with previous experimental results, β‐hairpin formation is initially driven by the bending propensity of the turn segment. Nevertheless, the fine organization of the turn region appears to be a late event in the folding process. Copyright © 2004 European Peptide Society and John Wiley & Sons, Ltd.     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. I. Importance of hydrophobic interactions in stabilization of β‐hairpin structure

We previously studied a 16‐amino acid‐residue fragment of the C‐terminal β‐hairpin of the B3 domain (residues 46–61), [IG(46–61)] of the immunoglobulin binding protein G from Streptoccocus, and found that hydrophobic interactions and the turn region play an important role in stabilizing the structure. Based on these results, we carried out systematic structural studies of peptides derived from the sequence of IG (46–61) by systematically shortening the peptide by one residue at a time from both the C‐ and the N‐terminus. To determine the structure and stability of two resulting 12‐ and 14‐amino acid‐residue peptides, IG(48–59) and IG(47–60), respectively, we carried out circular dichroism, NMR, and calorimetric studies of these peptides in pure water. Our results show that IG(48–59) possesses organized three‐dimensional structure stabilized by hydrophobic interactions (Tyr50–Phe57 and Trp48–Val59) at T = 283 and 305 K. At T = 313 K, the structure breaks down because of increased chain entropy, but the turn region is preserved in the same position observed for the structure of the whole protein. The breakdown of structure occurs near the melting temperature of this peptide (T_m = 310 K) measured by differential scanning calorimetry (DSC). The melting temperature of IG(47–60) determined by DSC is T_m = 330 K and its structure is similar to that of the native β‐hairpin at all (lower) temperatures examined (283–313 K). Both of these truncated sequences are conserved in all known amino acid sequences of the B domains of the immunoglobulin binding protein G from bacteria. Thus, this study contributes to an understanding of the mechanism of folding of this whole family of proteins, and provides information about the mechanism of formation and stabilization of a β‐hairpin structural element. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

β‐strand interactions at the domain interface critical for the stability of human lens γD‐crystallin

Human age‐onset cataracts are believed to be caused by the aggregation of partially unfolded or covalently damaged lens crystallin proteins; however, the exact molecular mechanism remains largely unknown. We have used microseconds of molecular dynamics simulations with explicit solvent to investigate the unfolding process of human lens γD‐crystallin protein and its isolated domains. A partially unfolded folding intermediate of γD‐crystallin is detected in simulations with its C‐terminal domain (C‐td) folded and N‐terminal domain (N‐td) unstructured, in excellent agreement with biochemical experiments. Our simulations strongly indicate that the stability and the folding mechanism of the N‐td are regulated by the interdomain interactions, consistent with experimental observations. A hydrophobic folding core was identified within the C‐td that is comprised of a and b strands from the Greek key motif 4, the one near the domain interface. Detailed analyses reveal a surprising non‐native surface salt‐bridge between Glu135 and Arg142 located at the end of the ab folded hairpin turn playing a critical role in stabilizing the folding core. On the other hand, an in silico single E135A substitution that disrupts this non‐native Glu135‐Arg142 salt‐bridge causes significant destabilization to the folding core of the isolated C‐td, which, in turn, induces unfolding of the N‐td interface. These findings indicate that certain highly conserved charged residues, that is, Glu135 and Arg142, of γD‐crystallin are crucial for stabilizing its hydrophobic domain interface in native conformation, and disruption of charges on the γD‐crystallin surface might lead to unfolding and subsequent aggregation.     

Exploring the free energy landscape of a model β‐hairpin peptide and its isoform

Secondary structural transitions from α‐helix to β‐sheet conformations are observed in several misfolding diseases including Alzheimer's and Parkinson's. Determining factors contributing favorably to the formation of each of these secondary structures is therefore essential to better understand these disease states. β‐hairpin peptides form basic components of anti‐parallel β‐sheets and are suitable model systems for characterizing the fundamental forces stabilizing β‐sheets in fibrillar structures. In this study, we explore the free energy landscape of the model β‐hairpin peptide GB1 and its E2 isoform that preferentially adopts α‐helical conformations at ambient conditions. Umbrella sampling simulations using all‐atom models and explicit solvent are performed over a large range of end‐to‐end distances. Our results show the strong preference of GB1 and the E2 isoform for β‐hairpin and α‐helical conformations, respectively, consistent with previous studies. We show that the unfolded states of GB1 are largely populated by misfolded β‐hairpin structures which differ from each other in the position of the β‐turn. We discuss the energetic factors contributing favorably to the formation of α‐helix and β‐hairpin conformations in these peptides and highlight the energetic role of hydrogen bonds and non‐bonded interactions. Proteins 2014; 82:2394–2402. © 2014 Wiley Periodicals, Inc.     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin‐binding protein G from Streptococcus. IV. Implication for the mechanism of folding of the parent protein

A 34‐residue α/β peptide [IG(28–61)], derived from the C‐terminal part of the B3 domain of the immunoglobulin binding protein G from Streptoccocus, was studied using CD and NMR spectroscopy at various temperatures and by differential scanning calorimetry. It was found that the C‐terminal part (a 16‐residue‐long fragment) of this peptide, which corresponds to the sequence of the β‐hairpin in the native structure, forms structure similar to the β‐hairpin only at T = 313 K, and the structure is stabilized by non‐native long‐range hydrophobic interactions (Val47–Val59). On the other hand, the N‐terminal part of IG(28–61), which corresponds to the middle α‐helix in the native structure, is unstructured at low temperature (283 K) and forms an α‐helix‐like structure at 305 K, and only one helical turn is observed at 313 K. At all temperatures at which NMR experiments were performed (283, 305, and 313 K), we do not observe any long‐range connectivities which would have supported packing between the C‐terminal (β‐hairpin) and the N‐terminal (α‐helix) parts of the sequence. Such interactions are absent, in contrast to the folding pathway of the B domain of protein G, proposed recently by Kmiecik and Kolinski (Biophys J 2008, 94, 726–736), based on Monte‐Carlo dynamics studies. Alternative folding mechanisms are proposed and discussed. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 469–480, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Mechanism of formation of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. II. Interplay of local backbone conformational dynamics and long‐range hydrophobic interactions in hairpin formation

Two peptides, corresponding to the turn region of the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus, consisting of residues 51–56 [IG(51–56)] and 50–57 [IG(50–57)], respectively, were studied by circular dichroism and NMR spectroscopy at various temperatures and by differential scanning calorimetry. Our results show that the part of the sequence corresponding to the β‐turn in the native structure (DDATKT) of the B3 domain forms bent conformations similar to those observed in the native protein. The formation of a turn is observed for both peptides in a broad range of temperatures (T = 283–323 K), which confirms the conclusion drawn from our previous studies of longer sequences from the C‐terminal β‐hairpin of the B3 domain of the immunoglobulin binding protein G (16, 14, and 12 residues), that the DDATKT sequence forms a nucleation site for formation of the β‐hairpin structure of peptides corresponding to the C‐terminal part of all the B domains of the immunoglobulin binding protein G. We also show and discuss the role of long‐range hydrophobic interactions as well as local conformational properties of polypeptide chains in the mechanism of formation of the β‐hairpin structure. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

A molecular dynamics study of the 41-56 beta-hairpin from B1 domain of protein G.

The structural and dynamical behavior of the 41-56 beta-hairpin from the protein G B1 domain (GB1) has been studied at different temperatures using molecular dynamics (MD) simulations in an aqueous environment. The purpose of these simulations is to establish the stability of this hairpin in view of its possible role as a nucleation site for protein folding. The conformation of the peptide in the crystallographic structure of the protein GB1 (native conformation) was lost in all simulations. The new equilibrium conformations are stable for several nanoseconds at 300K (>10 ns), 350 K (>6.5 ns), and even at 450 K (up to 2.5 ns). The new structures have very similar hairpin-like conformations with properties in agreement with available experimental nuclear Overhauser effect (NOE) data. The stability of the structure in the hydrophobic core region during the simulations is consistent with the experimental data and provides further evidence for the role played by hydrophobic interactions in hairpin structures. Essential dynamics analysis shows that the dynamics of the peptide at different temperatures spans basically the same essential subspace. The main equilibrium motions in this subspace involve large fluctuations of the residues in the turn and ends regions. Of the six interchain hydrogen bonds, the inner four remain stable during the simulations. The space spanned by the first two eigenvectors, as sampled at 450 K, includes almost all of the 47 different hairpin structures found in the database. Finally, analysis of the hydration of the 300 K average conformations shows that the hydration sites observed in the native conformation are still well hydrated in the equilibrium MD ensemble.     

Basic conformers in β‐peptides

The conformation of oligomers of β‐amino acids of the general type Ac‐[β‐Xaa]_n‐NHMe (β‐Xaa = β‐Ala, β‐Aib, and β‐Abu; n = 1–4) was systematically examined at different levels of ab initio molecular orbital theory (HF/6‐31G*, HF/3‐21G). The solvent influence was considered employing two quantum‐mechanical self‐consistent reaction field models. The results show a wide variety of possibilities for the formation of characteristic elements of secondary structure in β‐peptides. Most of them can be derived from the monomer units of blocked β‐peptides with n = 1. The stability and geometries of the β‐peptide structures are considerably influenced by the side‐chain positions, by the configurations at the C^α‐ and C^β‐atoms of the β‐amino acid constituents, and especially by environmental effects. Structure peculiarities of β‐peptides, in particular those of various helix alternatives, are discussed in relation to typical elements of secondary structure in α‐peptides. © 1999 John Wiley & Sons, Inc. Biopoly 50: 167–184, 1999     

Redesigning the type II' β‐turn in green fluorescent protein to type I': Implications for folding kinetics and stability

Both Type I' and Type II' β‐turns have the same sense of the β‐turn twist that is compatible with the β‐sheet twist. They occur predominantly in two residue β‐hairpins, but the occurrence of Type I' β‐turns is two times higher than Type II' β‐turns. This suggests that Type I' β‐turns may be more stable than Type II' β‐turns, and Type I' β‐turn sequence and structure can be more favorable for protein folding than Type II' β‐turns. Here, we redesigned the native Type II' β‐turn in GFP to Type I' β‐turn, and investigated its effect on protein folding and stability. The Type I' β‐turns were designed based on the statistical analysis of residues in natural Type I' β‐turns. The substitution of the native “GD” sequence of i+1 and i+2 residues with Type I' preferred “(N/D)G” sequence motif increased the folding rate by 50% and slightly improved the thermodynamic stability. Despite the enhancement of in vitro refolding kinetics and stability of the redesigned mutants, they showed poor soluble expression level compared to wild type. To overcome this problem, i and i + 3 residues of the designed Type I' β‐turn were further engineered. The mutation of Thr to Lys at i + 3 could restore the in vivo soluble expression of the Type I' mutant. This study indicates that Type II' β‐turns in natural β‐hairpins can be further optimized by converting the sequence to Type I'. Proteins 2014; 82:2812–2822. © 2014 Wiley Periodicals, Inc.     

Acyclic peptides incorporating the d‐Phe‐2‐Abz turn motif: Investigations on antimicrobial activity and propensity to adopt β‐hairpin conformations

Three linear peptides incorporating d ‐Phe‐2‐Abz as the turn motif are reported. Peptide 1 , a hydrophobic β‐hairpin, served as a proof of principle for the design strategy with both NMR and CD spectra strongly suggesting a β‐hairpin conformation. Peptides 2 and 3, designed as amphipathic antimicrobials, exhibited broad spectrum antimicrobial activity, with potency in the nanomolar range against Staphylococcus aureus. Both compounds possess a high degree of selectivity, proving non‐haemolytic at concentrations 500 to 800 times higher than their respective minimal inhibitory concentrations (MICs) against S. aureus. Peptide 2 induced cell membrane and cell wall disintegration in both S. aureus and Pseudomonas aeruginosa as observed by transmission electron microscopy. Peptide 2 also demonstrated moderate antifungal activity against Candida albicans with an MIC of 50 μM. Synergism was observed with sub‐MIC levels of amphotericin B (AmB), leading to nanomolar MICs against C. albicans for peptide 2 . Based on circular dichroism spectra, both peptides 2 and 3 appear to exist as a mixture of conformers with the β‐hairpin as a minor conformer in aqueous solution, and a slight increase in hairpin population in 50% trifluoroethanol, which was more pronounced for peptide 3 . NMR spectra of peptide 2 in a 1:1 CD₃CN/H₂O mixture and 30 mM deuterated sodium dodecyl sulfate showed evidence of an extended backbone conformation of the β‐strand residues. However, inter‐strand rotating frame Overhauser effects (ROE) could not be detected and a loosely defined divergent hairpin structure resulted from ROE structure calculation in CD₃CN/H₂O. The loosely defined hairpin conformation is most likely a result of the electrostatic repulsions between cationic strand residues which also probably contribute towards maintaining low haemolytic activity.     

The ensemble folding kinetics of the FBP28 WW domain revealed by an all-atom Monte Carlo simulation in a knowledge-based potential

In this work, we apply a detailed all‐atom model with a transferable knowledge‐based potential to study the folding kinetics of Formin‐Binding protein, FBP28, which is a canonical three‐stranded β‐sheet WW domain. Replica exchange Monte Carlo simulations starting from random coils find native‐like (Cα RMSD of 2.68 Å) lowest energy structure. We also study the folding kinetics of FBP28 WW domain by performing a large number of ab initio Monte Carlo folding simulations. Using these trajectories, we examine the order of formation of two β‐hairpins, the folding mechanism of each individual β‐hairpin, and transition state ensemble (TSE) of FBP28 WW domain and compare our results with experimental data and previous computational studies. To obtain detailed structural information on the folding dynamics viewed as an ensemble process, we perform a clustering analysis procedure based on graph theory. Further, a rigorous P_fold analysis is used to obtain representative samples of the TSEs showing good quantitative agreement between experimental and simulated Φ values. Our analysis shows that the turn structure between first and second β strands is a partially stable structural motif that gets formed before entering the TSE in FBP28 WW domain and there exist two major pathways for the folding of FBP28 WW domain, which differ in the order and mechanism of hairpin formation. Proteins 2011. © 2011 Wiley‐Liss, Inc.     

The role of proline‐containing peptide triads in β‐sheet formation: A kinetic study

The design of biomimetic materials through molecular self‐assembly is a growing area of modern nanotechnology. With problems of protein folding, self‐assembly, and sequence–structure relationships as essential in nanotechnology as in biology, the effect of the nucleation of β‐hairpin formation by proline on the folding process has been investigated in model studies. Previously such studies were limited to investigations of the influence of proline on the formation of turns in short peptide sequences. The effect of proline‐based triads on the folding of an 11‐kDa amyloidogenic peptide GH₆[(GA)₃GY(GA)₃GE]₈GAH₆ ( YE8 ) was investigated by selective substitution of the proline‐substituted triads at the γ‐turn sites. The folding and fibrillation of the singly proline‐substituted polypeptides, e.g., GH6? [(GA)3GY(GA)3GE]7(GA)3GY(GA)3PD? GAH6 ( 8PD ), and doubly proline‐substituted polypeptides, e.g., GH6? [(GA)3GY(GA)3GE]3(GA)3GY(GA)3PD[(GA)3GY(GA)3GE]3(GA)3GY(GA)3PD? GAH6 ( 4,8PD ), were directly monitored by circular dichroism and deep UV resonance Raman and fluorescence spectroscopies. These findings were used to identify the essential folding domains, i.e., the minimum number of β‐strands necessary for stable folding. These experimental findings may be especially useful in the design and construction of peptidic materials for a wide range of applications as well as in understanding the mechanisms of folding critical to fibril formation. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 339–350, 2015.     

β‐Amino acids containing peptides and click‐cyclized peptide as β‐turn mimics: a comparative study with ‘conventional’ lactam‐ and disulfide‐bridged hexapeptides

The increasing interest in click chemistry and its use to stabilize turn structures led us to compare the propensity for β‐turn stabilization of different analogs designed as mimics of the β‐turn structure found in tendamistat. The β‐turn conformation of linear β‐amino acid‐containing peptides and triazole‐cyclized analogs were compared to ‘conventional’ lactam‐ and disulfide‐bridged hexapeptide analogs. Their 3D structures and their propensity to fold in β‐turns in solution, and for those not structured in solution in the presence of α‐amylase, were analyzed by NMR spectroscopy and by restrained molecular dynamics with energy minimization. The linear tetrapeptide Ac‐Ser‐Trp‐Arg‐Tyr‐NH₂ and both the amide bond‐cyclized, c[Pro‐Ser‐Trp‐Arg‐Tyr‐D ‐Ala] and the disulfide‐bridged, Ac‐c[Cys‐Ser‐Trp‐Arg‐Tyr‐Cys]‐NH₂ hexapeptides adopt dominantly in solution a β‐turn conformation closely related to the one observed in tendamistat. On the contrary, the β‐amino acid‐containing peptides such as Ac‐(R)‐β³‐hSer‐(S)‐Trp‐(S)‐β³‐hArg‐(S)‐β³‐hTyr‐NH₂, and the triazole cyclic peptide, c[Lys‐Ser‐Trp‐Arg‐Tyr‐βtA]‐NH₂, both specifically designed to mimic this β‐turn, do not adopt stable structures in solution and do not show any characteristics of β‐turn conformation. However, these unstructured peptides specifically interact in the active site of α‐amylase, as shown by TrNOESY and saturation transfer difference NMR experiments performed in the presence of the enzyme, and are displaced by acarbose, a specific α‐amylase inhibitor. Thus, in contrast to amide‐cyclized or disulfide‐bridged hexapeptides, β‐amino acid‐containing peptides and click‐cyclized peptides may not be regarded as β‐turn stabilizers, but can be considered as potential β‐turn inducers. Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

Propensities of peptides containing the Asn‐Gly segment to form β‐turn and β‐hairpin structures

The propensities of peptides that contain the Asn‐Gly segment to form β‐turn and β‐hairpin structures were explored using the density functional methods and the implicit solvation model in CH₂Cl₂ and water. The populations of preferred β‐turn structures varied depending on the sequence and solvent polarity. In solution, β‐hairpin structures with βI′ turn motifs were most preferred for the heptapeptides containing the Asn‐Gly segment regardless of the sequence of the strands. These preferences in solution are consistent with the corresponding X‐ray structures. The sequence, H‐bond strengths, solvent polarity, and conformational flexibility appeared to interact to determine the preferred β‐hairpin structure of each heptapeptide, although the β‐turn segments played a role in promoting the formation of β‐hairpin structures and the β‐hairpin propensity varied. In the heptapeptides containing the Asn‐Gly segment, the β‐hairpin formation was enthalpically favored and entropically disfavored at 25°C in water. The calculated results for β‐turns and β‐hairpins containing the Asn‐Gly segment imply that these structural preferences may be useful for the design of bioactive macrocyclic peptides containing β‐hairpin mimics and the design of binding epitopes for protein–protein and protein–nucleic acid recognitions. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 653–664, 2016.     

Aromatic interactions with naphthylalanine in a β‐hairpin peptide

Stable peptides have been explored as epitope mimics for protein–protein and protein–nucleic acid interactions; however, presentation of a regular structure is critical. Aromatic interactions are ubiquitous and are competent at stabilizing a β‐hairpin fold. The greatest stabilization has been reported from pairs of tryptophan side chains. Naphthylalanine residues are often used as tryptophan replacements, but it is not clear if 1‐naphthylalanine or 2‐naphthylalanine is adequate at replicating the geometry and stability observed with tryptophan aromatic interactions. Herein, a 12‐residue peptide has been constructed with laterally disposed aromatic amino acids. A direct comparison is made between tryptophan and other bicyclic, unnatural amino acids. Significant stabilization is gained from all bicyclic amino acids; however, geometric analysis shows that only 1‐naphthylalanine adopts a similar edge to face geometry as tryptophan, whereas the 2‐naphthylalanine appears most similar to a substituted phenylalanine. Copyright © 2013 European Peptide Society and John Wiley & Sons, Ltd.     

Protein prosthesis: β‐peptides as reverse‐turn surrogates

The introduction of non‐natural modules could provide unprecedented control over folding/unfolding behavior, conformational stability, and biological function of proteins. Success requires the interrogation of candidate modules in natural contexts. Here, expressed protein ligation is used to replace a reverse turn in bovine pancreatic ribonuclease (RNase A) with a synthetic β‐dipeptide: β²‐homoalanine–β³‐homoalanine. This segment is known to adopt an unnatural reverse‐turn conformation that contains a 10‐membered ring hydrogen bond, but one with a donor–acceptor pattern opposite to that in the 10‐membered rings of natural reverse turns. The RNase A variant has intact enzymatic activity, but unfolds more quickly and has diminished conformational stability relative to native RNase A. These data indicate that hydrogen‐bonding pattern merits careful consideration in the selection of beneficial reverse‐turn surrogates.     

Optimization of a β‐sheet‐cap for long loop closure

Protein loops make up a large portion of the secondary structure in nature. But very little is known concerning loop closure dynamics and the effects of loop composition on fold stability. We have designed a small system with stable β‐sheet structures, including features that allow us to probe these questions. Using paired Trp residues that form aromatic clusters on folding, we are able to stabilize two β‐strands connected by varying loop lengths and composition (an example sequence: R W ITVTI – loop – KKIRV W E). Using NMR and CD, both fold stability and folding dynamics can be investigated for these systems. With the 16 residue loop peptide (sequence: R W ITVTI‐(GGGGKK)₂GGGG‐KKIRV W E) remaining folded (ΔG_U = 1.6 kJ/mol at 295K). To increase stability and extend the series to longer loops, we added an additional Trp/Trp pair in the loop flanking position. With this addition to the strands, the 16 residue loop (sequence: R W ITVRI W ‐(GGGGKK)₂GGGG‐ W KTIRV W E) supports a remarkably stable β‐sheet (ΔG_U = 6.3 kJ/mol at 295 K, T_m = ～55°C). Given the abundance of loops in binding motifs and between secondary structures, these constructs can be powerful tools for peptide chemists to study loop effects; with the Trp/Trp pair providing spectroscopic probes for assessing both stability and dynamics by NMR.     

Residue-specific α-helix propensities from molecular simulation

Formation of α-helices is a fundamental process in protein folding and assembly. By studying helix formation in molecular simulations of a series of alanine-based peptides, we obtain the temperature-dependent α-helix propensities of all 20 naturally occurring residues with two recent additive force fields, Amber ff03w and Amber ff99SB¹. Encouragingly, we find that the overall helix propensity of many residues is captured well by both energy functions, with Amber ff99SB¹ being more accurate. Nonetheless, there are some residues that deviate considerably from experiment, which can be attributed to two aspects of the energy function: i), variations of the charge model used to determine the atomic partial charges, with residues whose backbone charges differ most from alanine tending to have the largest error; ii), side-chain torsion potentials, as illustrated by the effect of modifications to the torsion angles of I, L, D, N. We find that constrained refitting of residue charges for charged residues in Amber ff99SB¹ significantly improves their helix propensity. The resulting parameters should more faithfully reproduce helix propensities in simulations of protein folding and disordered proteins.