Determination of protein fold class from Raman or Raman optical activity spectra using random forests

Knowledge of the fold class of a protein is valuable because fold class gives an indication of protein function and evolution. Fold class can be accurately determined from a crystal structure or NMR structure, though these methods are expensive, time-consuming, and inapplicable to all proteins. In contrast, vibrational spectra [infra-red, Raman, or Raman optical activity (ROA)] are rapidly obtained for proteins under wide range of biological molecules under diverse experimental and physiological conditions. Here, we show that the fold class of a protein can be determined from Raman or ROA spectra by converting a spectrum into data of 10 cm⁻¹ bin widths and applying the random forest machine learning algorithm. Spectral data from 605 and 1785 cm⁻¹ were analyzed, as well as the amide I, II, and III regions in isolation and in combination. ROA amide II and III data gave the best performance, with 33 of 44 proteins assigned to one of the correct four top-level structural classification of proteins (SCOP) fold class (all α, all β, α and β, and disordered). The method also shows which spectral regions are most valuable in assigning fold class.     

Engineering of betabellin-15D: a 64 residue beta sheet protein that forms long narrow multimeric fibrils.

The betabellin target structure is a beta-sandwich protein consisting of two 32 residue beta-sheets packed against one another by interaction of their hydrophobic faces. The 32 residue chain of betabellin-15S (HSLTAKIpkLTFSIAphTYTCAV pkYTAKVSH, where p=DPro, k=DLys, and h=DHis) did not fold in water at pH 6.5. Air oxidation of betabellin-15S provided betabellin-15D, the 64 residue disulfide bridged two-chain molecule, which also remained unfolded in water at pH 6.5. By circular dichroic spectropolarimetry, the extent of beta structure observed for betabellin-15D increased with the pH and ionic strength of the solution and the betabellin-15D concentration. By electron microscopy, in 5.0 mM MOPS and 0.25 M NaCl at pH 6.9, betabellin-15D formed long narrow multimeric fibrils. A molecular model was constructed to show that the dimensions of these betabellin-15D fibrils are consistent with a single row of beta-sandwich molecules joined by multiple intersheet H-bonds.     

Measuring cooperativity in the formation of a three-stranded beta sheet (double hairpin)

Recently validated chemical shift measures of hairpin structuring have been applied to a series of turn mutants of the Schenck-Gellman three-strand beta-sheet model with the aim of measuring the entropic advantage associated with aligning an additional strand onto an existing hairpin versus aligning the same two strands in an initial hairpin formation. In a four-state analysis (unfolded, 2 single hairpins, and the double hairpin fold in equilibrium) a cooperativity index can be defined as the factor by which the equilibrium constant for hairpin formation is improved when one strand is prestructured. This cooperativity index is 2.7 +/- 0.7 for hairpin formation about a stable D-Pro-Gly turn locus and increases to 7.6 +/- 1.2 for an Asn-Gly turn locus. The latter corresponds to a cooperativity induced DeltaDeltaG increment of 4.9 kJ/mol for the folding of a hairpin. Although larger than previous experimental measures of folding cooperativity in three-stranded sheets, the magnitude of this effect (which is considerably less than the TDeltaDeltaS expectation for prestructuring three or more beta-strand residue sites) likely reflects the intrinsic preference of these designed sequences for extended conformations. If similar or larger effects apply to protein beta-sheet folding, it is not surprising that particularly favorable hairpin alignments serve as nucleation sites in protein folding pathways.     

The importance of surface loops for stabilizing an eightfold beta alpha barrel protein.

An important step in understanding how a protein folds is to determine those regions of the sequence that are critical to both its stability and its folding pathway. We chose phosphoribosyl anthranilate isomerase from Escherichia coli, which is a monomeric representative of the (beta alpha)8 barrel family of proteins, to construct a variant that carries an internal tandem duplication of the fifth beta alpha module. This (beta alpha)9 variant was enzymically active and therefore must have a wild-type (beta alpha)8 core. It had a choice a priori to fold to three different folding frames, which are distinguished by carrying the duplicated segment as an insert into one out of three different loops. Steady-state kinetic constants, the fluorescence properties of a crucial tryptophan residue, and limited proteolysis showed that the stable (beta alpha)9 variant carries the insertion between beta-strand 5 and alpha-helix 5. This preference can be explained by the important role of loops between alpha helices and beta strands in stabilizing the structure of the enzyme.     

