Multidomain organization of eukaryotic guanine nucleotide exchange translation initiation factor eIF-2B subunits revealed by analysis of conserved sequence motifs.

Computer-assisted analysis of amino acid sequences using methods for database screening with individual sequences and with multiple alignment blocks reveals a complex multidomain organization of yeast proteins GCD6 and GCD1, and mammalian homolog of GCD6-subunits of the eukaryotic translation initiation factor eIF-2B involved in GDP/GTP exchange on eIF-2. It is shown that these proteins contain a putative nucleotide-binding domain related to a variety of nucleotidyltransferases, most of which are involved in nucleoside diphosphate-sugar formation in bacteria. Three conserved motifs, one of which appears to be a variant of the phosphate-binding site (P-loop) and another that may be considered a specific version of the Mg(2+)-binding site of NTP-utilizing enzymes, were identified in the nucleotidyltransferase-related domain. Together with the third unique motif adjacent to the the P-loop, these motifs comprise the signature of a new superfamily of nucleotide-binding domains. A domain consisting of hexapeptide amino acid repeats with a periodic distribution of bulky hydrophobic residues (isoleucine patch), which previously have been identified in bacterial acetyltransferases, is located toward the C-terminus from the nucleotidyltransferase-related domain. Finally, at the very C-termini of GCD6, eIF-2B epsilon, and two other eukaryotic translation initiation factors, eIF-4 gamma and eIF-5, there is a previously undetected, conserved domain. It is hypothesized that the nucleotidyltransferase-related domain is directly involved in the GDP/GTP exchange, whereas the C-terminal conserved domain may be involved in the interaction of eIF-2B, eIF-4 gamma, and eIF-5 with eIF-2.     

Construction of a dictionary of sequence motifs that characterize groups of related proteins.

An automatic procedure is proposed to identify, from the protein sequence database, conserved amino acid patterns (or sequence motifs) that are exclusive to a group of functionally related proteins. This procedure is applied to the PIR database and a dictionary of sequence motifs that relate to specific superfamilies constructed. The motifs have a practical relevance in identifying the membership of specific superfamilies without the need to perform sequence database searches in 20% of newly determined sequences. The sequence motifs identified represent functionally important sites on protein molecules. When multiple blocks exist in a single motif they are often close together in the 3-D structure. Furthermore, occasionally these motif blocks were found to be split by introns when the correlation with exon structures was examined.     

The PAS fold. A redefinition of the PAS domain based upon structural prediction.

In the postgenomic era it is essential that protein sequences are annotated correctly in order to help in the assignment of their putative functions. Over 1300 proteins in current protein sequence databases are predicted to contain a PAS domain based upon amino acid sequence alignments. One of the problems with the current annotation of the PAS domain is that this domain exhibits limited similarity at the amino acid sequence level. It is therefore essential, when using proteins with low-sequence similarities, to apply profile hidden Markov model searches for the PAS domain-containing proteins, as for the PFAM database. From recent 3D X-ray and NMR structures, however, PAS domains appear to have a conserved 3D fold as shown here by structural alignment of the six representative 3D-structures from the PDB database. Large-scale modelling of the PAS sequences from the PFAM database against the 3D-structures of these six structural prototypes was performed. All 3D models generated (> 5700) were evaluated using prosaii. We conclude from our large-scale modelling studies that the PAS and PAC motifs (which are separately defined in the PFAM database) are directly linked and that these two motifs form the PAS fold. The existing subdivision in PAS and PAC motifs, as used by the PFAM and SMART databases, appears to be caused by major differences in sequences in the region connecting these two motifs. This region, as has been shown by Gardner and coworkers for human PAS kinase (Amezcua, C.A., Harper, S.M., Rutter, J. & Gardner, K.H. (2002) Structure 10, 1349-1361, [1]), is very flexible and adopts different conformations depending on the bound ligand. Some PAS sequences present in the PFAM database did not produce a good structural model, even after realignment using a structure-based alignment method, suggesting that these representatives are unlikely to have a fold resembling any of the structural prototypes of the PAS domain superfamily.     

Sequence motifs specific for cytosine methyltransferases

