Comprehensive sequence analysis of the 182 predicted open reading frames of yeast chromosome III

20 years of the SMART protein domain annotation resource

SMART (Simple Modular Architecture Research Tool) is a web resource (http://smart.embl.de) for the identification and annotation of protein domains and the analysis of protein domain architectures. SMART version 8 contains manually curated models for more than 1300 protein domains, with approximately 100 new models added since our last update article (1). The underlying protein databases were synchronized with UniProt (2), Ensembl (3) and STRING (4), doubling the total number of annotated domains and other protein features to more than 200 million. In its 20th year, the SMART analysis results pages have been streamlined again and its information sources have been updated. SMART’s vector based display engine has been extended to all protein schematics in SMART and rewritten to use the latest web technologies. The internal full text search engine has been redesigned and updated, resulting in greatly increased search speed.

Requirement of POL3 and POL4 on non-homologous and microhomology-mediated end joining in rad50/xrs2 mutants of Saccharomyces cerevisiae

Non-homologous end joining (NHEJ) directly joins two broken DNA ends without sequence homology. A distinct pathway called microhomology-mediated end joining (MMEJ) relies on a few base pairs of homology between the recombined DNA. The majority of DNA double-strand breaks caused by endogenous oxygen species or ionizing radiation contain damaged bases that hinder direct religation. End processing is required to remove mismatched nucleotides and fill in gaps during end joining of incompatible ends. POL3 in Saccharomyces cerevisiae encodes polymerase δ that is required for DNA replication and other DNA repair processes. Our previous results have shown that POL3 is involved in gap filling at 3' overhangs in POL4-independent NHEJ. Here, we studied the epistatic interaction between POL3, RAD50, XRS2 and POL4 in NHEJ using a plasmid-based endjoining assay in yeast. We demonstrated that either rad50 or xrs2 mutation is epistatic for end joining of compatible ends in the rad50 pol3-t or xrs2 pol3-t double mutants. However, the pol3-t and rad50 or pol3-t and xrs2 mutants caused an additive decrease in the end-joining efficiency of incompatible ends, suggesting that POL3 and RAD50 or POL3 and XRS2 exhibit independent functions in NHEJ. In the rad50 pol4 mutant, end joining of incompatible ends was not detected. In the rad50 or xrs2 mutants, NHEJ events did not contain any microhomology at the rejoined junctions. The pol3-t mutation restored MMEJ in the rad50 or xrs2 mutant backgrounds. Moreover, we demonstrated that NHEJ of incompatible ends required RAD50 and POL4 more than POL3. In conclusion, POL3 and POL4 have differential functions in NHEJ, independent of the RAD50-mediated repair pathway.

Characterization of the neuropsychological phenotype of glycine N-methyltransferase-/- mice and evaluation of its responses to clozapine and sarcosine treatments

Glycine N-methyltransferase (GNMT) affects cellular methylation capacity through regulating the ratio between S-adenosylmethionine (SAM) and S-adenosylhomocysteine (SAH). The product of its enzymatic reaction-sarcosine has antipsychotic effect in patients with schizophrenia. In this study, through RT-PCR and immunohistochemical staining, we demonstrated that GNMT expressed in various neurons located in the cerebral cortex, hippocampus, substantia nigra and cerebellum. Compared to the wild-type mice, Gnmt-/- mice had significantly lower level of sarcosine in the cerebral cortex. Real-time PCR identified genes involved in the methionine metabolism (Dnmt1 and Dnmt3a), ErbB (Nrg1 and ErbB4) and mTOR (Akt2, S6, S6k1 and S6k2) signaling pathways were dysregulated significantly in the cortex of Gnmt-/- mice. Acoustic startle reflex test demonstrated that Gnmt-/- mice had significantly lower level of prepulse inhibition and the deficit was ameliorated through clozapine or sarcosine treatment. Furthermore, liver-specific-human-GNMT transgenic with Gnmt-/- (Tg-GNMT/Gnmt-/-) mice were used to rule out that the phenotype was due to abnormal liver function. In summary, the neuropsychological abnormalities found in Gnmt-/- mice may represent an endophenotype of schizophrenia. GNMT plays an important role in maintaining normal physiological function of brain and Tg-GNMT/Gnmt-/- mice are useful models for development of therapeutics for patients with schizophrenia.

Solving the Problem: Genome Annotation Standards before the Data Deluge

The promise of genome sequencing was that the vast undiscovered country would be mapped out by comparison of the multitude of sequences available and would aid researchers in deciphering the role of each gene in every organism. Researchers recognize that there is a need for high quality data. However, different annotation procedures, numerous databases, and a diminishing percentage of experimentally determined gene functions have resulted in a spectrum of annotation quality. NCBI in collaboration with sequencing centers, archival databases, and researchers, has developed the first international annotation standards, a fundamental step in ensuring that high quality complete prokaryotic genomes are available as gold standard references. Highlights include the development of annotation assessment tools, community acceptance of protein naming standards, comparison of annotation resources to provide consistent annotation, and improved tracking of the evidence used to generate a particular annotation. The development of a set of minimal standards, including the requirement for annotated complete prokaryotic genomes to contain a full set of ribosomal RNAs, transfer RNAs, and proteins encoding core conserved functions, is an historic milestone. The use of these standards in existing genomes and future submissions will increase the quality of databases, enabling researchers to make accurate biological discoveries.

Detecting conserved secondary structures in RNA molecules using constrained structural alignment

Constrained sequence alignment has been studied extensively in the past. Different forms of constraints have been investigated, where a constraint can be a subsequence, a regular expression, or a probability matrix of symbols and positions. However, constrained structural alignment has been investigated to a much lesser extent. In this paper, we present an efficient method for constrained structural alignment and apply the method to detecting conserved secondary structures, or structural motifs, in a set of RNA molecules. The proposed method combines both sequence and structural information of RNAs to find an optimal local alignment between two RNA secondary structures, one of which is a query and the other is a subject structure in the given set. The method allows a biologist to annotate conserved regions, or constraints, in the query RNA structure and incorporates these regions into the alignment process to obtain biologically more meaningful alignment scores. A statistical measure is developed to assess the significance of the scores. Experimental results based on detecting internal ribosome entry sites in the RNA molecules of hepatitis C virus and Trypanosoma brucei demonstrate the effectiveness of the proposed method and its superiority over existing techniques.

Hymenoic acid, a novel specific inhibitor of human DNA polymerase lambda from a fungus of Hymenochaetaceae sp

Hymenoic acid (1) is a natural compound isolated from cultures of a fungus, Hymenochaetaceae sp., and this structure was determined by spectroscopic analyses. Compound 1 is a novel sesquiterpene, trans-4-[(1'E,5'S)-5'-carboxy-1'-methyl-1'-hexenyl]cyclohexanecarboxylic acid. This compound selectively inhibited the activity of human DNA polymerase lambda (pol lambda) in vitro, and 50% inhibition was observed at a concentration of 91.7microM. Compound 1 did not influence the activities of the other seven mammalian pols (i.e., pols alpha, gamma, delta, epsilon, eta, iota, and kappa), but also showed no effect even on the activity of pol beta, which is thought to have a very similar three-dimensional structure to the pol beta-like region of pol lambda. This compound also did not inhibit the activities of prokaryotic pols and other DNA metabolic enzymes tested. These results suggested that compound 1 could be a selective inhibitor of eukaryotic pol lambda. This compound had no inhibitory activities against two N-terminal truncated pol lambda, del-1 pol lambda (lacking nuclear localization signal (NLS), BRCA1 C-terminus (BRCT) domain [residues 133-575]), and del-2 pol lambda (lacking NLS, BRCT, domain and proline-rich region [residues 245-575]). The compound 1-induced inhibition of intact pol lambda activity was non-competitive with respect to both the DNA template-primer and the dNTP substrate. On the basis of these results, the pol lambda inhibitory mechanism of compound 1 is discussed.

Protein function space: viewing the limits or limited by our view?

Given that the number of protein functions on earth is finite, the rapid expansion of biological knowledge and the concomitant exponential increase in the number of protein sequences should, at some point, enable the estimation of the limits of protein function space. The functional coverage of protein sequences can be investigated using computational methods, especially given the massive amount of data being generated by large-scale environmental sequencing (metagenomics). In completely sequenced genomes, the fraction of proteins to which at least some functional features can be assigned has recently risen to as much as approximately 85%. Although this fraction is more uncertain in metagenomics surveys, because of environmental complexities and differences in analysis protocols, our global knowledge of protein functions still appears to be considerable. However, when we consider protein families, continued sequencing seems to yield an ever-increasing number of novel families. Until we reconcile these two views, the limits of protein space will remain obscured.

The gene identification problem: an overview for developers

The gene identification problem is the problem of interpreting nucleotide sequences by computer, in order to provide tentative annotation on the location, structure, and functional class of protein-coding genes. This problem is of self-evident importance, and is far from being fully solved, particularly for higher eukaryotes. Thus it is not surprising that the number of algorithm and software developers working in the area is rapidly increasing. The present paper is an overview of the field, with an emphasis on eukaryotes, for such developers.

Beta-sitosterol-3-O-beta-D-glucopyranoside: a eukaryotic DNA polymerase lambda inhibitor

Beta-sitosterol-3-O-beta-D-glucopyranoside (compound 1), a steroidal glycoside isolated from onion (Allium cepa L.) selectively inhibited the activity of mammalian DNA polymerase lambda (pol lambda) in vitro. The compound did not influence the activities of replicative DNA polymerases such as alpha, delta and epsilon, but also showed no effect even on the activity of pol beta which is thought to have a very similar three-dimensional structure to the pol beta-like region of pol lambda. Since parts of compound 1 such as beta-sitosterol (compound 2) and D-glucose (compound 3) did not influence the activities of any enzymes tested, the converted structure of compounds 2 and 3 might be important for pol lambda inhibition. The inhibitory effect of compound 1 on both intact pol lambda (i.e. residues 1-575) and a truncated pol lambda lacking the N-terminal BRCA1 C-terminus (BRCT) domain (133-575, del-1 pol lambda) was dose-dependent, and 50% inhibition was observed at a concentration of 9.1 and 5.4 microM, respectively. The compound 1-induced inhibition of del-1 pol lambda activity was non-competitive with respect to both the DNA template-primer and the dNTP substrate. On the basis of these results, the pol lambda inhibitory mechanism of compound 1 is discussed.

DNA polymerase 4 of Saccharomyces cerevisiae is important for accurate repair of methyl-methanesulfonate-induced DNA damage

The DNA polymerase 4 protein (Pol4) of Saccharomyces cerevisiae is a member of the X family of DNA polymerases whose closest human relative appears to be DNA polymerase lambda. Results from previous genetic studies conflict over the role of Pol4 in vivo. Here we show that deletion of Pol4 in a diploid strain of the SK1 genetic background results in sensitivity to methyl methanesulfonate (MMS). However, deletion of Pol4 in other strain backgrounds and in haploid strains does not yield an observable phenotype. The MMS sensitivity of a Pol4-deficient strain can be rescued by deletion of YKu70. We also show that deletion of Pol4 results in a 6- to 14-fold increase in the MMS-induced mutation frequency and in a significant increase in AT-to-TA transversions. Our studies suggest that Pol4 is critical for accurate repair of DNA lesions induced by MMS.

Characterization of reduced expression of glycine N-methyltransferase in cancerous hepatic tissues using two newly developed monoclonal antibodies

Glycine N-methyltransferase (GNMT) is a protein with multiple functions. Recently, two Italian siblings who had hepatomegaly and chronic elevation of serum transaminases were diagnosed to have GNMT deficiency caused by inherited compound heterozygosity of the GNMT gene with missence mutations. To evaluate the expression of GNMT in cell lines and tissues from hepatocellular carcinoma (HCC) patients, we produced two monoclonal antibodies (mAbs) 4-17 and 14-1 using two recombinant GNMT fusion proteins. M13 phage peptide display showed that the reactive epitopes of mAbs 4-17 and 14-1 were amino acid residues 11-15 and 272-276 of human GNMT, respectively. The dissociation constants of the binding between GNMT and mAbs were 1.7 x 10(-8) M for mAb 4-17 and 1.8 x 10(-9) M for mAb 14-1. Both mAbs can identify GNMT present in normal human and mouse liver tissues using Western blotting (WB) and immunohistochemical staining assay (IHC). In addition, WB with both mAbs showed that none of 2 hepatoblastoma and 5 HCC cell lines expressed GNMT. IHC demonstrated that 50% (13/26) of nontumorous liver tissues and 96% (24/25) of HCC tissues did not express GNMT. Therefore, the expression of GNMT was downregulated in human HCC.

Global perspectives on proteins: comparing genomes in terms of folds, pathways and beyond

The sequencing of complete genomes provides us with a global view of all the proteins in an organism. Proteomic analysis can be done on a purely sequence-based level, with a focus on finding homologues and grouping them into families and clusters of orthologs. However, incorporating protein structure into this analysis provides valuable simplification; it allows one to collect together very distantly related sequences, thus condensing the proteome into a minimal number of 'parts.' We describe issues related to surveying proteomes in terms of structural parts, including methods for fold assignment and formats for comparisons (eg top-10 lists and whole-genome trees), and show how biases in the databases and in sampling can affect these surveys. We illustrate our main points through a case study on the unique protein properties evident in many thermophile genomes (eg more salt bridges). Finally, we discuss metabolic pathways as an even greater simplification of genomes. In comparison to folds these allow the organization of many more genes into coherent systems, yet can nevertheless be understood in many of the same terms.

Over-expression of human DNA polymerase lambda in E. coli and characterization of the recombinant enzyme

BACKGROUND

DNA polymerase lambda (Pol lambda) was recently identified as a new member of the family X of DNA polymerases in eukaryotic cells. Pol lambda contains a nuclear localization signal (NLS), a BRCA1-C terminal (BRCT) domain, a proline-rich region, helix-hairpin-helix (HhH) and pol X motifs. Since the amino acid sequence for Pol lambda shares a high degree of homology to Pol beta, Pol lambda is considered to have a similar enzymatic nature to Pol beta.

RESULTS

Recombinant human Pol lambda was shown to possess template-directed DNA polymerase activity in its monomeric form. Pol lambda required either Mn2+ or Mg2+ as a metal co-factor to catalyse this activity, and optimal activity was detected at pH 8.5-9.0. Pol lambda was insensitive to aphidicolin, but was sensitive to dideoxynucleoside triphosphates or N-ethylmaleimide. By constructing the truncated Pol lambda, the proline rich region was shown to act in a suppression of its polymerization activity. A chimeric enzyme comprised of the Pol lambda N-terminal region and Pol beta also showed a reduced Pol beta activity. Proliferating cell nuclear antigen (PCNA) directly interacts with Pol lambda through its Pol beta like region in vitro.

CONCLUSIONS

Pol lambda possesses similar enzymatic nature to Pol beta; requirements of cations and optimal conditions for pH and NaCl concentration, aside from sensitivity to N-ethylmaleimide and template preference. The proline rich region of Pol lambda functions as a suppressor domain for its polymerization activity (SDPA). Pol lambda interacts directly with PCNA through its Pol beta like region. The functional consequence of this interaction is the negative regulation of Pol lambda activity.

The past, present and future of genome-wide re-annotation

Annotation, the process by which structural or functional information is inferred for genes or proteins, is crucial for obtaining value from genome sequences. We define the process of annotating a previously annotated genome sequence as 're-annotation', and examine the strengths and weaknesses of current manual and automatic genome-wide re-annotation approaches.

Characterization of 16 novel human genes showing high similarity to yeast sequences

The entire set of open reading frames (ORFs) of Saccharomyces cerevisiae has been used to perform systematic similarity searches against nucleic acid and protein databases: with the aim of identifying interesting homologies between yeast and mammalian genes. Many similarities were detected: mostly with known genes. However: several yeast ORFs were only found to match human partial sequence tags: indicating the presence of human transcripts still uncharacterized that have a homologous counterpart in yeast. About 30 such transcripts were further studied and named HUSSY (human sequence similar to yeast). The 16 most interesting are presented in this paper along with their sequencing and mapping data. As expected: most of these genes seem to be involved in basic metabolic and cellular functions (lipoic acid biosynthesis: ribulose-5-phosphate-3-epimerase: glycosyl transferase: beta-transducin: serine-threonine-kinase: ABC proteins: cation transporters). Genes related to RNA maturation were also found (homologues to DIM1: ROK1-RNA-elicase and NFS1). Furthermore: five novel human genes were detected (HUSSY-03: HUSSY-22: HUSSY-23: HUSSY-27: HUSSY-29) that appear to be homologous to yeast genes whose function is still undetermined. More information on this work can be obtained at the website http://grup.bio.unipd.it/hussy

Extreme halophiles synthesize betaine from glycine by methylation

Glycine betaine is a compatible solute, which is able to restore and maintain osmotic balance of living cells. It is synthesized and accumulated in response to abiotic stress. Betaine acts also as a methyl group donor and has a number of important applications including its use as a feed additive. The known biosynthetic pathways of betaine are universal and very well characterized. A number of enzymes catalyzing the two-step oxidation of choline to betaine have been isolated. In this work we have studied a novel betaine biosynthetic pathway in two phylogenically distant extreme halophiles, Actinopolyspora halophila and Ectothiorhodospira halochloris. We have identified a three-step series of methylation reactions from glycine to betaine, which is catalyzed by two methyltransferases, glycine sarcosine methyltransferase and sarcosine dimethylglycine methyltransferase, with partially overlapping substrate specificity. The methyltransferases from the two organisms show high sequence homology. E. halochloris methyltransferase genes were successfully expressed in Escherichia coli, and betaine accumulation and improved salt tolerance were demonstrated.

Genomic structure, expression, and chromosomal localization of the human glycine N-methyltransferase gene

The glycine N-methyltransferase (GNMT) gene encodes a protein that not only acts as an enzyme to regulate the ratio of S-adenosylmethionine to S-adenosylhomocysteine, but also participates in the detoxification pathway in liver cells. Previously, we reported that the expression level of GNMT was diminished in human hepatocellular carcinoma. In this study, the human GNMT gene was cloned and characterized. It contains six exons and spans about 10 kb. Instead of a TATA box, it has a transcriptional initiator located 801 bp upstream from the translation start codon. The gene was localized to chromosome 6p12 using fluorescence in situ hybridization. Northern blot analysis of 16 tissues from different human organs showed that GNMT was expressed only in liver, pancreas, and prostate.

Spb1p is a yeast nucleolar protein associated with Nop1p and Nop58p that is able to bind S-adenosyl-L-methionine in vitro

We present here the characterization of SPB1, an essential yeast gene that is required for ribosome synthesis. A cold-sensitive allele for that gene (referred to here as spb1-1) had been previously isolated as a suppressor of a mutation affecting the poly(A)-binding protein gene (PAB1) and a thermosensitive allele (referred to here as spb1-2) was isolated in a search for essential genes required for gene silencing in Saccharomyces cerevisiae. The two mutants are able to suppress the deletion of PAB1, and they both present a strong reduction in their 60S ribosomal subunit content. In an spb1-2 strain grown at the restrictive temperature, processing of the 27S pre-rRNA into mature 25S rRNA and 5.8S is completely abolished and production of mature 18S is reduced, while the abnormal 23S species is accumulated. Spb1p is a 96.5-kDa protein that is localized to the nucleolus. Coimmunoprecipitation experiments show that Spb1p is associated in vivo with the nucleolar proteins Nop1p and Nop5/58p. Protein sequence analysis reveals that Spb1p possesses a putative S-adenosyl-L-methionine (AdoMet)-binding domain, which is common to the AdoMet-dependent methyltransferases. We show here that Spb1p is able to bind [(3)H]AdoMet in vitro, suggesting that it is a novel methylase, whose possible substrates will be discussed.

Genome-based structural biology

Spectacular achievements in whole genome sequencing open up new possibilities for structural research. Protein structures can now be studied in their natural genomic context. On the other hand, structure prediction algorithms can be improved using species-specific tendencies in folding patterns. Finally, efficient strategies to select targets for structure determination can be devised. In this review we consider new computational approaches and results in protein structure analysis stemming from the availability of complete genomes.

Patterns of protein-fold usage in eight microbial genomes: a comprehensive structural census

Eight microbial genomes are compared in terms of protein structure. Specifically, yeast, H. influenzae, M. genitalium, M. jannaschii, Synechocystis, M. pneumoniae, H. pylori, and E. coli are compared in terms of patterns of fold usage-whether a given fold occurs in a particular organism. Of the approximately 340 soluble protein folds currently in the structure databank (PDB), 240 occur in at least one of the eight genomes, and 30 are shared amongst all eight. The shared folds are depleted in allhelical structure and enriched in mixed helix-sheet structure compared to the folds in the PDB. The top-10 most common of the shared 30 are enriched in superfolds, uniting many non-homologous sequence families, and are especially similar in overall architecture-eight having helices packed onto a central sheet. They are also very different from the common folds in the PBD, highlighting databank biases. Folds can be ranked in terms of expression as well as genome duplication. In yeast the top-10 most highly expressed folds are considerably different from the most highly duplicated folds. A tree can be constructed grouping genomes in terms of their shared folds. This has a remarkably similar topology to more conventional classifications, based on very different measures of relatedness. Finally, folds of membrane proteins can be analyzed through transmembrane-helix (TM) prediction. All the genomes appear to have similar usage patterns for these folds, with the occurrence of a particular fold falling off rapidly with increasing numbers of TM-elements, according to a "Zipf-like" law. This implies there are no marked preferences for proteins with particular numbers of TM-helices (e.g. 7-TM) in microbial genomes.

Comparing genomes in terms of protein structure: surveys of a finite parts list

We give an overview of the emerging field of structural genomics, describing how genomes can be compared in terms of protein structure. As the number of genes in a genome and the total number of protein folds are both quite limited, these comparisons take the form of surveys of a finite parts list, similar in respects to demographic censuses. Fold surveys have many similarities with other whole-genome characterizations, e.g., analyses of motifs or pathways. However, structure has a number of aspects that make it particularly suitable for comparing genomes, namely the way it allows for the precise definition of a basic protein module and the fact that it has a better defined relationship to sequence similarity than does protein function. An essential requirement for a structure survey is a library of folds, which groups the known structures into 'fold families.' This library can be built up automatically using a structure comparison program, and we described how important objective statistical measures are for assessing similarities within the library and between the library and genome sequences. After building the library, one can use it to count the number of folds in genomes, expressing the results in the form of Venn diagrams and 'top-10' statistics for shared and common folds. Depending on the counting methodology employed, these statistics can reflect different aspects of the genome, such as the amount of internal duplication or gene expression. Previous analyses have shown that the common folds shared between very different microorganisms, i.e., in different kingdoms, have a remarkably similar structure, being comprised of repeated strand-helix-strand super-secondary structure units. A major difficulty with this sort of 'fold-counting' is that only a small subset of the structures in a complete genome are currently known and this subset is prone to sampling bias. One way of overcoming biases is through structure prediction, which can be applied uniformly and comprehensively to a whole genome. Various investigators have, in fact, already applied many of the existing techniques for predicting secondary structure and transmembrane (TM) helices to the recently sequenced genomes. The results have been consistent: microbial genomes have similar fractions of strands and helices even though they have significantly different amino acid composition. The fraction of membrane proteins with a given number of TM helices falls off rapidly with more TM elements, approximately according to a Zipf law. This latter finding indicates that there is no preference for the highly studied 7-TM proteins in microbial genomes. Continuously updated tables and further information pertinent to this review are available over the web at http://bioinfo.mbb.yale.edu/genome.

Fowlpox virus encodes nonessential homologs of cellular alpha-SNAP, PC-1, and an orphan human homolog of a secreted nematode protein

The genome of fowlpox virus (FWPV), type species of the Avipoxviridae, is considerably rearranged compared with that of vaccinia virus (the prototypic poxvirus and type species of the Orthopoxviridae) and is 30% larger. It is likely that the genome of FWPV contains genes in addition to those found in vaccinia virus, probably involved with its replication and survival in the chicken. A 7,470-bp segment of the FWPV genome has five open reading frames (ORFs), two of which encode ankyrin repeat proteins, many examples of which have been found in poxviruses. The remaining ORFs encode homologs of cellular genes not reported in any other virus. ORF-2 encodes a homolog of the yeast Sec17p and mammalian SNAP proteins, crucial to vesicular transport in the exocytic pathway. ORF-3 encodes a homolog of an orphan human protein, R31240_2, encoded on 19p13.2. ORF-3 is also homologous to three proteins (YLS2, YMV6, and C07B5.5) from the free-living nematode Caenorhabditis elegans and to a 43-kDa antigen from the parasitic nematode Trichinella spiralis. ORF-5 encodes a homolog of the mammalian plasma cell antigen PC-1, a type II glycoprotein with exophosphodiesterase activity. The ORFs are present in the virulent precursor, HP1, of the sequenced attenuated virus (FP9) and are conserved in other strains of FWPV. They were shown, by deletion mutagenesis, to be nonessential to virus replication in tissue culture. RNA encoding the viral homolog of PC-1 is expressed strongly early and late in infection, but RNAs encoding the homologs of SNAP and R31240_2 are expressed weakly and late.

Biosynthesis of 2-aceto-2-hydroxy acids: acetolactate synthases and acetohydroxyacid synthases

Two groups of enzymes are classified as acetolactate synthase (EC 4. 1.3.18). This review deals chiefly with the FAD-dependent, biosynthetic enzymes which readily catalyze the formation of acetohydroxybutyrate from pyruvate and 2-oxobutyrate, as well as of acetolactate from two molecules of pyruvate (the ALS/AHAS group). These enzymes are generally susceptible to inhibition by one or more of the branched-chain amino acids which are ultimate products of the acetohydroxyacids, as well as by several classes of herbicides (sulfonylureas, imidazolinones and others). Some ALS/AHASs also catalyze the (non-physiological) oxidative decarboxylation of pyruvate, leading to peracetic acid; the possible relationship of this process to oxygen toxicity is considered. The bacterial ALS/AHAS which have been well characterized consist of catalytic subunits (around 60 kDa) and smaller regulatory subunits in an alpha2beta2 structure. In the case of Escherichia coli isozyme III, assembly and dissociation of the holoenzyme has been studied. The quaternary structure of the eukaryotic enzymes is less clear and in plants and yeast only catalytic polypeptides (homologous to those of bacteria) have been clearly identified. The presence of regulatory polypeptides in these organisms cannot be ruled out, however, and genes which encode putative ALS/AHAS regulatory subunits have been identified in some cases. A consensus sequence can be constructed from the 21 sequences which have been shown experimentally to represent ALS/AHAS catalytic polypeptides. Many other sequences fit this consensus, but some genes identified as putative 'acetolactate synthase genes' are almost certainly not ALS/AHAS. The solution of the crystal structures of several thiamin diphosphate (ThDP)-dependent enzymes which are homologous to ALS/AHAS, together with the availability of many amino acid sequences for the latter enzymes, has made it possible for two laboratories to propose similar, reasonable models for a dimer of catalytic subunits of an ALS/AHAS. A number of characteristics of these enzymes can now be better understood on the basis of such models: the nature of the herbicide binding site, the structural role of FAD and the binding of ThDP-Mg2+. The models are also guides for experimental testing of ideas concerning structure-function relationships in these enzymes, e.g. the nature of the substrate recognition site. Among the important remaining questions is how the enzyme suppresses alternative reactions of the intrinsically reactive hydroxyethylThDP enamine formed by the decarboxylation of the first substrate molecule and specifically promotes its condensation with 2-oxobutyrate or pyruvate.

The Pol beta-14 dominant negative rat DNA polymerase beta mutator mutant commits errors during the gap-filling step of base excision repair in Saccharomyces cerevisiae

We demonstrated recently that dominant negative mutants of rat DNA polymerase beta (Pol beta) interfere with repair of alkylation damage in Saccharomyces cerevisiae. To identify the alkylation repair pathway that is disrupted by the Pol beta dominant negative mutants, we studied the epistatic relationship of the dominant negative Pol beta mutants to genes known to be involved in repair of DNA alkylation damage in S. cerevisiae. We demonstrate that the rat Pol beta mutants interfere with the base excision repair pathway in S. cerevisiae. In addition, expression of one of the Pol beta dominant negative mutants, Pol beta-14, increases the spontaneous mutation rate of S. cerevisiae whereas expression of another Pol beta dominant negative mutant, Pol beta-TR, does not. Expression of the Pol beta-14 mutant in cells lacking APN1 activity does not result in an increase in the spontaneous mutation rate. These results suggest that gaps are required for mutagenesis to occur in the presence of Pol beta-14 but that it is not merely the presence of a gap that results in mutagenesis. Our results suggest that mutagenesis can occur during the gap-filling step of base excision repair in vivo.

Chain termination codons and polymerase-induced frameshift mutations

The consensus sequence for single-base deletions in non-reiterated runs during in vitro DNA-dependent DNA polymerisation is refined using data available in the literature. This leads to the observation that chain termination codons are hotspots for single-base deletions. The evolutionary implications are discussed in two models which differ in whether polymerases evolved while the genetic code emerged or after the genetic code was fixed. A possible answer to the question 'Why are stop codons just what they are?' is suggested.

A structural census of genomes: comparing bacterial, eukaryotic, and archaeal genomes in terms of protein structure

Representative genomes from each of the three kingdoms of life are compared in terms of protein structure, in particular, those of Haemophilus influenzae (a bacteria), Methanococcus jannaschii (an archaeon), and yeast (a eukaryote). The comparison is in the form of a census (or comprehensive accounting) of the relative occurrence of secondary and tertiary structures in the genomes, which particular emphasis on patterns of supersecondary structure. Comparison of secondary structure shows that the three genomes have nearly the same overall secondary-structure content, although they differ markedly in amino acid composition. Comparison of super-secondary structure, using a novel "frequent-words" approach, shows that yeast has a preponderance of consecutive strands (e.g. beta-beta-beta patterns), Haemophilus, consecutive helices (alpha-alpha-alpha), and Methanococcus, alternating helix-strand structures (beta-alpha-beta). Yeast also has significantly more helical membrane proteins than the other two genomes, with most of the differences concentrated in proteins containing two transmembrane segments. Comparison of tertiary structure (by sequence matching and domain-level clustering) highlights the substantial duplication in each genome (approximately 30% to 50%), with the degree of duplication following similar patterns in all three. Many sequence families are shared among the genomes, with the degree of overlap between any two genomes being roughly similar. In total, the three genomes contain 148 of the approximately 300 known protein folds. Forty-five of these 148 that are present in all three genomes are especially enriched in mixed super-secondary structures (alpha/beta). Moreover, the five most common of these 45 (the "top-5") have a remarkably similar super-secondary structure architecture, containing a central sheet of parallel strands with helices packed onto at least one face and beta-alpha-beta connections between adjacent strands. These most basic molecular parts, which, presumably, were present in the last common ancestor to the three Kingdoms, include the TIM-barrel, Rossmann, flavodoxin, thiamin-binding, and P-loop-hydrolase folds.

Characterization of an African swine fever virus 20-kDa DNA polymerase involved in DNA repair

African swine fever virus (ASFV) encodes a novel DNA polymerase, constituted of only 174 amino acids, belonging to the polymerase (pol) X family of DNA polymerases. Biochemical analyses of the purified enzyme indicate that ASFV pol X is a monomeric DNA-directed DNA polymerase, highly distributive, lacking a proofreading 3'-5'-exonuclease, and with a poor discrimination against dideoxynucleotides. A multiple alignment of family X DNA polymerases, together with the extrapolation to the crystal structure of mammalian DNA polymerase beta (pol beta), showed the conservation in ASFV pol X of the most critical residues involved in DNA binding, nucleotide binding, and catalysis of the polymerization reaction. Therefore, the 20-kDa ASFV pol X most likely represents the minimal functional version of an evolutionarily conserved pol beta-type DNA polymerase core, constituted by only the "palm" and "thumb" subdomains. It is worth noting that such an "unfingered" DNA polymerase is able to handle templated DNA polymerization with a considerable high fidelity at the base discrimination level. Base excision repair is considered to be a cellular defense mechanism repairing modified bases in DNA. Interestingly, the fact that ASFV pol X is able to conduct filling of a single nucleotide gap points to a putative role in base excision repair during the ASFV life cycle.

Characterization of new proteins found by analysis of short open reading frames from the full yeast genome

We have analysed short open reading frames (between 150 and 300 base pairs long) of the yeast genome (Saccharomyces cerevisiae) with a two-step strategy. The first step selects a candidate set of open reading frames from the DNA sequence based on statistical evaluation of DNA and protein sequence properties. The second step filters the candidate set by selecting open reading frames with high similarity to other known sequences (from any organism). As a result, we report ten new predicted proteins not present in the current sequence databases. These include a new alcohol dehydrogenase, a protein probably related to the cell cycle, as well as a homolog of the prokaryotic ribosomal protein L36 likely to be a mitochondrial ribosomal protein coded in the nuclear genome. We conclude that the analysis of short open reading frames leads to biologically interesting discoveries, even though the quantitative yield of new proteins is relatively low.

Protein structural classes in five complete genomes

The predicted distribution of globular proteins over folding types in five complete genomes differs from the tendencies observed in known protein structures. The ratio between the number of predicted membrane and globular proteins is conserved.

FKHL15, a new human member of the forkhead gene family located on chromosome 9q22

FKHL15 was isolated from a cDNA library enriched for transcripts from 9q22. Isolation and sequencing of a 3.5-kb cDNA clone identified a putative 376-amino-acid protein with greater than 80% homology over a 100-amino-acid stretch to the forkhead DNA-binding domain. The FKHL15 gene contains a region rich in alanine residues, frequently associated with transcriptional repression. The forkhead genes are believed to play important roles in development and differentiation in many different organisms and have also been implicated in the development of some tumors. The map position of FKHL15 on 9q22 places the gene within the candidate regions for the cancer predisposition syndrome multiple self-healing squamous epitheliomata and the degenerative neurological disorder hereditary sensory neuropathy type I. This is a region frequently lost in squamous cell cancer.

The sequence of 32b on the left arm of yeast chromosome XII reveals six known genes, a new member of the seripauperins family and a new ABS transporter homologous to the human multidrug resistance protein

The analysis of a 32 kb DNA fragment from cosmid 2G12 on the left arm of chromosome XII identifies 14 open reading frames (ORFs) numbered L0948 to L1325, a new tRNA for proline, a delta remnant and two putative ARS. Six ORFs have been previously identified: HSP104, SSA2, SPA2, KNS1, DPS1/APS and SDC25. Three putative ORFs have significant homology with known proteins: L0968 is a new member of the very large 'seripauperins' family, comprising at least 20 yeast members; L1313 is a new ABC transporter highly homologous to the yeast cadmium resistance protein Ycf1p and to the human multidrug resistance protein hMRP1; the C-terminal part of L1325 present in our sequence is very homologous to the fruit fly abdominal segment formation protein Pumilio. Finally, two ORFs, L1201 and L1205, have weak homology with two yeast hypothetical proteins of unknown function identified by the yeast systematic sequencing genome. Since our nucleotide sequence overlaps by 11.6 kb the cosmid 2B18 sequenced by Miosga and Zimmerman (1996) on the right end, we have not reported here the analysis of the ORFs L1313, L1321 and L1325.

Complete genome sequences of cellular life forms: glimpses of theoretical evolutionary genomics

The availability of complete genome sequences of cellular life forms creates the opportunity to explore the functional content of the genomes and evolutionary relationships between them at a new qualitative level. With the advent of these sequences, the construction of a minimal gene set sufficient for sustaining cellular life and reconstruction of the genome of the last common ancestor of bacteria, eukaryotes, and archaea become realistic, albeit challenging, research projects. A version of the minimal gene set for modern-type cellular life derived by comparative analysis of two bacterial genomes, those of Haemophilus influenzae and Mycoplasma genitalium, consists of approximately 250 genes. A comparison of the protein sequences encoded in these genes with those of the proteins encoded in the complete yeast genome suggests that the last common ancestor of all extant life might have had an RNA genome.

Functional analysis of YCL09C: evidence for a role as the regulatory subunit of acetolactate synthase

We have analysed the function of the open reading frame (ORF) YCL09C. The deletion of this ORF from chromosome III does not affect the physiology of the corresponding yeast strain enough to give a distinct phenotype. Nevertheless a computational analysis reveals high homology between this ORF and the enterobacterial genes encoding the regulatory subunit of acetolactate synthase. We have therefore tested the possibility that yc109cp is the regulatory subunit of yeast acetolactate synthase by in vitro enzymatic analysis. The acetolactate synthase was previously shown to be retroinhibited by its final product valine. In Escherichia coli this retro-control is assured by the regulatory subunit. Using a yeast strain carrying a complete deletion of YCL09C, we have observed the loss of such retro-inhibition. These results together with the computational predictions show that YCL09C encodes the regulatory subunit of yeast acetolactate synthase.

Characterization of the human forkhead gene FREAC-4. Evidence for regulation by Wilms' tumor suppressor gene (WT-1) and p53

We describe the cloning and sequence analysis of a nearly full-length cDNA as well as a corresponding 5.2-kilobase pair genomic fragment encoding FREAC-4, a member of the forkhead family of transcription factors. The cDNA is collinear with respect to the coding region of the intronless genomic clone. The conceptual translation product predicts a protein of 465 amino acids with a hyperacidic amino-terminal end, a DNA binding forkhead domain and a carboxyl-terminal part that is rich in homopolymeric runs of prolines and alanines. The transcription start is identified using an RNase protection assay. A 2.7-kilobase pair genomic DNA fragment, located immediately upstream of the translation start, was fused to a luciferase reporter gene. Significant levels of luciferase activity were detected when this construct was transfected into two kidney-derived cell lines, 293 and COS-7 cells, whereas only background reporter gene expression was observed in a cell line of nonkidney origin. Cotransfections with plasmids expressing WT-1, WTAR (a mutated form of WT-1), p53, and a mutated form of p53 revealed a complex pattern of regulation with a 3-fold induction with WT-1, a 7-fold induction with mutated p53, and a 4-fold repression with wild-type p53. A 5'-promoter deletion series delimits a DNA fragment necessary for WT-1 inducibility in cotransfection experiments. This fragment is shown to contain at least one cis-element that is capable of interacting with recombinant WT-1.

PWP2, a member of the WD-repeat family of proteins, is an essential Saccharomyces cerevisiae gene involved in cell separation

WD-repeat proteins contain four to eight copies of a conserved motif that usually ends with a tryptophan-aspartate (WD) dipeptide. The Saccharomyces cerevisiae PWP2 gene, identified by sequencing of chromosome III, is predicted to contain eight so-called WD-repeats, flanked by nonhomologous extensions. This gene is expressed as a 3.2-kb mRNA in all cell types and encodes a protein of 104 kDa. The PWP2 gene is essential for growth because spores carrying the pwp2 delta 1::HIS3 disruption germinate before arresting growth with one or two large buds. The growth defect of pwp2 delta 1::HIS3 cells was rescued by expression of PWP2 or epitope-tagged HA-PWP2 using the galactose-inducible GALI promoter. In the absence of galactose, depletion of Pwp2p resulted in multibudded cells with defects in bud site selection, cytokinesis, and hydrolysis of the septal junction between mother and daughter cells. In cell fractionation studies, HA-Pwp2p was localized in the particulate component of cell lysates, from which it would be solubulized by high salt and alkaline buffer but not by nonionic detergents or urea. Indirect immunofluorescence microscopy indicated that HA-Pwp2p was clustered at multiple points in the cytoplasm. These results suggest that Pwp2p exists in a proteinaceous complex, possibly associated with the cytoskeleton, where it functions in control of cell growth and separation.

Five years on the wings of fork head

Since its discovery five years ago the conserved family of fork head/HNF-3-related transcription factors has gained increasing importance for the analysis of gene regulatory mechanisms during embryonic development and in differentiated cells. Different members of this family, which is defined by a conserved 110 amino acid residues encompassing DNA binding domain of winged helix structure, serve as regulatory keys in embryogenesis, in tumorigenesis or in the maintenance of differentiated cell states. The purpose of this review is to summarize the accumulating amount of data on structure, expression and function of fork head/HNF-3-related transcription factors.

Discrimination between fortuitous and biologically constrained open reading frames in DNA sequences of Saccharomyces cerevisiae

The systematic sequencing of the yeast genome has raised the problem of the biological significance of the open reading frames (ORFs) revealed: it is possible that some of these are fortuitous. To avoid the analysis of such fortuitous ORFs, a minimum length of 100 sense codons was adopted. Nevertheless, the presence of fortuitous ORFs of more than 100 codons cannot be excluded. Thus, in the context of functional analysis, a method for discrimination between fortuitous and biologically active ORFs may be useful. The discrimination method described here is based on multiple criteria: ORF length, codon bias, and both amino-acid and dipeptide composition of the corresponding polypeptide. The thresholds for each criterion are based on the comparison between two learning sets: one drawn from random DNA sequences and the second from known genes. The method was validated by two test sets (one random and one biological) and then applied to the ORFs of chromosomes I, II, III, V, VIII, IX and XI. This method predicts 123 fortuitous ORFs among the 1773 identified on these chromosomes.

Differential activation of lung-specific genes by two forkhead proteins, FREAC-1 and FREAC-2

We describe the cDNA sequences for two human transcription factors, Forkhead RElated ACtivator (FREAC)-1 and -2, that belong to the forkhead family of eukaryotic DNA binding proteins. FREAC-1 and -2 are encoded by distinct genes, are almost identical within their DNA binding domains and in the COOH termini, but are otherwise divergent. Cotransfections with a reporter carrying FREAC binding sites showed that both proteins are transcriptional activators, and deletions located the activation domains to the COOH-terminal side of the forkhead domains. Expression of FREAC-1 and FREAC-2 is restricted to lung and placenta. We show that the promoters of genes for lung-specific proteins such as pulmonary surfactant proteins A, B, and C (SPA, SPB, and SPC) and the Clara cell 10-kDa protein (CC10) contain potential binding sites for FREAC-1 and FREAC-2. DNaseI footprinting verified that FREAC proteins bind to the predicted sites in the CC10 and SPB promoters. While an SPB promoter construct could be transactivated by both FREAC-1 and FREAC-2, CC10 was only activated by FREAC-1. Efficient activation of CC10 by FREAC-1 is shown to be specific for a lung cell line with Clara cell characteristics (H441) and to involve a region of the FREAC-1 protein unable to activate in other cell types.

Dominant negative rat DNA polymerase beta mutants interfere with base excision repair in Saccharomyces cerevisiae

DNA polymerase beta is one of the smallest known eukaryotic DNA polymerases. This polymerase has been very well characterized in vitro, but its functional role in vivo has yet to be determined. Using a novel competition assay in Escherichia coli, we isolated two DNA polymerase beta dominant negative mutants. When we overexpressed the dominant negative mutant proteins in Saccharomyces cerevisiae, the cells became sensitive to methyl methanesulfonate. Interestingly, overexpression of the same polymerase beta mutant proteins did not confer sensitivity to UV damage, strongly suggesting that the mutant proteins interfere with the process of base excision repair but not nucleotide excision repair in S. cerevisiae. Our data implicate a role for polymerase IV, the S. cerevisiae polymerase beta homolog, in base excision repair in S. cerevisiae.

New protein functions in yeast chromosome VIII

The analysis of the 269 open reading frames of yeast chromosome VIII by computational methods has yielded 24 new significant sequence similarities to proteins of known function. The resulting predicted functions include three particularly interesting cases of translation-associated proteins: peptidyl-tRNA hydrolase, a ribosome recycling factor homologue, and a protein similar to cytochrome b translational activator CBS2. The methodological limits of the meaningful transfer of functional information between distant homologues are discussed.

A Drosophila hsp70 gene contains long, antiparallel, coupled open reading frames (LAC ORFs) conserved in homologous loci

A clone isolated from a Drosophila auraria heat-shock cDNA library presents two long, antiparallel, coupled (LAC) open reading frames (ORFs). One strand ORF is 1,929 nucleotides long and exhibits great identity (87.5% at the nucleotide level and 94% at the amino acid level) with the hsp70 gene copies of D. melanogaster, while the second strand ORF, in antiparallel in-frame register arrangement, is 1,839 nucleotides long and exhibits 32% identity with a putative, recently identified, NAD(+)-dependent glutamate dehydrogenase (NAD(+)-GDH). The overlap of the two ORFs is 1,824 nucleotides long. Computational analysis shows that this LAC ORF arrangement is conserved in other hsp70 loci in a wide range of organisms, raising questions about possible evolutionary benefits of such a peculiar genomic organization.

A 29.425 kb segment on the left arm of yeast chromosome XV contains more than twice as many unknown as known open reading frames

The nucleotide sequence of a 29.425 kb fragment localized on the left arm of chromosome XV from Saccharomyces cerevisiae has been determined. The sequence contains 13 open reading frames (ORFs) of which four encode the known genes ADH1, COQ3, MSH2 and RCF4. Predictions are made concerning the functions of the unknown ORFs. Some of the ORFs contain sequences similar to expressed sequence tags (EST) found in the database made available by TIGR. In particular, the highly expressed ADH1 gene is represented in this database by no less than 20 EST sequences. Two ARS sequences and a putative functional GCN4 motif have also been detected. One ORF (O0953) containing nine putative transmembrane segments is similar to a hypothetical membrane protein of Arabidopsis thaliana. Characteristic features of the other ORFs include ATP/GTP binding sites, a fungal Zn(2)-Cys(6) binuclear centre, an endoplasmic reticulum targeting sequence, a beta-transducin repeat signature and in two instances, good similarity to the prokaryotic lipoprotein signal peptide motif.

Exploring the Mycoplasma capricolum genome: a minimal cell reveals its physiology

We report on the analysis of 214kb of the parasitic eubacterium Mycoplasma capricolum sequenced by genomic walking techniques. The 287 putative proteins detected to date represent about half of the estimated total number of 500 predicted for this organism. A large fraction of these (75%) can be assigned a likely function as a result of similarity searches. Several important features of the functional organization of this small genome are already apparent. Among these are (i) the expected relatively large number of enzymes involved in metabolic transport and activation, for efficient use of host cell nutrients; (ii) the presence of anabolic enzymes; (iii) the unexpected diversity of enzymes involved in DNA replication and repair; and (iv) a sizeable number of orthologues (82 so far) in Escherichia coli. This survey is beginning to provide a detailed view of how M. capricolum manages to maintain essential cellular processes with a genome much smaller than that of its bacterial relatives.

Novel protein families in archaean genomes

In a quest for novel functions in archaea, all archaean hypothetical open reading frames (ORFs), as annotated in the Swiss-Prot protein sequence database, were used to search the latest databases for the identification of characterized homologues. Of the 95 hypothetical archaean ORFs, 25 were found to be homologous to another hypothetical archaean ORF, while 36 were homologous to non-archaean proteins, of which as many as 30 were homologous to a characterized protein family. Thus the level of sequence similarity in this set reaches 64%, while the level of function assignment is only 32%. Of the ORFs with predicted functions, 12 homologies are reported here for the first time and represent nine new functions and one gene duplication at an acetyl-coA synthetase locus. The novel functions include components of the transcriptional and translational apparatus, such as ribosomal proteins, modification enzymes and a translation initiation factor. In addition, new enzymes are identified in archaea, such as cobyric acid synthase, dCTP deaminase and the first archaean homologues of a new subclass of ATP binding proteins found in fungi. Finally, it is shown that the putative laminin receptor family of eukaryotes and an archaean homologue belong to the previously characterized ribosomal protein family S2 from eubacteria. From the present and previous work, the major implication is that archaea seem to have a mode of expression of genetic information rather similar to eukaryotes, while eubacteria may have proceeded into unique ways of transcription and translation. In addition, with the detection of proteins in various metabolic and genetic processes in archaea, we can further predict the presence of additional proteins involved in these processes.