Helix geometry in proteins

In this report we describe a general survey of all helices found in 57 of the known protein crystal structures, together with a detailed analysis of 48 alpha-helices found in 16 of the structures that are determined to high resolution. The survey of all helices reveals a total of 291 alpha-helices, 71 3(10)-helices and no examples of pi-helices. The conformations of the observed helices are significantly different from the "ideal" linear structures. The mean phi, psi angles for the alpha- and 3(10)-helices found in proteins are, respectively, (-62 degrees, -41 degrees) and (-71 degrees, -18 degrees). A computer program, HBEND, is used to characterize and to quantify the different types of helix distortion. alpha-Helices are classified as regular or irregular, linear, curved or kinked. Of the 48 alpha-helices analysed, only 15% are considered to be linear; 17% are kinked, and 58% are curved. The curvature of helices is caused by differences in the peptide hydrogen bonding on opposite faces of the helix, reflecting carbonyl-solvent/side-chain interactions for the exposed residues, and packing constraints for residues involved in the hydrophobic core. Kinked helices arise either as a result of included proline residues, or because of conflicting requirements for the optimal packing of the helix side-chains. In alpha-helices where there are kinks caused by proline residues, we show that the angle of kink is relatively constant (approximately 26 degrees), and that there is minimal disruption of the helix hydrogen bonding. The proline residues responsible for the kinks are highly conserved, suggesting that these distortions may be structurally/functionally important.     

Proline in alpha-helix: stability and conformation studied by dynamics simulation

Free-energy simulations have been used to estimate the change in the conformational stability of short polyalanine alpha-helices when one of the alanines is replaced by a proline residue. For substituting proline in the middle of the helix the change in free energy of folding (delta delta G degrees) was calculated as 14 kJ/mol (3.4 kcal/mol), in excellent agreement with the one available experimental value. The helix containing proline was found to be strongly kinked; the free energy for reducing the angle of the kink from 40 degrees to 15 degrees was calculated, and found to be small. A tendency to alternate hydrogen bonding schemes was observed in the proline-containing helix. These observations for the oligopeptide agree well with the observation of a range of kink angles (18-35 degrees) and variety of hydrogen bonding schemes, in the rare instances where proline occurs in helices in globular proteins. For substituting proline at the N-terminus of the helix the change in free energy of folding (delta delta G degrees) was calculated as -4 kJ/mol in the first helical position (N1) and +6 kJ/mol in the second helical position (N2). The observed frequent occurrence of proline in position N1 in alpha-helices in proteins therefore has its origin in stability differences of secondary structure. The conclusion reached here that proline may be a better helix former in position N1 than (even) alanine, and thus be a helix initiator may be testable experimentally by measurements of fraction helical conformation of individual residues in oligopeptides of appropriate sequence. The relevance of these results in regards to the frequent occurrence of proline-containing helices in certain membrane proteins is discussed.     

Proline-induced distortions of transmembrane helices

Proline residues in the transmembrane (TM) alpha-helices of integral membrane proteins have long been suspected to play a key role for helix packing and signal transduction by inducing regions of helix distortion and/or dynamic flexibility (hinges). In this study we try to characterise the effect of proline on the geometric properties of TM alpha-helices. We have examined 199 transmembrane alpha-helices from polytopic membrane proteins of known structure. After examining the location of proline residues within the amino acid sequences of TM helices, we estimated the helix axes either side of a hinge and hence identified a hinge residue. This enabled us to calculate helix kink and swivel angles. The results of this analysis show that proline residues occur with a significant concentration in the centre of sequences of TM alpha-helices. In this location, they may induce formation of molecular hinges, located on average about four residues N-terminal to the proline residue. A superposition of proline-containing TM helices structures shows that the distortion induced is anisotropic and favours certain relative orientations (defined by helix kink and swivel angles) of the two helix segments.     

High-resolution prediction of protein helix positions and orientations

We have developed a new method for predicting helix positions in globular proteins that is intended primarily for comparative modeling and other applications where high precision is required. Unlike helix packing algorithms designed for ab initio folding, we assume that knowledge is available about the qualitative placement of all helices. However, even among homologous proteins, the corresponding helices can demonstrate substantial differences in positions and orientations, and for this reason, improperly positioned helices can contribute significantly to the overall backbone root-mean-square deviation (RMSD) of comparative models. A helix packing algorithm for use in comparative modeling must obtain high precision to be useful, and for this reason we utilize an all-atom protein force field (OPLS) and a Generalized Born continuum solvent model. To reduce the computational expense associated with using a detailed, physics-based energy function, we have developed new hierarchical and multiscale algorithms for sampling the helices and flanking loops. We validate the method using a test suite of 33 cases, which are drawn from a diverse set of high-resolution crystal structures. The helix positions are reproduced with an average backbone RMSD of 0.6 A, while the average backbone RMSD of the complete loop-helix-loop region (i.e., the helix with the surrounding loops, which are also repredicted) is 1.3 A.     

