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.     

Coupling backbone flexibility and amino acid sequence selection in protein design.

Using a protein design algorithm that considers side-chain packing quantitatively, the effect of explicit backbone motion on the selection of amino acids in protein design was assessed in the core of the streptococcal protein G beta 1 domain (G beta 1). Concerted backbone motion was introduced by varying G beta 1's supersecondary structure parameter values. The stability and structural flexibility of seven of the redesigned proteins were determined experimentally and showed that core variants containing as many as 6 of 10 possible mutations retain native-like properties. This result demonstrates that backbone flexibility can be combined explicitly with amino acid side-chain selection and that the selection algorithm is sufficiently robust to tolerate perturbations as large as 15% of G beta 1's native supersecondary structure parameter values.     

Structural classification of αββ and ββα supersecondary structure units in proteins

We present a fully automatic structural classification of supersecondary structure units, consisting of two hydrogen-bonded β strands, preceded or followed by an α helix. The classification is performed on the spatial arrangement of the secondary structure elements, irrespective of the length and conformation of the intervening loops. The similarity of the arrangements is estimated by a structure alignment procedure that uses as similarity measure the root mean square deviation of superimposed backbone atoms. Applied to a set of 141 well-resolved nonhomologous protein structures, the classification yields 11 families of recurrent arrangements. In addition, fragments that are structurally intermediate between the families are found; they reveal the continuity of the classification. The analysis of the families shows that the α helix and β hairpin axes can adopt virtually all relative orientations, with, however, some preferable orientations; moreover, according to the orientation, preferences in the left/right handedness of the α–β connection are observed. These preferences can be explained by favorable side by side packing of the α helix and the β hairpin, local interactions in the region of the α–β connection or stabilizing environments in the parent protein. Furthermore, fold recognition procedures and structure prediction algorithms coupled to database-derived potentials suggest that the preferable nature of these arrangements does not imply their intrinsic stability. They usually accommodate a large number of sequences, of which only a subset is predicted to stabilize the motif. The motifs predicted as stable could correspond to nuclei formed at the very beginning of the folding process. Proteins 30:193–212, 1998. © 1998 Wiley-Liss, Inc.     

Contact patterns between helices and strands of sheet define protein folding patterns

Comparing and classifying protein folding patterns allows organizing the known structures and enumerating possible protein structural patterns including those not yet observed. We capture the essence of protein folding patterns in a concise tableau representation based on the order and contact patterns of secondary structures: helices and strands of sheet. The tableaux are intelligible to both humans and computers. They provide a database, derived from the Protein Data Bank, mineable in studies of protein architecture. Using this database, we have: (i) determined statistical properties of secondary structure contacts in an unbiased set of protein domains from ASTRAL, (ii) observed that in 98% of cases, the tableau is a faithful representation of the folding pattern as classified in SCOP, (iii) demonstrated that to a large extent the local structure of proteins indicates their complete folding topology, and (iv) studied the use of the representation for fold identification.     

The GD box: A widespread noncontiguous supersecondary structural element

Identification and characterization of recurrent supersecondary structural elements is central to understanding the rules governing protein tertiary structure. Here, we describe the GD box, a widespread noncontiguous supersecondary element, which we initially found in a group of topologically distinct but homologous β‐barrels—the cradle‐loop barrels. The GD box is similar both in sequence and structure and comprises two short unpaired β‐strands connected by an orthogonal type‐II β‐turn and a noncontiguous β‐strand forming hydrogen bonds with the β‐turn. Using structure‐based analysis, we have detected 518 instances of the GD box in a nonredundant subset of the SCOP database comprising 3771 domains. Apart from the cradle‐loop barrels, this motif is also found in a diverse set of nonhomologous folds including other topologically related β‐barrels. Since nonlocal interactions are fundamental in the formation of protein structure, systematic identification and characterization of other noncontiguous supersecondary structural elements is likely to prove valuable to protein structure modeling, validation, and prediction.     

Rules for connectivity of secondary structure elements in protein: Two–layer αβ sandwiches