The ensemble folding kinetics of the FBP28 WW domain revealed by an all-atom Monte Carlo simulation in a knowledge-based potential

In this work, we apply a detailed all‐atom model with a transferable knowledge‐based potential to study the folding kinetics of Formin‐Binding protein, FBP28, which is a canonical three‐stranded β‐sheet WW domain. Replica exchange Monte Carlo simulations starting from random coils find native‐like (Cα RMSD of 2.68 Å) lowest energy structure. We also study the folding kinetics of FBP28 WW domain by performing a large number of ab initio Monte Carlo folding simulations. Using these trajectories, we examine the order of formation of two β‐hairpins, the folding mechanism of each individual β‐hairpin, and transition state ensemble (TSE) of FBP28 WW domain and compare our results with experimental data and previous computational studies. To obtain detailed structural information on the folding dynamics viewed as an ensemble process, we perform a clustering analysis procedure based on graph theory. Further, a rigorous P_fold analysis is used to obtain representative samples of the TSEs showing good quantitative agreement between experimental and simulated Φ values. Our analysis shows that the turn structure between first and second β strands is a partially stable structural motif that gets formed before entering the TSE in FBP28 WW domain and there exist two major pathways for the folding of FBP28 WW domain, which differ in the order and mechanism of hairpin formation. Proteins 2011. © 2011 Wiley‐Liss, Inc.     

Effects of detergents on the secondary structures of prion protein peptides as studied by CD spectroscopy.

Pathogenic prion proteins (PrP(Sc)) are thought to be produced by alpha-helical to beta-sheet conformational changes in the normal cellular prion proteins (PrP(C)) located solely in the caveolar compartments. In order to inquire into the possible conformational changes due to the influences of hydrophobic environments within caveolae, the secondary structures of prion protein peptides were studied in various kinds of detergents by CD spectra. The peptides studied were PrP(129-154) and PrP(192-213); the former is supposed to assume beta-sheets and the latter alpha-helices, in PrP(Sc). The secondary structure analyses for the CD spectra revealed that in buffer solutions, both PrP(129-154) and PrP(192-213) mainly adopted random-coils (approximately 60%), followed by beta-sheets (30%-40%). PrP(129-154) showed no changes in the secondary structures even in various kinds of detergents such as octyl-beta-D-glucopyranoside (OG), octy-beta-D-maltopyranoside (OM). sodium dodecyl sulfate (SDS), Zwittergent 3-14 (ZW) and dodecylphosphocholine (DPC). In contrast, PrP(192-213) changed its secondary structure depending on the concentration of the detergents. SDS, ZW, OG and OM increased the alpha-helical content, and decreased the beta-sheet and random-coil contents. DPC also increased the alpha-helical content, but to a lesser extent than did SDS, ZW, OG or OM. These results indicate that PrP(129-154) has a propensity to adopt predominantly beta-sheets. On the other hand, PrP(192-213) has a rather fickle propensity and varies its secondary structure depending on the environmental conditions. It is considered that the hydrophobic environments provided by these detergents may mimic those provided by gangliosides in caveolae, the head groups of which consist of oligosaccharide chains containing sialic acids. It is concluded that PrP(C) could be converted into a nascent PrP(Sc) having a transient PrP(Sc) like structureunder the hydrophobic environments produced by gangliosides.     

Ranking the factors that contribute to protein beta-sheet folding

