Amino acid code of protein secondary structure

The calculation of protein three-dimensional structure from the amino acid sequence is a fundamental problem to be solved. This paper presents principles of the code theory of protein secondary structure, and their consequence--the amino acid code of protein secondary structure. The doublet code model of protein secondary structure, developed earlier by the author (Shestopalov, 1990), is part of this theory. The theory basis are: 1) the name secondary structure is assigned to the conformation, stabilized only by the nearest (intraresidual) and middle-range (at a distance no more than that between residues i and i + 5) interactions; 2) the secondary structure consists of regular (alpha-helical and beta-structural) and irregular (coil) segments; 3) the alpha-helices, beta-strands and coil segments are encoded, respectively, by residue pairs (i, i + 4), (i, i + 2), (i, i = 1), according to the numbers of residues per period, 3.6, 2, 1; 4) all such pairs in the amino acid sequence are codons for elementary structural elements, or structurons; 5) the codons are divided into 21 types depending on their strength, i.e. their encoding capability; 6) overlappings of structurons of one and the same structure generate the longer segments of this structure; 7) overlapping of structurons of different structures is forbidden, and therefore selection of codons is required, the codon selection is hierarchic; 8) the code theory of protein secondary structure generates six variants of the amino acid code of protein secondary structure. There are two possible kinds of model construction based on the theory: the physical one using physical properties of amino acid residues, and the statistical one using results of statistical analysis of a great body of structural data. Some evident consequences of the theory are: a) the theory can be used for calculating the secondary structure from the amino acid sequence as a partial solution of the problem of calculation of protein three-dimensional structure from the amino acid sequence, and the calculated secondary structure and codon strength distribution can be used for simulating the next step of protein folding; b) one can propose that the same secondary structures can be folded into different tertiary structures and, vice versa, different secondary structures can be folded into the same tertiary structures, provided codon distributions are considered also; c) codons can be considered as first elements of protein three-dimensional structure language.     

Cholera toxin A subunit: functional sites correlated with regions of secondary structure

The A subunit of cholera toxin contains the ADP-ribosyltransferase activity in its major constituent polypeptide A1 (Mr 23,000) which is responsible for the elevation of cAMP typically observed with most mammalian cell types after exposure to the toxin. The primary structure of the A subunit, recently established by sequence analyses, is presented and used as the basis for the secondary structure prediction according to the method of Chou and Fasman. The results indicated the presence of 27% alpha-helix, 25% beta-structure, 12% beta-turn, and 36% random coil. The majority of the beta-structure consisted of six strands located in the NH2-terminal portion of the molecule (residues 33-106) covering one-half of the region corresponding to the A1 polypeptide portion. The beta-sheet domain led immediately into the active site region characterized by the alternating structures of beta-pleated sheet and alpha-helix (residues 95-140) similar to that reported for other NAD+ binding proteins. The presence of this structural feature in the region was confirmed by the use of another predictive method (J. Garnier et al., J. Mol. Biol. 1978, 120, 97-120). In addition, two regions (residues 14-18 and 200-214), previously identified to contain binding sites for the B subunit as evidenced by chemical modification and monoclonal antibody studies, were found to be in alpha-helix configuration.     

A critical assessment of the secondary structure alpha-helices and their termini in proteins

Secondary structure prediction from amino acid sequence is a key component of protein structure prediction, with current accuracy at approximately 75%. We analysed two state-of-the-art secondary structure prediction methods, PHD and JPRED, comparing predictions with secondary structure assigned by the algorithms DSSP and STRIDE. The specific focus of our study was alpha-helix N-termini, as empirical free energy scales are available for residue preferences at N-terminal positions. Although these prediction methods perform well in general at predicting the alpha-helical locations and length distributions in proteins, they perform less well at predicting the correct helical termini. For example, although most predicted alpha-helices overlap a real alpha-helix (with relatively few completely missed or extra predicted helices), only one-third of JPRED and PHD predictions correctly identify the N-terminus. Analysis of neighbouring N-terminal sequences to predicted helical N-termini shows that the correct N-terminus is often within one or two residues. More importantly, the true N-terminal motif is, on average, more favourable as judged by our experimentally measured free energies. This suggests a simple, but powerful, strategy to improve secondary structure prediction using empirically derived energies to adjust the predicted output to a more favourable N-terminal sequence.     