Using a new alignment method, the sequences of 13 m5C methyltransferases (MTases) have been examined. Five extremely well-conserved blocks of sequence have been detected and have been used as fixed points for the alignment of the 13 sequences. Following this initial alignment, five further blocks of similarity have been identified to give a total of ten recognizable blocks of sequence homology that are all arranged in a common order. The structures of these MTases consist of a variable-length N-terminal arm followed by eight well-conserved blocks each separated by small variable-length regions. A large variable-length segment of 90 to 270 amino acids (aa) then follows. After this are two blocks, and a variable-length C-terminal segment completes the sequence. Within the final alignment, 20 aa in the protein sequences, and 86 nucleotides in the nucleotide sequences are invariant. The strongest conservation is found in proximity to a suspected functional site that contains the dipeptide proline-cysteine. Consensus patterns can be defined for the five best conserved blocks and, when used as search motifs, are able to clearly distinguish between the m5C MTases and all other identified proteins in the PIR database. This suggests they may be of use in identifying putative MTases among protein sequences of unknown function.     

Ancestral gene fusion in cellobiose dehydrogenases reflects a specific evolution of GMC oxidoreductases in fungi

Cellobiose dehydrogenases (CDHs) are extracellular hemoflavoenzymes that are thought to be involved in the degradation of two of the most abundant biopolymers in the biosphere, cellulose and lignin. To date, these enzymes, consisting of a cytochrome domain and a flavin domain, have been detected and sequenced exclusively in the kingdom of fungi. Independent phylogenetic analyses of two distinct domains of CDH genes reveal that they evolved in parallel as fused genes. Whereas the cytochrome domains are unique sequence motifs, the flavin domains clearly belong to the glucose-methanol-choline (GMC) oxidoreductase family--an evolution line of widespread flavoproteins extending from the Archae to higher eukaryotes. The most probable unrooted phylogenetic tree obtained from our analysis of 52 selected GMC members reveals five principal evolutionary branches: cellobiose dehydrogenase, cholesterol oxidase (COX), hydroxynitrile lyase, alcohol oxidase (AOX)/glucose oxidase (GOX)/choline dehydrogenase, and a branch of dehydrogenases with various specificities containing also an Archaeon open reading frame (ORF). Cellobiose dehydrogenases cluster with cholesterol oxidases and the clade of various specificities, whereas hydroxynitrile lyases are closely related to glucose oxidases, alcohol oxidases, and choline dehydrogenases. The results indicate that the evolutionary line from a primordial GMC flavoprotein to extant cellobiose dehydrogenases was augmented after an early acquisition of the cytochrome domain to form two distinct branches for basidiomycetes and ascomycetes. One ascomycetous evolutionary line of CDHs has acquired a carbohydrate-binding module (CBM) of type 1, the sequence of which is similar to that of corresponding domains in several glycosidases. This is the first attempt towards a comprehensive phylogenetic analysis of cellobiose dehydrogenases.     

Extracellular oxidative systems of the lignin-degrading Basidiomycete Phanerochaete chrysosporium

The US Department of Energy has assembled a high quality draft genome of Phanerochaete chrysosporium, a white rot Basidiomycete capable of completely degrading all major components of plant cell walls including cellulose, hemicellulose and lignin. Hundreds of sequences are predicted to encode extracellular enzymes including an impressive number of oxidative enzymes potentially involved in lignocellulose degradation. Herein, we summarize the number, organization, and expression of genes encoding peroxidases, copper radical oxidases, FAD-dependent oxidases, and multicopper oxidases. Possibly relevant to extracellular oxidative systems are genes involved in posttranslational processes and a large number of hypothetical proteins.     

Construction of a non-redundant human SH2 domain database

Domain database is essential for domain property research. Eliminating redundant information in database query is very important for database quality. Here we report the manual construction of a non-redundant human SH2 domain database. There are 119 human SH2 domains in 110 SH2-containing proteins. Human SH2s were aligned with ClustalX, and a homologous tree was generated. In this tree, proteins with similar known function were classified into the same group. Some proteins in the same group have been reported to have similar binding motifs experimentally. The tree might provide clues about possible functions of hypothetical proteins for further experimental verification.     