Role of stabilization centers in 4 helix bundle proteins

Stabilization centers (SCs) were shown to play an important role in preventing decay of three-dimensional protein structures. These residue clusters, stabilized by cooperative long-range interactions, were proposed to serve as anchoring points for arranging secondary structure elements. In all-alpha proteins, SC elements appear less frequently than in all-beta, alpha/beta, and alpha+beta proteins suggesting that tertiary structure formation of all-alpha proteins is governed by different principles than in other protein classes. Here we analyzed the relation between the formation of stabilization centers and the inter-axial angles (Omega) of alpha-helices in 4 helix bundle proteins. In the distance range, where dipoles have dominant effect on the helix pair arrangement, those helix pairs, where residues from both helices participate in SC elements, appear as parallel more frequently than those helices where no SC elements are present. For SC containing helix pairs, the energetic difference between the parallel and anti-parallel states decreases considerably from 1.1 kcal/mol to 0.4 kcal/mol. Although the observed effect is weak for more distant helices, a competition between the SC element formation and the optimal dipole-dipole interaction of alpha-helices is proposed as a mechanism for tertiary structure formation in 4 helix bundle proteins. The SC-forming potential of different arrangements as well as the pitfalls of the SC definition are also discussed.     

Interhelical hydrogen bonding drives strong interactions in membrane proteins

Polar residues in transmembrane alpha-helices may strongly influence the folding or association of integral membrane proteins. To test whether a motif that promotes helix association in a soluble protein could do the same within a membrane, we designed a model transmembrane helix based on the GCN4 leucine zipper. We found in both detergent micelles and biological membranes that helix association is driven strongly by asparagine, independent of the rest of the hydrophobic leucine and/or valine sequence. Hydrogen bonding between membrane helices gives stronger associations than the packing of surfaces in glycophorin A helices, creating an opportunity to stabilize structures, but also implying a danger that non-specific interactions might occur. Thus, membrane proteins may fold to avoid exposure of strongly hydrogen bonding groups at their lipid exposed surfaces.     

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.     

Identification of Helical Packing Motifs Common to Bacteriorhodopsin and G Protein-Coupled Receptors

A sequence analysis and comparison of transmembrane helices in bacteriorhodopsin (BR) and G protein-coupled receptors (GPCRs) is presented to identify potential regions of homology across protein families. The results show a common pattern of residues is conserved within the interhelical contact regions of BR that fit a knob-into-hole structural motif previously postulated for globular proteins and photosynthetic reaction centers. Based on an alignment of conserved prolines in transmembrane helices, it is inferred that analogous helix packing arrangements are possible in the rhodopsin-like GPCRs. Molecular models of GPCR helices V and VI indicate these interactions occur between aromatic and hydrophobic residues flanking the highly conserved prolines in these sequences. A similar packing arrangement is shown to occur in the X-ray structure of the melittin which also displays a unique pairing of proline-linked helices. The contact pattern identified is further applied to predict the packing of pairs of proline-containing helices in the pheromone-like and cAMP GPCRs. A potential role in stabilizing structure formation is also suggested for the contacts. The results and conclusions are supported by recent biophysical studies of zinc binding to kappa-opioid receptor mutants.     

Identification of Helical Packing Motifs Common to Bacteriorhodopsin and G Protein-Coupled Receptors

A sequence analysis and comparison of transmembrane helices in bacteriorhodopsin (BR) and G protein-coupled receptors (GPCRs) is presented to identify potential regions of homology across protein families. The results show a common pattern of residues is conserved within the interhelical contact regions of BR that fit a knob-into-hole structural motif previously postulated for globular proteins and photosynthetic reaction centers. Based on an alignment of conserved prolines in transmembrane helices, it is inferred that analogous helix packing arrangements are possible in the rhodopsin-like GPCRs. Molecular models of GPCR helices V and VI indicate these interactions occur between aromatic and hydrophobic residues flanking the highly conserved prolines in these sequences. A similar packing arrangement is shown to occur in the X-ray structure of the melittin which also displays a unique pairing of proline-linked helices. The contact pattern identified is further applied to predict the packing of pairs of proline-containing helices in the pheromone-like and cAMP GPCRs. A potential role in stabilizing structure formation is also suggested for the contacts. The results and conclusions are supported by recent biophysical studies of zinc binding to kappa-opioid receptor mutants.     

