A novel approach to predicting protein structural classes in a (20–1)-D amino acid composition space

The development of prediction methods based on statistical theory generally consists of two parts: one is focused on the exploration of new algorithms, and the other on the improvement of a training database. The current study is devoted to improving the prediction of protein structural classes from both of the two aspects. To explore a new algorithm, a method has been developed that makes allowance for taking into account the coupling effect among different amino acid components of a protein by a covariance matrix. To improve the training database, the selection of proteins is carried out so that they have (1) as many non-homologous structures as possible, and (2) a good quality of structure. Thus, 129 representative proteins are selected. They are classified into 30 α, 30 β, 30 α + β, 30 α/β, and 9 ζ (irregular) proteins according to a new criterion that better reflects the feature of the structural classes concerned. The average accuracy of prediction by the current method for the 4 × 30 regular proteins is 99.2%, and that for 64 independent testing proteins not included in the training database is 95.3%. To further validate its efficiency, a jackknife analysis has been performed for the current method as well as the previous ones, and the results are also much in favor of the current method. To complete the mathematical basis, a theorem is presented and proved in Appendix A that is instructive for understanding the novel method at a deeper level. © 1995 Wiley-Liss, Inc.     

Prediction of protein folding class from amino acid composition

An empirical relation between the amino acid composition and three-dimensional folding pattern of several classes of proteins has been determined. Computer simulated neural networks have been used to assign proteins to one of the following classes based on their amino acid composition and size: (1) 4α-helical bundles, (2) parallel (α/β)₈ barrels, (3) nucleotide binding fold, (4) immunoglobulin fold, or (5) none of these. Networks trained on the known crystal structures as well as sequences of closely related proteins are shown to correctly predict folding classes of proteins not represented in the training set with an average accuracy of 87%. Other folding motifs can easily be added to the prediction scheme once larger databases become available. Analysis of the neural network weights reveals that amino acids favoring prediction of a folding class are usually over represented in that class and amino acids with unfavorable weights are underrepresented in composition. The neural networks utilize combinations of these multiple small variations in amino acid composition in order to make a prediction. The favorably weighted amino acids in a given class also form the most intramolecular interactions with other residues in proteins of that class. A detailed examination of the contacts of these amino acids reveals some general patterns that may help stabilize each folding class. © 1993 Wiley-Liss, Inc.     

A weighting method for predicting protein structural class from amino acid composition.

A protein is generally classified into one of the following four structural classes: all alpha, all beta, alpha+beta and alpha/beta. In this paper, based on the weighting to the 20 constituent amino acids, a new method is proposed for predicting the structural class of a protein according to its amino acid composition. The 20 weighting parameters, which reflect the different properties of the 20 constituent amino acids, have been obtained from a training set of proteins through the linear-programming approach. The rate of correct prediction for a training set of proteins by means of the new method was 100%, whereas the highest rate of previous methods was 82.8%. Furthermore, the results showed that the more numerous training proteins, the more effective the new method.     

A new method to model membrane protein structure based on silent amino acid substitutions.

The importance of accurately modeling membrane proteins cannot be overstated, in lieu of the difficulties in solving their structures experimentally. Often, however, modeling procedures (e.g., global searching molecular dynamics) generate several possible candidates rather then pointing to a single model. Herein we present a new approach to select among candidate models based on the general hypothesis that silent amino acid substitutions, present in variants identified from evolutionary conservation data or mutagenesis analysis, do not affect the stability of a native structure but may destabilize the non-native structures also found. The proof of this hypothesis has been tested on the alpha-helical transmembrane domains of two homodimers, human glycophorin A and human CD3-zeta, a component of the T-cell receptor. For both proteins, only one structure was identified using all the variants. For glycophorin A, this structure is virtually identical to the structure determined experimentally by NMR. We present a model for the transmembrane domain of CD3-zeta that is consistent with predictions based on mutagenesis, homology modeling, and the presence of a disulfide bond. Our experiments suggest that this method allows the prediction of transmembrane domain structure based only on widely available evolutionary conservation data.     

Detection of gene duplication signals of Ig folds from their amino acid sequences

Protein folds may evolve from short peptide ancestors via gene duplication and fusion. For proteins with internal structural symmetry, this means that their sequences should be made up of identical repeats. However, many of these repeat signals can only be seen at the structural level yet. Motivated by the fact that proteins may have similar structures if their sequences have more than 25% identical amino acids, we suggest a method to detect the sequence repeats of proteins directly from their sequences. Using this method, we show that the internal repetitions of the immunoglobulin folds could be identified directly at the sequence level.     

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.     

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.     