Eukaryotic translation elongation factor 1γ contains a glutathione transferase domain—Study of a diverse,ancient protein super family using motif search and structural modeling

Using computer methods for multiple alignment, sequence motif search, and tertiary structure modeling, we show that eukaryotic translation elongation factor 1γ (EF1γ) contains an N-terminal domain related to class θ glutathione S-transferases (GST). GST-like proteins related to class θ comprise a large group including, in addition to typical GSTs and EF1γ, stress-induced proteins from bacteria and plants, bacterial reductive dehalogenases and β-etherases, and several uncharacterized proteins. These proteins share 2 conserved sequence motifs with GSTs of other classes (α, μ, and π). Tertiary structure modeling showed that in spite of the relatively low sequence similarity, the GST-related domain of EF1γ is likely to form a fold very similar to that in the known structures of class α, μ, and π GSTs. One of the conserved motifs is implicated in glutathione binding, whereas the other motif probably is involved in maintaining the proper conformation of the GST domain. We predict that the GST-like domain in EF1γ is enzymatically active and that to exhibit GST activity, EF1γ has to form homodimers. The GST activity may be involved in the regulation of the assembly of multisubunit complexes containing EF1 and aminoacyl-tRNA synthetases by shifting the balance between glutathione, disulfide glutathione, thiol groups of cysteines, and protein disulfide bonds. The GST domain is a widespread, conserved enzymatic module that may be covalently or noncovalently complexed with other proteins. Regulation of protein assembly and folding may be 1 of the functions of GST.     

Polyamine catabolism: target for antiproliferative therapies in animals and stress tolerance strategies in plants

Metabolism of polyamines spermidine and spermine, and their diamine precursor, putrescine, has been a target for antineoplastic therapy since these naturally occurring alkyl amines were found essential for normal mammalian cell growth. Intracellular polyamine concentrations are maintained at a cell type-specific set point through the coordinated and highly regulated interplay between biosynthesis, transport, and catabolism. A correlation between regulation of cell proliferation and polyamine metabolism is described. In particular, polyamine catabolism involves copper-containing amine oxidases and FAD-dependent polyamine oxidases. Several studies showed an important role of these enzymes in several developmental and disease-related processes in both animals and plants through a control on polyamine homeostasis in response to normal cellular signals, drug treatment, environmental and/or cellular stressors. The production of toxic aldehydes and reactive oxygen species, H(2)O(2) in particular, by these oxidases using extracellular and intracellular polyamines as substrates, suggests a mechanism by which the oxidases can be exploited as antineoplastic drug targets. This minireview summarizes recent advances on the physiological roles of polyamine catabolism in animals and plants in an attempt to highlight differences and similarities that may contribute to determine in detail the underlined mechanisms involved. This information could be useful in evaluating the possibility of this metabolic pathway as a target for new antiproliferative therapies in animals and stress tolerance strategies in plants.     

Protein structure mining using a structural alphabet

We present a comprehensive evaluation of a new structure mining method called PB-ALIGN. It is based on the encoding of protein structure as 1D sequence of a combination of 16 short structural motifs or protein blocks (PBs). PBs are short motifs capable of representing most of the local structural features of a protein backbone. Using derived PB substitution matrix and simple dynamic programming algorithm, PB sequences are aligned the same way amino acid sequences to yield structure alignment. PBs are short motifs capable of representing most of the local structural features of a protein backbone. Alignment of these local features as sequence of symbols enables fast detection of structural similarities between two proteins. Ability of the method to characterize and align regions beyond regular secondary structures, for example, N and C caps of helix and loops connecting regular structures, puts it a step ahead of existing methods, which strongly rely on secondary structure elements. PB-ALIGN achieved efficiency of 85% in extracting true fold from a large database of 7259 SCOP domains and was successful in 82% cases to identify true super-family members. On comparison to 13 existing structure comparison/mining methods, PB-ALIGN emerged as the best on general ability test dataset and was at par with methods like YAKUSA and CE on nontrivial test dataset. Furthermore, the proposed method performed well when compared to flexible structure alignment method like FATCAT and outperforms in processing speed (less than 45 s per database scan). This work also establishes a reliable cut-off value for the demarcation of similar folds. It finally shows that global alignment scores of unrelated structures using PBs follow an extreme value distribution. PB-ALIGN is freely available on web server called Protein Block Expert (PBE) at http://bioinformatics.univ-reunion.fr/PBE/.     

