Protein fold recognition using segmentation conditional random fields (SCRFs).

Protein fold recognition is an important step towards understanding protein three-dimensional structures and their functions. A conditional graphical model, i.e., segmentation conditional random fields (SCRFs), is proposed as an effective solution to this problem. In contrast to traditional graphical models, such as the hidden Markov model (HMM), SCRFs follow a discriminative approach. Therefore, it is flexible to include any features in the model, such as overlapping or long-range interaction features over the whole sequence. The model also employs a convex optimization function, which results in globally optimal solutions to the model parameters. On the other hand, the segmentation setting in SCRFs makes their graphical structures intuitively similar to the protein 3-D structures and more importantly provides a framework to model the long-range interactions between secondary structures directly. Our model is applied to predict the parallel beta-helix fold, an important fold in bacterial pathogenesis and carbohydrate binding/cleavage. The cross-family validation shows that SCRFs not only can score all known beta-helices higher than non-beta-helices in the Protein Data Bank (PDB), but also accurately locates rungs in known beta-helix proteins. Our method outperforms BetaWrap, a state-of-the-art algorithm for predicting beta-helix folds, and HMMER, a general motif detection algorithm based on HMM, and has the additional advantage of general application to other protein folds. Applying our prediction model to the Uniprot Database, we identify previously unknown potential beta-helices.     

Review: conformation and folding of novel beta-structural elements in viral fiber proteins: the triple beta-spiral and triple beta-helix

Apart from alpha-helical coiled coils and the collagen triple helices, fibrous proteins can contain beta-structure in various conformations. Elongated enzymes such as pectate lyase and the bacteriophage P22 tailspike protein contain single-stranded beta-helices. Virus and bacteriophage fibers, which are often trimeric, have been shown to contain novel triple-stranded beta-structures such as the triple beta-spiral and the triple beta-helix. The conformation and folding of viral fibers containing beta-structure are discussed.     

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.     

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.     

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.     

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.     

Beta-helix core packing within the triple-stranded oligomerization domain of the P22 tailspike

A right-handed parallel beta-helix of 400 residues in 13 tightly packed coils is a major motif of the chains forming the trimeric P22 tailspike adhesin. The beta-helix domains of three identical subunits are side-by-side in the trimer and make predominantly hydrophilic inter-subunit contacts (Steinbacher S et al., 1994, Science 265:383-386). After the 13th coil the three individual beta-helices terminate and the chains wrap around each other to form three interdigitated beta-sheets organized into the walls of a triangular prism. The beta-strands then separate and form antiparallel beta-sheets, but still defining a triangular prism in which each side is a beta-sheet from a different subunit (Seckler R, 1998, J Struct Biol 122:216-222). The subunit interfaces are buried in the triangular core of the prism, which is densely packed with hydrophobic side chains from the three beta-sheets. Examination of this structure reveals that its packed core maintains the same pattern of interior packing found in the left-handed beta-helix, a single-chain structure. This packing is maintained in both the interdigitated parallel region of the prism and the following antiparallel sheet section. This oligomerization motif for the tailspike beta-helices presumably contributes to the very high thermal and detergent stability that is a property of the native tailspike adhesin.     

Selection of a representative set of structures from Brookhaven Protein Data Bank.

Reliable structural and statistical analyses of three dimensional protein structures should be based on unbiased data. The Protein Data Bank is highly redundant, containing several entries for identical or very similar sequences. A technique was developed for clustering the known structures based on their sequences and contents of alpha- and beta-structures. First, sequences were aligned pairwise. A representative sample of sequences was then obtained by grouping similar sequences together, and selecting a typical representative from each group. The similarity significance threshold needed in the clustering method was found by analyzing similarities of random sequences. Because three dimensional structures for proteins of same structural class are generally more conserved than their sequences, the proteins were clustered also according to their contents of secondary structural elements. The results of these clusterings indicate conservation of alpha- and beta-structures even when sequence similarity is relatively low. An unbiased sample of 103 high resolution structures, representing a wide variety of proteins, was chosen based on the suggestions made by the clustering algorithm. The proteins were divided into structural classes according to their contents and ratios of secondary structural elements. Previous classifications have suffered from subjective view of secondary structures, whereas here the classification was based on backbone geometry. The concise view lead to reclassification of some structures. The representative set of structures facilitates unbiased analyses of relationships between protein sequence, function, and structure as well as of structural characteristics.     

The interdigitated beta-helix domain of the P22 tailspike protein acts as a molecular clamp in trimer stabilization

