Helix formation and the unfolded state of a 52-residue helical protein

A growing class of proteins in biological processes has been found to be unfolded on isolation under normal solution conditions. We have used NMR spectroscopy to characterize the structural and dynamic properties of the unfolded and partially folded states of a 52-residue alanine-rich protein (Ala-14) at temperatures from -5 degrees C to 40 degrees C. At 40 degrees C, alanine residues in Ala-14 adopt phi and psi angles, consistent with a significant ensemble population of polyproline II conformation. Analysis of relaxation rates in the protein reveals that a series of residues, Gln 35-Ala 36-Ala 37-Lys 38-Asp 39-Asp 40-Ala 41-Ala 42, displays slow motional dynamics at both -5 degrees C and 40 degrees C. Temperature-dependent chemical shift changes indicate that this region is the site of helix initiation. The remaining N-terminal residues become increasingly dynamic as they extend from the nucleation site. The C terminus remains dynamic and changes less with temperature, indicating it is relatively unstructured. Ala-14 provides a high-resolution portrait of the unfolded state and the process of helix nucleation and propagation in the absence of tertiary contacts, information that bears on early events in protein folding.     

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.     

Helix propensities of the amino acids measured in alanine-based peptides without helix-stabilizing side-chain interactions.

Helix propensities of the amino acids have been measured in alanine-based peptides in the absence of helix-stabilizing side-chain interactions. Fifty-eight peptides have been studied. A modified form of the Lifson-Roig theory for the helix-coil transition, which includes helix capping (Doig AJ, Chakrabartty A, Klingler TM, Baldwin RL, 1994, Biochemistry 33:3396-3403), was used to analyze the results. Substitutions were made at various positions of homologous helical peptides. Helix-capping interactions were found to contribute to helix stability, even when the substitution site was not at the end of the peptide. Analysis of our data with the original Lifson-Roig theory, which neglects capping effects, does not produce as good a fit to the experimental data as does analysis with the modified Lifson-Roig theory. At 0 degrees C, Ala is a strong helix former, Leu and Arg are helix-indifferent, and all other amino acids are helix breakers of varying severity. Because Ala has a small side chain that cannot interact significantly with other side chains, helix formation by Ala is stabilized predominantly by the backbone ("peptide H-bonds"). The implication for protein folding is that formation of peptide H-bonds can largely offset the unfavorable entropy change caused by fixing the peptide backbone. The helix propensities of most amino acids oppose folding; consequently, the majority of isolated helices derived from proteins are unstable, unless specific side-chain interactions stabilize them.     

α-Helix-forming propensities in peptides and proteins

