A comparison of three- and four-helix bundle TASP molecules.

We have designed, synthesized and characterized three- and four-helix bundle template-assembled synthetic proteins CTASPs). The TASPs were synthesized using disulphide bonds between the peptides and either the cyclotribenzylene (CTB) template, or the cavitand (BOWL) template, to form the three- and four helix bundles, respectively. The TASPs were constructed using peptides that were linked via their N-termini (peptide CGGGEELLKKXEELLKKG, where X = L, I, Nle or V), or via their C-termini (peptide GEELLKKLEELLKKGGGC). Each TASP was assayed for its structure, stability, 'native-like' characteristics and whether it was a monomer in solution. All TASPs were found to be highly helical, and highly resistant to chemical denaturation using guanidine hydrochloride (GnHCl). Analysis of the GnHCl-induced unfolding curves of the different TASPs demonstrated stability differences based on the number of helices in the bundle, the end of the helix that was attached to the template, and the identity of the core amino acid. The TASPs all had molten-globule structure, which is (generally) consistent with a degenerate sequence in the core. The four-helix bundle TASPs appeared to be monomers in solution, whereas there is some evidence that the three-helix bundle TASPs are weakly self-associating.     

Helix kinks are equally prevalent in soluble and membrane proteins

Helix kinks are a common feature of α‐helical membrane proteins, but are thought to be rare in soluble proteins. In this study we find that kinks are a feature of long α‐helices in both soluble and membrane proteins, rather than just transmembrane α‐helices. The apparent rarity of kinks in soluble proteins is due to the relative infrequency of long helices (≥20 residues) in these proteins. We compare length‐matched sets of soluble and membrane helices, and find that the frequency of kinks, the role of Proline, the patterns of other amino acid around kinks (allowing for the expected differences in amino acid distributions between the two types of protein), and the effects of hydrogen bonds are the same for the two types of helices. In both types of protein, helices that contain Proline in the second and subsequent turns are very frequently kinked. However, there are a sizeable proportion of kinked helices that do not contain a Proline in either their sequence or sequence homolog. Moreover, we observe that in soluble proteins, kinked helices have a structural preference in that they typically point into the solvent. Proteins 2014; 82:1960–1970. © 2014 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

Determination of alpha-helix N1 energies after addition of N1, N2, and N3 preferences to helix/coil theory

Surveys of protein crystal structures have revealed that amino acids show unique structural preferences for the N1, N2, and N3 positions in the first turn of the alpha-helix. We have therefore extended helix-coil theory to include statistical weights for these locations. The helix content of a peptide in this model is a function of N-cap, C-cap, N1, N2, N3, C1, and helix interior (N4 to C2) preferences. The partition function for the system is calculated using a matrix incorporating the weights of the fourth residue in a hexamer of amino acids and is implemented using a FORTRAN program. We have applied the model to calculate the N1 preferences of Gln, Val, Ile, Ala, Met, Pro, Leu, Thr, Gly, Ser, and Asn, using our previous data on helix contents of peptides Ac-XAKAAAAKAAGY-CONH2. We find that Ala has the highest preference for the N1 position. Asn is the most unfavorable, destabilizing a helix at N1 by at least 1.4 kcal mol(-1) compared to Ala. The remaining amino acids all have similar preferences, 0.5 kcal mol(-1) less than Ala. Gln, Asn, and Ser, therefore, do not stabilize the helix when at N1.     

Bihelix: Towards de novo structure prediction of an ensemble of G‐protein coupled receptor conformations