Low-ultraviolet circular dichroism spectroscopy of oligopeptides 1-95 and 96-168 derived from myelin basic protein of rabbit

Myelin basic protein (MBP) is a major protein constituent of the myelin sheath of the central nervous system, where it is believed to have functional alpha-helical segments. One element of the function of the protein might be "conformational adaptability" of specific regions of its amino acid sequence, since the purified protein appears to be largely devoid of ordered structure. To pursue this question, low-ultraviolet circular dichroism (CD) spectroscopy was conducted on the sequential thrombic peptides 1-95 and 96-168 of the protein in the presence of 0-92% trifluoroethanol (TFE), a solvent known to promote stable secondary structures in polypeptides. The series of CD spectra of the oligopeptides were subjected to a computerized best-fit analysis of four peptide conformations, the alpha-helix, beta-structure, beta-turn, and nonordered form. Agreement between experimental and best-fit composite spectra was achieved when standard CD curves of peptide conformations were derived from known theoretical spectra and experimental spectra of polypeptides. In dilute buffer alone, oligopeptides 1-95 and 96-168 evidence no alpha-helix but significant beta-structure (18% and 23%, respectively), as well as a predominant, extended nonordered conformation. However, the two parts of the protein differed in conformational adaptability. From 0% to 30% TFE, 96-168 exhibited concomitant transitions to 10% helix and 32% beta-structure from the nonordered form. In contrast, in 10-30% TFE, 1-95 underwent a transition to approximately 21% helix with partial loss of beta-structure as well as nonordered form; higher concentrations of TFE (40-75%) promoted additional transitions to both helix and beta-structure (totaling 33% and 25%, respectively).(ABSTRACT TRUNCATED AT 250 WORDS)     

Study and prediction of secondary structure for membrane proteins

In this paper we present a novel approach to membrane protein secondary structure prediction based on the statistical stepwise discriminant analysis method. A new aspect of our approach is the possibility to derive physical-chemical properties that may affect the formation of membrane protein secondary structure. The certain physical-chemical properties of protein chains can be used to clarify the formation of the secondary structure types under consideration. Another aspect of our approach is that the results of multiple sequence alignment, or the other kinds of sequence alignment, are not used in the frame of the method. Using our approach, we predicted the formation of three main secondary structure types (alpha-helix, beta-structure and coil) with high accuracy, that is Q(3) = 76%. Predicting the formation of alpha-helix and non-alpha-helix states we reached the accuracy which was measured as Q(2) = 86%. Also we have identified certain protein chain properties that affect the formation of membrane protein secondary structure. These protein properties include hydrophobic properties of amino acid residues, presence of Gly, Ala and Val amino acids, and the location of protein chain end.     

Secondary structure of human interleukin-3

1. The secondary structure of human interleuken-3 in solution was determined by circular dichroism spectroscopy. 2. The results were then compared with empirical secondary structure predictions based on primary sequence structure of the protein. 3. The two approaches are in extremently close agreement showing the protein to have 40% alpha-helix, 12% total beta-structure and 48% random coil content.     

Prediction of the Secondary Structure of Myelin Basic Protein

