Predicting the beta-helix fold from protein sequence data.

A method is presented that uses beta-strand interactions to predict the parallel right-handed beta-helix super-secondary structural motif in protein sequences. A program called BetaWrap implements this method and is shown to score known beta-helices above non-beta-helices in the Protein Data Bank in cross-validation. It is demonstrated that BetaWrap learns each of the seven known SCOP beta-helix families, when trained primarily on beta-structures that are not beta-helices, together with structural features of known beta-helices from outside the family. BetaWrap also predicts many bacterial proteins of unknown structure to be beta-helices; in particular, these proteins serve as virulence factors, adhesins, and toxins in bacterial pathogenesis and include cell surface proteins from Chlamydia and the intestinal bacterium Helicobacter pylori. The computational method used here may generalize to other beta-structures for which strand topology and profiles of residue accessibility are well conserved.     

Computational study of the fibril organization of polyglutamine repeats reveals a common motif identified in beta-helices

The formation of fibril aggregates by long polyglutamine sequences is assumed to play a major role in neurodegenerative diseases such as Huntington. Here, we model peptides rich in glutamine, through a series of molecular dynamics simulations. Starting from a rigid nanotube-like conformation, we have obtained a new conformational template that shares structural features of a tubular helix and of a beta-helix conformational organization. Our new model can be described as a super-helical arrangement of flat beta-sheet segments linked by planar turns or bends. Interestingly, our comprehensive analysis of the Protein Data Bank reveals that this is a common motif in beta-helices (termed beta-bend), although it has not been identified so far. The motif is based on the alternation of beta-sheet and helical conformation as the protein sequence is followed from the N to the C termini (beta-alpha(R)-beta-polyPro-beta). We further identify this motif in the ssNMR structure of the protofibril of the amyloidogenic peptide Abeta(1-40). The recurrence of the beta-bend suggests a general mode of connecting long parallel beta-sheet segments that would allow the growth of partially ordered fibril structures. The design allows the peptide backbone to change direction with a minimal loss of main chain hydrogen bonds. The identification of a coherent organization beyond that of the beta-sheet segments in different folds rich in parallel beta-sheets suggests a higher degree of ordered structure in protein fibrils, in agreement with their low solubility and dense molecular packing.     

Fold recognition and accurate sequence-structure alignment of sequences directing beta-sheet proteins

The ability to predict structure from sequence is particularly important for toxins, virulence factors, allergens, cytokines, and other proteins of public health importance. Many such functions are represented in the parallel beta-helix and beta-trefoil families. A method using pairwise beta-strand interaction probabilities coupled with evolutionary information represented by sequence profiles is developed to tackle these problems for the beta-helix and beta-trefoil folds. The algorithm BetaWrapPro employs a "wrapping" component that may capture folding processes with an initiation stage followed by processive interaction of the sequence with the already-formed motifs. BetaWrapPro outperforms all previous motif recognition programs for these folds, recognizing the beta-helix with 100% sensitivity and 99.7% specificity and the beta-trefoil with 100% sensitivity and 92.5% specificity, in crossvalidation on a database of all nonredundant known positive and negative examples of these fold classes in the PDB. It additionally aligns 88% of residues for the beta-helices and 86% for the beta-trefoils accurately (within four residues of the exact position) to the structural template, which is then used with the side-chain packing program SCWRL to produce 3D structure predictions. One striking result has been the prediction of an unexpected parallel beta-helix structure for a pollen allergen, and its recent confirmation through solution of its structure. A Web server running BetaWrapPro is available and outputs putative PDB-style coordinates for sequences predicted to form the target folds.     

Cold survival in freeze-intolerant insects: the structure and function of beta-helical antifreeze proteins.

Antifreeze proteins (AFPs) designate a class of proteins that are able to bind to and inhibit the growth of macromolecular ice. These proteins have been characterized from a variety of organisms. Recently, the structures of AFPs from the spruce budworm (Choristoneura fumiferana) and the yellow mealworm (Tenebrio molitor) have been determined by NMR and X-ray crystallography. Despite nonhomologous sequences, both proteins were shown to consist of beta-helices. We review the structures and dynamics data of these two insect AFPs to bring insight into the structure-function relationship and explore their beta-helical architecture. For the spruce budworm protein, the fold is a left-handed beta-helix with 15 residues per coil. The Tenebrio molitor protein consists of a right-handed beta-helix with 12 residues per coil. Mutagenesis and structural studies show that the insect AFPs present a highly rigid array of threonine residues and bound water molecules that can effectively mimic the ice lattice. Comparisons of the newly determined ryegrass and carrot AFP sequences have led to models suggesting that they might also consist of beta-helices, and indicate that the beta-helix might be used as an AFP structural motif in nonfish organisms.     

