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.     

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.     

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.     

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.     

Terminal sidechain packing of a designed β‐hairpin influences conformation and stability

While end capping in α‐helices is well understood, the concept of capping a β‐hairpin is a relatively recent development; to date, favorable Coulombic interactions are the only example of sidechains at the termini influencing the overall stability of a β‐hairpin. While cross‐strand hydrophobic residues generally provide hairpin stabilization, particular when flanking the turn region, those remote from this location appear to provide little stabilization. While probing for an optimal residue at a hydrogen bond position near the terminus of a designed β‐hairpin a conservative, hydrophobic, V → I mutation was observed to not only result in a significant change in fold population but also effected major changes in the structuring shifts at numerous sites in the peptide. Mutational studies reveal that there is an interaction between the sidechain at this H‐bonded site and the sidechain at the C‐terminal non‐H‐bonded site of the hairpin. This interaction, which appears to be hydrophobic in character, requires a highly twisted hairpin structure. Modifications at the C‐terminal site, for example an E → A mutation (ΔΔG_U = 6 kJ/mol), have profound affects on fold structure and stability. The data suggests that this may be a case of hairpin end capping by the formation of a hydrophobic cluster. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 557–564, 2009. 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     

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     

β‐Hairpin stabilization in a 28‐residue peptide derived from the β‐subunit sequence of human chorionic gonadotropin hormone

The β‐subunit of the human chorionic gonadotropin (hCG) hormone, which is believed to be related to certain types of cancer, contains three hairpin‐like fragments. To investigate the role of β‐hairpin formation in the early stages of the hCGβ folding, a 28‐residue peptide with the sequence RDVRFESIRLPGSPRGVNPVVSYAVALS, corresponding to the H3‐β hairpin fragment (residues 60–87) of the hCGβ subunit, was studied under various conditions using three optical spectroscopic methods: Fourier transform ir spectroscopy, electronic CD, and vibrational CD. Environmental conditions are critical factors for formation of secondary structure in this peptide. TFE : H₂O mixed solvents induced helical formation. Formation of β‐structure in this peptide, which may be related to the native β‐hairpin formation in the intact hormone, was found to be induced only under conditions such as high concentration, high temperature, and the presence of nonmicellar sodium dodecyl sulfate concentrations. These findings support a protein folding mechanism for the hCGβ subunit in which an initial hydrophobic collapse, which increases intermolecular interactions in hCGβ, is needed to induce the H3‐β hairpin formation. © 1999 John Wiley & Sons, Inc. Biopoly 50: 413–423, 1999     

β‐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.     

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.     

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.     

Computationally designed β‐turn foldamers of γ‐peptides based on 2‐(aminomethyl)cyclohexanecarboxylic acid

The γ‐peptide β‐turn structures have been designed computationally by the combination of chirospecific γ 2 , 3 ‐residues of 2‐(aminomethyl)cyclohexanecarboxylic acid (γAmc₆) with a cyclohexyl constraint on the C^α?C^β bond using density functional methods in water. The chirospecific γAmc₆ dipeptide with the (2S,3S)‐(2R,3R) configurations forms a stable turn structure in water, resembling a type II′ turn of α‐peptides, which can be used as a β‐turn motif in β‐hairpins of Ala‐based α‐peptides. The γAmc₆ dipeptide with homochiral (2S,3S)‐(2S,3S) configurations but different cyclohexyl puckerings shows the capability to be incorporated into one of two β‐turn motifs of gramicidin S. The overall structure of this gramicidin S analogue is quite similar to the native gramicidin S with the same patterns and geometries of hydrogen bonds. Our calculated results and the recently observed results may imply the wider applicability of chirospecific γ‐peptides with a cyclohexyl constraint on the backbone to form various peptide foldamers. © 2012 Wiley Periodicals, Inc. Biopolymers 97:1018–1025, 2012.     

β‐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.     

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.     

An amino acid code for β‐sheet packing structure

