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.     

Predicting helical segments in proteins by a helix-coil transition theory with parameters derived from a structural database of proteins

A novel helix-coil transition theory has been developed. This new theory contains more types of interactions than similar theories developed earlier. The parameters of the models were obtained from a database of 351 nonhomologous proteins. No manual adjustment of the parameters was performed. The interaction parameters obtained in this manner were found to be physically meaningful, consistent with current understanding of helix stabilizing/destabilizing interactions. Novel insights into helix stabilizing/destabilizing interactions have also emerged from this analysis. The theory developed here worked well in sorting out helical residues from amino acid sequences. If the theory was forced to make prediction on every residue of a given amino acid sequence, its performance was the best among ten other secondary structural prediction algorithms in distinguishing helical residues from nonhelical ones. The theory worked even better if one only required it to make prediction on residues that were “predictable” (identifiable by the theory); >90% predictive reliability could be achieved. The helical residues or segments identified by the helix-coil transition theory can be used as secondary structural contraints to speed up the prediction of the three-dimensional structure of a protein by reducing the dimension of a computational protein folding problem. Possible further improvements of this helix-coil transition theory are also discussed. Proteins 28:344–359, 1997. © 1997 Wiley-Liss, Inc.     

The complete amino acid sequence of the major component myoglobin of Amazon river dolphin (Inia geoffrensis).

The complete amino acid sequence of the major component myoglobin from Amazon River dolphin, Inia geoffrensis, was determined by specific cleavage of the protein to obtain large peptides which are readily degraded by the automatic sequencer. Three easily separable peptides were obtained by cleaving the protein with cyanogen bromide at the methionine residues and four peptides were obtained by cleaving the methyl-acetimidated protein with trypsin at the arginine residues. From these peptides over 85% of the sequence was completed. The remainder of the sequence was obtained by fragmentation of the large cyanogen bromide peptide with trypsin. This protein differs from that of the common porpoise, Phocoena phocoena, at seven positions, from that of the common dolphin, Delphinus delphis, at 11 positions, and from that of the sperm whale, Physeter catodon, at 15 positions. By comparison of this sequence with the three-dimensional structure of sperm whale myoglobin it appears that those residues close to the heme group are most conserved followed by those in nonhelical regions and lastly by those in the helical segments. All of the substitutions observed in this sequence fit easily into the three-dimensional structure of the sperm whale myoglobin.     

Structures and roles of the polymorphic forms of tobacco mosaic virus. V. Conservation of the abnormally titrating groups in tobacco mosaic virus

The two groups per subunit which titrate with an abnormal pK of about 7 in tobacco mosaic virus Vulgare have also been found in three other naturally occurring strains and two mutants. Possession of these groups, almost certainly carboxyl-carboxylate pairs, therefore appears to be a crucial feature of the virus protein structure, which has been conserved during the evolution of the strains. The residues responsible may be narrowed down to those carboxylic acid residues which are at the same position in all the tobacco mosaic virus variants titrated. Residues 115 and 116 make up one probable pair, while residue 145 and one other, as yet unidentified, make up the second probable pair.     

Amino acid sequences of alpha-helical segments from S-carbosymethylkerateine-A. Complete sequence of a type-I segment.

The amino acid sequence of a type-I helical segment from the low-sulphur protein (S-carboxymethylkerateine-A) of wool was determined by combining automatic and manual-sequencing data. Whereas in the type-II helical segment most of the cationic groups occur in pairs, 11 of the 22 anionic residues in the sequence of the type-I segment were situated next to a second anionic residue. This suggests possible interactions between type-I and type-II helical segments in alpha-keratin. As observed with the sequence of a type-II helical segment a model constructed on 3.6 residues per turn of helix shows a line of hydrophobic residues along the helix, thereby supporting the physicochemical evidence that the molecule is predominantly helical and forms part of a coiled-coil structure. Examination of the sequence data by predictive methods indicates the possibilty of extensive sections of alpha-helix interspersed with discontinuities. The molecule contains a number of regions with peptide sequences identical with those found by other workers after enzymic digestion of fractions from oxidized wool.     

On the recombinational origin of protein-sequence-subunit structure