Two-rung model of a left-handed beta-helix for prions explains species barrier and strain variation in transmissible spongiform encephalopathies

In this study, a new beta-helical model is proposed that explains the species barrier and strain variation in transmissible spongiform encephalopathies. The left-handed beta-helix serves as a structural model that can explain the seeded growth characteristics of beta-sheet structure in PrP(Sc) fibrils. Molecular dynamics simulations demonstrate that the left-handed beta-helix is structurally more stable than the right-handed beta-helix, with a higher beta-sheet content during the simulation and a better distributed network of inter-strand backbone-backbone hydrogen bonds between parallel beta-strands of different rungs. Multiple sequence alignments and homology modelling of prion sequences with different rungs of left-handed beta-helices illustrate that the PrP region with the highest beta-helical propensity (residues 105-143) can fold in just two rungs of a left-handed beta-helix. Even if no other flanking sequence participates in the beta-helix, the two rungs of a beta-helix can give the growing fibril enough elevation to accommodate the rest of the PrP protein in a tight packing at the periphery of a trimeric beta-helix. The folding of beta-helices is driven by backbone-backbone hydrogen bonding and stacking of side-chains in adjacent rungs. The sequence and structure of the last rung at the fibril end with unprotected beta-sheet edges selects the sequence of a complementary rung and dictates the folding of the new rung with optimal backbone hydrogen bonding and side-chain stacking. An important side-chain stack that facilitates the beta-helical folding is between methionine residues 109 and 129, which explains their importance in the species barrier of prions. Because the PrP sequence is not evolutionarily optimised to fold in a beta-helix, and because the beta-helical fold shows very little sequence preference, alternative alignments are possible that result in a different rung able to select for an alternative complementary rung. A different top rung results in a new strain with different growth characteristics. Hence, in the present model, sequence variation and alternative alignments clarify the basis of the species barrier and strain specificity in PrP-based diseases.     

SARS病毒S蛋白三维结构预测

蛋白质结构类型识别方法可以在没有序列同源性的蛋白质之间检测有没有结构相似性。利用蛋白质结构类型识别方法预测了SARS病毒S蛋白N端区域的结构。模建的SARS病毒S蛋白N端区域是一个全折叠的结构。     

Fold recognition and ab initio structure predictions using hidden markov models and β-strand pair potentials

Protein structure predictions were submitted for 9 of the target sequences in the competition that ran during 1994. Targets sequences were selected that had no known homology with any sequence of known structure and were members of a reasonably sized family of related but divergent sequences. The objective was either to recognize a compatible fold for the target sequence in the database of known structures or to predict ab initio its rough 3D topology. The main tools used were Hidden Markov models (HMM) for fold recognition, a β- strand pair potential to predict β-sheet topology, and the PHD server for secondary structure prediction. Compatible folds were correctly identified in a number of cases and the β-strand pair potential was shown to be a useful tool for ab initio topology prediction. © 1995 Wiley-Liss, Inc.     

Spruce budworm antifreeze protein: changes in structure and dynamics at low temperature

Antifreeze proteins (AFPs) prevent the growth of ice, and are used by some organisms that live in sub-zero environments for protection against freezing. All AFPs are thought to function by an adsorption inhibition process. In order to elucidate the ice-binding mechanism, the structures of several AFPs have been determined, and have been shown to consist of different folds. Recently, the first structures of the highly active insect AFPs have been characterized. These proteins have a beta-helix structure, which adds yet another fold to the AFP family. The 90-residue spruce budworm (Choristoneura fumiferana) AFP consists of a beta-helix with 15 residues per coil. The structure contains two ranks of aligned threonine residues (known as the TXT motif), which were shown by mutagenesis experiments to be located in the ice-binding face. In our previous NMR study of this AFP at 30 degrees C, we found that the TXT face was not optimally defined because of the broadening of NMR resonances potentially due to weak oligomerization. We present here a structure of spruce budworm AFP determined at 5 degrees C, where this broadening is reduced. In addition, the 1H-15N NMR dynamics of the protein were examined at 30 degrees C and 5 degrees C. The results show that the spruce budworm AFP is more structured at 5 degrees C, and support the general observation that AFPs become more rigid as the temperature is lowered.     