In protein structures, the fold is described according to the spatial arrangement of secondary structure elements (SSEs: α‐helices and β‐strands) and their connectivity. The connectivity or the pattern of links among SSEs is one of the most important factors for understanding the variety of protein folds. In this study, we introduced the connectivity strings that encode the connectivities by using the types, positions, and connections of SSEs, and computationally enumerated all the connectivities of two‐layer αβ sandwiches. The calculated connectivities were compared with those in natural proteins determined using MICAN, a nonsequential structure comparison method. For 2α‐4β, among 23,000 of all connectivities, only 48 were free from irregular connectivities such as loop crossing. Of these, only 20 were found in natural proteins and the superfamilies were biased toward certain types of connectivities. A similar disproportional distribution was confirmed for most of other spatial arrangements of SSEs in the two‐layer αβ sandwiches. We found two connectivity rules that explain the bias well: the abundances of interlayer connecting loops that bridge SSEs in the distinct layers; and nonlocal β‐strand pairs, two spatially adjacent β‐strands located at discontinuous positions in the amino acid sequence. A two‐dimensional plot of these two properties indicated that the two connectivity rules are not independent, which may be interpreted as a rule for the cooperativity of proteins.     

Folding pattern of the alpha-crystallin domain in alphaA-crystallin determined by site-directed spin labeling

The folding pattern of the alpha-crystallin domain, a conserved protein module encoding the molecular determinants of structure and function in the small heat-shock protein superfamily, was determined in the context of the lens protein alphaA-crystallin by systematic application of site-directed spin labeling. The sequence-specific secondary structure was assigned primarily from nitroxide scanning experiments in which the solvent accessibility and mobility of a nitroxide probe were measured as a function of residue number. Seven beta-strands were identified and their orientation relative to the aqueous solvent determined, thus defining the residues lining the hydrophobic core. The pairwise packing of adjacent strands in the primary structure was deduced from patterns of proximities in nitroxide pairs with one member on the exposed surface of each strand. In addition to identifying supersecondary structures, these proximities revealed that the seven strands are arranged in two beta-sheets. The overall packing of the two sheets was determined by application of the general rules of protein structure and from proximities in nitroxide pairs designed to distinguish between known all beta-sheet folds. Our data are consistent with an immunoglobulin-like fold consisting of two aligned beta-sheets. Comparison of this folding pattern to that of the evolutionary distant alpha-crystallin domain in Methanococcus jannaschii heat-shock protein 16.5 reveals a conserved core structure with the differences sequestered at one edge of the beta-sandwich. A beta-strand deletion in alphaA-crystallin disrupts a subunit interface and allows for a different dimerization motif. Putative substrate binding regions appear to include a buried loop and a buried turn, suggesting that the chaperone function involves a disassembly of the oligomer.     

Computer simulations aimed at structure prediction of supersecondary motifs in proteins

It is well established that protein structures are more conserved than protein sequences. One-third of all known protein structures can be classified into ten protein folds, which themselves are composed mainly of alpha-helical hairpin, beta hairpin, and betaalphabeta supersecondary structural elements. In this study, we explore the ability of a recent Monte Carlo-based procedure to generate the 3D structures of eight polypeptides that correspond to units of supersecondary structure and three-stranded antiparallel beta sheet. Starting from extended or misfolded compact conformations, all Monte Carlo simulations show significant success in predicting the native topology using a simplified chain representation and an energy model optimized on other structures. Preliminary results on model peptides from nucleotide binding proteins suggest that this simple protein folding model can help clarify the relation between sequence and topology.     

The N-terminal domain of MDM2 resembles calmodulin and its relatives.

It is shown here that the N-terminal domain of MDM2, which is not thought to bind calcium ions, otherwise bears a striking resemblance to a cluster of four EF-hand modules like those found in the calmodulin family. There are similarities in module arrangement, supersecondary structure and the main-chain to main-chain hydrogen-bonding pattern, especially in the vicinity of the short antiparallel beta-sheet, the two strands of which lie between the two E and F helices of tandem modules. Some conserved amino acid residues are identified that are associated with short side-chain to main-chain hydrogen-bonded motifs. Also, both types of domain bind a short, functionally important hydrophobic alpha-helix from another protein in a cavity between the two pairs of EF-hand, or EF-hand-like, modules.     