An investigation into the probable secondary structure of the myelin basic protein was carried out by the application of three procedures currently in use to predict the secondary structures of proteins from knowledge of their amino acid sequences. In order to increase the accuracy of the predictions, the amino acid substitutions that occur in the basic protein from different species were incorporated into the predictive algorithms. It was possible to locate regions of probable alpha-helix, beta-structure, beta-turn, and unordered conformation (coil) in the protein. One of the predictive methods introduces a bias into the algorithm to maximize or minimize the amounts of alpha-helix and/or beta-structure present; this made it possible to assess how conditions such as pH and protein concentration or the presence of anionic amphiphilic molecules could influence the protein's secondary structure. The predictions made by the three methods were in reasonably good agreement with one another. They were consistent with experimental data, provided that the stabilizing or destabilizing effects of the environment were taken into account. According to the predictions, the extent of possible alpha-helix and beta-structure formation in the protein s severely restricted by the low frequency and extensive scattering of hydrophobic residues, along with a high frequency and extensive scattering of residues that favor the formation of beta-turns and coils. Neither prolyl residues nor cationic residues per se are responsible for the low content of alpha-helix predicted in the protein. The principal ordered conformation predicted is the beta-turn. Many of the predicted beta-turns overlap extensively, involving in some cases up to 10 residues. In some of these structures it is possible for the peptide backbone to oscillate in a sinusoidal manner, generating a flat, pleated sheetlike structure. Cationic residues located in these structures would appear to be ideally oriented for interaction with lipid phosphate groups located at the cytoplasmic surface of the myelin membrane. An analysis of possible and probable conformations that the triproline sequence could assume questions the popular notion that this sequence produces a hairpin turn in the basic protein.     

On the Relation Between the Predicted Secondary Structure and the Protein Size

Accurately predicted protein secondary structure provides useful information for target selection, to analyze protein function and to predict higher dimensional structure. Existing research shows that more data + refined search = better prediction. We analyze relation between the prediction accuracy and another crucial factor, the protein size. Empirical tests performed with two secondary structure predictors on a large set of high-resolution, non-redundant proteins show that the average accuracies for small proteins (<100 residues) equal 73% and 54% for alpha-helices and beta-strands, respectively. The alpha-helix/beta-strand accuracies for very large proteins (>300 residues) equal 77%/68%, respectively. Similarly, the tests with three secondary structure content predictors show that the prediction errors for the small/very large proteins equal 0.13/0.09 and 0.09/0.06 for alpha-helix and beta-strand content, respectively. Our tests confirm that the secondary structure/content predictions for the very large proteins are characterized statistically significantly better quality than prediction for the small proteins. This is in contrast with the tertiary structure predictions in which higher accuracy is obtained for smaller proteins.     

Secondary structural changes of non-reduced and reduced ribonuclease A in solutions of urea, guanidine hydrochloride and sodium dodecyl sulfate

The secondary structures of ribonuclease A (RNAase A) before and after reduction of the disulfide bridges and blockage of the thiol groups with iodoacetamide were examined in solutions of urea, guanidine hydrochloride, and sodium dodecyl sulfate (SDS). The relative proportions of alpha-helix, beta-structure, and disordered structure were estimated by the curve-fitting method of circular dichroism (Chen, Y.H., Yang, J.T. and Chau, K.H. (1974) Biochemistry 13, 3350-3359). The native RNAase A, with the disulfide bridges intact, contained 19% helix and 38% beta-structure. Reduction of its disulfide bridges led to a decrease in the proportion of these structures to 9% for the alpha-helix and 17% for the beta-structure. The non-reduced RNAase A resisted unfolding in low concentrations of urea and guanidine hydrochloride. The beta-structure which remained after reduction appeared to be stable even in solutions of 6 M guanidine and 9 M urea. A considerable amount of the beta-structure in both the non-reduced and the reduced RNAase A remained unaffected by high concentrations of SDS.     

The synthesis and study of side-chain lactam-bridged peptides

