Folding mechanisms of individual beta-hairpins in a Go model of Pin1 WW domain by all-atom molecular dynamics simulations

This paper examines the folding mechanism of an individual beta-hairpin in the presence of other hairpins by using an off-lattice model of a small triple-stranded antiparallel beta-sheet protein, Pin1 WW domain. The turn zipper model and the hydrophobic collapse model originally developed for a single beta-hairpin in literature is confirmed to be useful in describing beta-hairpins in model Pin1 WW domain. We find that the mechanism for folding a specific hairpin is independent of whether it folds first or second, but the formation process are significantly dependent on temperature. More specifically, beta1-beta2 hairpin folds via the turn zipper model at a low temperature and the hydrophobic collapse model at a high temperature, while the folding of beta2-beta3 hairpin follows the turn zipper model at both temperatures. The change in folding mechanisms is interpreted by the interplay between contact stability (enthalpy) and loop lengths (entropy), the effect of which is temperature dependent.     

Alpha-helix and beta-hairpin Folding from experiment,analytical theory and molecular dynamics simulations

The alpha-helix and beta-hairpin are the minimal secondary structure elements of proteins. Identification of the factors governing the formation of these structures independently of the rest of the protein is important for understanding the determinants and rules driving the folding process to a unique native structure. It has been shown that some alpha-helices and beta-hairpins can fold autonomously into native-like structures, either in aqueous solution or in the presence of an organic co-solvent; possible mechanisms of these processes have been considered in literature. The characteristic times for folding of alpha and beta structures are estimated from experiments, simple analytical theories and more detailed computer models. Our aim is to review recent experimental and theoretical studies of folding of alpha and beta structures focusing much attention on beta-hairpins.     

Stability of cyclic beta-hairpins: asymmetric contributions from side chains of a hydrogen-bonded cross-strand residue pair

Amino acid structural propensities measured in "host-guest" model studies are often used in protein structure prediction or to choose appropriate residues in de novo protein design. While this concept has proven useful for helical structures, it is more difficult to apply successfully to beta-sheets. We have developed a cyclic beta-hairpin scaffold as a host for measurement of individual residue contributions to hairpin structural stability. Previously, we have characterized substitutions in non-backbone-hydrogen-bonded strand sites; relative stability differences measured in the cyclic host are highly predictive of changes in folding free energy for linear beta-hairpin peptides. Here, we examine the hydrogen-bonded strand positions of our host. Surprisingly, we find a large favorable contribution to stability from a valine (or isoleucine) substitution immediately preceding the C-terminal cysteine of the host peptide, but not at the cross-strand position of the host or in either strand of a folded linear beta-hairpin (trpzip peptide). Further substitutions in the peptides and NMR structural analysis indicate that the stabilizing effect of valine is general for CX(8)C cyclic hairpins and cannot be explained by particular side-chain-side-chain interactions. Instead, a localized decrease in twist of the peptide backbone on the N-terminal side of the cysteine allows the valine side chain to adopt a unique conformation that decreases the solvent accessibility of the peptide backbone. The conformation differs from the highly twisted (coiled) conformation of the trpzip hairpins and is more typical of conformations present in multistranded beta-sheets. This unexpected structural fine-tuning may explain why cyclic hairpins selected from phage-displayed libraries often have valine in the same position, preceding the C-terminal cysteine. It also emphasizes the diversity of structures accessible to beta-strands and the importance of considering not only "beta-propensity", but also hydrogen-bonding pattern and strand twist, when designing beta structures. Finally, we observe correlated, cooperative stabilization from side-chain substitutions on opposite faces of the hairpin. This suggests that cooperative folding in beta-hairpins and other small beta-structures is driven by cooperative strand-strand association.     

Site-specific relaxation kinetics of a tryptophan zipper hairpin peptide using temperature-jump IR spectroscopy and isotopic labeling