Distribution of charged residues stabilizes individual helices in myoglobins

The effect of the distribution of charged residues on stability of alpha helices in isolated peptides and in globular proteins exemplified by myoglobins from 62 different species is discussed. A highly simplified set of rules is used to account for the interaction of charged groups with the dipole of an alpha helix. Only the position and sign of a charge with respect to the center of the helix and its ability to participate in intrahelical salt bridges determine its effect. These rules lead to a linear correlation between the helicity in variant C-peptide helices from RNAse and the extent to which the charge distribution opposes the helix dipole. Of the sample of 496 helices in the myoglobins studied, 456 exhibit arrangements of charges which oppose the effective dipole moment of the helix according to this calculation. A number of variants occur which leave the backbone moment of helices A-D unchanged, or even add to it. However no such variants exist in the sequences of helices E-H. We suggest that the E, F, G and H helices in myoglobins which show the strongest reversal of the helix dipole participate in the structures of early intermediates in folding of the chain. Stable helix structures should be more likely to occur in these isolated sequences also, and introduction of charge alterations in helices E to H should affect the initial refolding rate of mutant myoglobins.     

Structure of the 1-36 N-terminal fragment of human phospholamban phosphorylated at Ser-16 and Thr-17

The structure of a 36-amino-acid-long N-terminal fragment of human phospholamban phosphorylated at Ser-16 and Thr-17 and Cys-36-->Ser mutated was determined from nuclear magnetic resonance data in aqueous solution containing 30% trifluoroethanol. The peptide assumes a conformation characterized by two alpha-helices connected by an irregular strand, which comprises the amino acids from Arg-13 to Pro-21. The proline is in a trans conformation. The two phosphate groups on Ser-16 and Thr-17 are shown to interact preferably with the side chains of Arg-14 and Arg-13, respectively. The helix comprising amino acids 22 to 35 is well determined (the rmsd for the backbone atoms, calculated for a family of 24 nuclear magnetic resonance structures is 0.69 +/- 0.28 A). The structures of phosphorylated and unphosphorylated phospholamban are compared, and the effect of the two phosphate groups on the relative spatial position of the two helices is examined. The packing parameters Omega (interhelical angle) and d (minimal interhelical distance) are calculated: in the case of the phosphorylated phospholamban, Omega = 100 +/- 35 degrees and d = 7.9 +/- 4.6 A, whereas for the unphosphorylated peptide the values are Omega = 80 +/- 20 degrees and d = 7.0 +/- 4.0 A. We conclude that 1) the phosphorylation does not affect the structure of the C terminus between residues 21 and 36 and 2) the phosphorylated phospholamban has more loose helical packing than the nonphosphorylated.     

Bundles of amphipathic transmembrane alpha-helices as a structural motif for ion-conducting channel proteins: studies on sodium channels and acetylcholine receptors

Channel proteins are transmembrane symmetric (or pseudosymmetric) oligomers organized around a central ionic pore. We present here a molecular model of the pore forming structures of two channel proteins with different primary structures and oligomeric size: the voltage-sensitive sodium channel and the nicotinic cholinergic receptor. We report low-energy arrangements of alpha-helical bundles calculated by semiempiricial potential energy functions and optimization routines and further refined using molecular dynamics. The ion-conducting pore is considered to be a symmetric or pseudosymmetric homooligomer of 3-5 amphipathic alpha-helices arranged such that the polar residues line a central hydrophilic pathway and the apolar residues face the hydrophobic bilayer interior. The channel lining exposes either charged (Asp, Glu, Arg, Lys) or polar-neutral (Ser, Thr) residues. A bundle of four parallel helices constrained to C4 symmetry, the helix axis aligned with the symmetry axis, and the helices constrained to idealized dihedral angles, produces a structure with a pore of the size inferred for the sodium channel protein (area approximately 16 A2). Similarly, a pentameric array optimized with constraints to maintain C5 symmetry and backbone torsions characteristic of alpha-helices adopts a structure that appears well suited to form the lining of the nicotinic cholinergic receptor (pore area approximately 46 A2). Thus, bundles of amphipathic alpha-helices satisfy the structural, energetic, and dynamic requirements to be the molecular structural motif underlying the function of ionic channels.     

Structural basis of ligand binding and release in insect pheromone-binding proteins: NMR structure of Antheraea polyphemus PBP1 at pH 4.5