Side-chain lactam bridges linking amino acid residues that are spaced several residues apart in the linear sequence offer a convenient and flexible method for introducing conformational constraints into a peptide structure. The availability of a variety of selectively cleavable protecting groups for amines and carboxylic acids allows for several approaches to the synthesis of monocyclic, dicyclic, and bicyclic lactam-bridged peptides by solid-phase methods. Multicyclic structures are also accessible, but segment-condensation syntheses with solution-phase cyclizations are most likely to provide the best synthetic approach to these more complex constrained peptides. Lactam bridges linking (i, i + 3)-, (i, i + 4), and (i, i + 7)-spaced residue pairs have all proven useful for stabilization of alpha helices, and (i, i + 3)-linked residues have also been demonstrated to stabilize beta-turns. These structures are finding an increasing number of applications in protein biology, including studies of protein folding, protein aggregation, peptide ligand-receptor recognition, and the development of more potent peptide therapeutics. Defining the functional roles of the amphiphilic alpha-helices in medium-sized peptide hormones, and studying helix propagation from rigid, alpha-helix initiating bicyclic peptides are among the most exciting developments currently underway in this field.     

Use of secondary structural information and C<Emphasis Type="Italic">α</Emphasis>-C<Emphasis Type="Italic">α</Emphasis> distance restraints to model protein structures with MODELLER

Protein secondary structure predictions and amino acid long range contact map predictions from primary sequence of proteins have been explored to aid in modelling protein tertiary structures. In order to evaluate the usefulness of secondary structure and 3D-residue contact prediction methods to model protein structures we have used the known Q3 (alpha-helix, beta-strands and irregular turns/loops) secondary structure information, along with residue-residue contact information as restraints for MODELLER. We present here results of our modelling studies on 30 best resolved single domain protein structures of varied lengths. The results shows that it is very difficult to obtain useful models even with 100% accurate secondary structure predictions and accurate residue contact predictions for up to 30% of residues in a sequence. The best models that we obtained for proteins of lengths 37, 70, 118, 136 and 193 amino acid residues are of RMSDs 4.17, 5.27, 9.12, 7.89 and 9.69, respectively. The results show that one can obtain better models for the proteins which have high percent of alpha-helix content. This analysis further shows that MODELLER restrain optimization program can be useful only if we have truly homologous structure(s) as a template where it derives numerous restraints, almost identical to the templates used. This analysis also clearly indicates that even if we satisfy several true residue-residue contact distances, up to 30% of their sequence length with fully known secondary structural information, we end up predicting model structures much distant from their corresponding native structures.     

DNA-induced secondary structure of the carboxyl-terminal domain of histone H1

We have studied the secondary structure of the carboxyl-terminal domains of linker histone H1 subtypes H1(0) (C-H1(0)) and H1t (C-H1t), free in solution and bound to DNA, by IR spectroscopy. The carboxyl-terminal domain has little structure in aqueous solution but becomes extensively folded upon interaction with DNA. The secondary structure elements present in the bound carboxyl-terminal domain include the alpha-helix, beta-structure, turns, and open loops. The structure of the bound domain shows a significant dependence on salt concentration. In low salt (10 mm NaCl), there is a residual amount of random coil, 7% in C-H1(0) and 12% in C-H1t. In physiological salt concentrations (140 mm NaCl), the carboxyl termini become fully structured. Under these conditions, C-H1(0) contained 24% alpha-helix, 25% beta-structure, 17% open loops, and 33% turns. The latter component could include a substantial proportion of the 3(10) helix. Despite their low sequence identity (approximately 30%), the representation of the different structural motifs in C-H1t was similar to that in C-H1(0). Examination of the changes in the amide I components in the 20-80 degrees C temperature interval showed that the secondary structure of the DNA-bound C-H1t is for the most part extremely stable. The H1 carboxyl-terminal domain appears to belong to the so-called disordered proteins, undergoing coupled binding and folding.     

Use of secondary structural information and C alpha-C alpha distance restraints to model protein structures with MODELLER