G‐Protein Coupled Receptors (GPCRs) play a critical role in cellular signal transduction pathways and are prominent therapeutic targets. Recently there has been major progress in obtaining experimental structures for a few GPCRs. Each GPCR, however, exhibits multiple conformations that play a role in their function and we have been developing methods aimed at predicting structures for all these conformations. Analysis of available structures shows that these conformations differ in relative helix tilts and rotations. The essential issue is, determining how to orient each of the seven helices about its axis since this determines how it interacts with the other six helices. Considering all possible helix rotations to ensure that no important packings are overlooked, and using rotation angle increments of 30° about the helical axis would still lead to 12⁷ or 35 million possible conformations each with optimal residue positions. We show in this paper how to accomplish this. The fundamental idea is to optimize the interactions between each pair of contacting helices while ignoring the other 5 and then to estimate the energies of all 35 million combinations using these pair‐wise interactions. This BiHelix approach dramatically reduces the effort to examine the complete set of conformations and correctly identifies the crystal packing for the experimental structures plus other near‐native packings we believe may play an important role in activation. This approach also enables a detailed structural analysis of functionally distinct conformations using helix‐helix interaction energy landscapes and should be useful for other helical transmembrane proteins as well. Proteins 2012. © 2011 Wiley Periodicals, Inc.     

Structure of a protein G helix variant suggests the importance of helix propensity and helix dipole interactions in protein design

Six helix surface positions of protein G (Gbeta1) were redesigned using a computational protein design algorithm, resulting in the five fold mutant Gbeta1m2. Gbeta1m2 is well folded with a circular dichroism spectrum nearly identical to that of Gbeta1, and a melting temperature of 91 degrees C, approximately 6 degrees C higher than that of Gbeta1. The crystal structure of Gbeta1m2 was solved to 2.0 A resolution by molecular replacement. The absence of hydrogen bond or salt bridge interactions between the designed residues in Gbeta1m2 suggests that the increased stability of Gbeta1m2 is due to increased helix propensity and more favorable helix dipole interactions.     

Solid-state NMR studies of collagen-based parchments and gelatin

Historical collagen-based parchments have been studied by solid-state NMR. In addition, new parchment (produced according to traditional methods) and gelatin from bovine skin were also studied. Wideline 1H and MAS 13C measurements were carried out directly on intact parchments. A simple approach is proposed for evaluation of the extent of parchment degradation based on the linewidth changes in the 13C CPMAS spectra relative to new parchment and gelatin. Structural (bound) water content was estimated from wideline 1H NMR lineshape and relaxation time measurements. It was found that the relative water content in parchments correlates linearly with 13C MAS linewidths. Its decrease on parchment degradation indicates that structural water molecules are of primary importance in stabilizing higher order collagen structures. Backbone and side chain dynamics of collagen in parchments were compared to those of gelatin based on the 13C dipolar-dephased experiments. Carbonyl 13C chemical shift anisotropies were measured to deduce the geometry of the collagen backbone motion. Unlike previous studies, we found that the collagen backbone motion is similar to that found in other proteins and occurs primarily via small-angle librations about internal bond directions.     

Stabilizing nonpolar/polar side-chain interactions in the alpha-helix.

A simplistic, yet often used, view of protein stability is that amino acids attract other amino acids with similar polarity, whereas nonpolar and polar side chains repel. Here we show that nonpolar/polar interactions, namely Val or Ile bonding to Lys or Arg in alpha-helices, can in fact be stabilizing. Residues spaced i, i + 4 in alpha-helices are on the same face of the helix, with potential to favorably interact and stabilize the structure. We observe that the nonpolar/polar pairs Ile-Lys, Ile-Arg, and Val-Lys occur in protein helices more often than expected when spaced i, i + 4. Partially helical peptides containing pairs of nonpolar/polar residues were synthesized. Controls with i, i + 5 spacing have the residues on opposite faces of the helix and are less helical than the test peptides with the i, i + 4 interactions. Experimental circular dichroism results were analyzed with helix-coil theory to calculate the free energy for the interactions. All three stabilize the helix with DeltaG between -0.14 and -0.32 kcal x mol(-1). The interactions are hydrophobic with contacts between Val or Ile and the alkyl groups in Arg or Lys. Side chains such as Lys and Arg can thus interact favorably with both polar and nonpolar residues.     

3(10)-helices in proteins are parahelices