The formation of beta-sheet domains in proteins involves five energetically important factors: the formation of networks of hydrogen bonds and hydrophobic faces, and the residue propensities, or preferences, to be found at the edges of the beta-sheet, to adopt the extended conformation, and to make contact with other residues. These relative energy contributions define a potential energy function. Here, we show how optimizing this potential energy function reveals the formation of hydrophobic faces as the utmost factor. The potential energy function was optimized to minimize the Z-scores of the native topologies among the exhaustive sets of over 400 different beta-sheets. These results corroborate with experimental data that showed the environment of a protein is an important modulator of beta-sheet folding. The contact propensities were found to be the least important, which could explain the poor predictive power of beta-strand alignment methods based on pair-wise contact matrices.     

Quantifying structure-function uncertainty: a graph theoretical exploration into the origins and limitations of protein annotation

Since the advent of investigations into structural genomics, research has focused on correctly identifying domain boundaries, as well as domain similarities and differences in the context of their evolutionary relationships. As the science of structural genomics ramps up adding more and more information into the databanks, questions about the accuracy and completeness of our classification and annotation systems appear on the forefront of this research. A central question of paramount importance is how structural similarity relates to functional similarity. Here, we begin to rigorously and quantitatively answer these questions by first exploring the consensus between the most common protein domain structure annotation databases CATH, SCOP and FSSP. Each of these databases explores the evolutionary relationships between protein domains using a combination of automatic and manual, structural and functional, continuous and discrete similarity measures. In order to examine the issue of consensus thoroughly, we build a generalized graph out of each of these databases and hierarchically cluster these graphs at interval thresholds. We then employ a distance measure to find regions of greatest overlap. Using this procedure we were able not only to enumerate the level of consensus between the different annotation systems, but also to define the graph-theoretical origins behind the annotation schema of class, family and superfamily by observing that the same thresholds that define the best consensus regions between FSSP, SCOP and CATH correspond to distinct, non-random phase-transitions in the structure comparison graph itself. To investigate the correspondence in divergence between structure and function further, we introduce a measure of functional entropy that calculates divergence in function space. First, we use this measure to calculate the general correlation between structural homology and functional proximity. We extend this analysis further by quantitatively calculating the average amount of functional information gained from our understanding of structural distance and the corollary inherent uncertainty that represents the theoretical limit of our ability to infer function from structural similarity. Finally we show how our measure of functional "entropy" translates into a more intuitive concept of functional annotation into similarity EC classes.     

An optimization approach to predicting protein structural class from amino acid composition.

Proteins are generally classified into four structural classes: all-alpha proteins, all-beta proteins, alpha + beta proteins, and alpha/beta proteins. In this article, a protein is expressed as a vector of 20-dimensional space, in which its 20 components are defined by the composition of its 20 amino acids. Based on this, a new method, the so-called maximum component coefficient method, is proposed for predicting the structural class of a protein according to its amino acid composition. In comparison with the existing methods, the new method yields a higher general accuracy of prediction. Especially for the all-alpha proteins, the rate of correct prediction obtained by the new method is much higher than that by any of the existing methods. For instance, for the 19 all-alpha proteins investigated previously by P.Y. Chou, the rate of correct prediction by means of his method was 84.2%, but the correct rate when predicted with the new method would be 100%! Furthermore, the new method is characterized by an explicable physical picture. This is reflected by the process in which the vector representing a protein to be predicted is decomposed into four component vectors, each of which corresponds to one of the norms of the four protein structural classes.     

A helix-turn-strand structural motif common in alpha-beta proteins.

By exhaustive structural comparisons, we have found that about one-third of the alpha-helix-turn-beta-strand polypeptides in alpha-beta barrel domains share a common structural motif. The chief characteristics of this motif are that first, the geometry of the turn between the alpha-helix and the beta-strand is somewhat constrained, and second, the beta-strand contains a hydrophobic patch that fits into a hydrophobic pocket on the alpha-helix. The geometry of the turn does not seem to be a major determinant of the alpha-beta unit, because the turns vary in length from four to six residues. However, the motif does not occur when there are few constraints on the geometry of the turn-for instance, when the turns between the alpha-helix and the beta-strands are very long. It also occurs much less frequently in flat-sheet alpha-beta proteins, where the topology is much less regular and the amount of twist on the sheet varies considerably more than in the barrel proteins. The motif may be one of the basic building blocks from which alpha-beta barrels are constructed.     

