Strategies and rationales for the de novo design of a helical hairpin peptide

Folding Simulations of an α-Helical Hairpin Motif αtα with Residue-Specific Force Fields

α-Helical hairpin (two-helix bundle) is a structure motif composed of two interacting helices connected by a turn or a short loop. It is an important model for protein folding studies, filling the gap between isolated α-helix and larger all-α domains. Here, we present, for the first time, successful folding simulations of an α-helical hairpin. Our RSFF1 and RSFF2 force fields give very similar predicted structures of this αtα peptide, which is in good agreement with its NMR structure. Our simulations also give site-specific stability of α-helix formation in good agreement with amide hydrogen exchange experiments. Combining the folding free energy landscapes and analyses of structures sampled in five different ranges of the fraction of native contacts (Q), a folding mechanism of αtα is proposed. The most stable sites of Q9-E15 in helix-1 and E24-A30 in helix-2 close to the loop region act as the folding initiation sites. The formation of interhelix side-chain contacts also initiates near the loop region, but some residues in the central parts of the two helices also form contacts quite early. The two termini fold at a final stage, and the loop region remains flexible during the whole folding process. This mechanism is similar to the "zipping out" pathway of β-hairpin folding.

Dynamic properties of pH-dependent structural organization of the amyloidogenic beta-protein (1-40)

The structural organization of the amyloidogenic beta-protein containing 40 amino acid residues (Abeta40) was studied by the high temperature molecular dynamics simulations in the acidic (pH approximately 3) and basic (pH approximately 8) pH regions. The obtained data suggest that the central Ala21-Gly29 segment of Abeta40 can adopt folded and partially unfolded structures. At the basic pH, this segment forms folded structures stabilized by electrostatic interactions and hydrogen bonds. At the acidic pH, it forms partially unfolded structures. Two other segments flanking to the central segment exhibit the propensity to adopt unstable interconverting alpha-helical, 3(10)-helical and turn-like structures. One of these segments is comprised of the Ala30-Val36 residues at both of the considered pHs. The second segment is comprised of the Glu11-Phe20 at the basic pH and of the Glu11-Val24 residues at the acidic pHs. The revealed pH-dependent structuration of the Abeta40 allowed us to suggest a possible scenario for initial Abeta aggregation. According to this scenario, the occurrence of the partially unfolded states of the Ala21-Gly29 segment plays main role in the Abeta oligomerization process.

Sampling of states for estimating the folding funnel entropy and energy landscape of a model alpha-helical hairpin peptide

Protein folding times are many orders of magnitude shorter than would occur if the peptide chain randomly sampled possible configurations, which implies that protein folding is a directed process. The detailed shape of protein's energy landscape determines the rate and reliability of folding to the native state, but the large number of structural degrees of freedom generates an energy landscape that is hard to visualize because of its high dimensionality. A commonly used picture is that of an energy funnel leading from high energy random coil state down to the low energy native state. As lattice computer models of protein dynamics become more realistic, the number of possible configurations becomes too large to count directly. Statistical mechanic and thermodynamic approaches allow us to count states in an approximate manner to quantify the entropy and energy of the energy landscape within a folding funnel for an alpha-helical protein. We also discuss the problems that arise in attempting to count the huge number of individual states of the random coil at the top of the funnel.

Fundamental processes of protein folding: measuring the energetic balance between helix formation and hydrophobic interactions

Theories of protein folding often consider contributions from three fundamental elements: loops, hydrophobic interactions, and secondary structures. The pathway of protein folding, the rate of folding, and the final folded structure should be predictable if the energetic contributions to folding of these fundamental factors were properly understood. alphatalpha is a helix-turn-helix peptide that was developed by de novo design to provide a model system for the study of these important elements of protein folding. Hydrogen exchange experiments were performed on selectively 15N-labeled alphatalpha and used to calculate the stability of hydrogen bonds within the peptide. The resulting pattern of hydrogen bond stability was analyzed using a version of Lifson-Roig model that was extended to include a statistical parameter for tertiary interactions. This parameter, x, represents the additional statistical weight conferred upon a helical state by a tertiary contact. The hydrogen exchange data is most closely fit by the XHC model with an x parameter of 9.25. Thus the statistical weight of a hydrophobic tertiary contact is approximately 5.8x the statistical weight for helix formation by alanine. The value for the x parameter derived from this study should provide a basis for the understanding of the relationship between hydrophobic cluster formation and secondary structure formation during the early stages of protein folding.