Functional and structural characterization of a thermostable acetyl esterase from Thermotoga maritima

TM0077 from Thermotoga maritima is a member of the carbohydrate esterase family 7 and is active on a variety of acetylated compounds, including cephalosporin C. TM0077 esterase activity is confined to short‐chain acyl esters (C2–C3), and is optimal around 100°C and pH 7.5. The positional specificity of TM0077 was investigated using 4‐nitrophenyl‐β‐D ‐xylopyranoside monoacetates as substrates in a β‐xylosidase‐coupled assay. TM0077 hydrolyzes acetate at positions 2, 3, and 4 with equal efficiency. No activity was detected on xylan or acetylated xylan, which implies that TM0077 is an acetyl esterase and not an acetyl xylan esterase as currently annotated. Selenomethionine‐substituted and native structures of TM0077 were determined at 2.1 and 2.5 Å resolution, respectively, revealing a classic α/β‐hydrolase fold. TM0077 assembles into a doughnut‐shaped hexamer with small tunnels on either side leading to an inner cavity, which contains the six catalytic centers. Structures of TM0077 with covalently bound phenylmethylsulfonyl fluoride and paraoxon were determined to 2.4 and 2.1 Å, respectively, and confirmed that both inhibitors bind covalently to the catalytic serine (Ser188). Upon binding of inhibitor, the catalytic serine adopts an altered conformation, as observed in other esterase and lipases, and supports a previously proposed catalytic mechanism in which Ser hydroxyl rotation prevents reversal of the reaction and allows access of a water molecule for completion of the reaction. Proteins 2012. © 2012 Wiley Periodicals, Inc.     

A new approach to predicting protein folding types

A new method is proposed for predicting the folding type of a protein according to its amino acid composition based on the following physical picture: (1) a protein is characterized as a vector of 20-dimensional space, in which its 20 components are defined by the compositions of its 20 amino acids; and (2) the similarity of two proteins is proportional to the mutual projection of their characterized vectors, and hence inversely proportional to the size of their correlation angle. Thus, the prediction is performed by calculating the correlation angles of the vector for the predicted protein with a set of standard vectors representing the norms of four protein folding types (i.e., alla, all ,a+, anda/). In comparison with the existing methods, the new method has the merits of yielding a higher rate of correct prediction, displaying a more intuitive physical picture, and being convenient in application. For instance, in predicting the 64 proteins in the development set based on which the standard vectors are derived, the average accuracy rate is 83.6%, which is higher than that obtained for the same set of proteins by any of the existing methods. The average accuracy predicted for an independent set of 35 proteins of known X-ray structure is 91.4%, which is significantly higher than any of the reported accuracies so far, implying that the new method is of great value in practical application. All of these have demonstrated that the new method as proposed in this paper is characterized by an improved feature in both self-consistency and extrapolating-effectiveness.On sabbatical leave from Department of Physics, Tianjin University, Tianjin, China.     

Probing structural requirements of fMLP receptor: on the size of the hydrophobic pocket corresponding to residue 2 of the tripeptide.

The conformationally constrained f-L-Met-Ac(n)c-L-Phe-OMe (n = 4,9-12) tripeptides, analogues of the chemoattractant f-L-Met-L-Leu-L-Phe-OH, were synthesized in solution by classical methods and fully characterized. These compounds and the published f-L-Met-Xxx-L-Phe-OMe (Xxx = Aib and Ac(n)c where n = 3, 5-8) analogues were compared to determine the combined effect of backbone preferred conformation and side-chain bulkiness at position 2 on the relation of 3D-structure to biological activity. A conformational study of all the analogues was performed in solution by FT-IR absorption and 1H-NMR techniques. In parallel, each peptide was tested for its ability to induce chemotaxis, superoxide anion production and lysozyme secretion from human neutrophils. The biological and conformational data are discussed in relation to the proposed model of the chemotactic receptor on neutrophils, in particular of the hydrophobic pocket accommodating residue 2 of the tripeptide.     

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.     

Backbone makes a significant contribution to the electrostatics of alpha/beta-barrel proteins.

