Assessing secondary structure assignment of protein structures by using pairwise sequence-alignment benchmarks

How to make an objective assignment of secondary structures based on a protein structure is an unsolved problem. Defining the boundaries between helix, sheet, and coil structures is arbitrary, and commonly accepted standard assignments do not exist. Here, we propose a criterion that assesses secondary structure assignment based on the similarity of the secondary structures assigned to pairwise sequence-alignment benchmarks, where these benchmarks are determined by prior structural alignments of the protein pairs. This criterion is used to rank six secondary structure assignment methods: STRIDE, DSSP, SECSTR, KAKSI, P-SEA, and SEGNO with three established sequence-alignment benchmarks (PREFAB, SABmark, and SALIGN). STRIDE and KAKSI achieve comparable success rates in assigning the same secondary structure elements to structurally aligned residues in the three benchmarks. Their success rates are between 1-4% higher than those of the other four methods. The consensus of STRIDE, KAKSI, SECSTR, and P-SEA, called SKSP, improves assignments over the best single method in each benchmark by an additional 1%. These results support the usefulness of the sequence-alignment benchmarks as a means to evaluate secondary structure assignment. The SKSP server and the benchmarks can be accessed at http://sparks.informatics.iupui.edu     

Estimation of the number of alpha-helical and beta-strand segments in proteins using circular dichroism spectroscopy

A simple approach to estimate the number of alpha-helical and beta-strand segments from protein circular dichroism spectra is described. The alpha-helix and beta-sheet conformations in globular protein structures, assigned by DSSP and STRIDE algorithms, were divided into regular and distorted fractions by considering a certain number of terminal residues in a given alpha-helix or beta-strand segment to be distorted. The resulting secondary structure fractions for 29 reference proteins were used in the analyses of circular dichroism spectra by the SELCON method. From the performance indices of the analyses, we determined that, on an average, four residues per alpha-helix and two residues per beta-strand may be considered distorted in proteins. The number of alpha-helical and beta-strand segments and their average length in a given protein were estimated from the fraction of distorted alpha-helix and beta-strand conformations determined from the analysis of circular dichroism spectra. The statistical test for the reference protein set shows the high reliability of such a classification of protein secondary structure. The method was used to analyze the circular dichroism spectra of four additional proteins and the predicted structural characteristics agree with the crystal structure data.     

Folding type-specific secondary structure propensities of amino acids,derived from α-helical, β-sheet, α/β, and α+β proteins of known structures

Folding type-specific secondary structure propensities of 20 naturally occurring amino acids have been derived from α-helical, β-sheet, α/β, and α+β proteins of known structures. These data show that each residue type of amino acids has intrinsic propensities in different regions of secondary structures for different folding types of proteins. Each of the folding types shows markedly different rank ordering, indicating folding type-specific effects on the secondary structure propensities of amino acids. Rigorous statistical tests have been made to validate the folding type-specific effects. It should be noted that α and β proteins have relatively small α-helices and β-strands forming propensities respectively compared with those of α+β and α/β proteins. This may suggest that, with more complex architectures than α and β proteins, α+β and α/β proteins require larger propensities to distinguish from interacting α-helices and β-strands. Our finding of folding type-specific secondary structure propensities suggests that sequence space accessible to each folding type may have differing features. Differing sequence space features might be constrained by topological requirement for each of the folding types. Almost all strong β-sheet forming residues are hydrophobic in character regardless of folding types, thus suggesting the hydrophobicities of side chains as a key determinant of β-sheet structures. In contrast, conformational entropy of side chains is a major determinant of the helical propensities of amino acids, although other interactions such as hydrophobicities and charged interactions cannot be neglected. These results will be helpful to protein design, class-based secondary structure prediction, and protein folding. © 1998 John Wiley & Sons, Inc. Biopoly 45: 35–49, 1998     

DSSPcont: Continuous secondary structure assignments for proteins