Removal of kinetic traps and enhanced protein folding by strategic substitution of amino acids in a model α-helical hairpin peptide

The presence of non-native kinetic traps in the free energy landscape of a protein may significantly lengthen the overall folding time so that the folding process becomes unreliable. We use a computational model alpha-helical hairpin peptide to calculate structural free energy landscapes and relate them to the kinetics of folding. We show how protein engineering through strategic changes in only a few amino acid residues along the primary sequence can greatly increase the speed and reliability of the folding process, as seen experimentally. These strategic substitutions also prevent the formation of long-lived misfolded configurations that can cause unwanted aggregations of peptides. These results support arguments that removal of kinetic traps, obligatory or nonobligatory, is crucial for fast folding.

Tuning the conformation properties of a peptide by glycosylation and phosphorylation

We have deployed the alpha-helical hairpin peptide (alpha-helix/turn/alpha-helix) and used it as a model system to explore how glycosylation and phosphorylation might affect the conformational properties of the peptide. The native conformations of the modified peptides in buffer solution have been compared with that of the wild-type peptide by nuclear magnetic resonance spectroscopy. Circular dichroism spectroscopy was used to probe the effects of an O-linked beta-GlcNAc and a phosphate group on the overall folding stability of the peptide. Finally, the rate of fibrillogenesis was used to infer the effects of these chemical modifications on the alpha-to-beta transition as well as the rate of nucleation of amyloidogenesis.

The Stabilizing Effects ofO-glycosylation on the Secondary Structural Integrity of the Designed α-loop-α motif by Molecular Dynamics Simulations

In this study, various 400 ps molecular dynamics simulations were conducted to determine the stabilizing effect of O-glycosylation on the secondary structural integrity of the design alpha-loop-alpha motif, which has the optimal loop length of 7 Gly residues (denoted as N-A16G7A16-C). In general, O-glycosylation stabilizes the structural integrity of the model peptide regardless of the length and position of glycosylation sites because it decreases the opportunity for water molecules to compete for the intramolecular hydrogen bonds. The designed peptide exhibits the highest helicity when residues 11 and 31 are replaced with Ser residues followed by O-linked with 3 galactose residues, representing the "face-to-face" glycosylation near the loop. In this case, the loop exhibits an extended conformation and several new hydrogen bonds are observed between the main chain of the loop and the galactose residues, resulting in decreasing the fluctuation and increasing the stability of the entire peptide. When the glycosylation are made close to the loop, the secondary structural integrity of the alpha-loop-alpha motif increases with the number of galactose residues. In addition, "face-to-face" glycosylation increases the structural integrity of this motif to a greater extent than "back-to-back" glycosylation. However, when the glycosylation are created away from the loop and near the N- and C-termini, no general rule is found for the stabilizing effect.

Molecular dynamics simulations to determine the optimal loop length in the helix-loop-helix motif

In this study, 100 ps molecular dynamics simulations were conducted to determine the optimal length of the Gly loop in the helix-loop-helix motif. The results revealed that the length of the loop significantly affected the stability of the two alpha-helices. In addition, the optimal loop length to maintain the highest helicity of the target peptide was found to be the one corresponding to 7 Gly residues. Longer loop length resulted in structural destabilization due to the entropic effect, while shorter loop length was not suitable to achieve intrachain packing of the two alpha-helices.

Application of the diffusion-collision model to the folding of three-helix bundle proteins