Two antiparallel beta-strands connected by a turn make beta-hairpins an ideal model system to analyze the interactions and dynamics of beta-sheets. Site-specific conformational dynamics were studied by temperature-jump IR spectroscopy and isotopic labeling in a model based on the tryptophan zipper peptide, Trpzip2, developed by Cochran et al. (Proc. Natl. Acad. Sci. U.S.A. 2001, 98, 5578). The modified Trpzip2C peptides have nearly identical equilibrium spectral behavior as Trpzip2 showing that they also form well-characterized beta-hairpin conformations in aqueous solution. Selective introduction of 13C=O groups on opposite strands lead to distinguishable cross-strand coupling of the labeled residues as monitored in the amide I' band. These frequency patterns reflect theoretical predictions, and the coupled 13C=O band loses intensity with increase in temperature and unfolding of the hairpin. Thermal relaxation kinetics were analyzed for unlabeled and cross-strand isotopically labeled variants. T-jumps of approximately 10 degrees C induce relaxation times of a few microseconds that decrease with increase of the peptide temperature. Differences in kinetic behavior for the loss of beta-strand and gain of disordered structure can be used to distinguish localized structure dynamics by comparison of nonlabeled and labeled amide I' components. Analysis of the data supports multistate dynamic and equilibrium behavior, but because of this process it is not possible to clearly define a folding and unfolding rate. Nonetheless, site-specific relaxation kinetics could be seen to be consistent with a hydrophobic collapse hypothesis for hairpin folding.     

A minimal peptide scaffold for beta-turn display: optimizing a strand position in disulfide-cyclized beta-hairpins.

Phage display of peptide libraries has become a powerful tool for the evolution of novel ligands that bind virtually any protein target. However, the rules governing conformational preferences in natural peptides are poorly understood, and consequently, structure-activity relationships in these molecules can be difficult to define. In an effort to simplify this process, we have investigated the structural stability of 10-residue, disulfide-constrained beta-hairpins and assessed their suitability as scaffolds for beta-turn display. Using disulfide formation as a probe, relative free energies of folding were measured for 19 peptides that differ at a one strand position. A tryptophan substitution promotes folding to a remarkable degree. NMR analysis confirms that the measured energies correlate well with the degree of beta-hairpin structure in the disulfide-cyclized peptides. Reexamination of a subset of the strand substitutions in peptides with different turn sequences reveals linear free energy relationships, indicating that turns and strand-strand interactions make independent, additive contributions to hairpin stability. Significantly, the tryptophan strand substitution is highly stabilizing with all turns tested, and peptides that display model turns or the less stable C'-C' ' turn of CD4 on this tryptophan "stem" are highly structured beta-hairpins in water. Thus, we have developed a small, structured beta-turn scaffold, containing only natural L-amino acids, that may be used to display peptide libraries of limited conformational diversity on phage.     

Molecular dynamics simulation of folding of a short helical peptide with many charged residues