The P22 tailspike adhesin is an elongated thermostable trimer resistant to protease digestion and to denaturation in sodium dodecyl sulfate. Monomeric, dimeric, and protrimeric folding and assembly intermediates lack this stability and are thermolabile. In the native trimer, three right-handed parallel beta-helices (residues 143-540), pack side-by-side around the three-fold axis. After residue 540, these single chain beta-helices terminate and residues 541-567 of the three polypeptide chains wrap around each other to form a three-stranded interdigitated beta-helix. Three mutants located in this region -- G546D, R563Q, and A575T -- blocked formation of native tailspike trimers, and accumulated soluble forms of the mutant polypeptide chains within cells. The substitutions R563Q and A575T appeared to prevent stable association of partially folded monomers. G546D, in the interdigitated region of the chain, blocked tailspike folding at the transition from the partially-folded protrimer to the native trimer. The protrimer-like species accumulating in the G546D mutant melted out at 42 degrees C and was trypsin and SDS sensitive. The G546D defect was not corrected by introduction of global suppressor mutations, which correct kinetic defects in beta-helix folding. The simplest interpretation of these results is that the very high thermostability (T(m) = 88 degrees C), protease and detergent resistance of the native tailspike acquired in the protrimer-to-trimer transition, depends on the formation of the three-stranded interdigitated region. This interdigitated beta-helix appears to function as a molecular clamp insuring thermostable subunit association in the native trimer.     

Wrap-and-Pack: a new paradigm for beta structural motif recognition with application to recognizing beta trefoils.

A method is presented that uses beta-strand interactions at both the sequence and the atomic level, to predict beta-structural motifs of protein sequences. A program called Wrap-and- Pack implements this method and is shown to recognize beta-trefoils, an important class of globular beta-structures, in the Protein Data Bank with 92% specificity and 92.3% sensitivity in cross-validation. It is demonstrated that Wrap-and-Pack learns each of the ten known SCOP beta-trefoil families, when trained primarily on beta-structures that are not beta-trefoils, together with three-dimensional structures of known beta-trefoils from outside the family. Wrap-and-Pack also predicts many proteins of unknown structure to be beta-trefoils. The computational method used here may generalize to other beta-structures for which strand topology and profiles of residue accessibility are well conserved.     

Molecular origin of polyglutamine aggregation in neurodegenerative diseases

Expansion of polyglutamine (polyQ) tracts in proteins results in protein aggregation and is associated with cell death in at least nine neurodegenerative diseases. Disease age of onset is correlated with the polyQ insert length above a critical value of 35-40 glutamines. The aggregation kinetics of isolated polyQ peptides in vitro also shows a similar critical-length dependence. While recent experimental work has provided considerable insights into polyQ aggregation, the molecular mechanism of aggregation is not well understood. Here, using computer simulations of isolated polyQ peptides, we show that a mechanism of aggregation is the conformational transition in a single polyQ peptide chain from random coil to a parallel beta-helix. This transition occurs selectively in peptides longer than 37 glutamines. In the beta-helices observed in simulations, all residues adopt beta-strand backbone dihedral angles, and the polypeptide chain coils around a central helical axis with 18.5 +/- 2 residues per turn. We also find that mutant polyQ peptides with proline-glycine inserts show formation of antiparallel beta-hairpins in their ground state, in agreement with experiments. The lower stability of mutant beta-helices explains their lower aggregation rates compared to wild type. Our results provide a molecular mechanism for polyQ-mediated aggregation.     

Do parallel beta-helix proteins have a unique fourier transform infrared spectrum?

Several polypeptides have been found to adopt an unusual domain structure known as the parallel beta-helix. These domains are characterized by parallel beta-strands, three of which form a single parallel beta-helix coil, and lead to long, extended beta-sheets. We have used ATR-FTIR (attenuated total reflectance-fourier transform infrared spectroscopy) to analyze the secondary structure of representative examples of this class of protein. Because the three-dimensional structures of parallel beta-helix proteins are unique, we initiated this study to determine if there was a corresponding unique FTIR signal associated with the parallel beta-helix conformation. Analysis of the amide I region, emanating from the carbonyl stretch vibration, reveals a strong absorbance band at 1638 cm(-1) in each of the parallel beta-helix proteins. This band is assigned to the parallel beta-sheet structure. However, components at this frequency are also commonly observed for beta-sheets in many classes of globular proteins. Thus we conclude that there is no unique infrared signature for parallel beta-helix structure. Additional contributions in the 1638 cm(-1) region, and at lower frequencies, were ascribed to hydrogen bonding between the coils in the loop/turn regions and amide side-chain interactions, respectively. A 13-residue peptide that forms fibrils and has been proposed to form beta-helical structure was also examined, and its FTIR spectrum was compared to that of the parallel beta-helix proteins.     

Beta-helix model for the filamentous haemagglutinin adhesin of Bordetella pertussis and related bacterial secretory proteins