家蚕醛氧化酶基因（BmAOXs）的鉴定与表达分析

醛氧化酶（aldehyde oxidases, AO, EC 1.2.3.1）是属于钼-黄素酶（molybdo-flavoenzymes, MFEs）家族的一类蛋白酶。为了探讨该酶在家蚕Bombyx mori中的功能, 本研究对家蚕醛氧化酶基因（BmAOXs）家族进行了鉴定和分析。以其他物种AO基因序列检索家蚕全基因组数据库, 获得了8个BmAOX候选基因, 均具有醛氧化酶保守的功能域。进化分析表明, BmAOX与其他昆虫AO聚为一簇。RT-PCR分析结果显示: BmAOX1, BmAOX2, BmAOX3, BmAOX5具有很强的组织特异性; 而BmAOX4, BmAOX6, BmAOX7, BmAOX8则在蛹和成虫的多个组织中均有表达, 提示它们在家蚕生理代谢活动中起重要作用。利用Native PAGE和活性染色方法, 对BmAOX编码的蛋白进行检测, 结果表明: 家蚕蛹中有5种有活性的醛氧化酶, 而成虫中有6种, 各组织中均有有活性的BmAOX, 只是种类和活性水平有所不同。本研究结果为深入探讨BmAOX家族的功能奠定了基础。     

PRINTS--a database of protein motif fingerprints.

PRINTS is a compendium of protein motif 'fingerprints'. A fingerprint is defined as a group of motifs excised from conserved regions of a sequence alignment, whose diagnostic power or potency is refined by iterative databasescanning (in this case the OWL composite sequence database). Generally, the motifs do not overlap, but are separated along a sequence, though they may be contiguous in 3D-space. The use of groups of independent, linearly- or spatially-distinct motifs allows protein folds and functionalities to be characterised more flexibly and powerfully than conventional single-component patterns or regular expressions. The current version of the database contains 200 entries (encoding 950 motifs), covering a wide range of globular and membrane proteins, modular polypeptides, and so on. The growth of the databaseis influenced by a number of factors; e.g. the use of multiple motifs; the maximisation of sequence information through iterative database scanning; and the fact that the database searched is a large composite. The information contained within PRINTS is distinct from, but complementary to the consensus expressions stored in the widely-used PROSITE dictionary of patterns.     

O-glycosylation in Toxoplasma gondii: identification and analysis of a family of UDP-GalNAc:polypeptide N-acetylgalactosaminyltransferases

The initiation of mucin-type O-glycosylation is catalysed by a family of UDP-GalNAc:polypeptide N-acetylgalactosaminyltransferases (EC 2.4.1.41). These enzymes are responsible for the transfer of N-acetylgalactosamine from the nucleotide sugar donor, UDP-GalNAc, to the hydroxyl group on specific serine or threonine residues in acceptor proteins. By screening a Toxoplasma gondii cDNA library, three distinct isoforms of the ppGalNAc-T gene family were cloned. Two additional isoforms were identified and partially cloned following analysis of the T. gondii genome sequence database. All of the cloned and identified ppGalNAc-T's are type II membrane proteins that share up to 50% amino acid sequence identity within the conserved catalytic domain. They each contain an N-terminal cytoplasmic domain, a hydrophobic transmembrane domain, and a lumenal domain; the latter consists of stem, catalytic, and lectin-like domains. Moreover, each of this ppGalNAc-T's contains important sequence motifs that are typical for this class of glycosyltransferases. These include a glycosyltransferase 1 motif containing the DXH sequence, a Gal/GalNAc-T motif, and the CLD and QXW sequence motifs located in alpha-, beta-, and gamma-repeats present within the lectin-like domain. The coding regions of T. gondii ppGalNAc-T1, -T2, and -T3 reside in multiple exons ranging in number from 6 to 10. Our results demonstrate that mucin-type O-glycosylation in T. gondii is catalysed by a multimember gene family, which is evolutionarily conserved from single-celled eukaryotes through nematodes and insects up to mammals. Taken together, this information creates the basis for future studies of the function of the ppGalNAc-T gene family in the pathobiology of this apicomplexan parasite.     