A molecular dynamics simulation of the folding of conantokin-T (con-T), a short helical peptide with 5 helical turns of 21 amino acids with 10 charged residues, was carried out to examine folding pathways for this peptide and to predict the folding rate. In the 18 trajectories run at 300 K, 16 trajectories folded, with an averaged folding time of approximately 50 ns. Two trajectories did not fold in up to 200 ns simulation. The folded structure in folded trajectories is in good agreement with experimental structure. An analysis of the trajectories showed that, at the beginning of a few nanoseconds, helix formation started from residues 5-9 with assistance of a hydrophobic clustering involving Tyr5, Met8, and Leu9. The peptide formed a U-shape mainly due to charge-charge interactions between charged residues at the N- and C-terminus segments. In the next approximately 10 ns, several nonnative charge-charge interactions were broken and nonnative Gla10-Lys18 (this denotes a salt bridge between Gal10 and Lys18) and/or Gla10-Lys19 interactions appeared more frequently in this folding step and the peptide became a fishhook J-shape. From this structure, the peptide folded to the folded state in 7 of all 16 folded trajectories in approximately 15 ns. Alternatively, in approximately 30 ns, the con-T went to a conformation in an L-shape with 4 helical turns and a kink at the Arg13 and Gla14 segment in the other 9 trajectories. Con-T in the L-shape then required another approximately 15 ns to fold into the folded state. In addition, in overall folding times, the former 7 trajectories folded faster with the total folding times all shorter than 45 ns, while the latter 9 trajectories folded at a time longer than 45 ns, resulting in an average folding time of approximately 50 ns. Two major folding intermediates found in 2 nonfolded trajectories are stabilized by charge clusters of 5 and 6 charged residues, respectively. With inclusion of friction and solvent-solvent interactions, which were ignored in the present GB/SA solvation model, the folding time obtained above should be multiplied by a factor of 1.25-1.7 according to a previous, similar simulation study. This results in a folding time of 65-105 ns, slightly shorter than the folding time of 127 ns for an alanine-based peptide of the same length. This suggests that the energy barrier of folding for this type of peptides with many charged residues is slightly lower than alanine-based helical peptides by less than 1 kcal/mol.     

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.     

The beta-peptide hairpin in solution: conformational study of a beta-hexapeptide in methanol by NMR spectroscopy and MD simulation.

The structural and thermodynamic properties of a 6-residue beta-peptide that was designed to form a hairpin conformation have been studied by NMR spectroscopy and MD simulation in methanol solution. The predicted hairpin would be characterized by a 10-membered hydrogen-bonded turn involving residues 3 and 4, and two extended antiparallel strands. The interproton distances and backbone torsional dihedral angles derived from the NMR experiments at room temperature are in general terms compatible with the hairpin conformation. Two trajectories of system configurations from 100-ns molecular-dynamics simulations of the peptide in solution at 298 and 340 K have been analyzed. In both simulations reversible folding to the hairpin conformation is observed. Interestingly, there is a significant conformational overlap between the unfolded state of the peptide at each of the temperatures. As already observed in previous studies of peptide folding, the unfolded state is composed of a (relatively) small number of predominant conformers and in this case lacks any type of secondary-structure element. The trajectories provide an excellent ground for the interpretation of the NMR-derived data in terms of ensemble averages and distributions as opposed to single-conformation interpretations. From this perspective, a relative population of the hairpin conformation of 20% to 30% would suffice to explain the NMR-derived data. Surprisingly, however, the ensemble of structures from the simulation at 340 K reproduces more accurately the NMR-derived data than the ensemble from the simulation at 298 K, a question that needs further investigation.     

A minimal proteinlike lattice model: an alpha-helix motif

A simple protein model of a four-helix bundle motif on a face-centered cubic lattice has been studied. Total energy of a conformation includes attractive interactions between hydrophobic residues, repulsive interactions between hydrophobic and polar residues, and a potential that favors helical turns. Using replica exchange Monte Carlo simulations we have estimated a set of parameters for which the native structure is a global minimum of conformational energy. Then we have shown that all the above types of interactions are necessary to guarantee the cooperativity of folding transition and to satisfy the thermodynamic hypothesis.     

Effects of turn stability on the kinetics of refolding of a hairpin in a beta-sheet

