NMR study and molecular dynamics simulations of optimized beta-hairpin fragments of protein G

The stability and structure of several beta-hairpin peptide variants derived from the C-terminus of the B1 domain of protein G were investigated by a number of experimental and computational techniques. Our analysis shows that the structure and stability of this hairpin can be greatly affected by one or a few simple mutations. For example, removing an unfavorable charge near the N-terminus of the peptide (Glu42 to Gln or Thr) or optimization of the N-terminal charge-charge interactions (Gly41 to Lys) both stabilize the peptide, even in water. Furthermore, a simple replacement of a charged residue in the turn (Asp47 to Ala) changes the beta-turn conformation. Finally, we show that the effects of combining these single mutations are additive, suggesting that independent stabilizing interactions can be isolated and evaluated in a simple model system. Our results indicate that the structure and stability of this beta-hairpin peptide can be modulated in numerous ways and thus contributes toward a more complete understanding of this important model beta-hairpin as well as to the folding and stability of larger peptides and proteins.     

Thermodynamic and kinetic characterization of a beta-hairpin peptide in solution: an extended phase space sampling by molecular dynamics simulations in explicit water

The folding of the amyloidogenic H1 peptide MKHMAGAAAAGAVV taken from the syrian hamster prion protein is explored in explicit aqueous solution at 300 K using long time scale all-atom molecular dynamics simulations for a total simulation time of 1.1 mus. The system, initially modeled as an alpha-helix, preferentially adopts a beta-hairpin structure and several unfolding/refolding events are observed, yielding a very short average beta-hairpin folding time of approximately 200 ns. The long time scale accessed by our simulations and the reversibility of the folding allow to properly explore the configurational space of the peptide in solution. The free energy profile, as a function of the principal components (essential eigenvectors) of motion, describing the main conformational transitions, shows the characteristic features of a funneled landscape, with a downhill surface toward the beta-hairpin folded basin. However, the analysis of the peptide thermodynamic stability, reveals that the beta-hairpin in solution is rather unstable. These results are in good agreement with several experimental evidences, according to which the isolated H1 peptide adopts very rapidly in water beta-sheet structure, leading to amyloid fibril precipitates [Nguyen et al., Biochemistry 1995;34:4186-4192; Inouye et al., J Struct Biol 1998;122:247-255]. Moreover, in this article we also characterize the diffusion behavior in conformational space, investigating its relations with folding/unfolding conditions.     

The folding landscape of an alpha-lytic protease variant reveals the role of a conserved beta-hairpin in the development of kinetic stability

Most secreted bacterial proteases, including alpha-lytic protease (alphaLP), are synthesized with covalently attached pro regions necessary for their folding. The alphaLP folding landscape revealed that its pro region, a potent folding catalyst, is required to circumvent an extremely large folding free energy of activation that appears to be a consequence of its unique unfolding transition. Remarkably, the alphaLP native state is thermodynamically unstable; a large unfolding free energy barrier is solely responsible for the persistence of its native state. Although alphaLP folding is well characterized, the structural origins of its remarkable folding mechanism remain unclear. A conserved beta-hairpin in the C-terminal domain was identified as a structural element whose formation and positioning may contribute to the large folding free energy barrier. In this article, we characterize the folding of an alphaLP variant with a more favorable beta-hairpin turn conformation (alphaLP(beta-turn)). Indeed, alphaLP(beta-turn) pro region-catalyzed folding is faster than that for alphaLP. However, instead of accelerating spontaneous folding, alphaLP(beta-turn) actually unfolds more slowly than alphaLP. Our data support a model where the beta-hairpin is formed early, but its packing with a loop in the N-terminal domain happens late in the folding reaction. This tight packing at the domain interface enhances the kinetic stability of alphaLP(beta-turn), to nearly the same degree as the change between alphaLP and a faster folding homolog. However, alphaLP(beta-turn) has impaired proteolytic activity that negates the beneficial folding properties of this variant. This study demonstrates the evolutionary limitations imposed by the simultaneous optimization of folding and functional properties.     

Complex folding pathways in a simple beta-hairpin

The determination of the folding mechanisms of proteins is critical to understand the topological change that can propagate Alzheimer and Creutzfeld-Jakobs diseases, among others. The computational community has paid considerable attention to this problem; however, the associated time scale, typically on the order of milliseconds or more, represents a formidable challenge. Ab initio protein folding from long molecular dynamics simulations or ensemble dynamics is not feasible with ordinary computing facilities and new techniques must be introduced. Here we present a detailed study of the folding of a 16-residue beta-hairpin, described by a generic energy model and using the activation-relaxation technique. From a total of 90 trajectories at 300 K, three folding pathways emerge. All involve a simultaneous optimization of the complete hydrophobic and hydrogen bonding interactions. The first two pathways follow closely those observed by previous theoretical studies (folding starting at the turn or by interactions between the termini). The third pathway, never observed by previous all-atom folding, unfolding, and equilibrium simulations, can be described as a reptation move of one strand of the beta-sheet with respect to the other. This reptation move indicates that non-native interactions can play a dominant role in the folding of secondary structures. Furthermore, such a mechanism mediated by non-native hydrogen bonds is not available for study by unfolding and Gō model simulations. The exact folding path followed by a given beta-hairpin is likely to be influenced by its sequence and the solvent conditions. Taken together, these results point to a more complex folding picture than expected for a simple beta-hairpin.     

