Secondary structure predictions and medium range interactions

Several authors have proposed that predictions of protein secondary structure derived from statistical information about the known structures can be improved when information about neighboring residues participating in short and medium range interactions is included. A substantial improvement shown here indicates that current methods of including this information are not more successful than methods that do not. Evaluations of the Chou and Fasman method (Adv. Enzymol. 47 (1978) 45-148), that does not include information about interactions (except in averaging), have shown it to be about 49% correct for three states (helix, beta-sheet and undefined). In comparison, the method of Garnier et al. (J. Mol. Biol. 120 (1978) 97-120), that explicitly includes information about neighboring residues, has an accuracy of 57% residues correct for three states. However, we have obtained an 8% improvement for predictions of secondary structure based on the algorithm by Chou and Fasman. The improvements are obtained by eliminating many rules and by choosing the best decision constants for structure assignments. The simplified method described here is 57% correct for three states using preference values calculated in 1978.     

Predicted secondary structures of four penicillin beta-lactamases and a comparison with two lysozymes.

We have predicted the secondary structures of four beta-lactamases (Bacillus cereus, Bacillus licheniformis, Staphylococcus aureus, and Escherichia coli R-TEM) by the statistical method of Chou & Fasman as well as by the information theory method of Garnier et al. The secondary structures of all four beta-lactamases are of the alpha/beta type (Levitt & Chothia's nomenclature), with helices at N- and C-termini. There are about eight short regions each of alpha-helical (30--50%) and beta-strand (10--20%) structure separated by about 20 reverse turns. The conformation of the Gram-positive and Gram-negative beta-lactamases are generally similar although a few differences are predicted between the S.aureus and E.coli structures. Surprisingly, the two bacilli structures differ significantly in three short regions. In all four enzymes the region near the catalytically-implicated tyrosine has similar secondary structure. The secondary structure of hen egg white lysozyme, a penicillin-binding enzyme, as well as T4 phage lysozyme, has similarities to the N-terminal half of the penicillin-destroying beta-lactamases.     

A simple method for predicting the secondary structure of globular proteins: implications and accuracy

A method is presented for predicting the secondary structureof globular proteins from their amino acid sequence. It is basedon a rigorous statistical exploitation of the well-known biologicalfact that the amino acid compositions of each secondary structureare different. We also propose an evaluation process that allowsus to estimate the capacity of a method to predict the secondarystructure of a new protein which does not have any homologousproteins whose structure is already known. This evaluation processshows that our method has a prediction accuracy of 58.7% overthree states for the 62 proteins of the Kabsch and Sander (1983a)data bank. This result is better than that obtained by the mostwidely used methods—Lim (1974), Chou and Fasman (1978)and Garnier et al. (1978)—and also than that obtainedby a recent method based on local homologies (Levin et al.,1986). Our prediction method is very simple and may be implementedon any microcomputer and even on programmable pocket calculators.A simple Pascal implementation of the method prediction algorithmis given. The interpretation of our results in terms of proteinfolding and directions for further work are discussed. Received on December 15, 1987; accepted on April 12, 1988     

Physicochemical characterization of bovine retinal arrestin

The native conformation of bovine retinal arrestin has been characterized by a variety of spectroscopic methods. The purified protein gives rise to a near uv absorption band centered at 279 nm which results from the absorbance of its 14 tyrosine and one tryptophan residue. The extinction coefficient for this absorption band was determined to be 38.64 mM-1, cm-1 using the tyrosinate-tyrosine difference spectrum method; this extinction coefficient is ca. 17% lower than the previously reported value, and provides estimates of protein concentration which are in good agreement with estimates from the Bradford colorimetric assay. When native arrestin is purified to homogeneity, it displays a fluorescence spectrum which is dominated by tyrosine emission with no discernible contribution from tryptophan. Observation of the tyrosine-like fluorescence is dependent on the purity and structural integrity of the protein. Denaturation of arrestin by guanidine hydrochloride results in a diminution of tyrosine fluorescence and the concomitant appearance of a second fluorescence maximum at ca. 340 nm, presumably due to the single tryptophan residue. Thermal denaturation of arrestin leads to a conformation characterized by a broad fluorescence band centered at ca. 325 nm. Study of the arrestin fluorescence spectrum as a function of temperature indicates that the thermal denaturation is well modeled as a two-state transition with a transition midpoint of 60 degrees C. Temperature-dependent far uv circular dichroism studies indicate that changes in secondary structure occur coincident with the change in fluorescence. Studies of the temperature dependence of arrestin binding to light-adapted phosphorylated rhodopsin shows a strong correlation between the fluorescence spectral features of arrestin and its ability to bind rhodopsin. These data suggest that the relative intensities of tyrosine and tryptophan fluorescence are sensitive to the structural integrity of the native (i.e., rhodopsin binding) state of arrestin, and can thus serve as useful markers of conformational transitions of this protein. The lack of tryptophan fluorescence for native arrestin suggests an unusual environment for this residue. Possible mechanisms for this tryptophan fluorescence quenching are discussed.     

Structural analysis of staphylococcal enterotoxins B and C1 using circular dichroism and fluorescence spectroscopy

Secondary and tertiary structural parameters of two functionally and serologically related proteins, staphylococcal enterotoxins B and C1, have been determined by using circular dichroism and fluorescence spectroscopy. The secondary structures derived from the respective far-UV circular dichroic spectra were 9.5% alpha-helix, 55.0% beta-pleated sheets, 16.5% beta-turns, and 19.0% random coils for enterotoxin B and 15.0% alpha-helix, 38.0% beta-pleated sheets, 25.5% beta-turns, and 21.5% random coils for staphylococcal enterotoxin C1. The values matched well with the secondary structures derived from the amino acid sequences (Chou and Fasman method). Seven antigenic sites have been predicted for both staphylococcal enterotoxins B and C1 by using the hydrophilicity and the secondary structure information. Three of these antigenic sites appear similar. Fluorescence quantum yield of the single tryptophan residue (Trp-197) of both the enterotoxins showed the tryptophan residue in staphylococcal enterotoxin B to be approximately 46% more fluorescent than in staphylococcal enterotoxin C1. Tryptophan fluorescence quenching by the surface quencher I- and the neutral quencher acrylamide revealed that the single tryptophan residue in each of the enterotoxins is buried in the protein matrix and is not accessible to the surface quencher I-. The tryptophan residue in staphylococcal enterotoxin C1 is 14% less accessible to acrylamide than in staphylococcal enterotoxin B. The data, in general, reflect several similarities and significant differences between the two related enterotoxins.     

Protein structure prediction

Current methods developed for predicting protein structure are reviewed. The most widely used algorithms of Chou and Fasman and Garnier et al for predicting secondary structure are compared to the most recent ones including sequence similarity methods, neural network, pattern recognition or joint prediction methods. The best of these methods correctly predict 63-65% of the residues in the database with cross-validation for 3 conformations, helix, beta strand and coli with a standard deviation of 6-8% per protein. However, when a homologous protein is already in the database, the accuracy of prediction by the similarity peptide method of Levin and Garnier reaches about 90%. Some conclusions can be drawn on the mechanism of protein folding. As all the prediction methods only use the local sequence for prediction (+/- 8 residues maximum) one can infer that 65% of the conformation of a residue is dictated on average by the local sequence, the rest is brought by the folding. The best predicted proteins or peptide segments are those for which the folding has less effect on the conformation. Presently, prediction of tertiary structure is only of practical use when the structure of a homologous protein is already known. Amino acid alignment to define residues of equivalent spatial position is critical for modelling of the protein. We showed for serine proteases that secondary structure prediction can help to define a better alignment. Non-homologous segments of the polypeptide chain, such as loops, libraries of known loops and/or energy minimization with various force fields, are used without yet giving satisfactory solutions. An example of modelling by homology, aided by secondary structure prediction on 2 regulatory proteins, Fnr and FixK is presented.     

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.     

Secondary structure of the variant surface glycoproteins of trypanosomes

The secondary structure of seven variant surface glycoproteins (VSGs) of trypanosomes has been determined by Raman spectroscopy. They are all predominantly alpha-helical, the alpha-helix content varying between 50 and 60%. The beta-strand content varies between 20 and 25%, and the content of beta-turn and nonregular structures is about 25%. For three VSGs the N-terminal domain obtained by proteolytic cleavage was found to have essentially the same secondary structure as the complete VSGs. For three VSGs a secondary structure prediction has been performed applying the rules of Chou and Fasman. In all cases, two long alpha-helices extending over about 50 residues or 80 A are predicted in agreement with the X-ray diffraction data of Freymann et al. [(1984) Nature 311, 167-169] and Metcalf et al. [(1987) Nature 325, 84-86]. The region between the two alpha-helical segments exhibits a high potential of beta-turns, suggesting that this segment may be exposed on the cell surface and carry major antigenic determinants.     

Ligand exchange during unfolding of cytochrome c.

The productive folding pathway of cytochrome c passes through an obligatory HW intermediate in which the heme is coordinated by a solvent water molecule and a native ligand, His-18, prior to the formation of the folded HM state with both the native His-18 and Met-80 heme coordination. Two off pathway intermediates, a five-coordinated state (5C) and a bis-histidine state (HH), were also identified during the folding reaction. In the present work, the thermodynamics and the kinetics of the unfolding reaction of cytochrome c were investigated with resonance Raman scattering, tryptophan fluorescence spectroscopy, and circular dichroism. The objective of these experiments was to determine if the protein opens up and diverges into the differing heme ligation states through a many pathway mechanism or if it passes through intermediate states analogous to those observed during the folding reaction. Equilibrium unfolding results indicate that, in contrast to 5C, the stability of HH with respect to HW decreases as the concentration of GdnHCl increases. The difference in their response to the denaturant indicates that the polypeptide structure of 5C is relatively loose as compared with HH in which the polypeptide is misfolded. Time-resolved resonance Raman measurements show that strikingly similar ligand exchange reactions occur during unfolding as were observed during folding. Combined with fluorescence data, a kinetic model is proposed in which local structural rearrangements controlled by heme ligand exchange reactions appear prior to the global relaxation of the polypeptide chain.     

A spectroscopic and conformational study of pertussis toxin

The conformation of native pertussis toxin has been investigated by secondary structure prediction and by circular dichroism, fluorescence and second-derivative ultraviolet absorption spectroscopy. The far-ultraviolet circular dichroic spectrum is characteristic of a protein of high beta-sheet and low alpha-helix content. This is also shown by an analysis of the circular dichroic spectrum with the Contin programme which indicates that the toxin possesses 53% beta-sheet, 10% alpha-helix and 37% beta-turn/loop secondary structure. Second-derivative ultraviolet absorption spectroscopy suggests that 34 tyrosine residues are solvent-exposed and quenching of tryptophan fluorescence emission has shown that 4 tryptophan residues are accessible to iodide ions. One of these tryptophans appears to be in close proximity to a positively charged side-chain, since only 3 tryptophans are accessible to caesium ion fluorescence quenching. When excited at 280 nm, the emission spectrum contains a significant contribution from tyrosine fluorescence, which may be a consequence of the high proportion (55%) of surface-exposed tyrosines. No changes in the circular dichroic spectra of the toxin were found in the presence of the substrate NAD. However, NAD did quench both tyrosine and tryptophan fluorescence emission but did not change the shape of the emission spectrum, or the accessibility of the tryptophans to either the ionic fluorescence quenchers or the neutral quencher acrylamide.     

Predicted secondary structure of calcitonin in relation to the biological activity.

The secondary structure of the natural analogues of calcitonin have been predicted using the methods of Burgess et al., Chou and Fasman, and Lim. The predicted structures were similar in the N-terminal part of the chain, but were variable in other regions of the molecules. The higher hypocalcaemic potency of ultimobranchial calcitonins as compared to thyroidal calcitonins may be related to the higher potential for helical structure in the former.     

p53 unfolding detected by CD but not by tryptophan fluorescence

Full-length human p53 protein was examined using tryptophan fluorescence and circular dichroism spectroscopy (CD) to monitor unfolding. No significant alteration in tryptophan fluorescence for the tetrameric protein was detectable over a wide range of either urea or guanidine hydrochloride concentrations, in contrast to results with the isolated DNA binding domain [Bullock et al. (1997) Proc. Natl. Acad. Sci. USA 94, 14338]. Under similar denaturant conditions, CD demonstrated significant protein unfolding for the full-length wild-type protein, with increased apparent structure loss compared to that detected during thermal denaturation [Nichols and Matthews (2001) Biochemistry 40, 3847]. Examination of X-ray structures containing two of the four tryptophan residues of a p53 monomer suggested local environments consistent with quenched fluorophores. Exploration of p53 fluorescence using potassium iodide as a quencher confirmed that these fluorophores are already substantially quenched in the native structure, and this quenching is not relieved during protein unfolding.     

Primary structure of Stoichactis helianthus cytolysin III

The primary structure of Stoichactis helianthus cytolysin III has been determined by automated Edman degradation of the intact protein and of peptides derived therefrom by hydrolysis with trypsin and staphylococcal protease and by chemical cleavage with cyanogen bromide and o-iodosobenzoic acid. As a result of these studies, the positions of all 153 amino acid residues of toxin III have been unambiguously determined. Most regions of sequence were determined two times in different types of digests of the protein. A number of highly hydrophobic regions of sequence, which may be functionally significant, have been identified, including a region rich in tyrosine and tryptophan (residues 86-98). The secondary structure of toxin III has been predicted by Chou-Fasman analysis (Chou, P.Y., and Fasman, G.D. (1978) Annu. Rev. Biochem. 47, 251-276) of the primary structure. The predicted secondary structure contains 16% alpha-helix and 31% beta-structure.     

The prediction of the conformation of membrane proteins from the sequence of amino acids.

The methods of Chou & Fasman [Biochemistry (1974) 13, 211-222, 222-245] and of Lim [J. Mol. Biol. (1974)88, 857-872, 873-894] for predicting secondary structure from amino acid sequence have been applied to five predominantly helical membrane-associated peptides. The predictions from the method of Lim (1974a,b) are consistent with the experimental observations, whereas those from Chou & Fasman (1974a,b), although not inconsistent with alpha-helix, favour a beta-structure for several very hydrophobic regions. The results may be rationalized in terms of the effect of the solvent on the conformation of a polypeptide.     

Tryptophan replacements in the trp aporepressor from Escherichia coli: probing the equilibrium and kinetic folding models.

Mutants of the dimeric Escherichia coli trp aporepressor are constructed by replacement of the two tryptophan residues in each subunit in order to assess the effects on equilibrium and kinetic fluorescence properties of the folding reaction. The three kinetic phases detected by intrinsic tryptophan fluorescence in refolding of the wild-type aporepressor are also observed in folding of both Trp 19 to Phe and Trp 99 to Phe single mutants, demonstrating that these phases correspond to global rather than local conformational changes. Comparison of equilibrium fluorescence (Royer, C.A., Mann, C.J., & Matthews, C.R., 1993, Protein Sci. 2, 1844-1852) and circular dichroism transition curves induced by urea shows that replacement of either Trp 19 or Trp 99 results in noncoincident behavior. Unlike the wild-type protein (Gittelman, M.S. & Matthews, C.R., 1990, Biochemistry 29, 7011-7020), tertiary and/or quaternary structures are disrupted at lower denaturant concentration than is secondary structure. The equilibrium results can be interpreted in terms of enhancement in the population of a monomeric folding intermediate in which the lone tryptophan residue is highly exposed to solvent, but in which substantial secondary structure is retained. The location of both mutations at the interface between the two subunits (Zhang, R.G., et al., 1987, Nature 327, 591-597) provides a simple explanation for this phenomenon.     

PEPPLOT, a protein secondary structure analysis program for the UWGCG sequence analysis software package.

We describe a program for the analysis of protein secondary structure that operates with the Sequence Analysis Software Package of the University of Wisconsin Genetics Computer Group (UWGCG). The program produces both graphic and printed output. Structure prediction using the Chou and Fasman and Robson et al methods, and hydropathy analysis by the method of Kyte and Doolittle are included along with a simplified method of hydrophobic moment analysis. The power of the program is the coordinated presentation of many different kinds of structural information on the same plot.     

The complete amino acid sequence of echinoidin, a lectin from the coelomic fluid of the sea urchin Anthocidaris crassispina. Homologies with mammalian and insect lectins

The complete amino acid sequence of echinoidin, the proposed name for a lectin from the coelomic fluid of the sea urchin Anthocidaris crassispina, has been determined by sequencing the peptides obtained from tryptic, Staphylococcus aureus V8 protease, chymotryptic, and thermolysin digestions. Echinoidin is a multimeric protein (Giga, Y., Sutoh, K., and Ikai, A. (1985) Biochemistry 24, 4461-4467) whose subunit consists of a total of 147 amino acid residues and one carbohydrate chain attached to Ser38. The molecular weight of the polypeptide without carbohydrate was calculated to be 16,671. Each polypeptide chain contains seven half-cystines, and six of them form three disulfide bonds in the single polypeptide chain (Cys3-Cys14, Cys31-Cys141, and Cys116-Cys132), while Cys2 is involved in an interpolypeptide disulfide linkage. From secondary structure prediction by the method of Chou and Fasman (Chou, P. Y., and Fasman, G. D. (1974) Biochemistry 13, 211-222) the protein appears to be rich in beta-sheet and beta-turn structures and poor in alpha-helical structure. The sequence of the COOH-terminal half of echinoidin is highly homologous to those of the COOH-terminal carbohydrate recognition portions of rat liver mannose-binding protein and several other hepatic lectins. This COOH-terminal region of echinoidin is also homologous to the central portion of the lectin from the flesh fly Sarcophaga peregrina. Moreover, echinoidin contains an Arg-Gly-Asp sequence which has been proposed to be a basic functional unit in cellular recognition proteins.     

Nanosecond dynamics of the single tryptophan reveals multi-state equilibrium unfolding of protein GB1

Unfolding of the immunoglobulin binding domain B1 of streptococcal protein G (GB1) was induced by guanidine hydrochloride (GdnHCl) and studied by circular dichroism, steady-state, and time-resolved fluorescence spectroscopy. The fluorescence methods employed the single tryptophan residue of GB1 as an intrinsic reporter. While the transitions monitored by circular dichroism and steady-state fluorescence coincided with each other, the transitions followed by dynamic fluorescence were markedly different. Specifically, fluorescence anisotropy data showed that a relaxation spectrum of tryptophan contained a slow motion with relaxation times of 9 ns in the native state and 4 ns in the unfolded state in 6 M GdnHCl. At intermediate GdnHCl concentrations of 3.8-4.2 M, however, the slow relaxation time increased to 18 ns. The fast nanosecond motion had an average time of 0.8 ns and showed no dependence on the formation of native structure. Overall, dynamic fluorescence revealed two preliminary stages in GB1 folding, which are equated with the formation of local structure in the beta(3)-strand hairpin and the initial collapse. Both stages exist without alpha-helix formation, i. e., before the appearance of any ordered secondary structure detectable by circular dichroism. Another stage in GB1 folding might exist at very low ( approximately 1 M) GdnHCl concentrations.     

Conformational aspects of the acid-induced fusion mechanism of influenza virus hemagglutinin. Circular dichroism and fluorescence studies

Circular dichroism and tryptophan fluorescence spectroscopy have been used to investigate the structures of the influenza virus membrane glycoprotein hemagglutinin, acid-treated hemagglutinin, and fragments of hemagglutinin derived by proteolysis. The conformational change in hemagglutinin which occurs at the pH of membrane fusion (pH 5-6) was associated with a significant change of the environment of tyrosine residues, a change in the environment of tryptophan residues, but no changes in secondary structure. Tryptic digestion of the hemagglutinin in its low pH conformation which releases one of the subunit polypeptides (HA1) caused minimal changes in tyrosine and tryptophan environments but a small secondary structural change in HA1. The secondary structure of the remainder of the molecule (HA2) was very similar to that predicted from the known x-ray crystallographic structure of the native molecule. However, fluorescence spectroscopy indicated a tertiary change in structure in the coiled coil of alpha-helices which form the fibrous central stem of the molecule. These results are consistent with a conformational change required for membrane fusion which involves a decrease of HA1/HA1, HA1/HA2 interactions and changes in tertiary structure not accompanied by changes in secondary structure.     

Structural characterization of myo-inositol monophosphatase from bovine brain by secondary structure prediction,fluorescence, circular dichroism and Raman spectroscopy

Structural aspects of myo-inositol monophosphatase were examined by spectroscopic techniques and empirical prediction methods. The enzyme belongs to the α/β class of proteins, with approx. 33% α-helix and 29% β-sheet, as shown by circular dichroism (CD), Raman spectroscopy and prediction based on the amino-acid sequence. The Raman spectrum also suggests that the three tryptophan residues in myo-inositol monophosphatase are not expose to solvent. This was confirmed by a blue shift of 25 nm in the fluorescence emission spectrum, as compared to tryptophan in water, and by quenching studies with acrylamide. The enzyme shows a transition temperature of 87°C for the CD signal at 222 nm. This remarkable heat stability is not due to the presence of disulfide bonds, since both the Raman spectrum and chemical modification studies clearly indicate that all six cysteine residues are in the reduced state.