As part of our continuing study of the effects of the turn sequence on the conformational stability as well as the mechanism of folding of a beta-sheet structure, we have undertaken a parallel investigation of the solution structure, conformational stability, and kinetics of refolding of the beta-sheet VFIVDGOTYTEV(D)PGOKILQ. The latter peptide is an analogue of the original Gellman beta-sheet VFITS(D)PGKTYTEV(D)PGOKILQ, wherein the TS(D)PGK turn sequence in the first hairpin has been replaced by VDGO. Thermodynamics studies revealed comparable conformational stability of the two peptides. However, unlike the Gellman peptide, which showed extremely rapid refolding of the first hairpin, early kinetic events associated with the refolding of the corresponding hairpin in the VDGO mutant were found to be significantly slower. A detailed study of the conformation of the modified peptide suggested that hydrophobic interactions might be contributing to its stability. Accordingly, we surmise that the early kinetic events are sensitive to whether the formation of the hairpin is nucleated at the turn or by sequestering of the hydrophobic residues across the strand, before structural rearrangements to produce the nativelike topology. Nucleation of the hairpin at the turn is expected to be intrinsically rapid for a strong turn. However, if the process must involve collapse of hydrophobic side chains, the nucleation should be slower as solvent molecules must be displaced to sequester the hydrophobic residues. These findings reflect the contribution of different forces toward nucleation of hairpins in the mechanism of folding of beta-sheets.     

Intramolecular Charge Interactions as a Tool to Control the Coiled‐Coil‐to‐Amyloid Transformation

Under the influence of a changed environment, amyloid‐forming proteins partially unfold and assemble into insoluble β‐sheet rich fibrils. Molecular‐level characterization of these assembly processes has been proven to be very challenging, and for this reason several simplified model systems have been developed over recent years. Herein, we present a series of three de novo designed model peptides that adopt different conformations and aggregate morphologies depending on concentration, pH value, and ionic strength. The design strictly follows the characteristic heptad repeat of the α‐helical coiled‐coil structural motif. In all peptides, three valine residues, known to prefer the β‐sheet conformation, have been incorporated at the solvent‐exposed b, c, and f positions to make the system prone to amyloid formation. Additionally, pH‐controllable intramolecular electrostatic repulsions between equally charged lysine (peptide A) or glutamate (peptide B) residues were introduced along one side of the helical cylinder. The conformational behavior was monitored by circular dichroism spectroscopic analysis and thioflavin T fluorescence, and the resulting aggregates were further characterized by transmission electron microscopy. Whereas uninterrupted α‐helical aggregates are found at neutral pH, Coulomb repulsions between lysine residues in peptide A destabilize the helical conformation at acidic pH values and trigger an assembly into amyloid‐like fibrils. Peptide B features a glutamate‐based switch functionality and exhibits opposite pH‐dependent folding behavior. In this case, α‐helical aggregates are found under acidic conditions, whereas amyloids are formed at neutral pH. To further validate the pH switch concept, peptide C was designed by including serine residues, thus resulting in an equal distribution of charged residues. Surprisingly, amyloid formation is observed at all pH values investigated for peptide C. The results of further investigations into the effect of different salts, however, strongly support the crucial role of intramolecular charge repulsions in the model system presented herein.     

Folding kinetics of a naturally occurring helical peptide: implication of the folding speed limit of helical proteins

The folding mechanism and dynamics of a helical protein may strongly depend on how quickly its constituent alpha-helices can fold independently. Thus, our understanding of the protein folding problem may be greatly enhanced by a systematic survey of the folding rates of individual alpha-helical segments derived from their parent proteins. As a first step, we have studied the relaxation kinetics of the central helix (L9:41-74) of the ribosomal protein L9 from the bacterium Bacillus stearothermophilus , in response to a temperature-jump ( T-jump) using infrared spectroscopy. L9:41-74 has been shown to exhibit unusually high helicity in aqueous solution due to a series of side chain-side chain interactions, most of which are electrostatic in nature, while still remaining monomeric over a wide concentration range. Thus, this peptide represents an excellent model system not only for examining how the folding rate of naturally occurring helices differs from that of the widely studied alanine-based peptides, but also for estimating the folding speed limit of (small) helical proteins. Our results show that the T-jump induced relaxation rate of L9:41-74 is significantly slower than that of alanine-based peptides. For example, at 11 degrees C its relaxation time constant is about 2 micros, roughly seven times slower than that of SPE(5), an alanine-rich peptide of similar chain length. In addition, our results show that the folding rate of a truncated version of L9:41-74 is even slower. Taken together, these results suggest that individual alpha-helical segments in proteins may fold on a time scale that is significantly slower than the folding time of alanine-based peptides. Furthermore, we argue that the relaxation rate of L9:41-74 measured between 8 and 45 degrees C provides a realistic estimate of the ultimate folding rate of (small) helical proteins over this temperature range.     