The electrostatic properties of seven alpha/beta-barrel enzymes selected from different evolutionary families were studied: triose phosphate isomerase, fructose-1,6-bisphosphate aldolase, pyruvate kinase, mandelate racemase, trimethylamine dehydrogenase, glycolate oxidase, and narbonin, a protein without any known enzymatic activity. The backbone of the alpha/beta-barrel has a distinct electrostatic field pattern, which is dipolar along the barrel axis. When the side chains are included in the calculations the general effect is to modulate the electrostatic pattern so that the electrostatic field is generally enhanced and is focused into a specific area near the active site. We use the electrostatic flux through a square surface near the active site to gauge the functionally relevant magnitude of the electrostatic field. The calculations reveal that in six out of the seven cases the backbone itself contributes greater than 45% of the total flux. The substantial electrostatic contribution of the backbone correlates with the known preference of alpha/beta-barrel enzymes for negatively charged substrates.     

Kinetic investigation of the effect of the amino acid side chains in the selective acidolysis of N-acyl-N,alpha,alpha-trialkyl glycine amides.

Accurate kinetic measurements of the rate constants for the acidolysis of five N-acetyl-N-(4-methoxybenzyl)-alpha,alpha-trialkyl glycine cyclohexyl amides in TFA were performed at 25 degrees C and the reactions monitored by HPLC. In all cases the results were consistent with a first order behaviour with respect to the substrate. No direct correlation was obtained with these data between the rate constant values and structure, but a good correlation coefficient was obtained when a multiple regression analysis was applied by taking advantage of a Taft equation using appropriate polar and steric substituent parameters. In a plot of the values observed for log k against those calculated by this equation all five points fell very close to the line of perfect correlation. The calculated sensitivity coefficients to polar and steric contributions were used to discuss the experimental results and showed that the acidolysis were comparatively less affected by steric effects than expected.     

Mannose-specific lectins bind alpha-2-macroglobulin and an unknown protein from human plasma.

GNA, the mannose-specific lectin from Galanthus nivalis was confirmed to bind alpha-2-macroglobulin (A2M) but another protein was copurified with A2M from total human plasma. A total of 23 other lectins with diverse specificities were tested for reaction with human A2M and with three other members of the A2M family. NPA, a mannose-specific lectin isolated from Narcissus pseudonarcissus bulbs, and RSA, the Rhizoctonia solani agglutinin, were selected for further testing. For isolation of A2M, immobilized NPA was superior to GNA because its binding capacity was an order of magnitude higher. The specificity of these lectins must be very similar however, because the same unknown plasma protein was also bound by NPA. A2M and the unknown protein must share a unique mannose carbohydrate structure not present in any other human plasma protein. The copurified protein subunit size of 185 kDa is very similar to that of A2M, but the native molecular mass of 350 kDa indicated a noncovalent homodimer structure. Together with the acid isoelectric point this is not typical for any known plasma protein nor for any unidentified spot on the two-dimensional map of human plasma proteins. No immunological reaction with available antisera was evident. A specific antiserum raised to the unknown protein demonstrated its presence in all human plasma samples examined. The N-terminal residue was blocked, whereas internal protein sequences obtained after CNBr fragmentation and proteolysis were not homologous to any known protein sequence. These data demonstrate that this protein is unknown and not a proteinase inhibitor of the A2M family.     

Reversible unfolding and refolding behavior of a monomeric aldolase from Staphylococcus aureus.

Thermal and GdmCl-induced unfolding transitions of aldolase from Staphylococcus aureus are reversible under a variety of solvent conditions. Analysis of the transitions reveals that no partially folded intermediates can be detected under equilibrium conditions. The stability of the enzyme is very low with a delta G0 value of -9 +/- 2 kJ/mol at 20 degrees C. The kinetics of unfolding and refolding of aldolase are complex and comprise at least one fast and two slow reactions. This complexity arises from prolyl isomerization reactions in the unfolded chain, which are kinetically coupled to the actual folding reaction. Comparison with model calculations shows that at least two prolyl peptide bonds give rise to the observed slow folding reactions of aldolase and that all of the involved bonds are presumably in the trans conformation in the native state. The rate constant of the actual folding reaction is fast with a relaxation time of about 15 s at the midpoint of the folding transition at 15 degrees C. The data presented on the folding and stability of aldolase are comparable to the properties of much smaller proteins. This might be connected with the simple and highly repetitive tertiary structure pattern of the enzyme, which belongs to the group of alpha/beta barrel proteins.     

Stimulation of hyaluronan synthesis by interleukin-1beta involves activation of protein kinase C betaII in fibroblasts from patients with Graves' ophthalmopathy