Protein secondary structure predictions and amino acid long range contact map predictions from primary sequence of proteins have been explored to aid in modelling protein tertiary structures. In order to evaluate the usefulness of secondary structure and 3D-residue contact prediction methods to model protein structures we have used the known Q3 (alpha-helix,beta-strands and irregular turns/loops) secondary structure information, along with residue-residue contact information as restraints for MODELLER. We present here results of our modelling studies on 30 best resolved single domain protein structures of varied lengths. The results shows that it is very difficult to obtain useful models even with 100% accurate secondary structure predictions and accurate residue contact predictions for up to 30% of residues in a sequence. The best models that we obtained for proteins of lengths 37, 70, 118, 136 and 193 amino acid residues are of RMSDs 4.17, 5.27, 9.12, 7.89 and 9.69,respectively. The results show that one can obtain better models for the proteins which have high percent of alpha-helix content. This analysis further shows that MODELLER restrain optimization program can be useful only if we have truly homologous structure(s) as a template where it derives numerous restraints, almost identical to the templates used. This analysis also clearly indicates that even if we satisfy several true residue-residue contact distances, up to 30%of their sequence length with fully known secondary structural information, we end up predicting model structures much distant from their corresponding native structures.     

Prediction of Protein Secondary Structure Using Feature Selection and Analysis Approach

The prediction of the secondary structure of a protein from its amino acid sequence is an important step towards the prediction of its three-dimensional structure. However, the accuracy of ab initio secondary structure prediction from sequence is about 80 % currently, which is still far from satisfactory. In this study, we proposed a novel method that uses binomial distribution to optimize tetrapeptide structural words and increment of diversity with quadratic discriminant to perform prediction for protein three-state secondary structure. A benchmark dataset including 2,640 proteins with sequence identity of less than 25 % was used to train and test the proposed method. The results indicate that overall accuracy of 87.8 % was achieved in secondary structure prediction by using ten-fold cross-validation. Moreover, the accuracy of predicted secondary structures ranges from 84 to 89 % at the level of residue. These results suggest that the feature selection technique can detect the optimized tetrapeptide structural words which affect the accuracy of predicted secondary structures.     

Structures of elicitin isoforms secreted by Phytophthora drechsleri.

Most of the phytopathogenic fungi Phytophthora secrete holoproteins (elicitins) responsible for the incompatible reaction and systemic leaf necroses on tobacco. We found that Phytophthora drechsleri produces several elicitin isoforms of various toxicity on tobacco. The CD spectra showed that their secondary structure was largely conserved, exhibiting ca 50% alpha-helix and little or no beta-structure. These 98 residue proteins were sequenced and compared with other known elicitins. Only one point mutation correlated with the differences in necrotic activities. This residue could be either an active or a regulatory site, involved in the interaction with a receptor responsible for necrosis induction.     

Limits on alpha-helix prediction with neural network models.

Using a backpropagation neural network model we have found a limit for secondary structure prediction from local sequence. By including only sequences from whole alpha-helix and non-alpha-helix structures in our training and test sets--sequences spanning boundaries between these two structures were excluded--it was possible to investigate directly the relationship between sequence and structure for alpha-helix. A group of non-alpha-helix sequences, that was disrupting overall prediction success, was indistinguishable to the network from alpha-helix sequences. These sequences were found to occur at regions adjacent to the termini of alpha-helices with statistical significance, suggesting that potentially longer alpha-helices are disrupted by global constraints. Some of these regions spanned more than 20 residues. On these whole structure sequences, 10 residues in length, a comparatively high prediction success of 78% with a correlation coefficient of 0.52 was achieved. In addition, the structure of the input space, the distribution of beta-sheet in this space, and the effect of segment length were also investigated.     

Purification of homogeneous gamma-cystathionase and study of its structure by circular dichroism

Rat liver gamma-cystathionase has been purified to homogeneity (verified by SDS electrophoresis and ultracentrifugation). The secondary and tertiary structures of the enzyme were studied by circular dichroism spectra. Our studies revealed that the holoenzyme molecule comprises approximately 22% of alpha-helices, 14% of beta-structure, 14% of beta-bends, and 50% of unordered structure. Conformational alterations of the enzyme molecule resulting from enzyme PLP elimination, reduction with sodium borohydride and irreversible inhibition by propargylglycine were examined. The enzyme's secondary structure was shown to be stable whereas the tertiary structure is labile. Saturation with PLP maintains the enzyme's optimal (catalytically active) tridimensional structure. Sodium dodecylsulfate alters its secondary (the amount of alpha-helix being raised to 34%) and tertiary structures.     

