Use of Helical Wheels to Represent the Structures of Proteins and to Identify Segments with Helical Potential

The three-dimensional structures of alpha-helices can be represented by two-dimensional projections which we call helical wheels. Initially, the wheels were employed as graphical restatements of the known structures determined by Kendrew, Perutz, Watson, and their colleagues at the University of Cambridge and by Phillips and his coworkers at The Royal Institution. The characteristics of the helices, discussed by Perutz et al. (1965), and Blake et al. (1965), can be readily visualized by examination of these wheels. For example, the projections for most helical segments of myoglobin, hemoglobin, and lysozyme have distinctive hydrophobic arcs. Moreover, the hydrophobic residues tend to be clustered in the n +/- 3, n, n +/- 4 positions of adjacent helical turns. Such hydrophobic arcs are not observed when the sequences of nonhelical segments are plotted on the wheels. Since the features of these projections are also distinctive, however, the wheels can be used to divide sequences into segments with either helical or nonhelical potential. The sequences of insulin, cytochrome c, ribonuclease A, chymotrypsinogen A, tobacco mosaic virus protein, and human growth hormone were chosen for application of the wheels for this purpose.     

Comparison of helix stability in wild-type and mutant LamB signal sequences

Previous studies of isolated peptides corresponding to the wild-type signal sequence of the LamB protein of Escherichia coli and to several export-impaired mutants demonstrated that a high tendency to adopt an alpha-helical conformation in low dielectric environments was a property of functional sequences. We have now used nuclear magnetic resonance to establish further characteristics of the helical conformation of these signal peptides in a solvent mixture (50% trifluoroethanol, by volume, in water) which mimics the conformational distribution of these peptides in lipid vesicles. The interactions of signal sequences in vivo may depend on the location of the helix in the sequence, on the length of the helical segment, and on the stability of the helix. We find that the hydrophobic core has the most persistent helix conformation and that the stability of this helix correlates with in vivo function of different mutants of the LamB signal sequence. In the family of signal peptides studied here, the length of the helix required for function appears to be less rigidly restricted since a signal peptide from a functional pseudorevertant with 4 residues deleted from the hydrophobic core takes up helix as stably as wild type but incorporates fewer residues in the helix.     

Recent advances in RNA-protein interaction studies

Conclusions The RNA binding sites for several small proteins have been characterised. These sites include double helical regions with hairpins, bulged bases and internal loops. As seen in Flock House virus structure, some proteins may recognise phosphate backbone of the canonical A-form helix not in a sequence-specific manner. If sequence-specific base contacts are to be made, then the A-helic major groove must be widened. This can be accomplished by introducing bulges, internal loops and hairpin loops into double helical regions. In these cases proteins may recognise both distorted backbone conformations and read out base sequences in a widened major groove. Crystallographic studies on complexes of aminoacyl-tRNA synthetase and tRNA showed that even RNAs with stable tertiary fold undergo substantial structural changes upon binding to the synthetases. The structural variability of RNA as well as the ability of RNA to distort upon protein binding may be crucial in RNA-protein interactions.     

Initial structural and dynamic characterization of the M2 protein transmembrane and amphipathic helices in lipid bilayers

Amphipathic helices in membrane proteins that interact with the hydrophobic/hydrophilic interface of the lipid bilayer have been difficult to structurally characterize. Here, the backbone structure and orientation of an amphipathic helix in the full-length M2 protein from influenza A virus has been characterized. The protein has been studied in hydrated DMPC/DMPG lipid bilayers above the gel to liquid-crystalline phase transition temperature by solid-state NMR spectroscopy. Characteristic PISA (Polar Index Slant Angle) wheels reflecting helical wheels have been observed in uniformly aligned bilayer preparations of both uniformly 15N labeled and amino acid specific labeled M2 samples. Hydrogen/deuterium exchange studies have shown the very slow exchange of some residues in the amphipathic helix and more rapid exchange for the transmembrane helix. These latter results clearly suggest the presence of an aqueous pore. A variation in exchange rate about the transmembrane helical axis provides additional support for this claim and suggests that motions occur about the helical axes in this tetramer to expose the entire backbone to the pore.     