Description and recognition of regular and distorted secondary structures in proteins using the automated protein structure analysis method

The Automated Protein Structure Analysis (APSA) method, which describes the protein backbone as a smooth line in three‐dimensional space and characterizes it by curvature κ and torsion τ as a function of arc length s, was applied on 77 proteins to determine all secondary structural units via specific κ(s) and τ(s) patterns. A total of 533 α‐helices and 644 β‐strands were recognized by APSA, whereas DSSP gives 536 and 651 units, respectively. Kinks and distortions were quantified and the boundaries (entry and exit) of secondary structures were classified. Similarity between proteins can be easily quantified using APSA, as was demonstrated for the roll architecture of proteins ubiquitin and spinach ferridoxin. A twenty‐by‐twenty comparison of all α domains showed that the curvature‐torsion patterns generated by APSA provide an accurate and meaningful similarity measurement for secondary, super secondary, and tertiary protein structure. APSA is shown to accurately reflect the conformation of the backbone effectively reducing three‐dimensional structure information to two‐dimensional representations that are easy to interpret and understand. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Conservation of polyproline II helices in homologous proteins: implications for structure prediction by model building.

Left-handed polyproline II (PPII) helices commonly occur in globular proteins in segments of 4-8 residues. This paper analyzes the structural conservation of PPII-helices in 3 protein families: serine proteinases, aspartic proteinases, and immunoglobulin constant domains. Calculations of the number of conserved segments based on structural alignment of homologous molecules yielded similar results for the PPII-helices, the alpha-helices, and the beta-strands. The PPII-helices are consistently conserved at the level of 100-80% in the proteins with sequence identity above 20% and RMS deviation of structure alignments below 3.0 A. The most structurally important PPII segments are conserved below this level of sequence identity. These results suggest that the PPII-helices, in addition to the other 2 secondary structure classes, should be identified as part of structurally conserved regions in proteins. This is supported by similar values for the local RMS deviations of the aligned segments for the structural classes of PPII-helices, alpha-helices, and beta-strands. The PPII-helices are shown to participate in supersecondary elements such as PPII-helix/alpha-helix. The conservation of PPII-helices depends on the conservation of a supersecondary element as a whole. PPII-helices also form links, possibly flexible, in the interdomain regions. The role of the PPII-helices in model building by homology is 2-fold; they serve as additional conserved elements in the structure allowing improvement of the accuracy of a model and provide correct chain geometry for modeling of the segments equivalenced to them in a target sequence. The improvement in model building is demonstrated in 2 test studies.     

Prediction of the tertiary structure of the alpha-subunit of tryptophan synthase

The tertiary structure of the alpha-subunit of tryptophan synthase was proposed using a combination of experimental data and computational methods. The vacuum-ultraviolet circular dichroism spectrum was used to assign the protein to the alpha/beta-class of supersecondary structures. The two-domain structure of the alpha-subunit (Miles et al.: Biochemistry 21:2586, 1982; Beasty and Matthews: Biochemistry 24:3547, 1985) eliminated consideration of a barrel structure and focused attention on a beta-sheet structure. An algorithm (Cohen et al.: Biochemistry 22:4894, 1983) was used to generate a secondary structure prediction that was consistent with the sequence data of the alpha-subunit from five species. Three potential secondary structures were then packed into tertiary structures using other algorithms. The assumption of nearest neighbors from second-site revertant data eliminated 97% of the possible tertiary structures; consideration of conserved hydrophobic packing regions on the beta-sheet eliminated all but one structure. The native structure is predicted to have a parallel beta-sheet flanked on both sides by alpha-helices, and is consistent with the available data on chemical cross-linking, chemical modification, and limited proteolysis. In addition, an active site region containing appropriate residues could be identified as well as an interface for beta 2-subunit association. The ability of experimental data to facilitate the prediction of protein structure is discussed.     

