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.     

Conversion of non-fibrillar beta-sheet oligomers into amyloid fibrils in Alzheimer's disease amyloid peptide aggregation

Abeta(1-40) is one of the main components of the fibrils found in amyloid plaques, a hallmark of brains affected by Alzheimer's disease. It is known that prior to the formation of amyloid fibrils in which the peptide adopts a well-ordered intermolecular beta-sheet structure, peptide monomers associate forming low and high molecular weight oligomers. These oligomers have been previously described in electron microscopy, AFM, and exclusion chromatography studies. Their specific secondary structures however, have not yet been well established. A major problem when comparing aggregation and secondary structure determinations in concentration-dependent processes such as amyloid aggregation is the different concentration range required in each type of experiment. In the present study we used the dye Thioflavin T (ThT), Fourier-transform infrared spectroscopy, and electron microscopy in order to structurally characterize the different aggregated species which form during the Abeta(1-40) fibril formation process. A unique sample containing 90microM peptide was used. The results show that oligomeric species which form during the lag phase of the aggregation kinetics are a mixture of unordered, helical, and intermolecular non-fibrillar beta-structures. The number of oligomers and the amount of non-fibrillar beta-structures grows throughout the lag phase and during the elongation phase these non-fibrillar beta-structures are transformed into fibrillar (amyloid) beta-structures, formed by association of high molecular weight intermediates.     

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.     

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.     

Synthesis and beta-conformation of copolypeptides with alternating hydrophilic and hydrophobic residues

Copolypeptides with alternating hydrophilic and hydrophobic residues were prepared, and their ability to form beta-structures in aqueous solutions was investigated by circular dichroism. Optically pure samples of poly (Lys-Leu-Lys-Leu) and poly (Leu-Glu-Leu-Glu), obtained via the 2-hydroxyphenyl esters, undergo a coil-to-beta transition in presence of salt. The beta-structures obtained under identical conditions with partially racemized samples of poly (Leu-Lys)Np and poly (Leu-Glu)Np, prepared by polycondensation of the corresponding dipeptide p-nitrophenyl esters, appear to be less regular. Non-alternating poly (Gly-Lys-Leu-Lys-Leu) does not form beta-structures in presence of NaCl as does alternating poly (Lys-Leu-Lys-Leu) indicating that the amino acid sequence can dramatically change the tendency to form beta-structures.     

Relation between the spatial structure and antigenic activity of trophoblast-specific beta1-glycoprotein

Temperature- and pH-dependence of spatial structure of a native trophoblast-specific beta-glycoprotein (TSG), its desialated and deglycosylated derivatives, as well as of a fragment obtained by partial acid hydrolysis on the temperature and pH variations has been studied using CD and UV spectroscopy. Within the range 45-50 degrees C a conformational transition of the protein moiety of TSG, leading to partially reversible alterations in tertiary and secondary structures of this glycoprotein after cooling the solution to 20 degrees C has been found out. The results of spectral studies of the spatial structure of the TSG protein component have been compared with the data on antigen activity of native, temperature- and pH-denaturated, desialated, and deglycosilated TSG. It has been concluded that the protein moiety of TSG consists mainly of beta-structures; the greater part of antigen determinants of TSG is topographic and belongs to the protein component of TSG, and only 15% of antigen determinants are not connected with TSG's spatial structure.     

The molecular structure and physical properties of elastin fibers as revealed by Raman microspectroscopy

Raman microspectroscopy has been used to investigate the structure of alpha-elastin and fibrous elastin from ligament and aorta, and to explore changes associated with mechanical strain and temperature. Although no vibrational modes associated with cross-linking of the fibers could be identified, the secondary structure of dehydrated fibrous elastin was significantly different from alpha-elastin. The former differed from previous experimental measurements, but was close to the theoretical predictions with 36% beta-structures, 46% unordered, and 18% alpha-helix. alpha-Elastin contained 29% beta-structures, 53% unordered, and 18% alpha-helix. In nuchal fibers the amide I mode was polarized, consistent with the peptide bond. Strains of up to 60% in ligament fiber bundles resulted in no significant shifts in peak position or in secondary structure. Polarization measurements revealed that the peptide bonds and several side chains re-orientated closer to the fiber axis. Heating nuchal fibers to 60 degrees C to increase the energetic component of the elasticity was associated with a 30% increase in the proportion of beta-structures in the amide I band, a 50% increase in the amide III band, and a 50% reduction in the signal from bound water. (c) 2008 Wiley Periodicals, Inc. Biopolymers 89: 931-940, 2008.This article was originally published online as an accepted preprint. The "Published Online" date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com.     