Bordetella pertussis establishes infection by attaching to epithelial cells of the respiratory tract. One of its adhesins is filamentous haemagglutinin (FHA), a 500-A-long secreted protein that is rich in beta-structure and contains two regions, R1 and R2, of tandem 19-residue repeats. Two models have been proposed in which the central shaft is (i) a hairpin made up of a pairing of two long antiparallel beta-sheets; or (ii) a beta-helix in which the polypeptide chain is coiled to form three long parallel beta-sheets. We have analysed a truncated variant of FHA by electron microscopy (negative staining, shadowing and scanning transmission electron microscopy of unstained specimens): these observations support the latter model. Further support comes from detailed sequence analysis and molecular modelling studies. We applied a profile search method to the sequences adjacent to and between R1 and R2 and found additional "covert" copies of the same motifs that may be recognized in overt form in the R1 and R2 sequence repeats. Their total number is sufficient to support the tenet of the beta-helix model that the shaft domain--a 350 A rod--should consist of a continuous run of these motifs, apart from loop inserts. The N-terminus, which does not contain such repeats, was found to be weakly homologous to cyclodextrin transferase, a protein of known immunoglobulin-like structure. Drawing on crystal structures of known beta-helical proteins, we developed structural models of the coil motifs putatively formed by the R1 and R2 repeats. Finally, we applied the same profile search method to the sequence database and found several other proteins--all large secreted proteins of bacterial provenance--that have similar repeats and probably also similar structures.     

Nanostructure design using protein building blocks enhanced by conformationally constrained synthetic residues

Increasing efforts are being invested in the construction of nanostructures with desired shapes and physical and chemical properties. Our strategy involves nanostructure design using naturally occurring protein building blocks. Inspection of the protein structural database (PDB) reveals the richness of the conformations, shapes, and chemistries of proteins and their building blocks. To increase the population of the native fold in the selected building block, we mutate natural residues by engineered, constrained residues that restrict the conformational freedom at the targeted site and have favorable interactions, geometry, and size. Here, as a model system, we construct nanotubes using building blocks from left-handed beta-helices which are commonly occurring repeat protein architectures. We pick two-turn beta-helical segments, duplicate and stack them, and using all-atom molecular dynamics simulations (MD) with explicit solvent probe the structural stability of these nanotubular structures as indicated by their capacity to retain the initial organization and their conformational dynamics. Comparison of the results for the wild-type and mutated sequences shows that the introduction of the conformationally restricted 1-aminocyclopropanecarboxylic acid (Ac3c) residue in loop regions greatly enhances the stability of beta-helix nanotubes. The Ac3c geometrical confinement effect is sequence-specific and position-specific. The achievement of high stability of nanotubular structures originates not only from the reduction of mobility at the mutation site induced by Ac3c but also from stabilizing association forces between building blocks such as hydrogen bonds and hydrophobic contacts. For the selected synthetic residue, similar size, hydrophobicity, and backbone conformational tendencies are desirable as in the Ac3c.     

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.     

An asymmetric ion channel derived from gramicidin A. Synthesis, function and NMR structure

The biological ion channel gramicidin A (gA) was modified by synthetic means to obtain the tail-to-tail linked asymmetric gA-derived dimer compound 3. Single-channel current measurements for 3 in planar lipid bilayers exhibit an Eisenman I ion selectivity for alkali cations. The structural asymmetry does not lead to an observable functional asymmetry. The structure of 3 in solution without and with Cs cations was investigated by 1H-NMR spectroscopy. In CDCl3/CD3OH (1 : 1, v/v), 3 forms a mixture of double-stranded beta-helices. Upon addition of excess CsCl, the double-stranded species are converted completely into one new conformer: the right-handed single-stranded beta-helix. A combination of DQF-COSY and TOCSY was used for the assignment of the 1H-NMR spectrum of the Cs-3 complex in CDCl3/CD3OH (1 : 1, v/v). A total of 69 backbone, 27 long-range, and 64 side-chain distance restraints were obtained from NOESY together with 25 phi and 14 chi1 torsion angles obtained from coupling constants. These data were used as input for structure calculation with dyana built in sybyl 6.8. A final set of 11 structures with an average rmsd for the backbone of 0.45 A was obtained (PDB: 1TKQ). The structure of the Cs-3 complex in solution is equivalent to the bioactive channel conformation in the membrane environment.     

Cloning and characterization of the gene encoding endo-beta-1,3-glucanase from Arthrobacter sp. NHB-10

The gluA gene, encoding an endo-beta-1,3-glucanase from Arthrobacter sp. (strain NHB-10), was cloned and analyzed. The deduced endo-beta-1,3-glucanase amino acid sequence was 750 amino acids long and contained a 42 amino acid signal peptide with a mature protein of 708 amino acids. There was no similarity to known endo-beta-1,3-glucanases, but GluA was partially similar to two fungal exo-beta-1,3-glucanases in glycoside hydrolase (GH) family 55. Of five possible residues for catalysis and two motifs in two beta-helix heads of GH family 55, three residues and one motif were conserved in GluA, suggesting that GluA is the first bacterial endo-beta-1,3-glucanase in GH family 55. Significant similarity was also found to two proteins of unknown function from Streptomyces coelicolor A3(2) and S. avermitilis.     