Analysis of interactive packing of secondary structural elements in alpha/beta units in proteins

An alpha-helix and a beta-strand are said to be interactively packed if at least one residue in each of the secondary structural elements loses 10% of its solvent accessible contact area on association with the other secondary structural element. An analysis of all such 5,975 nonidentical alpha/beta units in protein structures, defined at < or = 2.5 A resolution, shows that the interaxial distance between the alpha-helix and the beta-strand is linearly correlated with the residue-dependent function, log[(V/nda)/n-int], where V is the volume of amino acid residues in the packing interface, nda is the normalized difference in solvent accessible contact area of the residues in packed and unpacked secondary structural elements, and n-int is the number of residues in the packing interface. The beta-sheet unit (beta u), defined as a pair of adjacent parallel or antiparallel hydrogen-bonded beta-strands, packing with an alpha-helix shows a better correlation between the interaxial distance and log(V/nda) for the residues in the packing interface. This packing relationship is shown to be useful in the prediction of interaxial distances in alpha/beta units using the interacting residue information of equivalent alpha/beta units of homologous proteins. It is, therefore, of value in comparative modeling of protein structures.     

Crystal-state conformation of Calpha,alpha-dialkylated peptides containing chiral beta-homo-residues.

Secondary structure formation and stability are essential features in the knowledge of complex folding topology of biomolecules. To better understand the relationships between preferred conformations and functional properties of beta-homo-amino acids, the synthesis and conformational characterization by X-ray diffraction analysis of peptides containing conformationally constrained Calpha,alpha-dialkylated amino acid residues, such as alpha-aminoisobutyric acid or 1-aminocyclohexane-1-carboxylic acid and a single beta-homoamino acid, differently displaced along the peptide sequence have been carried out. The peptides investigated are: Boc-betaHLeu-(Ac6c)2-OMe, Boc-Ac6c-betaHLeu-(Ac6c)2-OMe and Boc-betaHVal-(Aib)5-OtBu, together with the C-protected beta-homo-residue HCl.H-betaHVal-OMe. The results indicate that the insertion of a betaH-residue at position 1 or 2 of peptides containing strong helix-inducing, bulky Calpha,alpha-disubstituted amino acid residues does not induce any specific conformational preferences. In the crystal state, most of the NH groups of beta-homo residues of tri- and tetrapeptides are not involved in intramolecular hydrogen bonds, thus failing to achieve helical structures similar to those of peptides exclusively constituted of Calpha,alpha-disubstituted amino acid residues. However, by repeating the structural motifs observed in the molecules investigated, a beta-pleated sheet secondary structure, and a new helical structure, named (14/15)-helix, were generated, corresponding to calculated minimum-energy conformations. Our findings, as well as literature data, strongly indicate that conformations of betaH-residues, with the micro torsion angle equal to -60 degrees, are very unlikely.     

Recurrent alpha beta loop structures in TIM barrel motifs show a distinct pattern of conserved structural features.

A systematic survey of seven parallel alpha/beta barrel protein domains, based on exhaustive structural comparisons, reveals that a sizable proportion of the alpha beta loops in these proteins--20 out of a total of 49--belong to either one of two loop types previously described by Thornton and co-workers. Six loops are of the alpha beta 1 type, with one residue between the alpha-helix and beta-strand, and 13 are of the alpha beta 3 type, with three residues between the helix and the strand. Protein fragments embedding the identified loops, and termed alpha beta connections since they contain parts of the flanking helix and strand, have been analyzed in detail revealing that each type of connection has a distinct set of conserved structural features. The orientation of the beta-strand relative to the helix and loop portions is different owing to a very localized difference in backbone conformation. In alpha beta 1 connections, the chain enters the beta-strand via a residue adopting an extended conformation, while in alpha beta 3 it does so via a residue in a near alpha-helical conformation. Other conserved structural features include distinct patterns of side chain orientation relative to the beta-sheet surface and of main chain H-bonds in the loop and the beta-strand moieties. Significant differences also occur in packing interactions of conserved hydrophobic residues situated in the last turn of the helix. Yet the alpha-helix surface of both types of connections adopts similar orientations relative to the barrel sheet surface. Our results suggest furthermore that conserved hydrophobic residues along the sequence of the connections, may be correlated more with specific patterns of interactions made with neighboring helices and sheet strands than with helix/strand packing within the connection itself. A number of intriguing observations are also made on the distribution of the identified alpha beta 1 and alpha beta 3 loops within the alpha/beta-barrel motifs. They often occur adjacent to each other; alpha beta 3 loops invariably involve even numbered beta-strands, while alpha beta 1 loops involve preferentially odd beta-strands; all the analyzed proteins contain at least one alpha beta 3 loop in the first half of the eightfold alpha/beta barrel. Possible origins of all these observations, and their relevance to the stability and folding of parallel alpha/beta barrel motifs are discussed.     