Structures of the transmembrane helices of the G-protein coupled receptor, rhodopsin.

An hypothesis is tested that individual peptides corresponding to the transmembrane helices of the membrane protein, rhodopsin, would form helices in solution similar to those in the native protein. Peptides containing the sequences of helices 1, 4 and 5 of rhodopsin were synthesized. Two peptides, with overlapping sequences at their termini, were synthesized to cover each of the helices. The peptides from helix 1 and helix 4 were helical throughout most of their length. The N- and C-termini of all the peptides were disordered and proline caused opening of the helical structure in both helix 1 and helix 4. The peptides from helix 5 were helical in the middle segment of each peptide, with larger disordered regions in the N- and C-termini than for helices 1 and 4. These observations show that there is a strong helical propensity in the amino acid sequences corresponding to the transmembrane domain of this G-protein coupled receptor. In the case of the peptides from helix 4, it was possible to superimpose the structures of the overlapping sequences to produce a construct covering the whole of the sequence of helix 4 of rhodopsin. As similar superposition for the peptides from helix 1 also produced a construct, but somewhat less successfully because of the disordering in the region of sequence overlap. This latter problem was more severe for helix 5 and therefore a single peptide was synthesized for the entire sequence of this helix, and its structure determined. It proved to be helical throughout. Comparison of all these structures with the recent crystal structure of rhodopsin revealed that the peptide structures mimicked the structures seen in the whole protein. Thus similar studies of peptides may provide useful information on the secondary structure of other transmembrane proteins built around helical bundles.     

Bulge loops used to measure the helical twist of RNA in solution.

Bulge loops are commonly found in helical segments of cellular RNAs. When incorporated into long double-stranded RNAs, they may introduce points of flexibility or permanent bend that can be detected by the altered electrophoretic gel mobility of the RNA. We find that a single An or Un bulge loop near the middle of a long RNA helix significantly retards the RNA during polyacrylamide gel electrophoresis if n greater than or equal to 2. The mobility of an RNA containing two A2 bulges various periodically with the number of base pairs between the bulges. We interpret this to mean that A2 bulges varies periodically with the number of base pairs between the bulges. We interpret this to mean that Z2 bulges form torsionally stiff bends in the helix; the gel mobility reaches a minimum when the total helical twist between the bulges rotates the arms of the molecule into a cis conformation. The gel mobilities are proportional to the predicted end-to-end distance of the RNA if the average RNA helical repeat is 11.8 +/- 0.2 bp/turn and there is no helical twist (3 +/- 9 degrees) associated with the bulge (data obtained in 0.15 M Na+). Other sizes and sequences of bulges have very different effects on RNA helix conformation and flexibility. U2 bulges bend the helix to a much smaller degree than A2 bulges, while longer A or U bulge sequences probably allow bends of 90 degrees or more; all of these may be fairly flexible joints.(ABSTRACT TRUNCATED AT 250 WORDS)     

Design principles for chlorophyll-binding sites in helical proteins

The cyclic tetrapyrroles, viz. chlorophylls (Chl), their bacterial analogs bacteriochlorophylls, and hemes are ubiquitous cofactors of biological catalysis that are involved in a multitude of reactions. One systematic approach for understanding how Nature achieves functional diversity with only this handful of cofactors is by designing de novo simple and robust protein scaffolds with heme and/or (bacterio)chlorophyll [(B)Chls]-binding sites. This strategy is currently mostly implemented for heme-binding proteins. To gain more insight into the factors that determine heme-/(B)Chl-binding selectivity, we explored the geometric parameters of (B)Chl-binding sites in a nonredundant subset of natural (B)Chl protein structures. Comparing our analysis to the study of a nonredundant database of heme-binding helical histidines by Negron et al. (Proteins 2009;74:400-416), we found a preference for the m-rotamer in (B)Chl-binding helical histidines, in contrast to the preferred t-rotamer in heme-binding helical histidines. This may be used for the design of specific heme- or (B)Chl-binding sites in water-soluble helical bundles, because the rotamer type defines the positioning of the bound cofactor with respect to the helix interface and thus the protein-binding site. Consensus sequences for (B)Chl binding were identified by combining a computational and database-derived approach and shown to be significantly different from the consensus sequences recommended by Negron et al. (Proteins 2009;74:400-416) for heme-binding helical proteins. The insights gained in this work on helix- (B)Chls-binding pockets provide useful guidelines for the construction of reasonable (B)Chl-binding protein templates that can be optimized by computational tools.     