Examination of the folding of a short alanine-based helical peptide with salt bridges using molecular dynamics simulation

A molecular dynamics simulation of the folding of a short alanine-based helical peptide of 17 residues with three Glu...Lys (i, i + 4) salt bridge pairs, referred to as the AEK17 peptide, was carried out. The simulation gave an estimated simulation folding time of 2.5 ns, shorter than 12 ns for an alanine-based peptide of 16 residues with three Lys residues only, referred to as the AK16 peptide, simulated previously. After folded, the AEK17 peptide had a helical content of 77%, in excellent agreement with the experimentally determined value of 80%. An examination of the folding pathways of AEK17 indicated that the peptide proceeded via three-turn helix conformations more than the helix-turn-helix conformation in the folding pathways. An analysis of interactions indicated that the formation of hydrogen bonds between Lys residue side chains and backbone carbonyls is a major factor in the abundant conformation of the three-turn helix intermediate. The substitution of three Ala with Glu residues reduces the extent of hydrophobic interaction in alanine-based AK peptides with the result that the breaking of the interactions of Lys epsilon-NH3+(side chain)...C=O(backbone) is a major activation action for the AEK17 to achieve a complete fold, in contrast to the AK16 peptide, in which breaking non-native hydrophobic interaction is the rate-determining step.     

Cyclic modular beta-sheets