Determination of the secondary structure of proteins from circular dichroism spectra. IV. Contribution of aromatic amino acid residues into circular dichroism spectra of proteins in the peptide region

A new method for determination of the secondary protein structure from the CD spectra taking into account the contribution of aromatic amino acid residues is proposed. New proteins reference CD spectra for five secondary structures (alpha-helices, antiparallel and parallel beta-structures, beta-bends and irregular form) without contribution of aromatic residues are obtained. By means of this new method the secondary structure of sixteen different proteins was analysed. There is a good correlation of these results with the X-ray data.     

Conformation and aggregation of M13 coat protein studied by molecular dynamics.

Molecular dynamics (MD) simulations are performed on M13 coat protein, a small membrane protein for which both alpha- and beta-structures have been suggested. The simulations are started from initial conformations that are either monomers or dimers of alpha-helices or U-shaped beta-sheets. The lipid bilayer is represented by a hydrophobic potential. The results are analyzed in terms of stability, energy and secondary structure. The U-shaped beta-structure changes from a planar to a twisted form with larger twist for the monomer than the dimer. The beta-sheet is much more flexible than the alpha-helix as monitored by the root mean square (rms) fluctuations of the C alpha atoms. A comparison of the energies after 100 ps MD simulation shows that of the monomers, the alpha-helix has the lowest energy. The energy difference between alpha- and beta-structures decreases from 266 kJ/mol to 148 kJ/mol, when going from monomers to dimers. It is expected that this difference will decrease with higher aggregation numbers.     

Dual nature of the infectious prion protein revealed by high pressure

Crude brain homogenates of terminally diseased hamsters infected with the 263K strain of scrapie (PrP(Sc)) and purified prion fibrils were heated or pressurized at 800 megapascals and 60 degrees C for 2 h in different buffers and in water. Prion proteins (PrP) were analyzed for their proteinase K resistance in immunoblots and for their infectivity in hamster bioassays. A notable decrease in the proteinase K resistance of unpurified prion proteins, probably because of pressure-induced changes in the protein conformation of native PrP(Sc) or the N-truncated PrP-(27-30), could be demonstrated when pressurized at initially neutral conditions in several buffers and in water but not in a slightly acidic pH. A subsequent 6-7 log(10) reduction of infectious units/g in phosphate-buffered saline buffer, pH 7.4, was found. The proteinase K-resistant core was also not detectable after purification of prions extracted from pressurized samples, confirming pressure effects at the level of the secondary structure of prion proteins. However, opposite results were found after pressurizing purified prions, arguing for the existence of pressure-sensitive beta-structures (PrP(Sc)(DeltaPsen)) and extremely pressure-resistant beta-structures (PrP(Sc)(DeltaPres)). Remarkably, after the first centrifugation step at 540,000 x g during isolation, prions remained proteinase K-resistant when pressurized in all tested buffers and in water. It is known that purified fibrils retain infectivity, but the isolated protein (full and N-truncated) behaved differently from native PrP(Sc) under pressure, suggesting a kind of semicrystalline polymer structure.     

PoPMuSiC, an algorithm for predicting protein mutant stability changes: application to prion proteins

A novel tool for computer-aided design of single-site mutations in proteins and peptides is presented. It proceeds by performing in silico all possible point mutations in a given protein or protein region and estimating the stability changes with linear combinations of database-derived potentials, whose coefficients depend on the solvent accessibility of the mutated residues. Upon completion, it yields a list of the most stabilizing, destabilizing or neutral mutations. This tool is applied to mouse, hamster and human prion proteins to identify the point mutations that are the most likely to stabilize their cellular form. The selected mutations are essentially located in the second helix, which presents an intrinsic preference to form beta-structures, with the best mutations being T183-->F, T192-->A and Q186-->A. The T183 mutation is predicted to be by far the most stabilizing one, but should be considered with care as it blocks the glycosylation of N181 and this blockade is known to favor the cellular to scrapie conversion. Furthermore, following the hypothesis that the first helix might induce the formation of hydrophilic beta-aggregates, several mutations that are neutral with respect to the structure's stability but improve the helix hydrophobicity are selected, among which is E146-->L. These mutations are intended as good candidates to undergo experimental tests.     