Molecular modeling of the amphipathic helices of the plasma apolipoproteins.

In this paper we propose a classification of the amphipathic helical repeats occurring in the plasma apolipoprotein sequences. It is based upon the calculation of the molecular hydrophobicity potential around the helical segments. The repeats were identified using a new autocorrelation matrix, based upon similarities of hydrophobic and hydrophilic properties of the amino acid residues within the apolipoprotein sequences. The helices were constructed by molecular modeling, the molecular hydrophobicity potential was calculated, and isopotential contour lines drawn around the helices yielded a three-dimensional visualization of the hydrophobicity potential. Two classes of apolipoproteins could be differentiated by comparing the hydrophobic angles obtained by projection of the isopotential contour lines on a plane perpendicular to the long axis of the helix. The isopotential contour lines around apo AI, AIV, and E are more hydrophilic than hydrophobic, whereas they are of similar intensity for apo AII, CI, and CIII. In both cases discoidal lipid-protein complexes are generated, with the amphipathic helices around the edge of the lipid core. The long axis of the helices is oriented parallel to the phospholipid acyl chains and the hydrophilic side of the helix toward the aqueous phase. As a result of the differences in hydrophobicity potential, the contact between the hydrophobic side of the helices and the phospholipid acyl chains is larger for apo AII, CI, and CIII than for the other apolipoproteins. This might account for the greater stability of the discoidal complexes generated between phospholipids and these apoproteins.     

The probable arrangement of the helices in G protein-coupled receptors.

G protein-coupled receptors form a large family of integral membrane proteins whose amino acid sequences have seven hydrophobic segments containing distinctive sequence patterns. Rhodopsin, a member of the family, is known to have transmembrane alpha-helices. The probable arrangement of the seven helices, in all receptors, was deduced from structural information extracted from a detailed analysis of the sequences. Constraints established include: (1) each helix must be positioned next to its neighbours in the sequence; (2) helices I, IV and V must be most exposed to the lipid surrounding the receptor and helix III least exposed. (1) is established from the lengths of the shortest loops. (2) is determined by considering: (i) sites of the most conserved residues; (ii) other sites where variability is restricted; (iii) sites that accommodate polar residues; (iv) sites of differences in sequence between pairs or within groups of closely related receptors. Most sites in the last category should be in unimportant positions and are most useful in determining the position and extent of lipid-facing surface in each helix. The structural constraints for the receptors are used to allocate particular helices to the peaks in the recently published projection map of rhodopsin and to propose a tentative three-dimensional arrangement of the helices in G protein-coupled receptors.     

Stabilization of alpha-helix structure by polar side-chain interactions: complex salt bridges,cation-pi interactions,and C-H em leader O H-bonds

It is generally understood that helical proteins are stabilized by a combination of hydrophobic and packing interactions, together with H-bonds and electrostatic interactions. Here we show that polar side-chain interactions on the surface can play an important role in helix formation and stability. We review studies on model helical peptides that reveal the effect of weak interactions between side chains on helix stability, focusing on some nonclassical side-chain-side-chain interactions: complex salt bridges, cation-pi, and C-H em leader O H-bonding interactions. Each of these can be shown to contribute to helix stability, and thus must be included in a comprehensive catalogue of helix stabilizing effects. The issue of the structure of the unfolded states of helical peptides is also discussed, in the light of recent experiments showing that these contain substantial amounts of polyproline II conformation.     