The 3(10)-helix is characterized by having at least two consecutive hydrogen bonds between the main-chain carbonyl oxygen of residue i and the main-chain amide hydrogen of residue i + 3. The helical parameters--pitch, residues per turn, radius, and root mean square deviation (rmsd) from the best-fit helix--were determined by using the HELFIT program. All 3(10)-helices were classified as regular or irregular based on rmsd/(N - 1)1/2 where N is the helix length. For both there are systematic, position-specific shifts in the backbone dihedral angles. The average phi, psi shift systematically from approximately -58 degrees, approximately -32 degrees to approximately -90 degrees, approximately -4 degrees for helices 5, 6, and 7 residues long. The same general pattern is seen for helices, N = 8 and 9; however, in N = 9, the trend is repeated with residues 6, 7, and 8 approximately repeating the phi, psi of residues 2, 3, and 4. The residues per turn and radius of regular 3(10)-helices decrease with increasing length of helix, while the helix pitch and rise per residue increase. That is, regular 3(10)-helices become thinner and longer as N increases from 5 to 8. The fraction of regular 3(10)-helices decreases linearly with helix length. All longer helices, N > or = 9 are irregular. Energy minimizations show that regular helices become less stable with increasing helix length. These findings indicate that the definition of 3(10)-helices in terms of average, uniform dihedral angles is not appropriate and that it is inherently unstable for a polypeptide to form an extended, regular 3(10)-helix. The 3(10)-helices observed in proteins are better referred to parahelices.     

Predicting the helix‐helix interactions from correlated residue mutations

Helix‐helix interactions are crucial in the structure assembly, stability and function of helix‐rich proteins including many membrane proteins. In spite of remarkable progresses over the past decades, the accuracy of predicting protein structures from their amino acid sequences is still far from satisfaction. In this work, we focused on a simpler problem, the prediction of helix‐helix interactions, the results of which could facilitate practical protein structure prediction by constraining the sampling space. Specifically, we started from the noisy 2D residue contact maps derived from correlated residue mutations, and utilized ridge detection to identify the characteristic residue contact patterns for helix‐helix interactions. The ridge information as well as a few additional features were then fed into a machine learning model HHConPred to predict interactions between helix pairs. In an independent test, our method achieved an F‐measure of ～60% for predicting helix‐helix interactions. Moreover, although the model was trained mainly using soluble proteins, it could be extended to membrane proteins with at least comparable performance relatively to previous approaches that were generated purely using membrane proteins. All data and source codes are available at http://166.111.152.91/Downloads.html or https://github.com/dpxiong/HHConPred .     

On the unyielding hydrophobic core of villin headpiece

Villin headpiece (HP67) is a small, autonomously-folding domain that has become a model system for understanding the fundamental tenets governing protein folding. In this communication, we explore the role that Leu61 plays in the structure and stability of the construct. Deletion of Leu61 results in a completely unfolded protein that cannot be expressed in Escherichia coli. Omission of only the aliphatic leucine side chain (HP67 L61G) perturbed neither the backbone conformation nor the orientation of local hydrophobic side chains. As a result, a large, solvent-exposed hydrophobic pocket, a negative replica of the leucine side-chain, was created on the surface. The loss of the hydrophobic interface between leucine 61 and the hydrophobic pocket destabilized the construct by ~3.3 kcal/mol. Insertion of a single glycine residue immediately before Leu61 (HP67 L61[GL]) was also highly destabilizing and had the effect of altering the backbone conformation (α-helix to π-helix) in order to precisely preserve the wild-type position and conformation of all hydrophobic residues, including Leu61. In addition to demonstrating that the hydrophobic side-chain of Leu61 is critically important for the stability of villin headpiece, our results are consistent with the notion that the precise interactions present within the hydrophobic core, rather than the hydrogen bonds that define the secondary structure, specify a protein's fold.     

Massive conformation change in the prion protein: Using dual-basin structure-based models to find misfolding pathways

We employ all-atom structure-based models with a force field with multiple energetic basins for the C-terminal (residues 166-226) of the mammalian prion protein. One basin represents the known alpha-helical (αH) structure while the other represents the same residues in a left-handed beta-helical (LHBH) conformation. The LHBH structure has been proposed to help describe one class of in vitro grown fibrils, as well as possibly self-templating the conversion of normal cellular prion protein to the infectious form. Yet, it is unclear how the protein may make this global rearrangement. Our results demonstrate that the conformation changes are not strongly limited by large-scale geometry modification and that there may exist an overall preference for the LHBH conformation. Furthermore, our model presents novel intermediate trapping conformations with twisted LHBH structure.     