Structure and Action of the Myxobacterial Chondrochloren Halogenase CndH: A New Variant of FAD-dependent Halogenases

The crystal structure of the FAD-dependent chondrochloren halogenase CndH has been established at 2.1 Å resolution. The enzyme contains the characteristic FAD-binding scaffold of the glutathione reductase superfamily. Except for its C-terminal domain, the chainfold of CndH is virtually identical with those of FAD-dependent aromatic hydroxylases. When compared to the structurally known FAD-dependent halogenases PrnA and RebH, CndH lacks a 45 residue segment near position 100 and deviates in the C-terminal domain. Both variations are near the active center and appear to reflect substrate differences. Whereas PrnA and RebH modify free tryptophan, CndH halogenates the tyrosyl group of a chondrochloren precursor that is most likely bound to a carrier protein. In contrast to PrnA and RebH, which enclose their small substrate completely, CndH has a large non-polar surface patch that may accommodate the putative carrier. Apart from the substrate binding site, the active center of CndH corresponds to those of PrnA and RebH. At the halogenation site, CndH has the characteristic lysine (Lys76) but lacks the required base Glu346 (PrnA). This base may be supplied by a residue of its C-terminal domain or by the carrier. These differences were corroborated by an overall sequence comparison between the known FAD-dependent halogenases, which revealed a split into a PrnA-RebH group and a CndH group. The two functionally established members of the CndH group use carrier-bound substrates, whereas three members of PrnA-RebH group are known to accept a free amino acid. Given the structural and functional distinction, we classify CndH as a new variant B of the FAD-dependent halogenases, adding a new feature to the structurally established variant A enzymes PrnA and RebH.     

Discovering motif pairs at interaction sites from protein sequences on a proteome-wide scale

MOTIVATION: Protein-protein interaction, mediated by protein interaction sites, is intrinsic to many functional processes in the cell. In this paper, we propose a novel method to discover patterns in protein interaction sites. We observed from protein interaction networks that there exist a kind of significant substructures called interacting protein group pairs, which exhibit an all-versus-all interaction between the two protein-sets in such a pair. The full-interaction between the pair indicates a common interaction mechanism shared by the proteins in the pair, which can be referred as an interaction type. Motif pairs at the interaction sites of the protein group pairs can be used to represent such interaction type, with each motif derived from the sequences of a protein group by standard motif discovery algorithms. The systematic discovery of all pairs of interacting protein groups from large protein interaction networks is a computationally challenging problem. By a careful and sophisticated problem transformation, the problem is solved using efficient algorithms for mining frequent patterns, a problem extensively studied in data mining. RESULTS: We found 5349 pairs of interacting protein groups from a yeast interaction dataset. The expected value of sequence identity within the groups is only 7.48%, indicating non-homology within these protein groups. We derived 5343 motif pairs from these group pairs, represented in the form of blocks. Comparing our motifs with domains in the BLOCKS and PRINTS databases, we found that our blocks could be mapped to an average of 3.08 correlated blocks in these two databases. The mapped blocks occur 4221 out of total 6794 domains (protein groups) in these two databases. Comparing our motif pairs with iPfam consisting of 3045 interacting domain pairs derived from PDB, we found 47 matches occurring in 105 distinct PDB complexes. Comparing with another putative domain interaction database InterDom, we found 203 matches. AVAILABILITY: http://research.i2r.a-star.edu.sg/BindingMotifPairs/resources. SUPPLEMENTARY INFORMATION: http://research.i2r.a-star.edu.sg/BindingMotifPairs and Bioinformatics online.     

Gene cluster of Rhodothermus marinus high-potential iron-sulfur Protein: oxygen oxidoreductase, a caa(3)-type oxidase belonging to the superfamily of heme-copper oxidases