Geometrical principles of homomeric β‐barrels and β‐helices: Application to modeling amyloid protofilaments

Examples of homomeric β‐helices and β‐barrels have recently emerged. Here we generalize the theory for the shear number in β‐barrels to encompass β‐helices and homomeric structures. We introduce the concept of the “β‐strip,” the set of parallel or antiparallel neighboring strands, from which the whole helix can be generated giving it n‐fold rotational symmetry. In this context, the shear number is interpreted as the sum around the helix of the fixed register shift between neighboring identical β‐strips. Using this approach, we have derived relationships between helical width, pitch, angle between strand direction and helical axis, mass per length, register shift, and number of strands. The validity and unifying power of the method is demonstrated with known structures including α‐hemolysin, T4 phage spike, cylindrin, and the HET‐s(218‐289) prion. From reported dimensions measured by X‐ray fiber diffraction on amyloid fibrils, the relationships can be used to predict the register shift and the number of strands within amyloid protofilaments. This was used to construct models of transthyretin and Alzheimer β(40) amyloid protofilaments that comprise a single strip of in‐register β‐strands folded into a “β‐strip helix.” Results suggest both stabilization of an individual β‐strip helix and growth by addition of further β‐strip helices can involve the same pair of sequence segments associating with β‐sheet hydrogen bonding at the same register shift. This process would be aided by a repeat sequence. Hence, understanding how the register shift (as the distance between repeat sequences) relates to helical dimensions will be useful for nanotube design.     

Alpha helix capping in synthetic model peptides by reciprocal side chain–main chain interactions: Evidence for an N terminal “capping box”

A significant fraction of the amino acids in proteins are alpha helical in conformation. Alpha helices in globular proteins are short, with an average length of about twelve residues, so that residues at the ends of helices make up an important fraction of all helical residues. In the middle of a helix, H-bonds connect the NH and CO groups of each residue to partners four residues along the chain. At the ends of a helix, the H-bond potential of the main chain remains unfulfilled, and helix capping interactions involving bonds from polar side chains to the NH or CO of the backbone have been proposed and detected. In a study of synthetic helical peptides, we have found that the sequence Ser-Glu-Asp-Glu stabilizes the alpha helix in a series of helical peptides with consensus sequences. Following the report by Harper and Rose, which identifies SerXaaXaaGlu as a member of a class of common motifs at the N termini of alpha helices in proteins that they refer to as “capping boxes,” we have reexamined the side chain–main chain interactions in a varient sequence using ¹H NMR, and find that the postulated reciprocal side chain-backbone bonding between the first Ser and last Glu side chains and their peptide NH partners can be resolved: Deletion of two residues N terminal to the Ser-Glu-Asp-Glu sequence in these peptides has no effect on the initiation of helical structure, as defined by two-dimensional (2D) NMR experiments on this variant. Thus the capping box sequence Ser-Glu-Asp-Glu inhibits N terminal fraying of the N terminus of alpha helix in these peptides, and shows the side chain–main chain interactions proposed by Harper and Rose. It thus acts as a helix initiating signal. Since normal a helix cannot propagate beyond the N terminus of this structure, the box acts as a termination signal in this direction as well. © 1994 John Wiley & Sons, Inc.     

Design of stable α‐helices using global sequence optimization