Chemical synthesis and evaluation of a backbone‐cyclized minimized 2‐helix Z‐domain

The Z‐molecule is a small, engineered IgG‐binding affinity protein derived from the immunoglobulin‐binding domain B of Staphylococcus aureus protein A. The Z‐domain consists of 58 amino acids forming a well‐defined antiparallel three‐helix structure. Two of the three helices are involved in ligand binding, whereas the third helix provides structural support to the three‐helix bundle. The small size and the stable three‐helix structure are two attractive properties comprised in the Z‐domain, but a further reduction in size of the protein is valuable for several reasons. Reduction in size facilitates synthetic production of any protein‐based molecule, which is beneficial from an economical viewpoint. In addition, a smaller protein is easier to manipulate through chemical modifications. By omitting the third stabilizing helix from the Z‐domain and joining the N‐ and C‐termini by a native peptide bond, the affinity protein obtains the advantageous properties of a smaller scaffold and in addition becomes resistant to exoproteases. We here demonstrate the synthesis and evaluation of a novel cyclic two‐helix Z‐domain. The molecule has retained affinity for its target protein, is resistant to heat treatment, and lacks both N‐ and C‐termini. Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

Insights into dimerization and four-helix bundle formation found by dissection of the dimer interface of the GrpE protein from Escherichia coli

The GrpE heat shock protein from Escherichia coli has a homodimeric structure. The dimer interface encompasses two long alpha-helices at the NH(2)-terminal end from each monomer (forming a "tail"), which lead into a small four-helix bundle from which each monomer contributes two short sequential alpha-helices in an antiparallel topological arrangement. We have created a number of different deletion mutants of GrpE that have portions of the dimer interface to investigate requirements for dimerization and to study four-helix bundle formation. Using chemical crosslinking and analytical ultracentrifugation techniques to probe for multimeric states, we find that a mutant containing only the long alpha-helical tail portion (GrpE1-88) is unable to form a dimer, most likely due to a decrease in alpha-helical content as determined by circular dichroism spectroscopy, thus one reason for a dimeric structure for the GrpE protein is to support the tail region. Mutants containing both of the short alpha-helices (GrpE1-138 and GrpE88-197) are able to form a dimer and presumably the four-helix bundle at the dimer interface. These two mutants have equilibrium constants for the monomer-dimer equilibrium that are very similar to the full-length protein suggesting that the tail region does not contribute significantly to the stability of the dimer. Interestingly, one mutant that contains just one of the short alpha-helices (GrpE1-112) exists as a tetrameric species, which presumably is forming a four-helix bundle structure. A proposed model is discussed for this mutant and its relevance for factors influencing four-helix bundle formation.     

De novo design of α,β-didehydrophenylalanine containing peptides: from models to applications

The de novo design of peptides and proteins has emerged as an approach for investigating protein structure and function. The success relies heavily on the ability to design relatively short peptides that can adopt stable secondary structures. To this end, substitution with α,β-dehydroamino acids, especially α,β-didehydrophenylalanine (ΔPhe or ΔF) has blossomed in manifold directions, providing a rich diversity of well-defined structural motifs. Introduction of α,β-didehydrophenylalanine induces β-bends in small and 3(10)-helices in longer peptide sequences. Most favorable conformation of ΔF residues are (φ,ψ) ～(60°, 30°), (-60°, -30°), (-60°, 150°), and (60°, -150°). These features have been exploited in designing helix-turn-helix, helical bundle arrangements, and glycine zipper type super secondary structural motifs. The unusual capability of α,β-didehydrophenylalanine ring to form a variety of multicentered interactions (N-H…O, C-H…O, C-H…π, and N-H…π) suggests its possible exploitation for future de novo design of supramolecular structures. This work has now been extended to the de novo design of peptides with antibiotic, antifibrillization activity, etc. More recently, self-assembling properties of small dehydropeptides have been explored. This review focuses primarily on the structural and functional behavior of α,β-didehydrophenylalanine containing peptides.