Since 1929 the concept that proteins are built from subunits of certain standard size (Svedberg 1929) has been revisited several times, each time with a new demonstration that, indeed, there are certain preferred protein sizes. According to recent estimates the overrepresented sizes are close to multiples of 125 amino acid (aa) residues for eukaryotes and 150 residues for prokaryotes. To explain these preferences, a hypothesis is suggested, and quantitatively developed, on the recombinational nature of this regularity. The protein-coding sequences are assumed to evolve at some early stage via recombinational events—insertions of DNA circles of a certain optimal size. The contour lengths of the protein-coding DNA circles had to be simultaneously divisible by three and, to minimize torsional constraint, by the DNA helical repeat. With these two conditions satisfied, the calculated contour lengths of the DNA circles, 250–500 base pairs (bp), turn out to correspond well to known optimal DNA circularization sizes and to the predicted range of the protein sequence subunit sizes: 80–170 as residues, which covers experimentally observed values. The subunit size is found to be strongly influenced by the helical repeat of DNA. The sizes 125 and 150 as are derived when the corresponding helical repeats of DNA are set within fractions of promilles from the 10.54 by/turn value. This fits to the experimentally estimated mean for natural mixed DNA sequences, 10.53–10.57 by/turn. The suggested recombinational mechanism thus not only gives a qualitative explanation for the observed underlying order in the protein sequences but also quantitatively links the observed protein sequence sizes with the optimal DNA circularization size and with the helical repeat of DNA. It also offers a versatile molecular model of early protein evolution by fusion and insertion of preexisting proteins of standard subunit sizes.     

Interacting helical surfaces of the transmembrane segments of subunits a and c' of the yeast V-ATPase defined by disulfide-mediated cross-linking