To understand the relationship between protein sequence and structure, this work extends the knob‐socket model in an investigation of β‐sheet packing. Over a comprehensive set of β‐sheet folds, the contacts between residues were used to identify packing cliques: sets of residues that all contact each other. These packing cliques were then classified based on size and contact order. From this analysis, the two types of four‐residue packing cliques necessary to describe β‐sheet packing were characterized. Both occur between two adjacent hydrogen bonded β‐strands. First, defining the secondary structure packing within β‐sheets, the combined socket or XY:HG pocket consists of four residues i, i+2 on one strand and j, j+2 on the other. Second, characterizing the tertiary packing between β‐sheets, the knob‐socket XY:H+B consists of a three‐residue XY:H socket (i, i+2 on one strand and j on the other) packed against a knob B residue (residue k distant in sequence). Depending on the packing depth of the knob B residue, two types of knob‐sockets are found: side‐chain and main‐chain sockets. The amino acid composition of the pockets and knob‐sockets reveal the sequence specificity of β‐sheet packing. For β‐sheet formation, the XY:HG pocket clearly shows sequence specificity of amino acids. For tertiary packing, the XY:H+B side‐chain and main‐chain sockets exhibit distinct amino acid preferences at each position. These relationships define an amino acid code for β‐sheet structure and provide an intuitive topological mapping of β‐sheet packing. Proteins 2014; 82:2128–2140. © 2014 Wiley Periodicals, Inc.     

Dissecting the behaviour of β‐amyloid peptide variants during oligomerization and fibrillation

The oligomerization and fibrillation of β‐amyloid (Aβ) peptides are important events in the pathogenesis of Alzheimer's disease. However, the motifs within the Aβ sequence that contribute to oligomerization and fibrillation and the complex interplay among these short motifs are unclear. In this study, the oligomerization and fibrillation abilities of the Aβ variants Aβ1–28, Aβ1–36, Aβ11–42, Aβ17–42, Aβ1–40 and Aβ1–42 were examined by thioflavin T fluorescence, western blotting and transmission electron microscopy. Compared with two C‐terminal‐truncated peptides (i.e. Aβ1–28 and Aβ1–36), Aβ11–42, Aβ17–42 and Aβ1–42 had stronger abilities to form oligomers. This indicated that amino acids 37–42 strengthen the β‐hairpin structure of Aβ. Both Aβ1–42 and Aβ1–40 could form fibres, but Aβ17–42 formed irregular fibres, suggesting that amino acids 1–17 were essential for Aβ fibre formation. Aβ1–28 and Aβ1–36 exhibited weak oligomerization and fibrillation, implying that they formed an unstable β‐hairpin structure owing to the incomplete C‐terminal region. Intermediate peptides were likely to form a stable structure, consistent with previous results. This work explains the roles and interplay among motifs within Aβ during oligomerization and fibrillation. Copyright © 2017 European Peptide Society and John Wiley & Sons, Ltd.     

VCD spectroscopic properties of the β‐hairpin forming miniprotein CLN025 in various solvents

Electronic and vibrational circular dichroism are often used to determine the secondary structure of proteins, because each secondary structure has a unique spectrum. Little is known about the vibrational circular dichroic spectroscopic features of the β‐hairpin. In this study, the VCD spectral features of a decapeptide, YYDPETGTWY (CLN025), which forms a stable β‐hairpin that is stabilized by intramolecular weakly polar interactions and hydrogen bonds were determined. Molecular dynamics simulations and ECD spectropolarimetry were used to confirm that CLN025 adopts a β‐hairpin in water, TFE, MeOH, and DMSO and to examine differences in the secondary structure, hydrogen bonds, and weakly polar interactions. CLN025 was synthesized by microwave‐assisted solid phase peptide synthesis with N^α‐Fmoc protected amino acids. The VCD spectra displayed a (?,+,?) pattern with bands at 1640 to 1656 cm^?1, 1667 to 1687 cm^?1, and 1679 to 1686 cm^?1 formed by the overlap of a lower frequency negative couplet and a higher frequency positive couplet. A maximum IR absorbance was observed at 1647 to 1663 cm^?1 with component bands at 1630 cm^?1, 1646 cm^?1, 1658 cm^?1, and 1675 to 1680 cm^?1 that are indicative of the β‐sheet, random meander, either random meander or loop and turn, respectively. These results are similar to the results of others, who examined the VCD spectra of β‐hairpins formed by ^DPro‐Xxx turns and indicated that observed pattern is typical of β‐hairpins. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 442–450, 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     