Much effort has been invested in seeking to understand the thermodynamic basis of helix stability in both peptides and proteins. Recently, several groups have measured the helix-forming propensities of individual residues (Lyu, P. C., Liff, M. I., Marky, L. A., Kallenbach, N. R. Science 250:669–673, 1990; O'Neil, K. T., DeGrado, W. F. Science 250:646–651, 1990; Padmanabhan, S., Marqusee, S., Ridgeway, T., Laue, T. M., Baldwin, R. L. Nature (London) 344:268–270, 1990). Using Monte Carlo computer simulations, we tested the hypothesis that these differences in measured helix-forming propensity are due primarily to loss of side chain conformational entropy upon helix formation (Creamer, T. P., Rose, G. D. Proc. Natl. Acad. Sci. U.S.A. 89:5937–5941, 1992). Our previous study employed a rigid helix backbone, which is here generalized to a completely flexible helix model in order to ensure that earlier results were not a methodological artifact. Using this flexible model, side chain rotamer distributions and entropy losses are calculated and shown to agree with those obtained earlier. We note that the side chain conformational entropy calculated for Trp in our previous study was in error; a corrected value is presented. Extending earlier work, calculated entropy losses are found to correlate strongly with recent helix propensity scales derived from substitutions made within protein helices (Horovitz, A., Matthews, J. M., Fersht, A. R. J. Mol. Biol. 227:560–568, 1992; Blaber, M., Zhang, X.-J., Matthews, B. M. Science 260:1637–1640, 1993). In contrast, little correlation is found between these helix propensity scales and the accessible surface area buried upon formation of a model polyalanyl α-helix. Taken in sum, our results indicate that loss of side chain entropy is a major determinant of the helix-forming tendency of residues in both peptide and protein helices. © 1994 Wiley-Liss, Inc.     

Effects of alanine substitutions in alpha-helices of sperm whale myoglobin on protein stability.

The peptide backbones in folded native proteins contain distinctive secondary structures, alpha-helices, beta-sheets, and turns, with significant frequency. One question that arises in folding is how the stability of this secondary structure relates to that of the protein as a whole. To address this question, we substituted the alpha-helix-stabilizing alanine side chain at 16 selected sites in the sequence of sperm whale myoglobin, 12 at helical sites on the surface of the protein, and 4 at obviously internal sites. Substitution of alanine for bulky side chains at internal sites destabilizes the protein, as expected if packing interactions are disrupted. Alanine substitutions do not uniformly stabilize the protein, either in capping positions near the ends of helices or at mid-helical sites near the surface of myoglobin. When corrected for the extent of exposure of each side chain replaced by alanine at a mid-helix position, alanine replacement still has no clear effect in stabilizing the native structure. Thus linkage between the stabilization of secondary structure and tertiary structure in myoglobin cannot be demonstrated, probably because of the relatively small free energy differences between side chains in stabilizing isolated helix. By contrast, about 80% of the variance in free energy observed can be accounted for by the loss in buried surface area of the native residue substituted by alanine. The differential free energy of helix stabilization does not account for any additional variation.     

Structural insights for designed alanine-rich helices: comparing NMR helicity measures and conformational ensembles from molecular dynamics simulation

The temperature dependence of helical propensities for the peptides Ac-ZGG-(KAAAA)(3)X-NH(2) (Z = Y or G, X = A, K, and D-Arg) were studied both experimentally and by MD simulations. Good agreement is observed in both the absolute helical propensities as well as relative helical content along the sequence; the global minimum on the calculated free energy landscape corresponds to a single alpha-helical conformation running from K4 to A18 with some terminal fraying, particularly at the C-terminus. Energy component analysis shows that the single helix state has favorable intramolecular electrostatic energy due to hydrogen bonds, and that less-favorable two-helix globular states have favorable solvation energy. The central lysine residues do not appear to increase helicity; however, both experimental and simulation studies show increasing helicity in the series X = Ala --> Lys --> D-Arg. This C-capping preference was also experimentally confirmed in Ac-(KAAAA)(3)X-GY-NH(2) and (KAAAA)(3)X-GY-NH(2) sequences. The roles of the C-capping groups, and of lysines throughout the sequence, in the MD-derived ensembles are analyzed in detail.     

Helix propensities of short peptides: molecular dynamics versus bioinformatics

Knowledge-based potential functions for protein structure prediction assume that the frequency of occurrence of a given structure or a contact in the protein database is a measure of its free energy. Here, we put this assumption to test by comparing the results obtained from sequence-structure cluster analysis with those obtained from long all-atom molecular dynamics simulations. Sixty-four eight-residue peptide sequences with varying degrees of similarity to the canonical sequence pattern for amphipathic helix were drawn from known protein structures, regardless of whether they were helical in the protein. Each was simulated using AMBER6.0 for at least 10 ns using explicit waters. The total simulation time was 1176 ns. The resulting trajectories were tested for reproducibility, and the helical content was measured. Natural peptides whose sequences matched the amphipathic helix motif with greater than 50% confidence were significantly more likely to form helix during the course of the simulation than peptides with lower confidence scores. The sequence pattern derived from the simulation data closely resembles the motif pattern derived from the database cluster analysis. The difficulties encountered in sampling conformational space and sequence space simultaneously are discussed.     

Retrospective analysis of a secondary structure prediction: the catalytic domain of matrix metalloproteinases.

Secondary structure prediction of the catalytic domain of matrix metalloproteinases is evaluated in the light of recently published experimentally determined structures. The prediction was made by combining conformational propensity, surface probability, and residue conservation calculated for an alignment of 19 sequences. The position of each observed secondary structure element was correctly predicted with a high degree of accuracy, with a single beta-strand falsely predicted. The domain fold was also anticipated from the prediction by analogy with the structural elements found in the distantly related metalloproteinases thermolysin, astacin, and adamalysin.     

Physical principles of protein structure and protein folding

Physical principles determining the protein structure and protein folding are reviewed: (i) the molecular theory of protein secondary structure and the method of its prediction based on this theory; (ii) the existence of a limited set of thermodynamically favourable folding patterns of α- and β-regions in a compact globule which does not depend on the details of the amino acid sequence; (iii) the moderns approaches to the prediction of the folding patterns of α- and β-regions in concrete proteins; (iv) experimental approaches to the mechanism of protein folding. The review reflects theoretical and experimental works of the author and his collaborators as well as those of other groups.     

A new perspective on analysis of helix-helix packing preferences in globular proteins

For many years, statistical analysis of protein databanks has led to the belief that the steric compatibility of helix interfaces may be the source of observed preferences for particular angles between neighboring helices. Several elegant models describing how side chains on helices can interdigitate without steric clashes were able to account quite reasonably for the observed distributions. However, it was later recognized that the 'bare' measured angle distribution should be corrected to avoid statistical bias.12 Disappointingly, the rescaled distributions dramatically lost their similarity with theoretical predictions, casting doubts on the validity of the geometrical assumptions and models. In this article, we elucidate a few points concerning the proper choice of a random reference distribution. In particular we demonstrate the need for corrections induced by unavoidable uncertainties in determining whether two helices are in face-to-face contact or not and their relative orientations. By using this new rescaling, we show that 'true' packing angle preferences are well described by regular packing models, thus proving that preferential angles between contacting helices do exist.     

Discovery of a significant,nontopological preference for antiparallel alignment of helices with parallel regions in sheets

To help elucidate the role of secondary structure packing preferences in protein folding, here we present an analysis of the packing geometry observed between alpha-helices and between alpha-helices and beta-sheets in 1316 diverse, nonredundant protein structures. Finite-length vectors were fit to the alpha-carbon atoms in each of the helices and strands, and the packing angle between the vectors, Omega, was determined at the closest point of approach within each helix-helix or helix-sheet pair. Helix-sheet interactions were found in 391 of the proteins, and the distributions of Omega values were calculated for all the helix-sheet and helix-helix interactions. The packing angle preferences for helix-helix interactions are similar to those previously observed. However, analysis of helix-strand packing preferences uncovered a remarkable tendency for helices to align antiparallel to parallel regions of beta-sheets, independent of the topological constraints or prevalence of beta-alpha-beta motifs in the proteins. This packing angle preference is significantly diminished in helix interactions involving mixed and antiparallel beta-sheets, suggesting a role for helix-sheet dipole alignment in guiding supersecondary structure formation in protein folding. This knowledge of preferred packing angles can be used to guide the engineering of stable protein modules.     

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.     

Conformational analysis of peptides corresponding to all the secondary structure elements of protein L B1 domain: secondary structure propensities are not conserved in proteins with the same fold.

The solution conformation of three peptides corresponding to the two beta-hairpins and the alpha-helix of the protein L B1 domain have been analyzed by circular dichroism (CD) and nuclear magnetic resonance spectroscopy (NMR). In aqueous solution, the three peptides show low populations of native and non-native locally folded structures, but no well-defined hairpin or helix structures are formed. In 30% aqueous trifluoroethanol (TFE), the peptide corresponding to the alpha-helix adopts a high populated helical conformation three residues longer than in the protein. The hairpin peptides aggregate in TFE, and no significant conformational change occurs in the NMR observable fraction of molecules. These results indicate that the helical peptide has a significant intrinsic tendency to adopt its native structure and that the hairpin sequences seem to be selected as non-helical. This suggests that these sequences favor the structure finally attained in the protein, but the contribution of the local interactions alone is not enough to drive the formation of a detectable population of native secondary structures. This pattern of secondary structure tendencies is different to those observed in two structurally related proteins: ubiquitin and the protein G B1 domain. The only common feature is a certain propensity of the helical segments to form the native structure. These results indicate that for a protein to fold, there is no need for large native-like secondary structure propensities, although a minimum tendency to avoid non-native structures and to favor native ones could be required.     

Fractals of Protein Backbone Structures

IIntroductionPYOt6lflgy8fTSCtslStfuCtLlf6Withgedffi6tTICMlf－Slffill8Yltyh8sberflStUdi6dfOTiT18flyy6grs,Includingproteinpolypeptldechainfractaldimensions（FD）andtheFDofproteinsurfacest。c－turell－7）．Inl％2,Havlinand匝n－Avraham〔8〕studiedthefractaldlmenslonalltyofcommonplymerchainsbasedonaself、a、idingwalks（SAW）medel,andpresentedamethedforcorn－pUt;。prpOfthepl、erCha;nStrUCtUr。,qu。t;一（1）～（2）．Allenet。l．（2）一putedp］teinFDbyandassumdfracta］inedelofpr…     

Prediction of the three-dimensional structure of human growth hormone

In recent years, the protein-folding problem has attracted the attention of molecular biologists. Efforts have focused on developing heuristic and energy-based algorithms to predict the three-dimensional structure of a protein from its amino acid sequence. We have applied a series of heuristic algorithms to the sequence of human growth hormone. A family of five structures which are generically right-handed fourfold alpha-helical bundles are found from an investigation of approximately 10(8) structures. A plausible receptor binding site is suggested. Independent crystallographic analysis confirms some aspects of these predictions. These methods only deal with the "core" structure, and conformations of many residues are not defined. Further work is required to identify a unique set of coordinates and to clarify the topological alternative available to alpha-helical proteins.     

Steric restrictions in protein folding: an alpha-helix cannot be followed by a contiguous beta-strand

Using only hard-sphere repulsion, we investigated short polyalanyl chains for the presence of sterically imposed conformational constraints beyond the dipeptide level. We found that a central residue in a helical peptide cannot adopt dihedral angles from strand regions without encountering a steric collision. Consequently, an alpha-helical segment followed by a beta-strand segment must be connected by an intervening linker. This restriction was validated both by simulations and by seeking violations within proteins of known structure. In fact, no violations were found within an extensive database of high-resolution X-ray structures. Nature's exclusion of alpha-beta hybrid segments, fashioned from an alpha-helix adjoined to a beta-strand, is built into proteins at the covalent level. This straightforward conformational constraint has far-reaching consequences in organizing unfolded proteins and limiting the number of possible protein domains.