Beta-hairpin formation in aqueous solution and in the presence of trifluoroethanol: a (1)H and (13)C nuclear magnetic resonance conformational study of designed peptides

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

Role of cation-pi interactions in single chain 'all-alpha' proteins

Cation–π interactions are known to be important contributors to protein stability and ligand–protein interactions. In this study, we have analyzed the influence of cation–π interactions in single chain ‘all-alpha’ proteins. We observed 135 cation–π interactions in a data set of 75 proteins. No significant correlation was observed between the total number of amino acid residues and number of cation–π interactions. These interactions are mainly formed by long-range contacts and there is preference of Arg over Lys in these interactions. Arg–Phe interactions are predominant among the various pairs analyzed. Despite the scarcity of interactions involving Trp, the average energy for Trp–cation interactions, was quite high. This information implies that the cation–π interactions involving Trp, maybe of high relevance to the proteins. Secondary structure analysis reveals that cation–π interactions are formed preferrably between residues, in which at least one of them, is in the secondary structure of alpha-helical segments. Among the various types of folds of ‘all-alpha’ proteins considered for the analysis, proteins belonging to alpha–alpha superhelix fold have the highest number of cation–π interaction forming residues.     

Free energy landscape and folding mechanism of a beta-hairpin in explicit water: a replica exchange molecular dynamics study

The free energy landscape and the folding mechanism of the C-terminal beta-hairpin of protein G is studied by extensive replica exchange molecular dynamics simulations (40 replicas and 340 ns total simulation time), using the GROMOS96 force field and the SPC explicit water solvent. The study reveals that the system preferentially adopts a beta-hairpin structure at biologically important temperatures, and that the helix content is low at all temperatures studied. Representing the free energy landscape as a function of several types of reaction coordinates, four local minima corresponding to the folded, partially folded, molten globule, and unfolded states are identified. The findings suggest that the folding of the beta-hairpin occurs as the sequence: collapse of hydrophobic core --> formation of H-bond --> formation of the turn. Identifying the folded and molten globule states as the main conformations, the free energy landscape of the beta-hairpin is consistent with a two-state behavior with a broad transition state. The temperature dependence of the folding-unfolding transition is investigated in some detail. The enthalpy and entropy jumps at the folding transition temperature are found to be about three times lower than the experimental estimates, indicating that the folding-unfolding transition in silico is less cooperative than its in vitro counterpart.     

Cation-pi interactions in protein-protein interfaces

Arginine is an abundant residue in protein-protein interfaces. The importance of this residue relates to the versatility of its side chain in intermolecular interactions. Different classes of protein-protein interfaces were surveyed for cation-pi interactions. Approximately half of the protein complexes and one-third of the homodimers analyzed were found to contain at least one intermolecular cation-pi pair. Interactions between arginine and tyrosine were found to be the most abundant. The electrostatic interaction energy was calculated to be approximately 3 kcal/mol, on average. A distance-based search of guanidinium:aromatic interactions was also performed using the Macromolecular Structure Database (MSD). This search revealed that half of the guanidinium:aromatic pairs pack in a coplanar manner. Furthermore, it was found that the cationic group of the cation-pi pair is frequently involved in intermolecular hydrogen bonds. In this manner the arginine side chain can participate in multiple interactions, providing a mechanism for inter-protein specificity. Thus, the cation-pi interaction is established as an important contributor to protein-protein interfaces.     

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.     

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

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

Investigation of the nature of the methionine-pi interaction in beta-hairpin peptide model systems

There are frequent contacts between aromatic rings and sulfur atoms in proteins. However, it is unclear to what degree this putative interaction is stabilizing and what the nature of the interaction is. We have investigated the aryl-sulfur interaction by placing a methionine residue diagonal to an aromatic ring on the same face of a beta-hairpin, which places the methionine side chain in close proximity to the aryl side chain. The methionine (Met)-aryl interaction was compared with an equivalent hydrophobic and cation-pi interaction in the context of the beta-hairpin. The interaction between phenylalanine (Phe), tryptophan (Trp), or cyclohexylalanine (Cha) and Met stabilized the beta-hairpin by -0.3 to -0.5 kcal mole(-1), as determined by double-mutant cycles. The peptides were subjected to thermal denaturations that suggest a hydrophobic driving force for the interactions between Met and Trp or Cha. The observed interaction of Met or norleucine (Nle) with Trp or Cha are quite similar, implying a hydrophobic driving force for the Met-pi interaction. However, the thermodynamic data suggest that there may be some differences between the interaction of Met with Trp and Phe and that there may be a small thermodynamic component to the Met...Phe interaction.     