The respiratory chain of the thermohalophilic bacterium Rhodothermus marinus contains an oxygen reductase, which uses HiPIP (high potential iron-sulfur protein) as an electron donor. The structural genes encoding the four subunits of this HiPIP:oxygen oxidoreductase were cloned and sequenced. The genes for subunits II, I, III, and IV (named rcoxA to rcoxD) are found in this order and seemed to be organized in an operon of at least five genes with a terminator structure a few nucleotides downstream of rcoxD. Examination of the amino acid sequence of the Rcox subunits shows that the subunits of the R. marinus enzyme have homology to the corresponding subunits of oxidases belonging to the superfamily of heme-copper oxidases. RcoxB has the conserved histidines involved in binding the binuclear center and the low-spin heme. All of the residues proposed to be involved in proton transfer channels are conserved, with the exception of the key glutamate residue of the D-channel (E(278), Paracoccus denitrificans numbering). Analysis of the homology-derived structural model of subunit I shows that the phenol group of a tyrosine (Y) residue and the hydroxyl group of the following serine (S) may functionally substitute the glutamate carboxyl in proton transfer. RcoxA has an additional sequence for heme C binding, after the Cu(A) domain, that is characteristic of caa(3) oxidases belonging to the superfamily. Homology modeling of the structure of this cytochrome domain of subunit II shows no marked electrostatic character, especially around the heme edge region, suggesting that the interaction with a redox partner is not of an electrostatic nature. This observation is analyzed in relation to the electron donor for this caa(3) oxidase, the HiPIP. In conclusion, it is shown that an oxidase, which uses an iron-sulfur protein as an electron donor, is structurally related to the caa(3) class of heme-copper cytochrome c oxidases. The data are discussed in the framework of the evolution of oxidases within the superfamily of heme-copper oxidases.     

The hydroxynitrile lyase from almond: a lyase that looks like an oxidoreductase

BACKGROUND: Cyanogenesis is a defense process of several thousand plant species. Hydroxynitrile lyase, a key enzyme of this process, cleaves a cyanohydrin into hydrocyanic acid and the corresponding aldehyde or ketone. The reverse reaction constitutes an important tool in biocatalysis. Different classes of hydroxynitrile lyases have convergently evolved from FAD-dependent oxidoreductases, alpha/beta hydrolases, and alcohol dehydrogenases. The FAD-dependent hydroxynitrile lyases (FAD-HNLs) carry a flavin cofactor whose redox properties appear to be unimportant for catalysis. RESULTS: We have determined the crystal structure of a 61 kDa hydroxynitrile lyase isoenzyme from Prunus amygdalus (PaHNL1) to 1.5 A resolution. Clear electron density originating from four glycosylation sites could be observed. As concerns the overall protein fold including the FAD cofactor, PaHNL1 belongs to the family of GMC oxidoreductases. The active site for the HNL reaction is probably at a very similar position as the active sites in homologous oxidases. CONCLUSIONS: There is strong evidence from the structure and the reaction product that FAD-dependent hydroxynitrile lyases have evolved from an aryl alcohol oxidizing precursor. Since key residues implicated in oxidoreductase activity are also present in PaHNL1, it is not obvious why this enzyme shows no oxidase activity. Similarly, features proposed to be relevant for hydroxy-nitrile lyase activity in other hydroxynitrile lyases, i.e., a general base and a positive charge to stabilize the cyanide, are not obviously present in the putative active site of PaHNL1. Therefore, the reason for its HNL activity is far from being well understood at this point.     

Local motifs involved in the canonical structure of the ligand-binding domain in the nuclear receptor superfamily

Structural and sequence alignment analyses have revealed the existence of class-dependent and -independent local motifs involved in the overall fold of the ligand-binding domain (LBD) in the nuclear receptor (NR) superfamily. Of these local motifs, three local motifs, i.e., AF-2 fixed motifs, were involved in the agonist conformation of the activation function-2 (AF-2) region of the LBD. Receptor–agonist interactions increased the stability of these AF-2 fixed motifs in the agonist conformation. In contrast, perturbation of the AF-2 fixed motifs by a ligand or another protein molecule led the AF-2 architecture to adopt an antagonist conformation. Knowledge of this process should provide us with novel insights into the ‘agonism’ and ‘antagonism’ of NRs.