Solution NMR structure of a D,L-alternating oligonorleucine as a model of beta-helix

beta-Helix structures are of particular interest due to their capacity to form transmembrane channels with different transport properties. However, the relatively large number of beta-helices configurations does not allow a direct conformational analysis of beta-helical oligopeptides. A synthetic alternating D,L-oligopeptide with twelve norleucines (XIIMe) has been used as a model to get insight in the conformational features of beta-helix structures. The spatial configuration of XIIMe in solution has been determined by NMR. An extensive set of distances (nuclear Overhauser effect) and dihedral (J coupling constants) constraints have been included in molecular dynamics calculations. The NMR experimental data and theoretical calculations clearly indicate that the XIIMe adopts a single beta(4.4)-helix-type conformation in nonpolar solvents.     

Structural Characteristics of Novel Protein Folds

Folds are the basic building blocks of protein structures. Understanding the emergence of novel protein folds is an important step towards understanding the rules governing the evolution of protein structure and function and for developing tools for protein structure modeling and design. We explored the frequency of occurrences of an exhaustively classified library of supersecondary structural elements (Smotifs), in protein structures, in order to identify features that would define a fold as novel compared to previously known structures. We found that a surprisingly small set of Smotifs is sufficient to describe all known folds. Furthermore, novel folds do not require novel Smotifs, but rather are a new combination of existing ones. Novel folds can be typified by the inclusion of a relatively higher number of rarely occurring Smotifs in their structures and, to a lesser extent, by a novel topological combination of commonly occurring Smotifs. When investigating the structural features of Smotifs, we found that the top 10% of most frequent ones have a higher fraction of internal contacts, while some of the most rare motifs are larger, and contain a longer loop region.     

Emergence of diverse biochemical activities in evolutionarily conserved structural scaffolds of proteins

Comparative analysis of numerous protein structures that have become available in the past few years, combined with genome comparison, has yielded new insights into the evolution of enzymes and their functions. In addition to the well-known diversification of substrate specificities, enzymes with several widespread catalytic folds, particularly the TIM barrel, the RRM-like domain and the double-stranded beta-helix (cupin) domain, have been extensively explored in 'reaction space', resulting in the evolution of numerous, diverse catalytic activities supported by the same structural scaffold. Common protein folds differ widely in the diversity of catalyzed reactions. The biochemical plasticity of a fold seems to hinge on the presence of a generic, symmetrical substrate-binding pocket as opposed to highly specialized binding sites.     

Computational exploration of the network of sequence flow between protein structures

We investigate small sequence adjustments (of one or a few amino acids) that induce large conformational transitions between distinct and stable folds of proteins. Such transitions are intriguing from evolutionary and protein‐design perspectives. They make it possible to search for ancient protein structures or to design protein switches that flip between folds and functions. A network of sequence flow between protein folds is computed for representative structures of the Protein Data Bank. The computed network is dense, on an average each structure is connected to tens of other folds. Proteins that attract sequences from a higher than expected number of neighboring folds are more likely to be enzymes and alpha/beta fold. The large number of connections between folds may reflect the need of enzymes to adjust their structures for alternative substrates. The network of the Cro family is discussed, and we speculate that capacity is an important factor (but not the only one) that determines protein evolution. The experimentally observed flip from all alpha to alpha + beta fold is examined by the network tools. A kinetic model for the transition of sequences between the folds (with only protein stability in mind) is proposed. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Constructing templates for protein structure prediction by simulation of protein folding pathways

How a one-dimensional protein sequence folds into a specific 3D structure remains a difficult challenge in structural biology. Many computational methods have been developed in an attempt to predict the tertiary structure of the protein; most of these employ approaches that are based on the accumulated knowledge of solved protein structures. Here we introduce a novel and fully automated approach for predicting the 3D structure of a protein that is based on the well accepted notion that protein folding is a hierarchical process. Our algorithm follows the hierarchical model by employing two stages: the first aims to find a match between the sequences of short independently-folding structural entities and parts of the target sequence and assigns the respective structures. The second assembles these local structural parts into a complete 3D structure, allowing for long-range interactions between them. We present the results of applying our method to a subset of the targets from CASP6 and CASP7. Our results indicate that for targets with a significant sequence similarity to known structures we are often able to provide predictions that are better than those achieved by two leading servers, and that the most significant improvements in comparison with these methods occur in regions of a gapped structural alignment between the native structure and the closest available structural template. We conclude that in addition to performing well for targets with known homologous structures, our method shows great promise for addressing the more general category of comparative modeling targets, which is our next goal.     