A comparative CD study of carbonic anhydrase isoenzymes with different number of tryptophans: impact on calculation of secondary structure content.

The CD spectra of human carbonic anhydrase I and II and bovine carbonic anhydrase III were recorded and analyzed. The 3D structures of these isoenzymes are known, showing very similar secondary structure and polypeptide-chain fold. The tryptophan content, however, differs between the isoenzymes, i.e., isoenzymes I, II, and III possess 6, 7, and 8 tryptophans, respectively. All of the tryptophans except the additional tryptophans in isoenzymes II and III, i.e., W245 and W47, are conserved. Despite the fact that X-ray structure determinations showed that the isoenzymes had highly similar secondary structure, the contents of alpha-helix and beta-sheet structure differed considerably when using different CD algorithms for estimation of the fractions of various secondary structural elements. This shows that aromatic amino acids also interfere in the wavelength region (far-UV) used to calculate the amount of secondary structure. Such interference is especially problematic when analyzing proteins like carbonic anhydrase, which consist mainly of beta-structure that gives rise to weak ellipticity bands, compared to the bands arising from alpha-helical structure.     

Studies on the role and mode of operation of the very-lysine-rich histones in eukaryote chromatin. The conformation of phi1 histones from marine invertebrate sperm.

Proton magnetic resonance, circular dichroism and infrared spectroscopy are used to investigate the secondary and tertiary structure of three very lysine-rich histones from marine invertebrate sperm. At high ionic strength both Arbacia lixula and Holothuria tubulosa histone phi 1 are observed to contain 25-30% alpha-helix, no beta-structure and to form specific folded structures. Both phi 1 proton magnetic resonance spectra have perturbed methyl resonances at chemical shifts close to those observed for calf thymus H1, suggesting analogies in tertiary structure. Mytilus edulis histone phi 1 however, shows no spectroscopic evidence of secondary and tertiary structure on salt addition.     

Significance of aromatic-backbone amide interactions in protein structure

Weakly polar interactions between aromatic rings of amino acids and hydrogens of backbone amides (Ar-HN) have been shown to support local structures in proteins. Their role in secondary structures, however, has not been elucidated. To investigate the relationship between Ar-HN interaction and the stability of local and secondary structures of polypeptides and to improve the prediction of this interaction based on amino acid sequence, the structures of 560 nonhomologous proteins, from the Protein Data Bank, were searched for Ar-HN interactions between the aromatic ring of each Phe, Tyr, and Trp residue at position i and the backbone amide group of any residue, except Pro, at the positions i, i - 1, i - 2, i - 3, i + 1, i + 2, and i + 3. Ar-HN interactions were identified by calculating the chemical shift of the amide hydrogen caused by the proximal aromatic ring. Ar(i)-HN(i + 1, i + 2 and i + 3) interactions were more common (7.10%, 2.08%, and 0.54%, respectively) than were Ar(i)-HN(i - 1, i - 2, and i - 3) interactions (0.66%, <0.1%, and 0.18%, respectively). The value of the chi(1) torsion angle of the aromatic residue in position i depended on the direction of the Ar-HN interaction. The position of the aromatic ring in Ar(i)-HN(i + 1, i + 2, and i + 3) interactions was mostly trans, in Ar(i)-HN(i - 1, i - 2, and i - 3) interactions mainly gauche(-), and in Ar(i)-HN(i) interactions mostly gauche(+). The analyses of the secondary structures of the protein fragments containing Ar-HN interactions showed that Ar-HN interactions were in all types of secondary structures. Search results suggest that Ar-HN interactions have a stabilizing effect on all types of secondary structures.