The rational design of peptide and protein helices is not only of practical importance for protein engineering but also is a useful approach in attempts to improve our understanding of protein folding. Recent modifications of theoretical models of helix‐coil transitions allow accurate predictions of the helix stability of monomeric peptides in water and provide new possibilities for protein design. We report here a new method for the design of α‐helices in peptides and proteins using AGADIR, the statistical mechanical theory for helix‐coil transitions in monomeric peptides and the tunneling algorithm of global optimization of multidimensional functions for optimization of amino acid sequences. CD measurements of helical content of peptides with optimized sequences indicate that the helical potential of protein amino acids is high enough to allow formation of stable α‐helices in peptides as short as of 10 residues in length. The results show the maximum achievable helix content (HC) of short peptides with fully optimized sequences at 5 °C is expected to be ～70–75%. Under certain conditions the method can be a powerful practical tool for protein engineering. Unlike traditional approaches that are often used to increase protein stability by adding a few favorable interactions to the protein structure, this method deals with all possible sequences of protein helices and selects the best one from them. Copyright © 2009 European Peptide Society and John Wiley & Sons, Ltd.     

Parameters of helix-coil transition theory for alanine-based peptides of varying chain lengths in water.

Thermal unfolding curves have been measured for a series of short alanine-based peptides that contain repeating sequences and varying chain lengths. Standard helix-coil theory successfully fits the observed transition curves, even for these short peptides. The results provide values for sigma, the helix nucleation constant, delta H0, the enthalpy change on helix formation, and for s (0 degree C), the average helix propagation parameter at 0 degree C. The enthalpy change agrees with the value determined calorimetrically. The success of helix-coil theory in describing the unfolding transitions of short peptides in water indicates that helical propensities, or s values, can be determined from substitution experiments in short alanine-based peptides.     

Helix packing moments reveal diversity and conservation in membrane protein structure

Helical membrane proteins are more tightly packed and the packing interactions are more diverse than those found in helical soluble proteins. Based on a linear correlation between amino acid packing values and interhelical propensity, we propose the concept of a helix packing moment to predict the orientation of helices in helical membrane proteins and membrane protein complexes. We show that the helix packing moment correlates with the helix interfaces of helix dimers of single pass membrane proteins of known structure. Helix packing moments are also shown to help identify the packing interfaces in membrane proteins with multiple transmembrane helices, where a single helix can have multiple contact surfaces. Analyses are described on class A G protein-coupled receptors (GPCRs) with seven transmembrane helices. We show that the helix packing moments are conserved across the class A family of GPCRs and correspond to key structural contacts in rhodopsin. These contacts are distinct from the highly conserved signature motifs of GPCRs and have not previously been recognized. The specific amino acid types involved in these contacts, however, are not necessarily conserved between subfamilies of GPCRs, indicating that the same protein architecture can be supported by a diverse set of interactions. In GPCRs, as well as membrane channels and transporters, amino acid residues with small side-chains (Gly, Ala, Ser, Cys) allow tight helix packing by mediating strong van der Waals interactions between helices. Closely packed helices, in turn, facilitate interhelical hydrogen bonding of both weakly polar (Ser, Thr, Cys) and strongly polar (Asn, Gln, Glu, Asp, His, Arg, Lys) amino acid residues. We propose the use of the helix packing moment as a complementary tool to the helical hydrophobic moment in the analysis of transmembrane sequences.     

Codification and evolution of experimentally observed specific recognition sites for restriction enzymes on DNA

The list of published restriction endonucleases along with their substrates provides an excellent data base for the evaluation of the evolution and codification of the key elements for specific recognition sites on the DNA. In this paper the considerations will be limited to palindromic tetramer-, pentamer-, and hexamer-sequences. It is basically assumed that each base pair within these sequences has to be recognized by directionally unique bidentate hydrogen bonds either within the plane of the base pair or by bridging the appropriate H-bond donor/acceptor groups of the neighbouring bases of the same strand. Thus sequence specificity is mediated by twelve (eight) H-bonds, originating from the protein recognition modules. Besides a pronounced preference for GC base pairs expressed by their high frequency in the most abundant sequences, serving the need of maximal thermodynamic stability of the double helical substrates, it can also be shown that the stacking of consecutive bases within the recognition site sequences plays a major role in shaping the particular DNA/protein interface. Finally it will be demonstrated that the full set of sequences discussed in this paper can readily be derived by stepwise expanding the vocabulary of three simple tetrameric sequences by inserting single base pairs into the centre of a minimal sequence, thus creating all the published pentameric restriction sites, or by inserting/adding two GC base pairs in a palindromic way, thus creating the known multiplicity of hexameric sites.     