A Fourier analysis of symmetry in protein structure

The score matrix from a structure comparison program (SAP) was used to search for repeated structures using a Fourier analysis. When tested with artificial data, a simple Fourier transform of the smoothed matrix provided a clear signal of the repeat periodicity that could be used to extract the repeating units with the SAP program. The strength of the Fourier signal was calibrated against the signal from model proteins. The most useful of these was the novel random-walk approach employed to generate realistic 'fake' structures. On the basis of these it was possible to conclude that only a small proportion of protein structures have an unexpected degree of symmetry. Artificially generated 'ideal' folds provided an upper limit on the strength of signal that could be expected from a 'perfectly' repeating compact structure. Unexpectedly, some of the very regular beta-propellor folds attained the same strength but the majority of symmetric structures lay below this region. When native proteins were ranked by the power of their spectrum a wide variety of fold types were seen to score highly. In the betaalpha class, these included the globular betaalpha proteins and the more repetitive leucine-rich betaalpha folds. In the all-beta class; beta-propellors, beta-prisms and beta-helices were found as well as the more globular gamma-crystalin domains. When this ranked list was filtered to remove proteins that contained detectable internal sequence similarity (using the program REPRO), the list became exclusively composed of just globular betaalpha class proteins and in the top 50 re-ranked proteins, only a single 4-fold propellor structure remained.     

A fast SCOP fold classification system using content-based E-Predict algorithm

Background  

Domain experts manually construct the Structural Classification of Protein (SCOP) database to categorize and compare protein structures. Even though using the SCOP database is believed to be more reliable than classification results from other methods, it is labor intensive. To mimic human classification processes, we develop an automatic SCOP fold classification system to assign possible known SCOP folds and recognize novel folds for newly-discovered proteins.     

Predicting protein folds with fold-specific PSSM libraries

Accurately assigning folds for divergent protein sequences is a major obstacle to structural studies. Herein, we outline an effective method for fold recognition using sets of PSSMs, each of which is constructed for different protein folds. Our analyses demonstrate that FSL (Fold-specific Position Specific Scoring Matrix Libraries) can predict/relate structures given only their amino acid sequences of highly divergent proteins. This ability to detect distant relationships is dependent on low-identity sequence alignments obtained from FSL. Results from our experiments demonstrate that FSL perform well in recognizing folds from the "twilight-zone" SABmark dataset. Further, this method is capable of accurate fold prediction in newly determined structures. We suggest that by building complete PSSM libraries for all unique folds within the Protein Database (PDB), FSL can be used to rapidly and reliably annotate a large subset of protein folds at proteomic level. The related programs and fold-specific PSSMs for our FSL are publicly available at: http://ccp.psu.edu/download/FSLv1.0/.     

Fibril formation by amyloid-beta proteins may involve beta-helical protofibrils.

We have proposed that amyloid fibrils contain subunits (protofibrils) that are formed from beta-strands wound into continuous 2-3 nm-diameter beta-helices. Subsequent lateral aggregation of the beta-helices to form the widely observed 5-12 nm-diameter fibrils could be promoted by hydrophobic residues on the exterior of the postulated beta-helix. A number of short peptide fragments of the amyloid-beta (A beta) proteins, such as A beta34-42 [LMVGGVVIA], the nine-residue, carboxyl-terminal portion of A beta1-42, can also form amyloid fibrils. In the present study, it was found that a beta-helix formed from A beta34-42 accounts for features suggested by published rotational resonance solid-state NMR data, including an anomalous conformation about the Gly-37-Gly-38 region and exaggerated pleating. An analogue of A beta34-42 was synthesized in which the hydrophobic groups on the exterior of the postulated beta-helix were replaced with glutamates, giving LEVGGVEIE. The analogue was completely soluble at pH 7, but at pH 2.5 it produced 2-2.5 nm-diameter fibrils which did not associate into larger-diameter bundles. The results of this study support the proposal that amyloid fibrils are formed from beta-helical subunits.