Role of loops connecting secondary structure elements in the stabilization of proteins isolated from thermophilic organisms

It has been recently discovered that the connection of secondary structure elements (ββ‐unit, βα‐ and αβ‐units) in proteins follows quite stringent principles regarding the chirality and the orientation of the structural units (Koga et al., Nature 2012;491:222–227). By exploiting these rules, a number of protein scaffolds endowed with a remarkable thermal stability have been designed (Koga et al., Nature 2012;491:222–227). By using structural databases of proteins isolated from either mesophilic or thermophilic organisms, we here investigate the influence of supersecondary associations on the thermal stability of natural proteins. Our results suggest that β‐hairpins of proteins from thermophilic organisms are very frequently characterized by shortenings of the loops. Interestingly, this shortening leads to states that display a very strong preference for the most common connectivity of the strands observed in native protein hairpins. The abundance of selective states in these proteins suggests that they may achieve a high stability by adopting a strategy aimed to reduce the possible conformations of the unfolded ensemble. In this scenario, our data indicate that the shortening is effective if it increases the adherence to these rules. We also show that this mechanism may operate in the stabilization of well‐known protein folds (thioredoxin and RNase A). These findings suggest that future investigations aimed at defining mechanism of protein stabilization should also consider these effects.     

Molecular Dynamics Simulation of Protein Folding with Supersecondary Structure Constraints

We integrate molecular dynamics simulation methods with a newly developed supersecondary structure prediction method and compute the structure of a protein molecule, crambin. The computed structure is similar to the crystal structure with an rms error of 3.94 Å.     

A novel fold recognition method using composite predicted secondary structures

In this work, we introduce a new method for fold recognition using composite secondary structures assembled from different secondary structure prediction servers for a given target sequence. An automatic, complete, and robust way of finding all possible combinations of predicted secondary structure segments (SSS) for the target sequence and clustering them into a few flexible clusters, each containing patterns with the same number of SSS, is developed. This program then takes two steps in choosing plausible homologues: (i) a SSS-based alignment excludes impossible templates whose SSS patterns are very different from any of those of the target; (ii) a residue-based alignment selects good structural templates based on sequence similarity and secondary structure similarity between the target and only those templates left in the first stage. The secondary structure of each residue in the target is selected from one of the predictions to find the best match with the template. Truncation is applied to a target where different predictions vary. In most cases, a target is also divided into N-terminal and C-terminal fragments, each of which is used as a separate subsequence. Our program was tested on the fold recognition targets from CASP3 with known PDB codes and some available targets from CASP4. The results are compared with a structural homologue list for each target produced by the CE program (Shindyalov and Bourne, Protein Eng 1998;11:739-747). The program successfully locates homologues with high Z-score and low root-mean-score deviation within the top 30-50 predictions in the overwhelming majority of cases.     

Discovery of a significant,nontopological preference for antiparallel alignment of helices with parallel regions in sheets

To help elucidate the role of secondary structure packing preferences in protein folding, here we present an analysis of the packing geometry observed between alpha-helices and between alpha-helices and beta-sheets in 1316 diverse, nonredundant protein structures. Finite-length vectors were fit to the alpha-carbon atoms in each of the helices and strands, and the packing angle between the vectors, Omega, was determined at the closest point of approach within each helix-helix or helix-sheet pair. Helix-sheet interactions were found in 391 of the proteins, and the distributions of Omega values were calculated for all the helix-sheet and helix-helix interactions. The packing angle preferences for helix-helix interactions are similar to those previously observed. However, analysis of helix-strand packing preferences uncovered a remarkable tendency for helices to align antiparallel to parallel regions of beta-sheets, independent of the topological constraints or prevalence of beta-alpha-beta motifs in the proteins. This packing angle preference is significantly diminished in helix interactions involving mixed and antiparallel beta-sheets, suggesting a role for helix-sheet dipole alignment in guiding supersecondary structure formation in protein folding. This knowledge of preferred packing angles can be used to guide the engineering of stable protein modules.     