The NMR structure of the Antheraea polyphemus pheromone-binding protein 1 at pH 4.5, ApolPBP1A, was determined at 20 degrees C. The structure consists of six alpha-helices, which are arranged in a globular fold that encapsulates a central helix alpha7 formed by the C-terminal polypeptide segment 131-142. The 3D arrangement of these helices is anchored by the three disulfide bonds 19-54, 50-108 and 97-117, which were identified by NMR. Superposition of the ApolPBP1A structure with the structure of the homologous pheromone-binding protein of Bombyx mori at pH 4.5, BmorPBPA, yielded an rmsd of 1.7 A calculated for the backbone heavy-atoms N, Calpha and C' of residues 10-142. In contrast, the present ApolPBP1A structure is different from a recently proposed molecular model for a low-pH form of ApolPBP1 that does not contain the central helix alpha7. ApolPBP1 exhibits a pH-dependent transition between two different globular conformations in slow exchange on the NMR chemical shift timescale similar to BmorPBP, suggesting that the two proteins use the same mechanism of ligand binding and ejection. The extensive sequence homology observed for pheromone-binding proteins from moth species further implies that the previously proposed mechanism of ligand ejection involving the insertion of a C-terminal helix into the pheromone-binding site is a general feature of pheromone signaling in moths.     

Helix packing in polytopic membrane proteins: role of glycine in transmembrane helix association.

The nature and distribution of amino acids in the helix interfaces of four polytopic membrane proteins (cytochrome c oxidase, bacteriorhodopsin, the photosynthetic reaction center of Rhodobacter sphaeroides, and the potassium channel of Streptomyces lividans) are studied to address the role of glycine in transmembrane helix packing. In contrast to soluble proteins where glycine is a noted helix breaker, the backbone dihedral angles of glycine in transmembrane helices largely fall in the standard alpha-helical region of a Ramachandran plot. An analysis of helix packing reveals that glycine residues in the transmembrane region of these proteins are predominantly oriented toward helix-helix interfaces and have a high occurrence at helix crossing points. Moreover, packing voids are generally not formed at the position of glycine in folded protein structures. This suggests that transmembrane glycine residues mediate helix-helix interactions in polytopic membrane proteins in a fashion similar to that seen in oligomers of membrane proteins with single membrane-spanning helices. The picture that emerges is one where glycine residues serve as molecular notches for orienting multiple helices in a folded protein complex.     

Is counterion delocalization responsible for collapse in RNA folding?

Although energetic and phylogenetic methods have been very successful for prediction of nucleic acid secondary structures, arrangement of these secondary structure elements into tertiary structure has remained a difficult problem. Here we explore the packing arrangements of DNA, RNA, and DNA/RNA hybrid molecules in crystals. In the conventional view, the highly charged double helix will be pushed toward isolation by favorable solvation effects; interactions with other like-charged stacks would be strongly disfavored. Contrary to this expectation, we find that most of the cases analyzed ( approximately 80%) exhibit specific, preferential packing between elements of secondary structure, which falls into three categories: (i) interlocking of major grooves of two helices, (ii) side-by-side parallel packing of helices, and (iii) placement of the ribose-phosphate backbone ridge of one helix into the major groove of another. The preponderance of parallel packing motifs is especially surprising. This category is expected to be maximally disfavored by charge repulsion. Instead, it comprises in excess of 50% of all packing interactions in crystals of A-form RNA and has also been observed in crystal structures of large RNA molecules. To explain this puzzle, we introduce a novel model for RNA folding. A simple calculation suggests that the entropy gained by a cloud of condensed cations surrounding the helices more than offsets the Coulombic repulsion of parallel arrangements. We propose that these condensed counterions are responsible for entropy-driven RNA collapse, analogous to the role of the hydrophobic effect in protein folding.     

Ca2+ -ATPase structure in the E1 and E2 conformations: mechanism,helix-helix and helix-lipid interactions