Proton translocation by the vacuolar (H+)-ATPase (or V-ATPase) has been shown by mutagenesis to be dependent upon charged residues present within transmembrane segments of subunit a as well as the three proteolipid subunits (c, c', and c"). Interaction between R735 in TM7 of subunit a and the glutamic acid residue in the middle of TM4 of subunits c and c' or TM2 of subunit c" has been proposed to be essential for proton release to the luminal compartment. In order to determine whether the helical face of TM7 of subunit a containing R735 is capable of interacting with the helical face of TM4 of subunit c' containing the essential glutamic acid residue (Glu-145), cysteine-mediated cross-linking between these subunits in yeast has been performed. Cys-less forms of subunits a and c' as well as forms containing unique cysteine residues were constructed, introduced together into a strain disrupted in both endogenous subunits, and tested for growth at neutral pH, for assembly competence and for cross-linking in the presence of cupric-phenanthroline by SDS-PAGE and Western blot analysis. Four different cysteine mutants of subunit a were each tested pairwise with ten different unique cysteine mutants of subunit c'. Strong cross-linking was observed for the pairs aS728C/c'I142C, aA731C/c'E145C, aA738C/c'F143C, aA738C/c'L147C, and aL739C/c'L147C. Partial cross-linking was observed for an additional 13 of 40 pairs analyzed. When arrayed on a helical wheel diagram, the results suggest that the helical face of TM7 of subunit a containing Arg-735 interacts with the helical face of TM4 of subunit c' centered on Val-146 and bounded by Glu-145 and Leu-147. The results are consistent with a possible rotational flexibility of one or both of these transmembrane segments as well as some flexibility of movement perpendicular to the membrane.     

The primary structure of component 8c-1, a subunit protein of intermediate filaments in wool keratin. Relationships with proteins from other intermediate filaments.

Component 8c-1, one of four highly homologous component-8 subunit proteins present in the microfibrils of wool, was isolated as its S-carboxymethyl derivative and its amino acid sequence was determined. Large peptides were isolated after cleaving the protein chemically or enzymically and the sequence of each was determined with an automatic Sequenator. The peptides were ordered by sequence overlaps and, in some instances, by homology with known sequences from other component-8 subunits. The C-terminal residues were identified by three procedures. Full details of the various procedures used have been deposited as Supplementary Publication SUP 50133 (4 pp.) at the British Library Lending Division, Boston Spa, Wetherby, West Yorkshire LS23 7BQ, U.K., from whom copies can be obtained on the terms indicated in Biochem. J. (1986) 233, 5. The result showed that the protein comprises 412 residues and has an Mr, including the N-terminal acetyl group, of 48,300. The sequence of residues 98-200 of component 8c-1 was found to correspond to the partial or complete sequences of four homologous type I helical segments previously isolated from helical fragments recovered from chymotryptic digests of microfibrillar proteins of wool [Crewther & Dowling (1971) Appl. Polym. Symp. 18, 1-20; Crewther, Gough, Inglis & McKern (1978) Text. Res. J. 48, 160-162; Gough, Inglis & Crewther (1978) Biochem. J. 173, 385]. Considered in relation to amino acid sequences of other intermediate-filament proteins, the sequence is in accord with the view that keratin filament proteins are of two types [Hanukoglu & Fuchs (1983) Cell (Cambridge, Mass.) 33, 915-924]. Filament proteins from non-keratinous tissues, such as desmin, vimentin, neurofilament proteins and the glial fibrillary acidic protein, which form monocomponent filaments, constitute a third type. It is suggested that as a whole the proteins from intermediate filaments be classed as filamentins, the three types at present identified forming subgroups of this class. The significant homologies between types I, II and III occur almost exclusively in segments of the chain that have been identified as having a coiled-coil structure together with the relatively short sections connecting these segments. The non-coiled-coil segments at the C- and N-termini show no significant homology between types, nor is homology in these segments apparent in all members of one type. Component 8c-1 does not show homology in its terminal segments with the known sequence of any other filamentin.(ABSTRACT TRUNCATED AT 400 WORDS)     

Stereochemical control of peptide folding

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

The differentially conserved residues of carbamoyl-phosphate synthetase

Carbamoyl-phosphate synthetase (CPS) from Escherichia coli is a heterodimeric protein. The larger of the two subunits (M(r) approximately 118,000) contains a pair of homologous domains of approximately 400 residues each that are approximately 40% identical in amino acid sequence. The carboxy phosphate (residues 1-400) and carbamoyl phosphate domains (residues 553-933) also contain approximately 79 differentially conserved residues. These are residues that are conserved throughout the bacterial evolution of CPS in one of these homologous domains but not the other. The role of these differentially conserved residues in the structural and catalytic properties of CPS was addressed by swapping segments of these residues from one domain to the other. Nine of these chimeric mutant enzymes were constructed, expressed, purified, and characterized. A majority of the mutants were unable to synthesize any carbamoyl phosphate and the rest were severely crippled. True tandem repeat chimeric proteins were constructed by the complete substitution of one homologous domain sequence for the other. Neither of the two possible chimeric proteins was structurally stable. These results have been interpreted to demonstrate that the two homologous domains in the large subunit of CPS are functionally and structurally nonequivalent. This nonequivalence is a direct result of the specific functions each of these domains must perform during the overall synthesis of carbamoyl phosphate in the wild type enzyme and the specific structural alterations imposed by the differentially conserved residues.     

Amino acid sequence of the b subunit of human factor XIII, a protein composed of ten repetitive segments

Factor XIII is a plasma protein that participates in the final stages of blood coagulation. The complete amino acid sequence of the b subunit of human factor XIII was determined by a combination of cDNA cloning and amino acid sequence analysis. A lambda gt11 cDNA library prepared from human liver mRNA was screened with an affinity-purified antibody against the b subunit of human factor XIII. Nine positive clones were isolated from 2 X 10(6) phage and plaque-purified. The largest cDNA insert was sequenced and shown to contain 2180 base pairs coding for a portion of the leader sequence (19 amino acids), the mature protein (641 amino acids), a stop codon (TGA), a 3' noncoding region (187 nucleotides), and a poly(A) tail. When the b subunit of human factor XIII was digested with cyanogen bromide, nine peptides were isolated by gel filtration and reverse-phase high-performance liquid chromatography. Amino acid sequence analyses of these peptides were performed with an automated sequenator, and 299 amino acid residues were identified. These amino acid sequences were in complete agreement with the amino acid sequence predicted from the cDNA. The b subunit of factor XIII contained 10 repetitive homologous segments, each composed of about 60 amino acids and 4 half-cystine residues. Each of these repeated segments is a member of a family of repeats present in human beta 2-glycoprotein I, complement factor B, and haptoglobin alpha 1 chain. Three potential Asn-linked carbohydrate attachment sites were also identified in the b subunit of factor XIII.     

Three-dimensional structure of poliovirus serotype 1 neutralizing determinants.

Antigenic mutants of poliovirus (Sabin strain, serotype 1) were isolated by the resistance of the virus to anti-Sabin neutralizing monoclonal antibodies. The amino acid replacements within the capsid protein sequence causing the altered antigenicity were identified for each of 63 isolates. The mutations cluster into distinct nonoverlapping peptide segments that group into three general immunological phenotypes on the basis of cross-neutralization analyses with 15 neutralizing anti-Sabin monoclonal antibodies. Location of the mutated amino acid residues within the three-dimensional structure of the virion indicates that the majority of these amino acid residues are highly exposed and located within prominent structural features of the viral surface. Those mutated amino acid residues that are less accessible to antibody interaction are often involved in hydrogen bonds or salt bridges that would stabilize the local tertiary structure of the antigenic site. The interactions of the peptide segments that form these neutralizing sites suggest specific models for the generation of neutralization-resistant variants and for the interaction between the viral surface and antibody.     

Free-energy profiles of membrane insertion of the M2 transmembrane peptide from influenza A virus

The insertion of the M2 transmembrane peptide from influenza A virus into a membrane has been studied with molecular-dynamics simulations. This system is modeled by an atomically detailed peptide interacting with a continuum representation of a membrane bilayer in aqueous solution. We performed replica-exchange molecular-dynamics simulations with umbrella-sampling techniques to characterize the probability distribution and conformation preference of the peptide in the solution, at the membrane interface, and in the membrane. The minimum in the calculated free-energy surface of peptide insertion corresponds to a fully inserted, helical peptide spanning the membrane. The free-energy profile also shows that there is a significant barrier for the peptide to enter into this minimum in a nonhelical conformation. The sequence of the peptide is such that hydrophilic amino acid residues at the ends of the otherwise primarily hydrophobic peptide create a trapped, U-shaped conformation with the hydrophilic residues associated with the aqueous phase and the hydrophobic residues embedded in the membrane. Analysis of the free energy shows that the barrier to insertion is largely enthalpic in nature, whereas the membrane-spanning global minimum is favored by entropy.     

Molecular dynamics simulations of helix folding: The effects of amino acid substitution

Molecular dynamics simulations were applied to helix folding of alanine-based synthetic peptides. A single alanine residue in the middle of the peptide was substituted with various nonpolar amino acids (leucine, isoleucine, valine, glycine or proline) to study the effect of the substitution. Unlike many other molecular dynamics simulations, nonhelical initial conformations were used in our simulations to study the folding process. An average solvent effect was included in the energy function to simplify the solvent calculation and to overcome the multiple minima problem. During the simulations, the peptides folded into helices in nanoseconds. Compact structures containing two helical segments were also observed. The calculated helical ratios of the peptides showed the same rank order as observed experimentally for the alanine-based peptides. Within a peptide, the helical ratio of each residue was calculated and a minimum was found near the center of the sequence for all peptides. The substitutions had different asymmetric effects on the helical ratios of the residues preceding and following the substitution site, indicating different helix capping preferences of the substituting amino acids. © 1997 John Wiley & Sons, Inc. Biopoly 42: 633–644, 1997     

Dynamic alpha-helices: conformations that do not conform

Structural transitions are important for the stability and function of proteins, but these phenomena are poorly understood. An extensive analysis of Protein Data Bank entries reveals 103 regions in proteins with a tendency to transform from helical to nonhelical conformation and vice versa. We find that these dynamic helices, unlike other helices, are depleted in hydrophobic residues. Furthermore, the dynamic helices have higher surface accessibility and conformational mobility (P-value = 3.35e-07) than the rigid helices. Contact analyses show that these transitions result from protein-ligand, protein-nucleic acid, and crystal-contacts. The immediate structural environment differs quantitatively (P-value = 0.003) as well as qualitatively in the two alternate conformations. Often, dynamic helix experiences more contacts in its helical conformation than in the nonhelical counterpart (P-value = 0.001). There is differential preference for the type of short contacts observed in two conformational states. We also demonstrate that the regions in protein that can undergo such large conformational transitions can be predicted with a reasonable accuracy using logistic regression model of supervised learning. Our findings have implications in understanding the molecular basis of structural transitions that are coupled with binding and are important for the function and stability of the protein. Based on our observations, we propose that several functionally relevant regions on the protein surface can switch over their conformation from coil to helix and vice-versa, to regulate the recognition and binding of their partner and hence these may work as "molecular switches" in the proteins to regulate certain biological process. Our results supports the idea that protein structure-function paradigm should transform from static to a highly dynamic one.     

Human interstitial retinoid-binding protein. Gene structure and primary structure

Interstitial retinoid-binding protein (IRBP) is synthesized and secreted by rod photoreceptor cells into the interphotoreceptor matrix and is known to bind retinoids and fatty acids. We have used cDNA clones encoding human IRBP to isolate a 15-kilobase genomic fragment that encompasses the complete human IRBP gene. The IRBP gene spans more than 11 kilobases and is interrupted by three introns, all of which are positioned near the 3'-end of the coding sequence. The 3741-base pair coding region of IRBP appears to have been generated by quadruplication of an approximately 900 base pair long ancestral gene. The deduced amino acid sequence predicts a mature protein of 1,230 residues (calculated molecular weight 133,000). The protein sequence can be aligned into four homologous segments, each consisting of about 300 residues. Sequence similarity between segments is as high as 60% when conservative substitutions are taken into account. Two putative N-linked glycosylation sites are located in highly conserved domains in the center of the first and second segment of IRBP. A domain consisting of 41 residues at the COOH-terminal end of the third segment has 15 matching residues (38%) with an intradiscal loop of rhodopsin, a retinal-binding protein in rod photoreceptors.     

Do the ends justify the mean? Proline mutations at the ends of the keratin coiled-coil rod segment are more disruptive than internal mutations

Intermediate filament (IF) assembly is remarkable, in that it appears to be self-driven by the primary sequence of IF proteins, a family (40-220 kd) with diverse sequences, but similar secondary structures. Each IF polypeptide has a central 310 amino acid residue alpha-helical rod domain, involved in coiled-coil dinner formation. Two short (approximately 10 amino acid residue) stretches at the ends of this rod are more highly conserved than the rest, although the molecular basis for this is unknown. In addition, the rod is segmented by three short nonhelical linkers of conserved location, but not sequence. To examine the degree to which different conserved helical and nonhelical rod sequences contribute to dimer, tetramer, and higher ordered interactions, we introduced proline mutations in residues throughout the rod of a type I keratin, and we removed existing proline residues from the linker regions. To further probe the role of the rod ends, we introduced more subtle mutations near the COOH-terminus. We examined the consequences of these mutations on (a) IF network formation in vivo, and (b) 10-nm filament assembly in vitro. Surprisingly, all proline mutations located deep in the coiled-coil rod segment showed rather modest effects on filament network formation and 10-nm filament assembly. In addition, removing the existing proline residues was without apparent effect in vivo, and in vitro, these mutants assembled into 10-nm filaments with a tendency to aggregate, but with otherwise normal appearance. The most striking effects on filament network formation and IF assembly were observed with mutations at the very ends of the rod. These data indicate that sequences throughout the rod are not equal with respect to their role in filament network formation and in 10-nm filament assembly. Specifically, while the internal rod segments seem able to tolerate considerable changes in alpha-helical conformation, the conserved ends seem to be essential for creating a very specific structure, in which even small perturbations can lead to loss of IF stability and disruption of normal cellular interactions. These findings have important implications for the disease Epidermolysis Bullosa Simplex, arising from point mutations in keratins K5 or K14.     

Importance of context in protein folding: secondary structural propensities versus tertiary contact-assisted secondary structure formation

Molecular dynamics simulations can be used to reveal the detailed conformational behaviors of peptides and proteins. By comparing fragment and full-length protein simulations, we can investigate the role of each peptide segment in the folding process. Here, we take advantage of information regarding the helix formation process from our previous simulations of barnase and protein A as well as new simulations of four helical fragments from these proteins at three different temperatures, starting with both helical and extended structures. Segments with high helical propensity began the folding process by tethering the chain through side chain interactions involving either polar interactions, such as salt bridges, or hydrophobic staples. These tethers were frequently nonnative (i.e., not i --> i + 4 spacing) and provided a scaffold for other residues, thereby limiting the conformational search. The helical structure then propagated on both sides of the tether. Segments with low stability and propensity formed later in the folding process and utilized contacts with other portions of the protein when folding. These helices formed via a tertiary contact-assisted mechanism, primarily via hydrophobic contacts between residues distant in sequence. Thus, segments with different helical propensities appear to play different roles during protein folding. Furthermore, the active role of nonlocal side chains in helix formation highlights why we must move beyond simple hierarchical models of protein folding.