ARED: human AU-rich element-containing mRNA database reveals an unexpectedly diverse functional repertoire of encoded proteins

The adenylate uridylate-rich elements (AREs) mediate the rapid turnover of mRNAs encoding proteins that regulate cellular growth and body response to exogenous agents such as microbes, inflammatory and environmental stimuli. However, the full repertoire of ARE-containing mRNAs is unknown. Here, we explore the distribution of AREs in human mRNA sequences. Computational derivation of a 13-bp ARE pattern was performed using multiple expectation maximization for motif elicitations (MEME) and consensus analyses. This pattern was statistically validated for the specificity towards the 3′-untranslated region and not coding region. The computationally derived ARE pattern is the basis of a database which contains non-redundant full-length ARE-mRNAs. The ARE-mRNA database (ARED; http://rc.kfshrc.edu.sa/ared) reveals that ARE-mRNAs encode a wide repertoire of functionally diverse proteins that belong to different biological processes and are important in several disease states. Cluster analysis was performed using the ARE sequences to demonstrate potential relationships between the type and number of ARE motifs, and the functional characteristics of the proteins.

CluSTr: a database of clusters of SWISS-PROT+TrEMBL proteins

The CluSTr (Clusters of SWISS-PROT and TrEMBL proteins) database offers an automatic classification of SWISS-PROT and TrEMBL proteins into groups of related proteins. The clustering is based on analysis of all pairwise comparisons between protein sequences. Analysis has been carried out for different levels of protein similarity, yielding a hierarchical organisation of clusters. The database provides links to InterPro, which integrates information on protein families, domains and functional sites from PROSITE, PRINTS, Pfam and ProDom. Links to the InterPro graphical interface allow users to see at a glance whether proteins from the cluster share particular functional sites. CluSTr also provides cross-references to HSSP and PDB. The database is available for querying and browsing at http://www.ebi.ac.uk/clustr.

The COG database: new developments in phylogenetic classification of proteins from complete genomes

The database of Clusters of Orthologous Groups of proteins (COGs), which represents an attempt on a phylogenetic classification of the proteins encoded in complete genomes, currently consists of 2791 COGs including 45 350 proteins from 30 genomes of bacteria, archaea and the yeast Saccharomyces cerevisiae (http://www.ncbi.nlm.nih.gov/COG). In addition, a supplement to the COGs is available, in which proteins encoded in the genomes of two multicellular eukaryotes, the nematode Caenorhabditis elegans and the fruit fly Drosophila melanogaster, and shared with bacteria and/or archaea were included. The new features added to the COG database include information pages with structural and functional details on each COG and literature references, improvements of the COGNITOR program that is used to fit new proteins into the COGs, and classification of genomes and COGs constructed by using principal component analysis.

Proteome Analysis Database: online application of InterPro and CluSTr for the functional classification of proteins in whole genomes

The SWISS-PROT group at EBI has developed the Proteome Analysis Database utilising existing resources and providing comparative analysis of the predicted protein coding sequences of the complete genomes of bacteria, archaea and eukaryotes (http://www.ebi.ac.uk/proteome/). The two main projects used, InterPro and CluSTr, give a new perspective on families, domains and sites and cover 31–67% (InterPro statistics) of the proteins from each of the complete genomes. CluSTr covers the three complete eukaryotic genomes and the incomplete human genome data. The Proteome Analysis Database is accompanied by a program that has been designed to carry out InterPro proteome comparisons for any one proteome against any other one or more of the proteomes in the database.

Digging for dead genes: an analysis of the characteristics of the pseudogene population in the Caenorhabditis elegans genome

Pseudogenes are non-functioning copies of genes in genomic DNA, which may either result from reverse transcription from an mRNA transcript (processed pseudogenes) or from gene duplication and subsequent disablement (non-processed pseudogenes). As pseudogenes are apparently ‘dead’, they usually have a variety of obvious disablements (e.g., insertions, deletions, frameshifts and truncations) relative to their functioning homologs. We have derived an initial estimate of the size, distribution and characteristics of the pseudogene population in the Caenorhabditis elegans genome, performing a survey in ‘molecular archaeology’. Corresponding to the 18 576 annotated proteins in the worm (i.e., in Wormpep18), we have found an estimated total of 2168 pseudogenes, about one for every eight genes. Few of these appear to be processed. Details of our pseudogene assignments are available from http://bioinfo.mbb.yale.edu/genome/worm/pseudogene. The population of pseudogenes differs significantly from that of genes in a number of respects: (i) pseudogenes are distributed unevenly across the genome relative to genes, with a disproportionate number on chromosome IV; (ii) the density of pseudogenes is higher on the arms of the chromosomes; (iii) the amino acid composition of pseudogenes is midway between that of genes and (translations of) random intergenic DNA, with enrichment of Phe, Ile, Leu and Lys, and depletion of Asp, Ala, Glu and Gly relative to the worm proteome; and (iv) the most common protein folds and families differ somewhat between genes and pseudogenes—whereas the most common fold found in the worm proteome is the immunoglobulin fold and the most common ‘pseudofold’ is the C-type lectin. In addition, the size of a gene family bears little overall relationship to the size of its corresponding pseudogene complement, indicating a highly dynamic genome. There are in fact a number of families associated with large populations of pseudogenes. For example, one family of seven-transmembrane receptors (represented by gene B0334.7) has one pseudogene for every four genes, and another uncharacterized family (represented by gene B0403.1) is approximately two-thirds pseudogenic. Furthermore, over a hundred apparent pseudogenic fragments do not have any obvious homologs in the worm.

Smad proteins function as co-modulators for MEF2 transcriptional regulatory proteins

An emerging theme in transforming growth factor-β (TGF-β) signalling is the association of the Smad proteins with diverse groups of transcriptional regulatory proteins. Several Smad cofactors have been identified to date but the diversity of TGF-β effects on gene transcription suggests that interactions with other co-regulators must occur. In these studies we addressed the possible interaction of Smad proteins with the myocyte enhancer-binding factor 2 (MEF2) transcriptional regulators. Our studies indicate that Smad2 and 4 (Smad2/4) complexes cooperate with MEF2 regulatory proteins in a GAL4-based one-hybrid reporter gene assay. We have also observed in vivo interactions between Smad2 and MEF2A using co-immunoprecipitation assays. This interaction is confirmed by glutathione S-transferase pull-down analysis. Immunofluorescence studies in C2C12 myotubes show that Smad2 and MEF2A co-localise in the nucleus of multinuclear myotubes during differentiation. Interestingly, phospho-acceptor site mutations of MEF2 that render it unresponsive to p38 MAP kinase signalling abrogate the cooperativity with the Smads suggesting that p38 MAP Kinase-catalysed phosphorylation of MEF2 is a prerequisite for the Smad–MEF2 interaction. Thus, the association between Smad2 and MEF2A may subserve a physical link between TGF-β signalling and a diverse array of genes controlled by the MEF2 cis element.

Combined transcriptome and proteome analysis as a powerful approach to study genes under glucose repression in Bacillus subtilis

We used 2D protein gel electrophoresis and DNA microarray technologies to systematically analyze genes under glucose repression in Bacillus subtilis. In particular, we focused on genes expressed after the shift from glycolytic to gluconeogenic at the middle logarithmic phase of growth in a nutrient sporulation medium, which remained repressed by the addition of glucose. We also examined whether or not glucose repression of these genes was mediated by CcpA, the catabolite control protein of this bacterium. The wild-type and ccpA1 cells were grown with and without glucose, and their proteomes and transcriptomes were compared. 2D gel electrophoresis allowed us to identify 11 proteins, the synthesis of which was under glucose repression. Of these proteins, the synthesis of four (IolA, I, S and PckA) was under CcpA-independent control. Microarray analysis enabled us to detect 66 glucose-repressive genes, 22 of which (glmS, acoA, C, yisS, speD, gapB, pckA, yvdR, yxeF, iolA, B, C, D, E, F, G, H, I, J, R, S and yxbF ) were at least partially under CcpA-independent control. Furthermore, we found that CcpA and IolR, a repressor of the iol divergon, were involved in the glucose repression of the synthesis of inositol dehydrogenase encoded by iolG included in the above list. The CcpA-independent glucose repression of the iol genes appeared to be explained by inducer exclusion.

Identification of protein oligomerization states by analysis of interface conservation

The discrimination of true oligomeric protein–protein contacts from nonspecific crystal contacts remains problematic. Criteria that have been used previously base the assignment of oligomeric state on consideration of the area of the interface and/or the results of scoring functions based on statistical potentials. Both techniques have a high success rate but fail in more than 10% of cases. More importantly, the oligomeric states of several proteins are incorrectly assigned by both methods. Here we test the hypothesis that true oligomeric contacts should be identifiable on the basis of an increased degree of conservation of the residues involved in the interface. By quantifying the degree of conservation of the interface and comparing it with that of the remainder of the protein surface, we develop a new criterion that provides a highly effective complement to existing methods.

A Family of Small Coiled-Coil–forming Proteins Functioning at the Late Endosome in Yeast

The multispanning membrane protein Ste6, a member of the ABC-transporter family, is transported to the yeast vacuole for degradation. To identify functions involved in the intracellular trafficking of polytopic membrane proteins, we looked for functions that block Ste6 transport to the vacuole upon overproduction. In our screen, we identified several known vacuolar protein sorting (VPS) genes (SNF7/VPS32, VPS4, and VPS35) and a previously uncharacterized open reading frame, which we named MOS10 (more of Ste6). Sequence analysis showed that Mos10 is a member of a small family of coiled-coil–forming proteins, which includes Snf7 and Vps20. Deletion mutants of all three genes stabilize Ste6 and show a “class E vps phenotype.” Maturation of the vacuolar hydrolase carboxypeptidase Y was affected in the mutants and the endocytic tracer FM4-64 and Ste6 accumulated in a dot or ring-like structure next to the vacuole. Differential centrifugation experiments demonstrated that about half of the hydrophilic proteins Mos10 and Vps20 was membrane associated. The intracellular distribution was further analyzed for Mos10. On sucrose gradients, membrane-associated Mos10 cofractionated with the endosomal t-SNARE Pep12, pointing to an endosomal localization of Mos10. The growth phenotypes of the mutants suggest that the “Snf7-family” members are involved in a cargo-specific event.

Identification, retinoid binding, and x-ray analysis of a human retinol-binding protein

Two cellular retinol-binding proteins (CRBP I and II) with distinct tissue distributions and retinoid-binding properties have been recognized thus far in mammals. Here, we report the identification of a human retinol-binding protein resembling type I (55.6% identity) and type II (49.6% identity) CRBPs, but with a unique H residue in the retinoid-binding site and a distinctively different tissue distribution. Additionally, this binding protein (CRBP III) exhibits a remarkable sequence identity (62.2%) with the recently identified ι-crystallin/CRBP of the diurnal gecko Lygodactylus picturatus [Werten, P. J. L., Röll, B., van Alten, D. M. F. & de Jong, W. W. (2000) Proc. Natl. Acad. Sci. USA 97, 3282–3287 (First Published March 21, 2000; 10.1073/pnas.050500597)]. CRBP III and all-trans-retinol form a complex (Kd ≈ 60 nM), the absorption spectrum of which is characterized by the peculiar fine structure typical of the spectra of holo-CRBP I and II. As revealed by a 2.3-Å x-ray molecular model of apo-CRBP III, the amino acid residues that line the retinol-binding site in CRBP I and II are positioned nearly identically in the structure of CRBP III. At variance with the human CRBP I and II mRNAs, which are most abundant in ovary and intestine, respectively, the CRBP III mRNA is expressed at the highest levels in kidney and liver thus suggesting a prominent role for human CRBP III as an intracellular mediator of retinol metabolism in these tissues.

Phylogenetic analysis of tmRNA genes within a bacterial subgroup reveals a specific structural signature

Bacterial tmRNA mediates a trans-translation reaction, which permits the recycling of stalled ribosomes and probably also contributes to the regulated expression of a subset of genes. Its action results in the addition of a small number of C-terminal amino acids to protein whose synthesis had stalled and these constitute a proteolytic recognition tag for the degradation of these incompletely synthesized proteins. Previous work has identified pseudoknots and stem–loops that are widely conserved in divergent bacteria. In the present work an alignment of tmRNA gene sequences within 13 β-proteobacteria reveals an additional sub-structure specific for this bacterial group. This sub-structure is in pseudoknot Pk2, and consists of one to two additional stem–loop(s) capped by stable GNRA tetraloop(s). Three-dimensional models of tmRNA pseudoknot 2 (Pk2) containing various topological versions of the additional sub-structure suggest that the sub-structures likely point away from the core of the RNA, containing both the tRNA and the mRNA domains. A putative tertiary interaction has also been identified.

A comprehensive two-hybrid analysis to explore the yeast protein interactome

Protein–protein interactions play crucial roles in the execution of various biological functions. Accordingly, their comprehensive description would contribute considerably to the functional interpretation of fully sequenced genomes, which are flooded with novel genes of unpredictable functions. We previously developed a system to examine two-hybrid interactions in all possible combinations between the ≈6,000 proteins of the budding yeast Saccharomyces cerevisiae. Here we have completed the comprehensive analysis using this system to identify 4,549 two-hybrid interactions among 3,278 proteins. Unexpectedly, these data do not largely overlap with those obtained by the other project [Uetz, P., et al. (2000) Nature (London) 403, 623–627] and hence have substantially expanded our knowledge on the protein interaction space or interactome of the yeast. Cumulative connection of these binary interactions generates a single huge network linking the vast majority of the proteins. Bioinformatics-aided selection of biologically relevant interactions highlights various intriguing subnetworks. They include, for instance, the one that had successfully foreseen the involvement of a novel protein in spindle pole body function as well as the one that may uncover a hitherto unidentified multiprotein complex potentially participating in the process of vesicular transport. Our data would thus significantly expand and improve the protein interaction map for the exploration of genome functions that eventually leads to thorough understanding of the cell as a molecular system.

Proteomic analysis of the bacterial cell cycle

A global approach was used to analyze protein synthesis and stability during the cell cycle of the bacterium Caulobacter crescentus. Approximately one-fourth (979) of the estimated C. crescentus gene products were detected by two-dimensional gel electrophoresis, 144 of which showed differential cell cycle expression patterns. Eighty-one of these proteins were identified by mass spectrometry and were assigned to a wide variety of functional groups. Pattern analysis revealed that coexpression groups were functionally clustered. A total of 48 proteins were rapidly degraded in the course of one cell cycle. More than half of these unstable proteins were also found to be synthesized in a cell cycle-dependent manner, establishing a strong correlation between rapid protein turnover and the periodicity of the bacterial cell cycle. This is, to our knowledge, the first evidence for a global role of proteolysis in bacterial cell cycle control.

Fusicoccin, 14-3-3 Proteins, and Defense Responses in Tomato Plants1

Fusicoccin (FC) is a fungal toxin that activates the plant plasma membrane H+-ATPase by binding with 14-3-3 proteins, causing membrane hyperpolarization. Here we report on the effect of FC on a gene-for-gene pathogen-resistance response and show that FC application induces the expression of several genes involved in plant responses to pathogens. Ten members of the FC-binding 14-3-3 protein gene family were isolated from tomato (Lycopersicon esculentum) to characterize their role in defense responses. Sequence analysis is suggestive of common biochemical functions for these tomato 14-3-3 proteins, but their genes showed different expression patterns in leaves after challenges. Different specific subsets of 14-3-3 genes were induced after treatment with FC and during a gene-for-gene resistance response. Possible roles for the H+-ATPase and 14-3-3 proteins in responses to pathogens are discussed.

Characterization of Two Novel Type I Ribosome-Inactivating Proteins from the Storage Roots of the Andean Crop Mirabilis expansa1

Two novel type I ribosome-inactivating proteins (RIPs) were found in the storage roots of Mirabilis expansa, an underutilized Andean root crop. The two RIPs, named ME1 and ME2, were purified to homogeneity by ammonium sulfate precipitation, cation-exchange perfusion chromatography, and C4 reverse-phase chromatography. The two proteins were found to be similar in size (27 and 27.5 kD) by sodium dodecyl sulfate-polyacrylamide gel electrophoresis, and their isoelectric points were determined to be greater than pH 10.0. Amino acid N-terminal sequencing revealed that both ME1 and ME2 had conserved residues characteristic of RIPs. Amino acid composition and western-blot analysis further suggested a structural similarity between ME1 and ME2. ME2 showed high similarity to the Mirabilis jalapa antiviral protein, a type I RIP. Depurination of yeast 26S rRNA by ME1 and ME2 demonstrated their ribosome-inactivating activity. Because these two proteins were isolated from roots, their antimicrobial activity was tested against root-rot microorganisms, among others. ME1 and ME2 were active against several fungi, including Pythium irregulare, Fusarium oxysporum solani, Alternaria solani, Trichoderma reesei, and Trichoderma harzianum, and an additive antifungal effect of ME1 and ME2 was observed. Antibacterial activity of both ME1 and ME2 was observed against Pseudomonas syringae, Agrobacterium tumefaciens, Agrobacterium radiobacter, and others.