The prion protein unstructured N-terminal region is a broad-spectrum molecular sensor with diverse and contrasting potential functions

The physiological function of the prion protein (PrP(C) ) and its conversion into its infectious form (PrP(Sc) ) are central issues to understanding the pathogenesis of prion diseases. The N-terminal moiety of PrP(C) (NH(2) -PrP(C) ) is an unstructured region with the characteristic of interacting with a broad range of partners. These interactions endow PrP(C) with multifunctional and sometimes contrasting capabilities, including neuroprotection and neurotoxicity. Recently, binding of β-sheet rich conformers to NH(2) -PrP(C) demonstrated a probable neurotoxic function for PrP(C) in Alzheimer's disease. NH(2) -PrP(C) also enhances the propagation of prions in vivo and is the target of the most potent antiprion compounds. Another level of complexity is provided by endoproteolysis and release of most of NH(2) -PrP(C) into the extracellular space. Further studies will be necessary to understand how NH(2) -PrP(C) regulates the physiological function of PrP(C) and how it is involved in the corruption of its normal function in diseases.     

Analysis of hydrophobicity in the alpha and beta chemokine families and its relevance to dimerization.

The chemokine family of chemotactic cytokines plays a key role in orchestrating the immune response. The family has been divided into 2 subfamilies, alpha and beta, based on the spacing of the first 2 cysteine residues, function, and chromosomal location. Members within each subfamily have 25-70% sequence identity, whereas the amino acid identity between members of the 2 subfamilies ranges from 20 to 40%. A quantitative analysis of the hydrophobic properties of 11 alpha and 9 beta chemokine sequences, based on the coordinates of the prototypic alpha and beta chemokines, interleukin-8 (IL-8), and human macrophage inflammatory protein-1 beta (hMIP-1 beta), respectively, is presented. The monomers of the alpha and beta chemokines have their strongest core hydrophobic cluster at equivalent positions, consistent with their similar tertiary structures. In contrast, the pattern of monomer surface hydrophobicity between the alpha and beta chemokines differs in a manner that is fully consistent with the observed differences in quaternary structure. The most hydrophobic surface clusters on the monomer subunits are located in very different regions of the alpha and beta chemokines and comprise in each case the amino acids that are buried at the interface of their respective dimers. The theoretical analysis of hydrophobicity strongly supports the hypothesis that the distinct dimers observed for IL-8 and hMIP-1 beta are preserved for all the alpha and beta chemokines, respectively. This provides a rational explanation for the lack of receptor crossbinding and reactivity between the alpha and beta chemokine subfamilies.     

Disruption of the beta-sheet structure of a protected pentapeptide, related to the beta-amyloid sequence 17-21, induced by a single, helicogenic C(alpha)-tetrasubstituted alpha-amino acid.

Fifteen years ago it was shown that an alpha-aminoisobutyric acid (Aib) residue is significantly more effective than an L-Pro or a D-amino acid residue in inducing beta-sheet disruption in short model peptides. As this secondary structure element is known to play a crucial role in the neuropathology of Alzheimer's disease, it was decided to check the effect of Aib (and other selected, helix inducer, C(alpha)-tetrasubstituted alpha-amino acids) on the beta-sheet conformation adopted by a protected pentapeptide related to the sequence 17-21 of the beta-amyloid peptide. By use of FT-IR absorption and 1H NMR techniques it was found that the strong self-association characterizing the pentapeptide molecules in weakly polar organic solvents is completely abolished by replacing a single residue with Aib or one of its congeners.     