Hyaluronan accumulation in the retroorbital connective tissue is one of the pathological features of Graves' ophthalmopathy. Interleukin-1beta (IL-1beta) is known to stimulate hyaluronan synthesis in orbital fibroblasts. In the present study, the intracellular signal transduction pathways involved in this stimulatory effect were investigated in cultured human retroorbital fibroblasts from patients with Graves' ophthalmopathy. IL-1beta-induced hyaluronan synthesis was significantly inhibited by pretreatment of the cells with two protein kinase C (PKC) inhibitors, chlerythrine chloride and H-7. In addition, treatment with phorbol 12-myristate 13-acetate (PMA), a direct PKC activator, also resulted in increased hyaluronan production. IL-1beta- or PMA-stimulated hyaluronan synthesis was blocked by the protein synthesis inhibitor, cycloheximide. Moreover, the intracellular Ca(2+) concentration of the orbital fibroblasts was also involved in the IL-1beta induced transduction pathway, the effect being completely inhibited by BAPTA, an internal calcium chelator. In addition, A23187, a calcium ionophore, increased hyaluronan synthesis in unstimulated cells. These results suggest that the Ca(2+)-dependent PKC signal transduction pathway plays an important role in the IL-1beta-induced hyaluronan synthesis. Moreover, IL-1beta treatment resulted in increased PKC activity and the rapid translocation of PKC betaII from the cytoplasm to the plasma membrane. These results indicate that cytosolic Ca(2+) and PKC betaII are involved in IL-1beta-induced hyaluronan synthesis in cultured orbital fibroblasts from patients with Graves' ophthalmopathy.     

Beta-turn propensities as paradigms for the analysis of structural motifs to engineer protein stability.

The thermodynamic stability of a protein provides an experimental metric for the relationship of protein sequence and native structure. We have investigated an approach based on an analysis of the structural database for stability engineering of an immunoglobulin variable domain. The most frequently occurring residues in specific positions of beta-turn motifs were predicted to increase the folding stability of mutants that were constructed by site-directed mutagenesis. Even in positions in which different residues are conserved in immunoglobulin sequences, the predictions were confirmed. Frequently, mutants with increased beta-turn propensities display increased folding cooperativities, suggesting pronounced effects on the unfolded state independent of the expected effect on conformational entropy. We conclude that structural motifs with predominantly local interactions can serve as templates with which patterns of sequence preferences can be extracted from the database of protein structures. Such preferences can predict the stability effects of mutations for protein engineering and design.     

Hematopoietic contribution to skeletal muscle regeneration in acid alpha-glucosidase knockout mice.

Recent studies have shown that cells from bone marrow (BM) can give rise to differentiated skeletal muscle fibers. However, the mechanisms and identities of the cell types involved remain unknown. We performed BM transplantation in acid alpha-glucosidase (GAA) knockout mice, a model of glycogen storage disease type II, and our observations suggested that the BM cells contribute to skeletal muscle fiber formation. Furthermore, we showed that most CD45+:Sca1+ cells have a donor character in regenerating muscle of recipient mice. Based on these findings, CD45+:Sca1+ cells were sorted from regenerating muscles. The cell number was increased with granulocyte colony-stimulating factor after cardiotoxin injury, and the cells were transplanted directly into the tibialis anterior (TA) muscles of GAA knockout mice. Sections of the TA muscles stained with anti-laminin-alpha2 antibody showed that the number of CD45+:Sca1+ cells contributing to muscle fiber formation and glycogen levels were decreased in transplanted muscles. Our results indicated that hematopoietic stem cells, such as CD45+:Sca1+ cells, are involved in skeletal muscle regeneration.     

Solution structure of the hypothetical protein Mth677 from Methanobacterium thermoautotrophicum: a novel alpha+beta fold

The structure of Mth677, a hypothetical protein from Methanobacterium thermoautotrophicum (Mth), has been determined by using heteronuclear nuclear magnetic resonance (NMR) methods on a double-labeled (15)N-(13)C sample. Mth677 adopts a novel alpha+beta fold, consisting of two alpha-helices (one N terminal and one C terminal) packed on the same side of a central beta-hairpin. This structure is likely shared by its three orthologs, detected in three other Archaebacteria. There are no clear features in the sequences of these proteins or in the genome organization of Mth to make a reliable functional assignment to this protein. However, the structural similarity to Escherichia coli MinE, the protein which controls that division occurs at the midcell site, lends support to the proposal that Mth677 might be, in Mth, the counterpart of the topological specificity domain of MinE in E. coli.