The DSSP program automatically assigns the secondary structure for each residue from the three-dimensional co-ordinates of a protein structure to one of eight states. However, discrete assignments are incomplete in that they cannot capture the continuum of thermal fluctuations. Therefore, DSSPcont (http://cubic.bioc.columbia.edu/services/DSSPcont) introduces a continuous assignment of secondary structure that replaces 'static' by 'dynamic' states. Technically, the continuum results from calculating weighted averages over 10 discrete DSSP assignments with different hydrogen bond thresholds. A DSSPcont assignment for a particular residue is a percentage likelihood of eight secondary structure states, derived from a weighted average of the ten DSSP assignments. The continuous assignments have two important features: (i) they reflect the structural variations due to thermal fluctuations as detected by NMR spectroscopy; and (ii) they reproduce the structural variation between many NMR models from one single model. Therefore, functionally important variation can be extracted from a single X-ray structure using the continuous assignment procedure.     

Continuum secondary structure captures protein flexibility

The DSSP program assigns protein secondary structure to one of eight states. This discrete assignment cannot describe the continuum of thermal fluctuations. Hence, a continuous assignment is proposed. Technically, the continuum results from averaging over ten discrete DSSP assignments with different hydrogen bond thresholds. The final continuous assignment for a single NMR model successfully reflected the structural variations observed between all NMR models in the ensemble. The structural variations between NMR models were verified to correlate with thermal motion; these variations were captured by the continuous assignments. Because the continuous assignment reproduces the structural variation between many NMR models from one single model, functionally important variation can be extracted from a single X-ray structure. Thus, continuous assignments of secondary structure may affect future protein structure analysis, comparison, and prediction.     

Molecular dynamics simulations of beta-turn forming tetra- and hexapeptides

It was previously shown that the structural ensemble of model peptides DDKG and GKDG (H. Ishii et al. Biopolymers 24, 2045-2056, 1985), DEKS (A. Otter et al. J. Biomol. Struct. Dyn. 7, 455-476, 1989) NPGQ (F. R. Carbone et al. Int. J. Pept. Protein. Res. 26, 498-508, 1985), SALN (H. Santa et al. J. Biomol. Struct. Dyn. 16, 1033-1041, 1999), SYPFDV and SYPYDV (J. Yao et al. J. Mol. Biol. 243, 736-753, 1994), VP(D)AH and VP(D)SH (B. Imperiali et al. J. Am. Chem. Soc. 114, 3182-3188, 1992) in solution contains a significant - or in some cases dominant - proportion of beta-turn conformation. In this study, a protein database was searched for the above, unprotected sequences which incorporate only L-amino acid residues. Simulated annealing and 25 ns MD simulations of structures were also performed. The DSSP and STRIDE secondary structure-assigning algorithms and clustering were used to analyze trajectories and i, i+3 hydrogen bonds were also sought. The DSSP analysis showed a fluctuation between beta-turn and random meander structure, although bend structures were not detected because of the insufficient length of peptide chains. This alternating trend was confirmed when the STRIDE algorithm was used to analyze trajectories, but STRIDE assigned more turn structures. The population of the strongest clusters was above 40% and the middle structures adopted beta-turn structure for most sequences. These results are in good agreement with previous experimental results and support the idea of the ultra-marginal stability of turns in the absence of stabilizing long-range interactions of the neighboring segments of a polypeptide chain. However, interactions between the side-chains in tetrapeptides could also contribute to turn stability and result in unusual stability in some cases. Our observations suggest that such interactions are the consequence rather than the driving force of turn formation.     

Molecular Dynamics Simulations of β-turn Forming Tetra- and Hexapeptides

Abstract

It was previously shown that the structural ensemble of model peptides DDKG and GKDG (H. Ishii et al. Biopolymers 24, 2045–2056, 1985), DEKS (A. Otter et al. J. Biomol. Struct. Dyn. 7, 455–476, 1989) NPGQ (F. R. Carbone et al. Int. J. Pept. Protein. Res. 26, 498–508, 1985), SALN (H. Santa et al. J. Biomol. Struct. Dyn. 16, 1033–1041, 1999), SYPFDV and SYPYDV (J. Yao et al. J. Mol. Biol. 243, 736–753, 1994), VP^DAH and VP^DSH (B. Imperiali et al. J. Am. Chem. Soc. 114, 3182–3188, 1992) in solution contains a significant—or in some cases dominant—proportion of β-turn conformation. In this study, a protein database was searched for the above, unprotected sequences which incorporate only L-amino acid residues. Simulated annealing and 25 ns MD simulations of structures were also performed. The DSSP and STRIDE secondary structure-assigning algorithms and clustering were used to analyze trajectories and i, i+3 hydrogen bonds were also sought. The DSSP analysis showed a fluctuation between β-turn and random meander structure, although bend structures were not detected because of the insufficient length of peptide chains. This alternating trend was confirmed when the STRIDE algorithm was used to analyze trajectories, but STRIDE assigned more turn structures. The population of the strongest clusters was above 40% and the middle structures adopted β-turn structure for most sequences. These results are in good agreement with previous experimental results and support the idea of the ultra-marginal stability of turns in the absence of stabilizing long-range interactions of the neighboring segments of a polypeptide chain. However, interactions between the side-chains in tetrapeptides could also contribute to turn stability and result in unusual stability in some cases. Our observations suggest that such interactions are the consequence rather than the driving force of turn formation.     

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.     

Chemical cross-linking between lysine groups in vimentin oligomers is dependent on local peptide conformations

Previous studies have shown that cytoplasmic intermediate filaments, other than the keratins, are each constructed from a single type of polypeptide chain. Studies involving chemical crosslinking between lysine groups have shown that assembly of the filaments begins with the formation of dimers in which the peptide chains are parallel and in exact register, and that these dimers further associate in antiparallel patterns having specific degrees of overlap. In the present study, molecular modeling of the conformations of vimentin molecules indicated that lysine side chains in identical positions in regions of α-helix in parallel chains might be unable to be linked because they are on opposite sides of the coiled coil hydrophobic core. Examination of published data on chemical crosslinking of lysines in vimentin confirmed that there were no instances of linkage within dimers between the nine pairs of identical lysines that lie more than one position within α-helical regions in parallel chains. Even among linkages that apparently were between dimers, only one of the 11 linkage products identified involved lysines that were both within an α-helical region. In 10 of the 11 identified linkages between dimers, one or both of the linked lysines were in regions of random coil conformation. These results of molecular modeling indicate that relative motion between polypeptide chains in oligomers of intermediate filament proteins is not sufficient to overcome an orientation of lysine groups that is unfavorable for their chemical linkage. This finding supports the interpretations of keratin cross-linking data indicating that parallel homodimers are the basis for keratin intermediate filament assembly. © 1996 Wiley-Liss, Inc.     

A new catalog of protein beta-sheets

Systematic protein folding studies depend on protein three-dimensional structure annotation, the assignment of amino acid structural types from atomic coordinates. Significant stabilizing factors between adjacent beta-sheet peptide chains have recently been characterized and were not considered during the development of previously published annotation methods. To produce an accurate beta-sheet domain catalog and to encompass the full beta-sheet spectacle, we developed a method, beta-Spider, which evaluates a packing energy between adjacent peptide chains in accordance with the newly discovered stabilizing factors. While considering important energetic factors, our approach also minimizes the use of subjective criteria, such as (phi,psi) boundaries and sets of H-bonding motifs that are used in other existing methods. As a result of the application of beta-Spider to a set of available high-resolution X-ray crystal structures, we present here a new beta-sheet catalog that differs considerably from the one produced by the most acclaimed DSSP method. The catalog includes new H-bonding motifs that were never reported.     

Oligoglycine surface structures: molecular dynamics simulation

The full-atomic molecular dynamics (MD) simulation of adsorption mode for diantennary oligoglycines [H-Gly4-NH(CH2)5]2 onto graphite and mica surface is described. The resulting structure of adsorption layers is analyzed. The peptide second structure motives have been studied by both STRIDE (structural identification) and DSSP (dictionary of secondary structure of proteins) methods. The obtained results confirm the possibility of polyglycine II (PGII) structure formation in diantennary oligoglycine (DAOG) monolayers deposited onto graphite surface, which was earlier estimated based on atomic-force microscopy measurements.     

Conformational analysis of protein structures derived from NMR data

A study is presented of the conformational characteristics of NMR-derived protein structures in the Protein Data Bank compared to X-ray structures. Both ensemble and energy-minimized average structures are analyzed. We have addressed the problem using the methods developed for crystal structures by examining the distribution of ?, Ψ, and χ angles as indicators of global conformational irregularity. All these features in NMR structures occur to varying degrees in multiple conformational states. Some measures of local geometry are very tightly constrained by the methods used to generate the structure, e.g., proline ? angles, α-helix ?, Ψ angles, ω angles, and C_α chirality. The more lightly restrained torsion angles do show increasead clustering as the number of overall experimental observations increases. ?, Ψ, and χ₁ angle conformational heterogeneity is strongly correlated with accessibility but shows additional differences which reflect the differing number of observations possible in NMR for the various side chains (e.g., many for Trp, few for Ser). In general, we find that the core is defined to a notional resolution of 2.0 to 2.3 Å. Of real interest is the behavior of surface residues and in particular the side chains where multiple rotameric states in different structures can vary from 10% to 88%. Later generation structures show a much tighter definition which correlates with increasing use of J-coupling information, stereospecific assignments, and heteronumclear techniques. A suite of programs is being developed to address the special needs of NMR-derived structures which will take into account the existence of increased mobility in solution. © 1993 Wiley-Liss, Inc.     

Surface structures of oligoglycines: A molecular dynamics simulation

An atomistic molecular dynamics (MD) simulation of the adsorption of biantennary oligoglycine [H-Gly₄-NH(CH₂)₅]₂ onto a graphite and mica surface is described. The structure of the resultant adsorption layers is analyzed. The secondary structure motifs of peptide blocks are studied by the STRIDE (structural identification) and DSSP (dictionary of the secondary structure of proteins) methods. The results of the study confirm the possibility of forming a polyglycine-II (PGII) structure in the monolayers of biantennary oligoglycines (BAOG) on a graphite surface previously supposed from the data of atomic force microscopy.     

Estimation of protein secondary structure from circular dichroism spectra: comparison of CONTIN, SELCON, and CDSSTR methods with an expanded reference set

We have expanded the reference set of proteins used in SELCON3 by including 11 additional proteins (selected from the reference sets of Yang and co-workers and Keiderling and co-workers). Depending on the wavelength range and whether or not denatured proteins are included in the reference set, five reference sets were constructed with the number of reference proteins varying from 29 to 48. The performance of three popular methods for estimating protein secondary structure fractions from CD spectra (implemented in software packages CONTIN, SELCON3, and CDSSTR) and a variant of CONTIN, CONTIN/LL, that incorporates the variable selection method in the locally linearized model in CONTIN, were examined using the five reference sets described here, and a 22-protein reference set. Secondary structure assignments from DSSP were used in the analysis. The performances of all three methods were comparable, in spite of the differences in the algorithms used in the three software packages. While CDSSTR performed the best with a smaller reference set and larger wavelength range, and CONTIN/LL performed the best with a larger reference set and smaller wavelength range, the performances for individual secondary structures were mixed. Analyzing protein CD spectra using all three methods should improve the reliability of predicted secondary structural fractions. The three programs are provided in CDPro software package and have been modified for easier use with the different reference sets described in this paper. CDPro software is available at the website: http://lamar.colostate.edu/ approximately sreeram/CDPro.     

Molecular modeling indicates that homodimers form the basis for intermediate filament assembly from human and mouse epidermal keratins

Mammalian epidermal keratin molecules adopt rod-shaped conformations that aggregate to form cytoplasmic intermediate filaments. To investigate these keratin conformations and the basis for their patterns of molecular association, graphical methods were developed to relate known amino acid sequences to probable spacial configurations. The results support the predominantly α-helical conformation of keratin chains, interrupted by short non-α-helical linkages. However, it was found that many of the linkages have amino acid sequences typical of β-strand conformations. Space-filling atomic models revealed that the β-strand sequences would permit the formation of 2-chain and 4-chain cylindrical β-helices, fully shielding the hydrophobic amino acid chains that alternate with hydrophilic residues in these sequences. Because of the locations of the β-helical regions in human and mouse stratum corneum keratin chains, only homodimers of the keratins could interact efficiently to form 2-chain and 4-chain β-helices. Tetramers having the directions and degrees of overlap of constituent dimers that have been identified by previous investigators are also predicted from the interactions of β-helical motifs. Heterotetramers formed from dissimilar homodimers could combine, through additional β-helical structures, to form higher oligomers having the dimensions seen in electron microscopic studies. Previous results from chemical crosslinking studies can be interpreted to support the concept of homodimers rather than heterodimers as the basis for keratin filament assembly. © 1995 Wiley-Liss, Inc.     

A consensus prediction of the secondary structure for the 6-phospho-β-D-galactosidase superfamily

Two separate unrefined models for the secondary structure of two subfamilies of the 6-phospho-β-D -galactosidase superfamily were independently constructed by examining patterns of variation and conservation within homologous protein sequences, assigning surface, interior, parsing, and active site residues to positions in the alignment, and identifying periodicities in these. A consensus model for the secondary structure of the entire superfamily was then built. The prediction tests the limits of an unrefined prediction made using this approach in a large protein with substantial functional and sequence divergence within the family. The protein belongs to the (α–β class), with the core β strands aligned parallel. The supersecondary structural elements that are readily identified in this model is a parallel β sheet built by strands C, D, and E, with helices 2 and 3 connecting strands (C + D) and (D + E), respectively, and an analogous α–β unit (strand G and helix 7) toward the end of the sequence. The resemblance of the supersecondary model to the tertiary structure formed by 8-fold α–β barrel proteins is almost certainly not coincidental. © 1995 Wiley-Liss, Inc.     

Effect of polyamines on the structure, thermal stability and 2,2,2-trifluoroethanol-induced aggregation of alpha-chymotrypsin

Naturally occurring polyamines are known to interact with a variety of biomolecules and critically involve in some important physiological processes. They have also been shown to influence protein aggregation in vitro in some cases. The aim of the present study was to investigate how polyamines may influence the structure and thermal stability of alpha-chymotrypsin and modulate alcohol-induced aggregation of this protein. Various techniques, including turbidity measurements, tensiometry, DSC, intrinsic fluorescence and far- and near-UV circular dichroism spectroscopy were used to examine the effect of putrescine and spermidine on alpha-chymotrypsin. While slight changes in the secondary and tertiary structure of the protein was observed, a clear stabilizing effect against its thermal unfolding was achieved. Moreover, the polyamines were found to inhibit TFE-induced aggregation at 32% TFE and promote formation of non-native alpha-helices in the protein structure. Based on the observed increase in surface tension induced by polyamines, it is suggested that their effects on enhancing thermal stability and alcohol-induced alpha-helices formation may be due to their kosmotropic properties.