The development of peptide beta-hairpins is problematic, because folding depends on the amino acid sequence and changes to the sequence can significantly decrease folding. Robust beta-hairpins that can tolerate such changes are attractive tools for studying interactions involving protein beta-sheets and developing inhibitors of these interactions. This paper introduces a new class of peptide models of protein beta-sheets that addresses the problem of separating folding from the sequence. These model beta-sheets are macrocyclic peptides that fold in water to present a pentapeptide beta-strand along one edge; the other edge contains the tripeptide beta-strand mimic Hao [JACS 2000, 122, 7654] and two additional amino acids. The pentapeptide and Hao-containing peptide strands are connected by two delta-linked ornithine (deltaOrn) turns [JACS 2003, 125, 876]. Each deltaOrn turn contains a free alpha-amino group that permits the linking of individual modules to form divalent beta-sheets. These "cyclic modular beta-sheets" are synthesized by standard solid-phase peptide synthesis of a linear precursor followed by solution-phase cyclization. Eight cyclic modular beta-sheets 1a-1h containing sequences based on beta-amyloid and macrophage inflammatory protein 2 were synthesized and characterized by 1H NMR. Linked cyclic modular beta-sheet 2, which contains two modules of 1b, was also synthesized and characterized. 1H NMR studies show downfield alpha-proton chemical shifts, deltaOrn delta-proton magnetic anisotropy, and NOE cross-peaks that establish all compounds but 1c and 1g to be moderately or well folded into a conformation that resembles a beta-sheet. Pulsed-field gradient NMR diffusion experiments show little or no self-association at low (相似文献   

Molecular Dynamics Simulations of a β-Hairpin Fragment of Protein G by Means of Atom-Bond Electronegativity Equalization Method Fused into Molecular Mechanics (ABEEMσπ/MM)

There are some controversial opinions about the origin of folding β‐hairpin stability in aqueous solution. In this study, the structural and dynamic behavior of a 16‐residue β‐hairpin from B1 domain of protein G has been investigated at 280, 300, 350 and 450 K using molecular dynamics (MD) simulations by means of Atom‐Bond Electronegativity Equalization Method Fused into Molecular Mechanics i.e., ABEEMδπ/MM and the explicit ABEEM‐7P water solvent model. In addition, a 300 K simulation of one mutant having the aromatic residues substituted with alanines has been performed. The hydrophobic surface area, hydrophilic surface area and some structural properties have been used to measure the role of the hydrophobic interactions. It is found that the aromatic residues substituted with alanines have shown an evident destabilization of the structure and unfolding started after 1.5 ns. It is also found that the number of the main chain hydrogen bonds have different distributions through three different simulations. All above demonstrate that the hydrophobic interactions and the main chain hydrogen bonds play an important role in the stability of the folding structure of β‐hairpin in solution. Furthermore, through the structural analyses of the β‐hairpin structures from four temperature simulations and the comparison with other MD simulations of β‐hairpin peptides, the new ABEEMδπ force field can reproduce the structural data in good agreement with the experimental data.     

Folding of a miniprotein with mixed fold

Using the 28 residue betabetaalpha protein FSD-EY as a target system, we examine correction terms for the ECEPP/3 force field. We find an increased probability of formation of the native state at low temperatures resulting from a reduced propensity to form alpha helices and increased formation of beta sheets. Our analysis of the observed folding events suggests that the C-terminal helix of FSD-EY is much more stable than the N-terminal beta hairpin and forms first. The hydrophobic groups of the helix provide a template which promotes the formation of the beta hairpin that is never observed to form without the helix.     

A photocontrolled beta-hairpin peptide

Beta-hairpins constitute the smallest beta-type structures in peptides and proteins. The development of highly stable, yet monomeric beta-hairpins based on the tryptophan zipper motif was therefore a remarkable success [A. G. Cochran, N. J. Skelton, M. A. Starovasnik, Proc. Natl. Acad. Sci USA 2001, 98, 5578-5583]. We have been able to design, synthesize and characterize a hairpin based on this motif which incorporates an azobenzene-based photoswitch, that allows for time-resolved folding studies of beta-structures with unprecedented time resolution. At room temperature the trans-azo isomer exhibits a mostly disordered structure; however, light-induced isomerization to the cis-azo form leads to a predominantly extended and parallel conformation of the two peptide parts, which are linked by the novel photoswitch, [3-(3-aminomethyl)phenylazo]phenylacetic acid (AMPP). While in the original sequence the dipeptide Asn-Gly forms a type I' beta-turn which connects the two strands of the hairpin, this role is adopted by the AMPP chromophore in our photoresponsive beta-hairpin that can apparently act as a beta I'-turn mimetic. The beta-hairpin structure was determined and confirmed by NMR spectroscopy, but the folding process can be monitored by pronounced changes in the CD, IR and fluorescence spectra. Finally, incorporation of the structurally and functionally important beta-hairpin motif into proteins by chemical ligation might allow for the photocontrol of protein structures and/or functions.     

beta-Hairpins,alpha-helices,and the intermediates among the secondary structures in the energy landscape of a peptide from a distal beta-hairpin of SH3 domain

Energy landscape of a peptide, extracted from a distal beta-hairpin of src SH3 domain, in explicit water was obtained with the multicanonical molecular dynamics. A variety of beta-hairpins with various strand-strand hydrogen bonds were found in the energy landscape at 300 K. There was no energy barrier between random-coil and hairpins. Thus, the peptide conformation can easily change from the random-coil to the hairpins in the thermal fluctuations at 300 K. The landscape also included two clusters of alpha-helices, among which an energy barrier existed, and besides, these helix clusters were separated from the other conformations. Thus, the free-energy barrier exists among the helices and the other conformations. Intermediate clusters were found between the helix and the hairpin clusters. The current study showed that the isolated state of this peptide in water fluctuates among random-coil, beta-hairpin, and alpha-helix. In SH3 domain, which has a topology of mainly beta-protein, the whole-protein folding may proceed when the segment is folded in the beta-hairpin and the other parts of the protein are coupled with the beta-hairpin in an energetically or kinetically favorite way.