The determination of the crystal structure of the Ca(2+)-ATPase of sarcoplasmic reticulum (SR) in its Ca(2+)-bound [Nature 405 (2000) 647] and Ca(2+)-free forms [Nature 418 (2002) 605] gives the opportunity for an analysis of conformational changes on the Ca(2+)-ATPase and of helix-helix and helix-lipid interactions in the transmembrane (TM) region of the ATPase. The locations of the ends of the TM alpha-helices on the cytoplasmic side of the membrane are reasonably well defined by the location of Trp residues and by the location of Lys-262 that snorkels up to the surface. The locations of the lumenal ends of the helices are less clear. The position of Lys-972 on the lumenal side of helix M9 suggests that the hydrophobic thickness of the protein is only about 21 A, rather than the normal 30 A. The experimentally determined TM alpha-helices do not agree well with those predicted theoretically. Charged headgroups are required for strong interaction of lipids with the ATPase, consistent with the large number of charged residues located close to the lipid-water interface. Helix packing appears to be rather irregular. Packing of helices M8 and M10 is of the 3-4 ridges-into-grooves or knobs-into-holes types. Packing of helices M5 and M7 involves two Gly residues in M7 and one Gly residue in M5. Packing of the other helices generally involves just one or two residues on each helix at the crossing point. The irregular packing of the TM alpha-helices in the Ca(2+)-ATPase, combined with the diffuse structure of the ATPase on the lumenal side of the membrane, is suggested to lead to a relative low activation energy for changing the packing of the TM alpha-helices, with changes in TM alpha-helical packing being important in the process of transfer of Ca(2+) ions across the membrane. The inhibitor thapsigargin binds in a cleft between TM alpha-helices M3, M5 and M7. It is suggested that this and other similar clefts provide binding sites for a variety of hydrophobic molecules affecting the activity of the Ca(2+)-ATPase.     

A segment-based approach to protein secondary structure prediction.

Amino acid sequence patterns have been used to identify the location of turns in globular proteins [Cohen et al. (1986) Biochemistry 25, 266-275]. We have developed sequence patterns that facilitate the prediction of helices in all helical proteins. Regular expression patterns recognize the component parts of a helix: the amino terminus (N-cap), the core of the helix (core), and the carboxy terminus (C-cap). These patterns recognize the core features of helices with a 95% success rate and the N- and C-capping features with success rates of 56% and 48%, respectively. A metapattern language, ALPPS, coordinates the recognition of turns and helical components in a scheme that predicts the location and extent of alpha-helices. On the basis of raw residue scoring, a 71% success rate is observed. By focusing on the recognition of core helical features, we achieve a 78% success rate. Amended scoring procedures are presented and discussed, and comparisons are made to other predictive schemes.     

Seven-helix bundles: molecular modeling via restrained molecular dynamics.

Simulated annealing via restrained molecular dynamics (SA/MD) has been used to model compact bundles of seven approximately (anti)parallel alpha-helices. Seven such helix bundles occur, e.g., in bacteriorhodopsin, in rhodopsin, and in the channel-forming N-terminal domain of Bacillus thuringiensis delta-endotoxin. Two classes of model are considered: (a) those consisting of seven Ala20 peptide chains; and (b) those containing a single polypeptide chain, made up of seven Ala20 helices linked by GlyN interhelix loops (where N = 5 or 10). Three different starting C alpha templates for SA/MD are used, in which the seven helices are arranged (a) on a left-handed circular template, (b) on a bacteriorhodopsin-like template, or (c) on a zig-zag template. The ensembles of models generated by SA/MD are analyzed in terms of their geometry and energetics, and the most stable structures from each ensemble are examined in greater detail. Structures resembling bacteriorhodopsin and structures resembling delta-endotoxin are both represented among the most stable structures. delta-Endotoxin-like structures arise from both circular and bacteriorhodopsin-like C alpha templates. A third helix-packing mode occurs several times among the stable structures, regardless of the C alpha template and of the presence or absence of interhelix loops. It is characterized by a "4 + 1" core, in which four helices form a distorted left-handed supercoil around a central, buried helix. The remaining two helices pack onto the outside of the core. This packing mode is comparable with that proposed for rhodopsin on the basis of two-dimensional electron crystallographic and sequence analysis studies.     

A comparative analysis of interhelical polar interactions of various alpha-helix packings in proteins

One hundred twenty globular proteins and forty five "leucine zippers" representing all types of packing of long alpha-helices were studied in terms of revealing and comparing their interhelical hydrogen and salt bonds. Many previous studies of "leucine zippers" and their analogs showed that interhelical interactions between polar groups could impart specificity to packing of an alpha-helix. The current comparison demonstrated that basically, globular proteins and "leucine zippers" had similar interhelical polar interactions with presumably a similar structural role. However, depending on packing of alpha-helices, the networks of interhelical polar bonds were shown to be distinct and determined both by physicochemical properties of involved amino acid residues and by the relative positions of hydrophobic and hydrophilic residues on the surface of alpha-helices. The revealed distinction is probably crucial for selecting the unique packing of an alpha-helix.