The diffusion-collision model has been successful in explaining many features of protein folding kinetics, particularly for helical proteins. In the model the folding reaction is described in terms of coupled chemical kinetic (Master) equations of coarse grained entities, called microdomains. Here, the diffusion-collision model is applied to compute the folding kinetics of four three-helix bundle proteins, all of which fold on a time scale of tens of microseconds and appear to have two-state folding. The native structure and the stability of the helical microdomains are used to determine the parameters of the model. The formulation allows computation of the overall rate and determination of the importance of kinetic intermediates. The proteins considered are the B domain of protein A (1BDC), the Engrailed Homeodomain (1ENH), the peripheral sub-unit-binding domain (1EBD C-chain) and the villin headpiece subdomain (1VII). The results for the folding time of protein A, the Engrailed Homeodomain, and 1EBD C-chain are in agreement with experiment, while 1VII is not stable in the present model. In the three proteins that are stable, two-state folding is predicted by the diffusion-collision model. This disagrees with published assertions that multistate kinetics would be obtained from the model. The contact order prediction agrees with experiment for protein A, but yields values that are a factor of 40, 30 and 15 too slow for 1ENH, 1EBD C-chain and 1VII. The effect of mutants on folding is described for protein A and it is demonstrated that significant intermediate concentrations (i.e. deviation from two-state folding) can occur if the stability of some of the helical microdomains is increased. A linear relationship between folding time and the length of the loop between helices B and C in protein A is demonstrated; this is not evident in the contact order description.

Evaluation of the energetic contribution of interhelical Coulombic interactions for coiled coil helix orientation specificity

Coiled coils are formed by two or more alpha-helices that align in a parallel or an antiparallel relative orientation. The factors that determine a preference for a given relative helix orientation are incompletely understood. The helix orientation preference for the designed coiled coil, Acid-a1-Base-a1, was measured previously. This model system therefore provides a means for the experimental determination of the energetic contribution of a variety of interactions to helix orientation specificity. The antiparallel preference for Acid-a1-Base-a1 is imparted by a single buried polar interaction. Interhelical Coulombic interactions between residues at the e and g positions have been proposed to influence helix orientation preference. In the Acid-a1-Base-a1 heterodimer, potentially attractive Coulombic interactions are expected in both orientations. To determine the energetic consequences of Coulombic interactions for helix orientation preference, we have positioned a single charged residue in each peptide such that exclusively favorable interhelical Coulombic interactions can occur only in the parallel orientation. In contrast, two potentially repulsive interactions are expected in the antiparallel orientation. Because the buried polar interaction can occur only in the antiparallel orientation, interhelical Coulombic interactions favor the parallel orientation and the potential to form a buried polar interaction favors the antiparallel orientation. We find no clear preference for an antiparallel orientation in the resulting heterodimer, Acid-Ke-Base-Eg, suggesting that interhelical Coulombic interactions and a buried polar interaction are of approximately equal importance for helix orientation specificity. Stability measurements indicate that maintenance of all favorable electrostatic interactions and/or avoidance of two potentially repulsive interactions contributes approximately 2.1 kcal/mol to helix orientation preference.

Uniquely folded mini-protein motifs

Mini-proteins containing fewer than 40 amino acids provide simple model systems for studying protein folding and stability as well as serving as scaffolds for the rational design of new functional motifs. This article reviews current progress on the design and characterization of discretely folded mini-protein motifs.

Stereochemical control of peptide folding

Stereochemically constrained amino acid residues that strongly favour specific backbone conformations may be used to nucleate and stabilize specific secondary structures in designed peptides. An overview of the use of alphaalpha-dialkyl amino acids in stabilizing helical structures in synthetic peptides is presented, with an emphasis on work carried out in the authors laboratory. Alpha-aminoisobutyric acid (Aib) and related achiral homologs facilitate stable helix formation in oligopeptides as exemplified by a large number of crystal structure determinations in the solid state. The ability to design conformationally rigid helical modules has been exploited in attempts to design structurally well characterized helix-linker helix, using potential nonhelical linking segments. Beta-hairpin design has been approached by exploiting the tendency of 'prime turns' to nucleate hairpin formation. The use of nucleating (D)Pro-Gly segments has resulted in the generation of several well characterized beta-hairpin structures, including the crystallographic observation of beta-hairpin in a synthetic apolar octapeptide. Extensions of this approach to three stranded beta-sheets and larger structures containing multiple (D)Pro-Gly segments appear readily possible.