HLA‐DQ7 β1 and β2 derived peptides as immunomodulators

Modulation of protein–protein interactions involved in the immune system by using small molecular mimics of the contact interfaces may lead to the blockage of the autoimmune response and the development of drugs for immunotherapy. The nonpolymorphic β‐regions, exposed to the microenvironment, of the modeled HLA‐DQ7, which is genetically linked to autoimmune diseases, were determined. Peptides 132–141 and 58–67, located at the β₁ and β₂ domains of HLA‐DQ7, respectively, were tested for their involvement in the interactions with CD4⁺ T lymphocytes. Linear, cyclic, and dimeric analogs that mimic the exposed surfaces of HLA‐DQ7 were designed and synthesized. Their immunosuppressory activities, found in the secondary, humoral immune response to sheep erythrocytes (SRBC) in mice in vitro, ranged from 11% to 53%. The significance of the total charge of the peptides, the pattern of the hydrogen bonding, and the presence of secondary structure were investigated in relation to the immunomodulatory effect of the peptides. Two dimeric analogs of the HLA‐DQ7 58–67 fragment, consisting of the two monomers covalently linked by a polyethylene glycol (PEG) spacer, able to mimic the superdimers, were also synthesized and studied. As the 58–67 segment is located at the β₁ region of HLA‐DQ7, close to the major histocompatibility complex (MHC) groove, one may assume that the 58–67 peptide could accommodate the association between T‐cell receptor (TCR) and human leukocyte antigen (HLA) by activating a co‐stimulatory molecule of the TCR/HLA interaction. This hypothesis is supported by the confocal laser image of the fluorescein‐labeled 58–67 peptide and by the fact that it is an immunostimulator at low concentration. Copyright © 2009 European Peptide Society and John Wiley & Sons, Ltd.     

Folding thermodynamics of β‐hairpins studied by replica‐exchange molecular dynamics simulations

We study the differences in folding stability of β‐hairpin peptides, including GB1 hairpin and a point mutant GB1 K10G, as well as tryptophan zippers (TrpZips): TrpZip1, TrpZip2, TrpZip3‐1, and TrpZip4. By performing replica‐exchange molecular dynamics simulations with Amber03* force field (a modified version of Amber ff03) in explicit solvent, we observe ab initio folding of all the peptides except TrpZip3‐1, which is experimentally known to be the least stable among the peptides studied here. By calculating the free energies of unfolding of the peptides at room temperature and folding midpoint temperatures for thermal unfolding of peptides, we find that TrpZip4 and GB1 K10G peptides are the most stable β‐hairpins followed by TrpZip1, GB1, and TrpZip2 in the given order. Hence, the proposed K10G mutation of GB1 peptide results in enhanced stability compared to wild‐type GB1. An important goal of our study is to test whether simulations with Amber 03* model can reproduce experimentally predicted folding stability differences between these peptides. While the stabilities of GB1 and TrpZip1 yield close agreement with experiment, TrpZip2 is found to be less stable than predicted by experiment. However, as heterogenous folding of TrpZip2 may yield divergent thermodynamic parameters by different spectroscopic methods, mismatching of results with previous experimental values are not conclusive of model shortcomings. For most of the cases, molecular simulations with Amber03* can successfully reproduce experimentally known differences between the mutated peptides, further highlighting the predictive capabilities of current state‐of‐the‐art all‐atom protein force fields. Proteins 2015; 83:1307–1315. © 2015 Wiley Periodicals, Inc.     

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.