Non-intertwined binary patterns of hydrophobic/nonhydrophobic amino acids are considerably better markers of regular secondary structures than nonconstrained patterns

Patterns of hydrophobic and hydrophilic residues (binary patterns) play an important role in protein architecture and can be roughly categorized into two classes regarding their preferential participation in α‐helices or β‐strands. However, a single binary pattern can be embedded into different longer patterns carrying opposite structural information and thus cannot be as much informative as expected. Here, we consider conditional binary patterns, or hydrophobic clusters, whose existence is conditioned by the presence of a minimum number of nonhydrophobic residues, called the connectivity distance, that separate two hydrophobic amino acids assumed to belong to two distinct patterns. Conditional binary patterns are distinct from simple ones in that they are not intertwined, i.e., they can not include or be included in other conditional patterns and therefore carry a much more differentiated information, in particular being dramatically better correlated with regular secondary structures (especially β ones). The distribution of these nonintertwined binary patterns in natural proteins was assessed relative to randomness, evidencing the structural bricks that are favored and disfavored by evolutionary selection. Several connectivity distances as well as several hydrophobic alphabets were tested, evidencing the clear superiority of a connectivity distance of 4, which mimics the minimum current length of loops in globular domains, and of the VILFMYW alphabet, selected from structural data (secondary structure propension and Voronoï tesselation), in highlighting fundamental properties of protein folds. Proteins 2003;51:236–244. © 2003 Wiley‐Liss, Inc.     

The major amino acid transporter superfamily has a similar core structure as Na+-galactose and Na+-leucine transporters

The sodium solute symporters (SSS) and neurotransmitter sodium symporters (NSS) are two families of secondary transporters that are not related in amino acid sequence. Nonetheless, recent crystal structures showed that the Na⁺/galactose (SSS) and Na⁺/leucine (NSS) transporters have similar core structures. The structural relatedness highlights the need for classification methods for membrane protein structures based on other criteria than amino acid similarity. Here, we demonstrate that a method based on hydropathy profile alignments convincingly identifies structural similarity between the NSS and SSS families. Most importantly, the method shows that one of the largest transporter families for which a crystal structure is elusive (the amino acid/polyamine/organocation or APC superfamily), also shares the similar core structure observed for the Na⁺/galactose and Na⁺/leucine transporters. The APC superfamily contains the major amino acid transporter families that are found throughout life. Insight into their structure will significantly facilitate the studies of this important group of transporters.     

Prediction of protein secondary structure content for the twilight zone sequences

Secondary protein structure carries information about local structural arrangements, which include three major conformations: alpha-helices, beta-strands, and coils. Significant majority of successful methods for prediction of the secondary structure is based on multiple sequence alignment. However, multiple alignment fails to provide accurate results when a sequence comes from the twilight zone, that is, it is characterized by low (<30%) homology. To this end, we propose a novel method for prediction of secondary structure content through comprehensive sequence representation, called PSSC-core. The method uses a multiple linear regression model and introduces a comprehensive feature-based sequence representation to predict amount of helices and strands for sequences from the twilight zone. The PSSC-core method was tested and compared with two other state-of-the-art prediction methods on a set of 2187 twilight zone sequences. The results indicate that our method provides better predictions for both helix and strand content. The PSSC-core is shown to provide statistically significantly better results when compared with the competing methods, reducing the prediction error by 5-7% for helix and 7-9% for strand content predictions. The proposed feature-based sequence representation uses a comprehensive set of physicochemical properties that are custom-designed for each of the helix and strand content predictions. It includes composition and composition moment vectors, frequency of tetra-peptides associated with helical and strand conformations, various property-based groups like exchange groups, chemical groups of the side chains and hydrophobic group, auto-correlations based on hydrophobicity, side-chain masses, hydropathy, and conformational patterns for beta-sheets. The PSSC-core method provides an alternative for predicting the secondary structure content that can be used to validate and constrain results of other structure prediction methods. At the same time, it also provides useful insight into design of successful protein sequence representations that can be used in developing new methods related to prediction of different aspects of the secondary protein structure.