Free energy surfaces of beta-hairpin and alpha-helical peptides generated by replica exchange molecular dynamics with the AGBNP implicit solvent model

We have studied the potential of mean force of two peptides, one known to adopt a beta-hairpin and the other an alpha-helical conformation in solution. These peptides are, respectively, residues 41-56 of the C-terminus (GEWTYDDATKTFTVTE) of the B1 domain of protein G and the 13 residue C-peptide (KETAAAKFERQHM) of ribonuclease A. Extensive canonical ensemble sampling has been performed using a parallel replica exchange method. The effective potential employed in this work consists of the OPLS all-atom force field (OPLS-AA) and an analytical generalized Born (AGB) implicit solvent model including a novel nonpolar solvation free energy estimator (NP). An additional dielectric screening parameter has been incorporated into the AGBNP model. In the case of the beta-hairpin, the nonpolar solvation free energy estimator provides the necessary effective interactions for the collapse of the hydrophobic core (W43, Y45, F52, and V54), which the more commonly used surface-area-dependent nonpolar model does not provide. For both the beta-hairpin and the alpha-helix, increased dielectric screening reduces the stability of incorrectly formed salt bridges, which tend to disrupt the formation of the hairpin and helix, respectively. The fraction of beta-hairpin and alpha-helix content we obtained using the AGBNP model agrees well with experimental results. The thermodynamic stability of the beta-hairpin from protein G and the alpha-helical C-peptide from ribonuclease A as modeled with the OPLS-AA/AGBNP effective potential reflects the balance between the nonpolar effective potential terms, which drive compaction, and the polar and hydrogen bonding terms, which promote secondary structure formation.     

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

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

Signature of n→π* interactions in α-helices