Solution structure of alpha t alpha, a helical hairpin peptide of de novo design

alpha t alpha is a 38-residue peptide designed to adopt a helical hairpin conformation in solution (Fezoui Y, Weaver DL Osterhout JJ, 1995, Protein Sci 4:286-295). A previous study of the carboxylate form of alpha t alpha by CD and two-dimensional NMR indicated that the peptide was highly helical and that the helices associated in approximately the intended orientation (Fezoui Y, Weaver DL, Osterhout JJ, 1994, Proc Natl Acad Sci USA 91:3675-3679). Here, the solution structure of alpha t alpha as determined by two-dimensional NMR is reported. A total of 266 experimentally derived distance restraints and 20 dihedral angle restraints derived from J-couplings were used. One-hundred initial structures were generated by distance geometry and refined by dynamical simulated annealing. Twenty-three of the lowest-energy structures consistent with the experimental restraints were analyzed. The results presented here show that alpha t alpha is comprised of two associating helices connected by a turn region.

Prediction of folding stability and degradability of the de novo designed protein MB-1 in cow rumen

The authors have recently reported on the design of a protein (MB-1) enriched in methionine, threonine, lysine, and leucine. The protein is intended to be produced by rumen bacteria, in a way that would provide high producing lactating cows with limiting amino acids. In this report, MB-1 stability in the rumen is assessed, i.e., where the protein might be found after cell lysis or after being secreted by rumen bacteria. Current in vitro methods used to predict proteolytic degradability in the rumen were used for MB-1, as well as other natural proteins for comparison. MB-1 was found to be more susceptible to degradation than cytochrome c and ribonuclease A. Data indicate that MB-1 will be rapidly degraded if exposed to the rumen environment without protection. The contribution of folding stability to proteolytic stability was also examined. Rumen liquor components were selected to formulate a solution compatible with constraints of thermal denaturation studies. Denaturation curves show that the natural proteins were folded at rumen temperature. The MB-1 denaturation curves indicated that MB-1 does not unfold in a cooperative transition when heated from 20 to 70 degrees C. This suggests that MB-1 structure may be progressively modified as temperature increases, and that a continuum of conformations are available to MB-1. At 39 degrees C, a significant (50%) portion of MB-1 molecules had their tertiary structure unfolded, contributing to proteolytic degradability. Despite the unusual constraints used in MB-1 design (i.e., a maximized content in selected essential amino acids), results show that MB-1 has structural properties similar to previously reported de novo designed proteins.

A designed beta-hairpin peptide in crystals

Beta-hairpin structures have been crystallographically characterized only in very short acyclic peptides, in contrast to helices. The structure of the designed beta-hairpin, t-butoxycarbonyl-Leu-Val-Val-D-Pro-Gly-Leu-Val-Val-OMe in crystals is described. The two independent molecules of the octapeptide fold into almost ideal beta-hairpin conformations with the central D-Pro-Gly segment adopting a Type II' beta-turn conformation. The definitive characterization of a beta-hairpin has implications for de novo peptide and protein design, particularly for the development of three- and four-stranded beta-sheets.

Theoretical approaches to designing novel sequences to fit a given fold

One of the major goals of molecular biology is to understand how protein chains fold into a unique three-dimensional structure. Given this knowledge, perhaps the most exciting prospect will be the possibility of designing new proteins to perform designated tasks. The eventual pinnacle of protein engineering will be the fully automated design of a protein with novel structure and function. Achievement of this aim lies far in the future, although some early progress has been made recently.

Native-like and structurally characterized designed alpha-helical bundles

A number of coiled coils and alpha-helical bundles have recently been designed, and many have now been structurally characterized by X-ray crystallography. Others have not been as well characterized structurally but exhibit native-like properties in aqueous solution. Both areas of investigation have contributed greatly to our understanding of the nature of specificity in this class of molecules.