The structure of the Haemophilus influenzae HMW1 pro-piece reveals a structural domain essential for bacterial two-partner secretion

In pathogenic Gram-negative bacteria, many virulence factors are secreted via the two-partner secretion pathway, which consists of an exoprotein called TpsA and a cognate outer membrane translocator called TpsB. The HMW1 and HMW2 adhesins are major virulence factors in nontypeable Haemophilus influenzae and are prototype two-partner secretion pathway exoproteins. A key step in the delivery of HMW1 and HMW2 to the bacterial surface involves targeting to the HMW1B and HMW2B outer membrane translocators by an N-terminal region called the secretion domain. Here we present the crystal structure at 1.92 A of the HMW1 pro-piece (HMW1-PP), a region that contains the HMW1 secretion domain and is cleaved and released during HMW1 secretion. Structural analysis of HMW1-PP revealed a right-handed beta-helix fold containing 12 complete parallel coils and one large extra-helical domain. Comparison of HMW1-PP and the Bordetella pertussis FHA secretion domain (Fha30) reveals limited amino acid homology but shared structural features, suggesting that diverse TpsA proteins have a common structural domain required for targeting to cognate TpsB proteins. Further comparison of HMW1-PP and Fha30 structures may provide insights into the keen specificity of TpsA-TpsB interactions.     

Biophysical characterisation reveals structural disorder in the nucleolar protein, Dribble

Dribble (DBE) is an essential protein in Drosophila that belongs to the evolutionarily conserved Krr1p protein family. Proteins in this family are localised in the cell nucleolus and are important for the processing of ribosomal RNAs. However, little is known about their structural and biophysical properties. We have expressed and purified full-length DBE protein from Escherichia coli. Consistent with the native role of DBE in RNA processing, recombinant DBE was shown to bind RNA homo-polymers in vitro. By bioinformatics, size-exclusion chromatography, equilibrium sedimentation analysis, controlled proteolysis, and a variety of spectroscopic techniques, we have found that DBE is a monomeric protein in solution containing both alpha- and beta-structures. Moreover, the structure of DBE is expanded and significantly disordered (approximately 45% disordered). Natively disordered proteins are thought to provide a disproportionately large surface area and structural plasticity for nucleic acid binding. We therefore propose that the presence of structural disorder is an important feature of DBE that facilitates the protein to interact with RNAs in the nucleolus.     

Plasticity and steric strain in a parallel beta-helix: rational mutations in the P22 tailspike protein

By means of genetic screens, a great number of mutations that affect the folding and stability of the tailspike protein from Salmonella phage P22 have been identified. Temperature-sensitive folding (tsf) mutations decrease folding yields at high temperature, but hardly affect thermal stability of the native trimeric structure when assembled at low temperature. Global suppressor (su) mutations mitigate this phenotype. Virtually all of these mutations are located in the central domain of tailspike, a large parallel beta-helix. We modified tailspike by rational single amino acid replacements at three sites in order to investigate the influence of mutations of two types: (1) mutations expected to cause a tsf phenotype by increasing the side-chain volume of a core residue, and (2) mutations in a similar structural context as two of the four known su mutations, which have been suggested to stabilize folding intermediates and the native structure by the release of backbone strain, an effect well known for residues that are primarily evolved for function and not for stability or folding of the protein. Analysis of folding yields, refolding kinetics and thermal denaturation kinetics in vitro show that the tsf phenotype can indeed be produced rationally by increasing the volume of side chains in the beta-helix core. The high-resolution crystal structure of mutant T326F proves that structural rearrangements only take place in the remarkably plastic lumen of the beta-helix, leaving the arrangement of the hydrogen-bonded backbone and thus the surface of the protein unaffected. This supports the notion that changes in the stability of an intermediate, in which the beta-helix domain is largely formed, are the essential mechanism by which tsf mutations affect tailspike folding. A rational design of su mutants, on the other hand, appears to be more difficult. The exchange of two residues in the active site expected to lead to a drastic release of steric strain neither enhanced the folding properties nor the stability of tailspike. Apparently, side-chain interactions in these cases overcompensate for backbone strain, illustrating the extreme optimization of the tailspike protein for conformational stability. The result exemplifies the view arising from the statistical analysis of the distribution of backbone dihedral angles in known three-dimensional protein structures that the adoption of straight phi/psi angles other than the most favorable ones is often caused by side-chain interactions. Proteins 2000;39:89-101.