Symmetry, homology, and phrasing in the recognition of helical regulatory sequences in DNA.

Regulatory regions in DNA which have been sequenced have generally been found to contain one or more axes of two-fold rotational symmetry. If this symmetry is to be maintained in the helical sequence, the axis of rotation must be aligned with one of the two dyad axes of the helix. This is equivalent to saying that the rotational symmetry of the sequence can only be seen from certain viewing points in a circuit about the helix. More surprising is the fact that new symmetrical sequence arrangements can be seen at +/- 36 degrees, +/- 72 degrees, +/- 108 degrees, and +/- 144 degrees relative to the point at which the rotational symmetry is seen. This "amplification" of symmetry suggests a three-dimensional approach to sequence analysis. A specific reading frame, suggested by the geometry of the helix, is examined with regard to its elucidation of intra- and inter-sequence homologies. Two sequences are thus identified as being recurrent in a number of different regulatory sequences.     

A unified model for the origin of DNA sequence-directed curvature

The fine structure of the DNA double helix and a number of its physical properties depend upon nucleotide sequence. This includes minor groove width, the propensity to undergo the B-form to A-form transition, sequence-directed curvature, and cation localization. Despite the multitude of studies conducted on DNA, it is still difficult to appreciate how these fundamental properties are linked to each other at the level of nucleotide sequence. We demonstrate that several sequence-dependent properties of DNA can be attributed, at least in part, to the sequence-specific localization of cations in the major and minor grooves. We also show that effects of cation localization on DNA structure are easier to understand if we divide all DNA sequences into three principal groups: A-tracts, G-tracts, and generic DNA. The A-tract group of sequences has a peculiar helical structure (i.e., B*-form) with an unusually narrow minor groove and high base-pair propeller twist. Both experimental and theoretical studies have provided evidence that the B*-form helical structure of A-tracts requires cations to be localized in the minor groove. G-tracts, on the other hand, have a propensity to undergo the B-form to A-form transition with increasing ionic strength. This property of G-tracts is directly connected to the observation that cations are preferentially localized in the major groove of G-tract sequences. Generic DNA, which represents the vast majority of DNA sequences, has a more balanced occupation of the major and minor grooves by cations than A-tracts or G-tracts and is thereby stabilized in the canonical B-form helix. Thus, DNA secondary structure can be viewed as a tug of war between the major and minor grooves for cations, with A-tracts and G-tracts each having one groove that dominates the other for cation localization. Finally, the sequence-directed curvature caused by A-tracts and G-tracts can, in both cases, be explained by the cation-dependent mismatch of A-tract and G-tract helical structures with the canonical B-form helix of generic DNA (i.e., a cation-dependent junction model).     

CD and NMR studies of prion protein (PrP) helix 1. Novel implications for its role in the PrPC-->PrPSc conversion process

The conversion of prion helix 1 from an alpha-helical into an extended conformation is generally assumed to be an essential step in the conversion of the cellular isoform PrPC of the prion protein to the pathogenic isoform PrPSc. Peptides encompassing helix 1 and flanking sequences were analyzed by nuclear magnetic resonance and circular dichroism. Our results indicate a remarkably high instrinsic helix propensity of the helix 1 region. In particular, these peptides retain significant helicity under a wide range of conditions, such as high salt, pH variation, and presence of organic co-solvents. As evidenced by a data base search, the pattern of charged residues present in helix 1 generally favors helical structures over alternative conformations. Because of its high stability against environmental changes, helix 1 is unlikely to be involved in the initial steps of the pathogenic conformational change. Our results implicate that interconversion of helix 1 is rather representing a barrier than a nucleus for the PrPC-->PrPSc conversion.