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.

MethDB—a public database for DNA methylation data

Methylation of cytosine in the 5 position of the pyrimidine ring is a major modification of the DNA in most organisms. In eukaryotes, the distribution and number of 5-methylcytosines (5mC) along the DNA is heritable but can also change with the developmental state of the cell and as a response to modifications of the environment. While DNA methylation probably has a number of functions, scientific interest has recently focused on the gene silencing effect methylation can have in eukaryotic cells. In particular, the discovery of changes in the methylation level during cancer development has increased the interest in this field. In the past, a vast amount of data has been generated with different levels of resolution ranging from 5mC content of total DNA to the methylation status of single nucleotides. We present here a database for DNA methylation data that attempts to unify these results in a common resource. The database is accessible via WWW (http://www.methdb.de). It stores information about the origin of the investigated sample and the experimental procedure, and contains the DNA methylation data. Query masks allow for searching for 5mC content, species, tissue, gene, sex, phenotype, sequence ID and DNA type. The output lists all available information including the relative gene expression level. DNA methylation patterns and methylation profiles are shown both as a graphical representation and as G/A/T/C/5mC-sequences or tables with sequence positions and methylation levels, respectively.

Estimation of divergence times from multiprotein sequences for a few mammalian species and several distantly related organisms

When many protein sequences are available for estimating the time of divergence between two species, it is customary to estimate the time for each protein separately and then use the average for all proteins as the final estimate. However, it can be shown that this estimate generally has an upward bias, and that an unbiased estimate is obtained by using distances based on concatenated sequences. We have shown that two concatenation-based distances, i.e., average gamma distance weighted with sequence length (d2) and multiprotein gamma distance (d3), generally give more satisfactory results than other concatenation-based distances. Using these two distance measures for 104 protein sequences, we estimated the time of divergence between mice and rats to be approximately 33 million years ago. Similarly, the time of divergence between humans and rodents was estimated to be approximately 96 million years ago. We also investigated the dependency of time estimates on statistical methods and various assumptions made by using sequence data from eubacteria, protists, plants, fungi, and animals. Our best estimates of the times of divergence between eubacteria and eukaryotes, between protists and other eukaryotes, and between plants, fungi, and animals were 3, 1.7, and 1.3 billion years ago, respectively. However, estimates of ancient divergence times are subject to a substantial amount of error caused by uncertainty of the molecular clock, horizontal gene transfer, errors in sequence alignments, etc.

Crystal structure of cytochrome P450 14α-sterol demethylase (CYP51) from Mycobacterium tuberculosis in complex with azole inhibitors

Cytochrome P450 14α-sterol demethylases (CYP51) are essential enzymes in sterol biosynthesis in eukaryotes. CYP51 removes the 14α-methyl group from sterol precursors such as lanosterol, obtusifoliol, dihydrolanosterol, and 24(28)-methylene-24,25-dihydrolanosterol. Inhibitors of CYP51 include triazole antifungal agents fluconazole and itraconazole, drugs used in treatment of topical and systemic mycoses. The 2.1- and 2.2-Å crystal structures reported here for 4-phenylimidazole- and fluconazole-bound CYP51 from Mycobacterium tuberculosis (MTCYP51) are the first structures of an authentic P450 drug target. MTCYP51 exhibits the P450 fold with the exception of two striking differences—a bent I helix and an open conformation of BC loop—that define an active site-access channel running along the heme plane perpendicular to the direction observed for the substrate entry in P450BM3. Although a channel analogous to that in P450BM3 is evident also in MTCYP51, it is not open at the surface. The presence of two different channels, with one being open to the surface, suggests the possibility of conformationally regulated substrate-in/product-out openings in CYP51. Mapping mutations identified in Candida albicans azole-resistant isolates indicates that azole resistance in fungi develops in protein regions involved in orchestrating passage of CYP51 through different conformational stages along the catalytic cycle rather than in residues directly contacting fluconazole. These new structures provide a basis for rational design of new, more efficacious antifungal agents as well as insight into the molecular mechanism of P450 catalysis.

ClpA mediates directional translocation of substrate proteins into the ClpP protease

The intracellular degradation of many proteins is mediated in an ATP-dependent manner by large assemblies comprising a chaperone ring complex associated coaxially with a proteolytic cylinder, e.g., ClpAP, ClpXP, and HslUV in prokaryotes, and the 26S proteasome in eukaryotes. Recent studies of the chaperone ClpA indicate that it mediates ATP-dependent unfolding of substrate proteins and directs their ATP-dependent translocation into the ClpP protease. Because the axial passageway into the proteolytic chamber is narrow, it seems likely that unfolded substrate proteins are threaded from the chaperone into the protease, suggesting that translocation could be directional. We have investigated directionality in the ClpA/ClpP-mediated reaction by using two substrate proteins bearing the COOH-terminal ssrA recognition element, each labeled near the NH2 or COOH terminus with fluorescent probes. Time-dependent changes in both fluorescence anisotropy and fluorescence resonance energy transfer between donor fluorophores in the ClpP cavity and the substrate probes as acceptors were measured to monitor translocation of the substrates from ClpA into ClpP. We observed for both substrates that energy transfer occurs 2–4 s sooner with the COOH-terminally labeled molecules than with the NH2-terminally labeled ones, indicating that translocation is indeed directional, with the COOH terminus of the substrate protein entering ClpP first.

The length of telomeric G-rich strand 3′-overhang measured by oligonucleotide ligation assay

A typical G-rich telomeric DNA strand, which runs 5′→3′ toward the chromosome ends, protrudes by several nucleotides in lower eukaryotes. In human chromosomes long G-rich 3′-overhangs have been found. Apart from the standard G-rich tail, several non-canonical terminal structures have been proposed. However, the mechanism of long-tail formation, the presence and the role of these structures in telomere maintenance or shortening are not completely understood. In a search for a simple method to accurately measure the 3′-overhang we have established a protocol based on the ligation of telomeric oligonucleotide hybridized to non-denatured DNA under stringent conditions (oligonucleotide ligation assay with telomeric repeat oligonucleotide). This method enabled us to detect a large proportion of G-rich single-stranded telomeric DNA that was as short as 24 nt. Nevertheless, we showed G-tails longer than 400 nt. In all tested cells the lengths ranging from 108 to 270 nt represented only 37% of the whole molecule population, while 56–62% were <90 nt. Our protocol provides a simple and sensitive method for measuring the length of naturally occurring unpaired repeated DNA.

Substrate-specific inhibition of RecQ helicase

The RecQ helicases constitute a small but highly conserved helicase family. Proteins in this family are of particular interest because they are critical to maintenance of genomic stability in prokaryotes and eukaryotes. Eukaryotic RecQ helicase family members have been shown to unwind not only DNA duplexes but also DNAs with alternative structures, including structures stabilized by G quartets (G4 DNAs). We report that Escherichia coli RecQ can also unwind G4 DNAs, and that unwinding requires ATP and divalent cation. RecQ helicase is comparably active on duplex and G4 DNA substrates, as measured by direct comparison of protein activity and by competition assays. The porphyrin derivative, N-methyl mesoporphyrin IX (NMM), is a highly specific inhibitor of RecQ unwinding activity on G4 DNA but not duplex DNA: the inhibition constant (Ki) for NMM inhibition of G4 DNA unwinding is 1.7 µM, approximately two orders of magnitude below the Ki for inhibition of duplex DNA unwinding (>100 µM). NMM may therefore prove to be a valuable compound for substrate-specific inhibition of other RecQ family helicases in vitro and in vivo.

A Selaginella lepidophylla Trehalose-6-Phosphate Synthase Complements Growth and Stress-Tolerance Defects in a Yeast tps1 Mutant1

The accumulation of the disaccharide trehalose in anhydrobiotic organisms allows them to survive severe environmental stress. A plant cDNA, SlTPS1, encoding a 109-kD protein, was isolated from the resurrection plant Selaginella lepidophylla, which accumulates high levels of trehalose. Protein-sequence comparison showed that SlTPS1 shares high similarity to trehalose-6-phosphate synthase genes from prokaryotes and eukaryotes. SlTPS1 mRNA was constitutively expressed in S. lepidophylla. DNA gel-blot analysis indicated that SlTPS1 is present as a single-copy gene. Transformation of a Saccharomyces cerevisiae tps1Δ mutant disrupted in the ScTPS1 gene with S. lepidophylla SlTPS1 restored growth on fermentable sugars and the synthesis of trehalose at high levels. Moreover, the SlTPS1 gene introduced into the tps1Δ mutant was able to complement both deficiencies: sensitivity to sublethal heat treatment at 39°C and induced thermotolerance at 50°C. The osmosensitive phenotype of the yeast tps1Δ mutant grown in NaCl and sorbitol was also restored by the SlTPS1 gene. Thus, SlTPS1 protein is a functional plant homolog capable of sustaining trehalose biosynthesis and could play a major role in stress tolerance in S. lepidophylla.

Mitochondrial Single-stranded DNA-binding Protein Is Required for Mitochondrial DNA Replication and Development in Drosophila melanogaster

The discovery that several inherited human diseases are caused by mtDNA depletion has led to an increased interest in the replication and maintenance of mtDNA. We have isolated a new mutant in the lopo (low power) gene from Drosophila melanogaster affecting the mitochondrial single-stranded DNA-binding protein (mtSSB), which is one of the key components in mtDNA replication and maintenance. lopo1 mutants die late in the third instar before completion of metamorphosis because of a failure in cell proliferation. Molecular, histochemical, and physiological experiments show a drastic decrease in mtDNA content that is coupled with the loss of respiration in these mutants. However, the number and morphology of mitochondria are not greatly affected. Immunocytochemical analysis shows that mtSSB is expressed in all tissues but is highly enriched in proliferating tissues and in the developing oocyte. lopo1 is the first mtSSB mutant in higher eukaryotes, and its analysis demonstrates the essential function of this gene in development, providing an excellent model to study mitochondrial biogenesis in animals.

Dynamics of Fos-Jun-NFAT1 complexes

Transcription initiation in eukaryotes is controlled by nucleoprotein complexes formed through cooperative interactions among multiple transcription regulatory proteins. These complexes may be assembled via stochastic collisions or defined pathways. We investigated the dynamics of Fos-Jun-NFAT1 complexes by using a multicolor fluorescence resonance energy transfer assay. Fos-Jun heterodimers can bind to AP-1 sites in two opposite orientations, only one of which is populated in mature Fos-Jun-NFAT1 complexes. We studied the reversal of Fos-Jun binding orientation in response to NFAT1 by measuring the efficiencies of energy transfer from donor fluorophores linked to opposite ends of an oligonucleotide to an acceptor fluorophore linked to one subunit of the heterodimer. The reorientation of Fos-Jun by NFAT1 was not inhibited by competitor oligonucleotides or heterodimers. The rate of Fos-Jun reorientation was faster than the rate of heterodimer dissociation at some binding sites. The facilitated reorientation of Fos-Jun heterodimers therefore can enhance the efficiency of Fos-Jun-NFAT1 complex formation. We also examined the influence of the preferred orientation of Fos-Jun binding on the stability and transcriptional activity of Fos-Jun-NFAT1 complexes. Complexes formed at sites where Fos-Jun favored the same binding orientation in the presence and absence of NFAT1 exhibited an 8-fold slower dissociation rate than complexes formed at sites where Fos-Jun favored the opposite binding orientation. Fos-Jun-NFAT1 complexes also exhibited greater transcription activation at promoter elements that favored the same orientation of Fos-Jun binding in the presence and absence of NFAT1. Thus, the orientation of heterodimer binding can influence both the dynamics and promoter selectivity of multiprotein transcription regulatory complexes.

The HMG-domain protein BAP111 is important for the function of the BRM chromatin-remodeling complex in vivo

The Drosophila trithorax group gene brahma (brm) encodes the ATPase subunit of a SWI/SNF-like chromatin-remodeling complex. A key question about chromatin-remodeling complexes is how they interact with DNA, particularly in the large genomes of higher eukaryotes. Here, we report the characterization of BAP111, a BRM-associated protein that contains a high mobility group (HMG) domain predicted to bind distorted or bent DNA. The presence of an HMG domain in BAP111 suggests that it may modulate interactions between the BRM complex and chromatin. BAP111 is an abundant nuclear protein that is present in all cells throughout development. By using gel filtration chromatography and immunoprecipitation assays, we found that the majority of BAP111 protein in embryos is associated with the BRM complex. Furthermore, heterozygosity for BAP111 enhanced the phenotypes resulting from a partial loss of brm function. These data demonstrate that the BAP111 subunit is important for BRM complex function in vivo.

Two enzymes of diacylglyceryl-O-4′-(N,N,N,-trimethyl)homoserine biosynthesis are encoded by btaA and btaB in the purple bacterium Rhodobacter sphaeroides

Betaine lipids are ether-linked, nonphosphorous glycerolipids that resemble the more commonly known phosphatidylcholine in overall structure. Betaine lipids are abundant in many eukaryotes such as nonseed plants, algae, fungi, and amoeba. Some of these organisms are entirely devoid of phosphatidylcholine and, instead, contain a betaine lipid such as diacylglyceryl-O-4′-(N,N,N,-trimethyl)homoserine. Recently, this lipid also was discovered in the photosynthetic purple bacterium Rhodobacter sphaeroides where it seems to replace phosphatidylcholine under phosphate-limiting growth conditions. This discovery provided the opportunity to study the biosynthesis of betaine lipids in a bacterial model system. Mutants of R. sphaeroides deficient in the biosynthesis of the betaine lipid were isolated, and two genes essential for this process, btaA and btaB, were identified. It is proposed that btaA encodes an S-adenosylmethionine:diacylglycerol 3-amino-3-carboxypropyl transferase and btaB an S-adenosylmethionine-dependent N-methyltransferase. Both enzymatic activities can account for all reactions of betaine lipid head group biosynthesis. Because the equivalent reactions have been proposed for different eukaryotes, it seems likely that orthologs of btaA/btaB may be present in other betaine lipid-containing organisms.

RNA polymerase II transcription initiation: A structural view

In eukaryotes, RNA polymerase II transcribes messenger RNAs and several small nuclear RNAs. Like RNA polymerases I and III, polymerase II cannot act alone. Instead, general initiation factors [transcription factor (TF) IIB, TFIID, TFIIE, TFIIF, and TFIIH] assemble on promoter DNA with polymerase II, creating a large multiprotein–DNA complex that supports accurate initiation. Another group of accessory factors, transcriptional activators and coactivators, regulate the rate of RNA synthesis from each gene in response to various developmental and environmental signals. Our current knowledge of this complex macromolecular machinery is reviewed in detail, with particular emphasis on insights gained from structural studies of transcription factors.

Origin of Genes

We discuss two tests of the hypothesis that the first genes were assembled from exons. The hypothesis of exon shuffling in the progenote predicts that intron phases will be correlated so that exons will be an integer number of codons and predicts that the exons will be correlated with compact regions of polypeptide chain. These predictions have been tested on ancient conserved proteins (proteins without introns in prokaryotes but with introns in eukaryotes) and hold with high statistical significance. We conclude that introns are correlated with compact features of proteins 15-, 22-, or 30-amino acid residues long, as was predicted by “The Exon Theory of Genes.”