The oxygen of a peptide bond has two lone pairs of electrons. One of these lone pairs is poised to interact with the electron-deficient carbon of the subsequent peptide bond in the chain. Any partial covalency that results from this n→π* interaction should induce pyramidalization of the carbon (C'(i)) toward the oxygen (O(i-1)). We searched for such pyramidalization in 14 peptides that contain both α- and β-amino acid residues and that assume a helical structure. We found that the α-amino acid residues, which adopt the main chain dihedral angles of an α-helix, display dramatic pyramidalization but the β-amino acid residues do not. Thus, we conclude that O(i-1) and C'(i) are linked by a partial covalent bond in α-helices. This finding has important ramifications for the folding and conformational stability of α-helices in isolation and in proteins.     

Determination of alpha-helix and beta-sheet stability in the solid state: a solid-state NMR investigation of poly(L-alanine)

The relative stability of alpha-helix and beta-sheet secondary structure in the solid state was investigated using poly(L-alanine) (PLA) as a model system. Protein folding and stability has been well studied in solution, but little is known about solid-state environments, such as the core of a folded protein, where peptide packing interactions are the dominant factor in determining structural stability. (13)C cross-polarization with magic angle spinning (CPMAS) NMR spectroscopy was used to determine the backbone conformation of solid powder samples of 15-kDa and 21.4-kDa PLA before and after various sample treatments. Reprecipitation from helix-inducing solvents traps the alpha-helical conformation of PLA, although the method of reprecipitation also affects the conformational distribution. Grinding converts the secondary structure of PLA to a final steady-state mixture of 55% beta-sheet and 45% alpha-helix at room temperature regardless of the initial secondary structure. Grinding PLA at liquid nitrogen temperatures leads to a similar steady-state mixture with 60% beta-sheet and 40% alpha-helix, indicating that mechanical shear force is sufficient to induce secondary structure interconversion. Cooling the sample in liquid nitrogen or subjecting it to high pressure has no effect on secondary structure. Heating the sample without grinding results in equilibration of secondary structure to 50% alpha-helix/50% beta-sheet at 100 degrees C when starting from a mostly alpha-helical state. No change was observed upon heating a beta-sheet sample, perhaps due to kinetic effects and the different heating rate used in the experiments. These results are consistent with beta-sheet approximately 260 J/mol more stable than alpha-helix in solid-state PLA.     

Turn stability in beta-hairpin peptides: Investigation of peptides containing 3:5 type I G1 bulge turns

The turn-forming ability of a series of three-residue sequences was investigated by substituting them into a well-characterized beta-hairpin peptide. The starting scaffold, bhpW, is a disulfide-cyclized 10-residue peptide that folds into a stable beta-hairpin with two antiparallel strands connected by a two-residue reverse turn. Substitution of the central two residues with the three-residue test sequences leads to less stable hairpins, as judged by thiol-disulfide equilibrium measurements. However, analysis of NMR parameters indicated that each molecule retains a significant folded population, and that the type of turn adopted by the three-residue sequence is the same in all cases. The solution structure of a selected peptide with a PDG turn contained an antiparallel beta-hairpin with a 3:5 type I + G1 bulge turn. Analysis of the energetic contributions of individual turn residues in the series of peptides indicates that substitution effects have significant context dependence, limiting the predictive power of individual amino acid propensities for turn formation. The most stable and least stable sequences were also substituted into a more stable disulfide-cyclized scaffold and a linear beta-hairpin scaffold. The relative stabilities remained the same, suggesting that experimental measurements in the bhpW context are a useful way to evaluate turn stability for use in protein design projects. Moreover, these scaffolds are capable of displaying a diverse set of turns, which can be exploited for the mimicry of protein loops or for generating libraries of reverse turns.     

The role of a beta-bulge in the folding of the beta-hairpin structure in ubiquitin

It is known that the peptide corresponding to the N-terminal beta-hairpin of ubiquitin, U(1-17), can populate the monomeric beta-hairpin conformation in aqueous solution. In this study, we show that the Gly-10 that forms the bulge of the beta-turn in this hairpin is very important to the stability of the hairpin. The deletion of this residue to desG10(1-16) unfolds the structure of the peptide in water. Even under denaturing conditions, this bulge appears to be important in maintaining the residual structure of ubiquitin, which involves tertiary interactions within the sequence 1 to 34 in the denatured state. We surmise that this residual structure functions as one of the nucleation centers in the folding process and is important in stabilizing the transition state. In accordance with this idea, deleting Gly-10 slows down the refolding and unfolding rate by about one half.     

One-dimensional self-assembly of a rational designed beta-structure peptide

Fabricating various nanostructures based on the self-assembly of diverse biological molecules is now of great interest to the field of bionanotechnology. In this study, we report a de novo designed peptide (T1) with a preferential beta-hairpin forming property that can spontaneously assemble into nanofibrils in ultrapure water. The nanofibrils assembled by T1 could grow up to tens of microns in length with a left-handed helical twist and an average height of 4.9 +/- 0.9 nm. Moreover, protofilaments and nucleus structures both with a similar height of 1.4 +/- 0.2 nm were observed during fibrilization as well as via sonication of the mature nanofibrils. A typical conformational transition from random coil to beta-structure was observed in association with the fibrilization. Molecular modeling of T1 assemblies displayed that the beta-hairpin molecules organize in a parallel fashion in which the beta-strands align in an antiparallel fashion and each adjoining beta-strand runs left-handed twist at about 2.9 degrees with respect to the one located before it along the fibrillar axis. It also revealed that the maximum thickness of the assembly intermediate, the helical tape structure, is about 1.4 nm and four tapes can further assemble into a fibril with a diameter of about 4.1 nm. Taken together the results obtained by AFM, CD, and molecular modeling, T1 fibrilization probably undergoes a hierarchy approach, in which the aromatic stacking and the electrostatic interactions between the assembled structures are most likely the two major factors directing the one-dimensional self-assembly. Based on these studies, we propose T1 can be used as a model peptide to investigate the beta-sheet based self-assembly process and could be a potential bioorganic template to develop functional materials.     

Stability of monomeric Cro variants: Isoenergetic transformation of a type I' to a type II' beta-hairpin by single amino acid replacements

The thermodynamic stabilities of three monomeric variants of the bacteriophage lambda Cro repressor that differ only in the sequence of two amino acids at the apex of an engineered beta-hairpin have been determined. The sequences of the turns are EVK-XX-EVK, where the two central residues are DG, GG, and GT, respectively. Standard-state unfolding free energies, determined from circular dichroism measurements as a function of urea concentration, range from 2.4 to 2.7 kcal/mole, while those determined from guanidine hydrochloride range from 2.8 to 3.3 kcal/mole for the three proteins. Thermal denaturation yields van't Hoff unfolding enthalpies of 36 to 40 kcal /mole at midpoint temperatures in the range of 53 to 58 degrees C. Extrapolation of the thermal denaturation free energies with heat capacities of 400 to 600 cal/mole deg gives good agreement with the parameters determined in denaturant titrations. As predicted from statistical surveys of amino acid replacements in beta-hairpins, energetic barriers to transformation from a type I' turn (DG) to a type II' turn (GT) can be quite small.