Conformational intermediate of the amyloidogenic protein beta 2-microglobulin at neutral pH.

Aggregation and fibrillation of beta(2)-microglobulin are hallmarks of dialysis-related amyloidosis. We characterize perturbations of the native conformation of beta(2)-microglobulin that may precede fibril formation. For a beta(2)-microglobulin variant cleaved at lysine 58, we show using capillary electrophoresis that two conformers spontaneously exist in aqueous buffers at neutral pH. Upon treatment of wild-type beta(2)-microglobulin with acetonitrile or trifluoroethanol, two conformations were also observed. These conformations were in equilibrium dependent on the sample temperature and the percentage of organic solvent present. Circular dichroism showed a loss of beta-structures and gain of alpha-helices. Reversal to the native conformation occurred when removing the organics. Affinity capillary electrophoresis experiments showed increased specific interactions of the nonnative beta(2)-microglobulin conformation with the dyes 8-anilino-1-naphthalene sulfonic acid and Congo red. The observations may relate to early folding events prior to amyloid fibrillation and facilitate the development of methods to detect and inhibit pro-amyloid protein and peptide conformations.     

Purification and structural study of the beta form of human cAMP-dependent protein kinase inhibitor.

The beta form of human cAMP-dependent protein kinase inhibitor (human PKIbeta), a novel heat-stable protein, was isolated with high yield using a bacterial expression system. Assays of PKI activity demonstrated that purified PKIbeta inhibits the catalytic subunit of cAMP-dependent protein kinase. FTIR, Raman spectroscopy and CD experiments implied that human PKIbeta contained only small amounts of alpha-helix and beta-structures, but large amounts of random coil and turn structures, which may explain its high thermostability. The details of its conformational changes in response to heat were studied by CD experiments for the first time, revealing that the protein unfolded at high temperature and refolded when decreased to room temperature.     

β-structures of polypeptides with L- and D-residues

The possible formation of beta-structures from polypeptide chains with L-and D-Residues randomly distributed was statistically analyzed within the frame of two hypotheses. Firstly, only those segments containing residues of identical chirality can associate to form antiparallel beta-structures, and secondly these segments must have a minimum length. The influence of different factors was examined: initial ratio of the L-and D-monomer, minimum length required for the segments to be incorporated into beta-sheets, average length of the peptide molecules, and stereoselectivity in the course of the polymerization process. The results show that in all cases nuclei of beta-sheets surrounded by random coil segments are formed, the optical activity of which very increases to purity when the initial ratio of monomers deviates from the racemic mixture. This suggests experiments to enrich the system in one enantiomer. Comparison is made with the corresponding behavior and properties of the alpha-helical structure.     

Determination of the secondary prothrombin structure under different degrees of citraconylation by means of ring dichroism]

A process of spontaneous activation of bovine prothrombin under acylation of cytacone anhydride is studied by means of ring dichroism and disc electrophoresis in sodium dodecylsulphate. Ring dichroism data have shown that native prothrombin contains 14% of alpha-helical structures, 36% of beta-structures and 50% of random coil. Cytraconylation results in a fragmentation of native prothrombin chain and is accompanies by the decrease of alpha-helical regions in the fragments formed. Poly-L-lysine, modified by citracone anydride, does not form alpha-helical regions.     

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.     

A method of calculating the discrete secondary structures of globular proteins

The model of formation of alpha-helices and beta-structures determined by joint action of the three elements: N-terminal, internal and C-terminal fragments are presented. Algorithm for calculation of their localization in a given amino acid sequence was constructed on the base of this model. The preference of the fragments of the amino acid sequence to a definite type of the secondary structure was estimated on the base of corresponding average values of linear discriminant functions dsk (s = alpha, beta, k = N, in, C). The latter were constructed in the previous paper on the base of the revealed significant characteristics. These integral characteristics are used for calculating the localisation of discrete secondary structures. The total prediction for 3 states (alpha, beta, c) given 71% correctly predicted residues (for 4 states alpha, beta, c, t) 62% for the training set, consisting of 72 proteins. For the control set (15 proteins) the accuracy of prediction is about 65%. The essential advantages of this method are: 1) the possibility to localize the discrete secondary structures; 2) the high accuracy of prediction of long secondary structures (for alpha-helices approximately 90%, for beta-structures approximately 80%), which is important for the determination of the protein folding. The influence of mutation on the secondary structure of proteins was investigated. The anormally high stability of the secondary structures of immunoglobulins to mutations was revealed. This probably results from the selection during evolution of such variants of amino acid sequences, which are able to provide the functional variability of antigenic determinants, but keep invariant the tertially structure of protein.     