Conformational analysis of endomorphin-2 by molecular dynamics methods

Endomorphin-2 (EM2, H-Tyr-Pro-Phe-Phe-NH(2)) is a highly potent and selective mu-opioid receptor agonist. A conformational analysis of EM2 was carried out by simulated annealing (SA) and molecular dynamics (MD) methods. Molecular modeling was conducted on both neutral (N-terminal NH(2)) and charged (N-terminal NH(3) (+)) molecules. Based on the results of NMR investigations showing an equilibrium mixture of cis and trans Tyr(1)-Pro(2) peptide bonds for EM2 in solution, simulations were performed with restrained cis-Pro and trans-Pro peptide bonds, too. A separate SA study with unrestrained Pro peptide bonds was also conducted. Preferred conformational states are presented in Ramachandran plots. The g(+), g(-), and trans populations of the aromatic amino acid residue side chains were determined in chi(1) space. The distances between the N-terminal N atom and the other backbone N and O atoms, and the distances between the centers of the aromatic rings and the Pro(2) ring, were determined. The energy distribution of the structures obtained by SA was calculated. The preferred secondary structural elements were different kinds of beta-turns, an inverse gamma-turn located in the N-terminal region, and regular and inverse gamma-turns located in the C-terminal region. These turns were stabilized by intramolecular H-bonds and bifurcated H-bonds.     

Identification,classification, and analysis of beta-bulges in proteins.

A beta-bulge is a region of irregularity in a beta-sheet involving two beta-strands. It usually involves two or more residues in the bulged strand opposite to a single residue on the adjacent strand. These irregularities in beta-sheets were identified and classified automatically, extending the definition of beta-bulges given by Richardson et al. (Richardson, J.S., Getzoff, E.D., & Richardson, D.C., 1978, Proc. Natl. Acad. Sci. USA 75, 2574-2578). A set of 182 protein chains (170 proteins) was used, and a total of 362 bulges were extracted. Five types of beta-bulges were found: classic, G1, wide, bent, and special. Their characteristic amino acid preferences were found for most classes of bulges. Basically, bulges occur frequently in proteins; on average there are more than two bulges per protein. In general, beta-bulges produce two main changes in the structure of a beta-sheet: (1) disrupt the normal alternation of side-chain direction; (2) accentuate the twist of the sheet, altering the direction of the surrounding strands.     

Dimerization of beta B2-crystallin: the role of the linker peptide and the N- and C-terminal extensions.

beta B2- and gamma B-crystallins of vertebrate eye lens are 2-domain proteins in which each domain consists of 2 Greek key motifs connected by a linker peptide. Although the folding topologies of beta B2- and gamma B-domains are very similar, gamma B-crystallin is always monomeric, whereas beta B2-crystallin associates to homodimers. It has been suggested that the linker or the protruding N- and C-terminal arms of beta B2-crystallin (not present in gamma B) are a necessary requirement for this association. In order to investigate the role of these segments for dimerization, we constructed two beta B2 mutants. In the first mutant, the linker peptide was replaced with the one from gamma B (beta B2 gamma L). In the second mutant, the N- and C-terminal arms of 15- and 12-residues length were deleted (beta B2 delta NC). The beta B2 gamma L mutant is monomeric, whereas the beta B2 delta NC mutant forms dimers and tetramers that cannot be interconverted without denaturation. The spectral properties of the beta B2 mutants, as well as their stabilities against denaturants, resemble those of wild-type beta B2-crystallin, thus indicating that the overall peptide fold of the subunits is not changed significantly. We conclude that the peptide linker in beta B2-crystallin is necessary for dimerization, whereas the N- and C-terminal arms appear to be involved in preventing the formation of higher homo-oligomers.