The structure of yeast enolase at 2.25-A resolution. An 8-fold beta + alpha-barrel with a novel beta beta alpha alpha (beta alpha)6 topology

The three-dimensional structure of yeast enolase has been determined by the multiple isomorphous replacement method followed by the solvent flattening technique. A polypeptide model, corresponding with the known amino acid sequence, has been fitted to the electron density map. Crystallographic restrained least-squares refinement of the model without solvent gave R = 20.0% for 6-2.25-A resolution with good geometry. A model with 182 water molecules and 1 sulfate which is still being refined has presently R = 17.0%. The molecule is a dimer with subunits related by 2-fold crystallographic symmetry. The subunit has dimensions 60 X 55 X 45 A and is built from two domains. The smaller N-terminal domain has an alpha + beta structure based on a three-stranded antiparallel meander and four helices. The main domain is an 8-fold beta + alpha-barrel. The enolase barrel is, however, different from the triose phosphate isomerase barrel; its topology is beta beta alpha alpha (beta alpha)6 rather than (beta alpha)8 as found in triose phosphate isomerase. The inner beta-barrel is not entirely parallel, the second strand is antiparallel to the other strands, and the direction of the first helix is also reversed with respect to the other helices. This supports the hypothesis that some enzymes evolved independently producing the stable structure of beta alpha barrels with either enolase or triose phosphate isomerase topology. The active site of enolase is located at the carboxylic end of the barrel. A fragment of the N-terminal domain and two long loops protruding from the barrel domain form a wide crevice leading to the active site region. Asp246, Glu295, and Asp320 are the ligands of the conformational cation. Other residues in the active site region are Glu168, Asp321, Lys345, and Lys396.     

Supramolecular amyloid-like assembly of the polypeptide sequence coded by exon 30 of human tropoelastin

Elastin is known to self-aggregate in twisted-rope filaments. However, an ultrastructural organization different from the fibrils typical of elastin, but rather similar to those shown by amyloid networks, is shown by the polypeptide sequence encoded by exon 30 of human tropoelastin. To better understand the molecular properties of this sequence to give amyloid fibers, we used CD, NMR, and FTIR (Fourier transform infrared spectroscopy) to identify the structural characteristics of the peptide. In this study, we have demonstrated, by FTIR, that antiparallel beta-sheet conformation is predominant in the exon 30 fibers. These physical-chemical studies were combined with transmission electron microscopy and atomic force microscopy to analyze the supramolecular structure of the self-assembled aggregate. These studies show the presence of fibrils that interact side-by-side probably originating from an extensive self-interaction of elemental cross beta-structures. Similar sequences, of the general type XGGZG(X, Z = V, L, A, I), are widely found in many proteins such as collagens IV and XVII, major prion protein precursor, amyloid beta A4 precursor protein-binding family, etc., thus suggesting that this sequence could be involved in contributing to the self-assembly of amyloid fibers even in other proteins.     

Cloning and sequence analysis of the human liver rhodanese: comparison with the bovine and chicken enzymes

The cDNA for the human rhodanese (thiosulfate: cyanide sulfurtransferase, EC 2.8.1.1), a nuclearly encoded protein of the mitochondrial matrix, was isolated from a human fetal liver cDNA library. Nucleotide sequence revealed an open reading frame coding for a polypeptide of 295 amino acids, which presented a 57% and 58% identity with the bovine and avian rhodanese, respectively. The analysis of the 5'-ends of the coding region gave no evidence for the presence of a cleavable signal sequence as found in other mitochondrial proteins. A comparison with two available amino acid sequences (cow and chicken) showed that sequence similarity is not restricted to the alpha-helices and beta-structures motifs which are remarkably superimposable in the two halves of bovine rhodanese, but extends to adjacent regions.