The B12-dependent ribonucleotide reductase from the archaebacterium Thermoplasma acidophila: an evolutionary solution to the ribonucleotide reductase conundrum

The case for an early biological origin of DNA

All life generates deoxyribonucleotides, the building blocks of DNA, via ribonucleotide reductases (RNRs). The complexity of this reaction suggests it did not evolve until well after the advent of templated protein synthesis, which in turn suggests DNA evolved later than both RNA and templated protein synthesis. However, deoxyribonucleotides may have first been synthesised via an alternative, chemically simpler route—the reversal of the deoxyriboaldolase (DERA) step in deoxyribonucleotide salvage. In light of recent work demonstrating that this reaction can drive synthesis of deoxyribonucleosides, we consider what pressures early adoption of this pathway would have placed on cell metabolism. This in turn provides a rationale for the replacement of DERA-dependent DNA production by RNR-dependent production.

Targeting the Large Subunit of Human Ribonucleotide Reductase for Cancer Chemotherapy

Ribonucleotide reductase (RR) is a crucial enzyme in de novo DNA synthesis, where it catalyses the rate determining step of dNTP synthesis. RRs consist of a large subunit called RR1 (α), that contains two allosteric sites and one catalytic site, and a small subunit called RR2 (β), which houses a tyrosyl free radical essential for initiating catalysis. The active form of mammalian RR is an α_nβ_m hetero oligomer. RR inhibitors are cytotoxic to proliferating cancer cells. In this brief review we will discuss the three classes of RR, the catalytic mechanism of RR, the regulation of the dNTP pool, the substrate selection, the allosteric activation, inactivation by ATP and dATP, and the nucleoside drugs that target RR. We will also discuss possible strategies for developing a new class of drugs that disrupts the RR assembly.

Molecular paleoscience: systems biology from the past

Experimental paleomolecular biology, paleobiochemistry, and paleogenetics are closely related emerging fields that infer the sequences of ancient genes and proteins from now-extinct organisms, and then resurrect them for study in the laboratory. The goal of paleogenetics is to use information from natural history to solve the conundrum of modern genomics: How can we understand deeply the function of biomolecular structures uncovered and described by modern chemical biology? Reviewed here are the first 20 cases where biomolecular resurrections have been achieved. These show how paleogenetics can lead to an understanding of the function of biomolecules, analyze changing function, and put meaning to genomic sequences, all in ways that are not possible with traditional molecular biological studies.

The cbiS gene of the archaeon Methanopyrus kandleri AV19 encodes a bifunctional enzyme with adenosylcobinamide amidohydrolase and alpha-ribazole-phosphate phosphatase activities

Here we report the initial biochemical characterization of the bifunctional alpha-ribazole-P (alpha-RP) phosphatase, adenosylcobinamide (AdoCbi) amidohydrolase CbiS enzyme from the hyperthermophilic methanogenic archaeon Methanopyrus kandleri AV19. The cbiS gene encodes a 39-kDa protein with two distinct segments, one of which is homologous to the AdoCbi amidohydrolase (CbiZ, EC 3.5.1.90) enzyme and the other of which is homologous to the recently discovered archaeal alpha-RP phosphatase (CobZ, EC 3.1.3.73) enzyme. CbiS function restored AdoCbi salvaging and alpha-RP phosphatase activity in strains of the bacterium Salmonella enterica where either step was blocked. The two halves of the cbiS genes retained their function in vivo when they were cloned separately. The CbiS enzyme was overproduced in Escherichia coli and was isolated to >95% homogeneity. High-performance liquid chromatography, UV-visible spectroscopy, and mass spectroscopy established alpha-ribazole and cobyric acid as the products of the phosphatase and amidohydrolase reactions, respectively. Reasons why the CbiZ and CobZ enzymes are fused in some archaea are discussed.

Studies of the CobA-type ATP:Co(I)rrinoid adenosyltransferase enzyme of Methanosarcina mazei strain Go1

Although methanogenic archaea use B(12) extensively as a methyl carrier for methanogenesis, little is known about B(12) metabolism in these prokaryotes or any other archaea. To improve our understanding of how B(12) metabolism differs between bacteria and archaea, the gene encoding the ATP:co(I)rrinoid adenosyltransferase in Methanosarcina mazei strain Gö1 (open reading frame MM3138, referred to as cobA(Mm) here) was cloned and used to restore coenzyme B(12) synthesis in a Salmonella enterica strain lacking the housekeeping CobA enzyme. cobA(Mm) protein was purified and its initial biochemical analysis performed. In vitro, the activity is enhanced 2.5-fold by the addition of Ca(2+) ions, but the activity was not enhanced by Mg(2+) and, unlike the S. enterica CobA enzyme, it was >50% inhibited by Mn(2+). The CobA(Mm) enzyme had a K(m)(ATP) of 3 microM and a K(m)(HOCbl) of 1 microM. Unlike the S. enterica enzyme, CobA(Mm) used cobalamin (Cbl) as a substrate better than cobinamide (Cbi; a Cbl precursor); the beta phosphate of ATP was required for binding to the enzyme. A striking difference between CobA(Se) and CobA(Mm) was the use of ADP as a substrate by CobA(Mm), suggesting an important role for the gamma phosphate of ATP in binding. The results from (31)P-nuclear magnetic resonance spectroscopy experiments showed that triphosphate (PPP(i)) is the reaction by-product; no cleavage of PPP(i) was observed, and the enzyme was only slightly inhibited by pyrophosphate (PP(i)). The data suggested substantial variations in ATP binding and probably corrinoid binding between CobA(Se) and CobA(Mm) enzymes.

Integrating protein structures and precomputed genealogies in the Magnum database: examples with cellular retinoid binding proteins

BACKGROUND

When accurate models for the divergent evolution of protein sequences are integrated with complementary biological information, such as folded protein structures, analyses of the combined data often lead to new hypotheses about molecular physiology. This represents an excellent example of how bioinformatics can be used to guide experimental research. However, progress in this direction has been slowed by the lack of a publicly available resource suitable for general use.

RESULTS

The precomputed Magnum database offers a solution to this problem for ca. 1,800 full-length protein families with at least one crystal structure. The Magnum deliverables include 1) multiple sequence alignments, 2) mapping of alignment sites to crystal structure sites, 3) phylogenetic trees, 4) inferred ancestral sequences at internal tree nodes, and 5) amino acid replacements along tree branches. Comprehensive evaluations revealed that the automated procedures used to construct Magnum produced accurate models of how proteins divergently evolve, or genealogies, and correctly integrated these with the structural data. To demonstrate Magnum's capabilities, we asked for amino acid replacements requiring three nucleotide substitutions, located at internal protein structure sites, and occurring on short phylogenetic tree branches. In the cellular retinoid binding protein family a site that potentially modulates ligand binding affinity was discovered. Recruitment of cellular retinol binding protein to function as a lens crystallin in the diurnal gecko afforded another opportunity to showcase the predictive value of a browsable database containing branch replacement patterns integrated with protein structures.

CONCLUSION

We integrated two areas of protein science, evolution and structure, on a large scale and created a precomputed database, known as Magnum, which is the first freely available resource of its kind. Magnum provides evolutionary and structural bioinformatics resources that are useful for identifying experimentally testable hypotheses about the molecular basis of protein behaviors and functions, as illustrated with the examples from the cellular retinoid binding proteins.

The last common ancestor: what's in a name?

Twenty completely sequenced cellular genomes from the three major domains were analyzed using twice one-way BLAST searches in order to define the set of the most conserved protein-encoding sequences to characterize the gene complement of the last common ancestor of extant life. The resulting set is dominated by different putative ATPases, and by molecules involved in gene expression and RNA metabolism. DEAD-type RNA helicase and enolase genes, which are known to be part of the RNA degradosome, are as conserved as many transcription and translation genes. This suggests the early evolution of a control mechanism for gene expression at the RNA level, providing additional support to the hypothesis that during early cellular evolution RNA molecules played a more prominent role. Conserved sequences related to biosynthetic pathways include those encoding putative phosphoribosyl pyrophosphate synthase and thioredoxin, which participate in nucleotide metabolism. Although the information contained in the available databases corresponds only to a minor portion of biological diversity, the sequences reported here are likely to be part of an essential and highly conserved pool of proteins domains common to all organisms.

ABC transporter for corrinoids in Halobacterium sp. strain NRC-1

We report evidence for the existence of a putative ABC transporter for corrinoid utilization in the extremely halophilic archaeon Halobacterium sp. strain NRC-1. Results from genetic and nutritional analyses of Halobacterium showed that mutants with lesions in open reading frames (ORFs) Vng1370G, Vng1371Gm, and Vng1369G required a 10(5)-fold higher concentration of cobalamin for growth than the wild-type or parent strain. The data support the conclusion that these ORFs encode orthologs of the bacterial cobalamin ABC transporter permease (btuC; Vng1370G), ATPase (btuD; Vng1371Gm), and substrate-binding protein (btuF; Vng1369G) components. Mutations in the Vng1370G, Vng1371Gm, and Vng1369G genes were epistatic, consistent with the hypothesis that their products work together to accomplish the same function. Extracts of btuF mutant strains grown in the presence of cobalamin did not contain any cobalamin molecules detectable by a sensitive bioassay, whereas btuCD mutant strain extracts did. The data are consistent with the hypothesis that the BtuF protein is exported to the extracellular side of the cell membrane, where it can bind cobalamin in the absence of BtuC and BtuD. Our data also provide evidence for the regulation of corrinoid transport and biosynthesis. Halobacterium synthesized cobalamin in a chemically defined medium lacking corrinoid precursors. To the best of our knowledge, this is the first genetic analysis of an archaeal corrinoid transport system.

Modern mRNA proofreading and repair: clues that the last universal common ancestor possessed an RNA genome?

RNA repair has now been demonstrated to be a genuine biological process and appears to be present in all three domains of life. In this article, we consider what this might mean for the transition from an early RNA-dominated world to modern cells possessing genetically encoded proteins and DNA. There are significant gaps in our understanding of how the modern protein-DNA world could have evolved from a simpler system, and it is currently uncertain whether DNA genomes evolved once or twice. Against this backdrop, the discovery of RNA repair in modern cells is timely food for thought and brings us conceptually one step closer to understanding how RNA genomes were replaced by DNA genomes. We have examined the available literature on multisubunit RNA polymerase structure and function and conclude that a strong case can be made that the Last Universal Common Ancestor (LUCA) possessed a repair-competent RNA polymerase, which would have been capable of acting on an RNA genome. However, while this lends credibility to the proposal that the LUCA had an RNA genome, the alternative, that LUCA had a DNA genome, cannot be completely ruled out.

A new pathway for salvaging the coenzyme B12 precursor cobinamide in archaea requires cobinamide-phosphate synthase (CbiB) enzyme activity

The ability of archaea to salvage cobinamide has been under question because archaeal genomes lack orthologs to the bacterial nucleoside triphosphate:5'-deoxycobinamide kinase enzyme (cobU in Salmonella enterica). The latter activity is required for cobinamide salvaging in bacteria. This paper reports evidence that archaea salvage cobinamide from the environment by using a pathway different from the one used by bacteria. These studies demanded the functional characterization of two genes whose putative function had been annotated based solely on their homology to the bacterial genes encoding adenosylcobyric acid and adenosylcobinamide-phosphate synthases (cbiP and cbiB, respectively) of S. enterica. A cbiP mutant strain of the archaeon Halobacterium sp. strain NRC-1 was auxotrophic for adenosylcobyric acid, a known intermediate of the de novo cobamide biosynthesis pathway, but efficiently salvaged cobinamide from the environment, suggesting the existence of a salvaging pathway in this archaeon. A cbiB mutant strain of Halobacterium was auxotrophic for adenosylcobinamide-GDP, a known de novo intermediate, and did not salvage cobinamide. The results of the nutritional analyses of the cbiP and cbiB mutants suggested that the entry point for cobinamide salvaging is adenosylcobyric acid. The data are consistent with a salvaging pathway for cobinamide in which an amidohydrolase enzyme cleaves off the aminopropanol moiety of adenosylcobinamide to yield adenosylcobyric acid, which is converted by the adenosylcobinamide-phosphate synthase enzyme to adenosylcobinamide-phosphate, a known intermediate of the de novo biosynthetic pathway. The existence of an adenosylcobinamide amidohydrolase enzyme would explain the lack of an adenosylcobinamide kinase in archaea.

The first experimental test of the hypothesis that enzymes have evolved to enhance hydrogen tunneling

The literature hypothesis that "the optimization of enzyme catalysis may entail the evolutionary implementation of chemical strategies that increase the probability of quantum-mechanical tunneling" is experimentally tested herein for the first time. The system employed is the key to being able to provide this first experimental test of the "enhanced hydrogen tunneling" hypothesis, one that requires a comparison of the three criteria diagnostic of tunneling (vide infra) for the same, or nearly the same, reaction with and without the enzyme. Specifically, studied herein are the adenosylcobalamin (AdoCbl, also known as coenzyme B(12))-dependent diol dehydratase model reactions of (i). H(D)(*) atom abstraction from ethylene glycol-d(0) and ethylene glycol-d(4) solvent by 5'-deoxyadenosyl radical (Ado(*)) and (ii.) the same H(*) abstraction reactions by the 8-methoxy-5'-deoxyadenosyl radical (8-MeOAdo(*)). The Ado(*) and 8-MeOAdo(*) radicals are generated by Co-C thermolysis of their respective precursors, AdoCbl and 8-MeOAdoCbl. Deuterium kinetic isotope effects (KIEs) of the H(*)(D(*)) abstraction reactions from ethylene glycol have been measured over a temperature range of 80-120 degrees C: KIE = 12.4 +/- 1.1 at 80 degrees C for Ado(*) and KIE = 12.5 +/- 0.9 at 80 degrees C for 8-MeOAdo(*) (values ca. 2-fold that of the predicted maximum primary times secondary ground-state zero-point energy (GS-ZPE) KIE of 6.4 at 80 degrees C). From the temperature dependence of the KIEs, zero-point activation energy differences ([E(D) - E(H)]) of 3.0 +/- 0.3 kcal mol(-)(1) for Ado(*) and 2.1 +/- 0.6 kcal mol(-)(1) for 8-MeOAdo(*) have been obtained, both of which are significantly larger than the nontunneling, zero-point energy only maximum of 1.2 kcal mol(-)(1). Pre-exponential factor ratios (A(H)/A(D)) of 0.16 +/- 0.07 for Ado(*) and 0.5 +/- 0.4 for 8-MeOAdo(*) are observed, both of which are significantly less than the 0.7 minimum for nontunneling behavior. The data provide strong evidence for the expected quantum mechanical tunneling in the Ado(*) and 8-MeOAdo(*)-mediated H(*) abstraction reactions from ethylene glycol. More importantly, a comparison of these enzyme-free tunneling data to the same KIE, (E(D) - E(H)) and A(H)/A(D) data for a closely related, Ado(*)-mediated H(*) abstraction reaction from a primary CH(3)- group in AdoCbl-dependent methylmalonyl-CoA mutase shows the enzymic and enzyme-free data sets are identical within experimental error. The Occam's Razor conclusion is that at least this adenosylcobalamin-dependent enzyme has not evolved to enhance quantum mechanical tunneling, at least within the present error bars. Instead, this B(12)-dependent enzyme simply exploits the identical level of quantum mechanical tunneling that is available in the enzyme-free, solution-based H(*) abstraction reaction. The results also require a similar, if not identical, barrier width and height within experimental error for the H(*) abstraction both within, and outside of, the enzyme.

The cobY gene of the archaeon Halobacterium sp. strain NRC-1 is required for de novo cobamide synthesis

Genetic and nutritional analyses of mutants of the extremely halophilic archaeon Halobacterium sp. strain NRC-1 showed that open reading frame (ORF) Vng1581C encodes a protein with nucleoside triphosphate:adenosylcobinamide-phosphate nucleotidyltransferase enzyme activity. This activity was previously associated with the cobY gene of the methanogenic archaeon Methanobacterium thermoautotrophicum strain DeltaH, but no evidence was obtained to demonstrate the direct involvement of this protein in cobamide biosynthesis in archaea. Computer analysis of the Halobacterium sp. strain NRC-1 ORF Vng1581C gene and the cobY gene of M. thermoautotrophicum strain DeltaH showed the primary amino acid sequence of the proteins encoded by these two genes to be 35% identical and 48% similar. A strain of Halobacterium sp. strain NRC-1 carrying a null allele of the cobY gene was auxotrophic for cobinamide-GDP, a known intermediate of the late steps of cobamide biosynthesis. The auxotrophic requirement for cobinamide-GDP was corrected when a wild-type allele of cobY was introduced into the mutant strain, demonstrating that the lack of cobY function was solely responsible for the observed block in cobamide biosynthesis in this archaeon. The data also show that Halobacterium sp. strain NRC-1 possesses a high-affinity transport system for corrinoids and that this archaeon can synthesize cobamides de novo under aerobic growth conditions. To the best of our knowledge this is the first genetic and nutritional analysis of cobalamin biosynthetic mutants in archaea.

Isolation of the gene for the B12-dependent ribonucleotide reductase from Anabaena sp. strain PCC 7120 and expression in Escherichia coli

The gene for ribonucleotide reductase from Anabaena sp. strain PCC 7120 was identified and expressed in Escherichia coli. This gene codes for a 1,172-amino-acid protein that contains a 407-amino-acid intein. The intein splices itself from the protein when it is expressed in E. coli, yielding an active ribonucleotide reductase of 765 residues. The mature enzyme was purified to homogeneity from E. coli extracts. Anabaena ribonucleotide reductase is a monomer with a molecular weight of approximately 88,000, as determined by sodium dodecyl sulfate-polyacrylamide gel electrophoresis and Superose 12 column chromatography. The enzyme reduces ribonucleotides at the triphosphate level and requires a divalent cation and a deoxyribonucleoside triphosphate effector. The enzyme is absolutely dependent on the addition of the cofactor, 5'-adenosylcobalamin. These properties are characteristic of the class II-type reductases. The cyanobacterial enzyme has limited sequence homology to other class II reductases; the greatest similarity (38%) is to the reductase from Lactobacillus leichmannii. In contrast, the Anabaena reductase shows over 90% sequence similarity to putative reductases found in genome sequences of other cyanobacteria, such as Nostoc punctiforme, Synechococcus sp. strain WH8102, and Prochlorococcus marinus MED4, suggesting that the cyanobacterial reductases form a closely related subset of the class II enzymes.

Molecular evolution before the origin of species

Amino acids at conserved sites in the residue sequence of 10 ancient proteins, from 844 phylogenetically diverse sources, were used to specify their time of origin in the interval before species divergence from the last common ancestor (LCA). The order of amino acid addition to the genetic code, based on biosynthesis path length and other molecular evidence, provided a reference for evaluating the 'code age' of each residue profile examined. Significantly earlier estimates were obtained for conserved amino acid residues in these proteins than non-conserved residues. Evidence from the primary structure of 'fossil' proteins thus corroborated the biosynthetic order of amino acid addition to the code.Low potential ferredoxin (Fdxn) had the earliest residue profile among the proteins in this study. A phylogenetic tree for 82 prokaryote Fdxn sequences was rooted midway between bacteria and archaea branches. LCA Fdxn had a 23-residue antecedent whose residue profile matched mid-expansion phase codon assignments and included an amide residue. It contained a highly acidic N-terminal region and a non-charged C-terminal region, with all four cysteine residues. This small protein apparently anchored a [4Fe-4S] cluster, ligated by C-terminal cysteines, to a positively charged mineral surface, consistent with mediating e(-) transfer in a primordial surface system before cells appeared. Its negatively charged N-terminal 'attachment site' was highly mutable during evolution of ancestral Fdxn for Bacteria and Archaea, consistent with a loss of function after cell formation. An initial glutamate to lysine substitution may link 'attachment site' removal to early post-expansion phase entry of basic amino acids to the code. As proteins evidently anchored non-charged amide residues initially, surface attachment of cofactors and other functional groups emerges as a general function of pre-cell proteins.A phylogenetic tree of 107 proteolipid (PL) helix-1 sequences from H(+)-ATPase of bacteria, archaea and eukaryotes had its root between prokaryote branches. LCA PL h1 residue profile optimally fit a late expansion phase codon array. Sequence repeats in transmembrane PL helices h1 and h2 indicated formation of the archetypal PL hairpin structure involved successive tandem duplications, initiated within the gene for an 11-residue (or 4-residue) hydrophobic peptide. Ancestral PL h1 lacked acidic residues, in a fundamental departure from the prototype pre-cell protein. By this stage, proteins with a hydrophobic domain had evolved. Its non-polar, late expansion phase residue profile point to ancestral PL being a component of an early permeable cell membrane. Other indicators of cell formation about this stage of code evolution include phospholipid biosynthesis path length, FtsZ residue profile, and late entry of basic amino acids into the genetic code. Estimates based on conserved residues in prokaryote cell septation protein, FtsZ, and proteins involved with synthesis, transcription and replication of DNA revealed FtsZ, ribonucleotide reductase, RNA polymerase core subunits and 5'-->3' flap exonuclease, FEN-1, originated soon after cells putatively evolved. While reverse transcriptase and topoisomerase I, Topo I, appeared late in the pre-divergence era, when the genetic code was essentially complete. The transition from RNA genes to a DNA genome seemingly proceeded via formation of a DNA-RNA heteroduplex. These results suggest formation of DNA awaited evolution of a catalyst with a hydrophobic domain, capable of sequestering radical bearing intermediates in its synthesis from ribonucleotide precursors. Late formation of topology altering protein, Topo I, further suggests consolidation of genes into chromosomes followed synthesis of comparatively thermostable DNA strands.

Divergent evolution of enzymatic function: mechanistically diverse superfamilies and functionally distinct suprafamilies

The protein sequence and structure databases are now sufficiently representative that strategies nature uses to evolve new catalytic functions can be identified. Groups of divergently related enzymes whose members catalyze different reactions but share a common partial reaction, intermediate, or transition state (mechanistically diverse superfamilies) have been discovered, including the enolase, amidohydrolase, thiyl radical, crotonase, vicinal-oxygen-chelate, and Fe-dependent oxidase superfamilies. Other groups of divergently related enzymes whose members catalyze different overall reactions that do not share a common mechanistic strategy (functionally distinct suprafamilies) have also been identified: (a) functionally distinct suprafamilies whose members catalyze successive transformations in the tryptophan and histidine biosynthetic pathways and (b) functionally distinct suprafamilies whose members catalyze different reactions in different metabolic pathways. An understanding of the structural bases for the catalytic diversity observed in super- and suprafamilies may provide the basis for discovering the functions of proteins and enzymes in new genomes as well as provide guidance for in vitro evolution/engineering of new enzymes.

Streptomyces spp. contain class Ia and class II ribonucleotide reductases: expression analysis of the genes in vegetative growth

Genes encoding two ribonucleotide reductases (RNRs) were identified in members of the genus Streptomyces. One gene, nrdJ, encoded an oligomeric protein comprising four identical subunits each with a molecular mass of approximately 108 kDa. The activity of this protein depended on the presence of 5'-deoxyadenosylcobalamine (coenzyme B12), establishing it as a class II RNR. The Streptomyces clavuligerus nrdJ gene was cloned, using internal peptide sequences from the purified protein, and was found to encode a polypeptide of 961 aa. Molecular phylogenetic analysis showed that the S. clavuligerus class II RNR shares significant similarity with most other bacterial and archaeal class II RNRs. Two other genes, nrdA and nrdB, were initially identified in the Streptomyces coelicolor genome database in unannotated ORFs as encoding a class Ia RNR. Southern analysis demonstrated that the nrdAB genes were present in different Streptomyces spp. The S. coelicolor nrdAB genes were cloned and expressed in Escherichia coli, and the recombinant proteins were shown to represent a class I RNR. It was shown, using quantitative real-time PCR, that the S. clavuligerus class Ia and class II RNR genes were differentially transcribed during vegetative growth. The copy number of the class II nrdJ transcripts was approximately constant throughout the exponential phase of vegetative growth (3-5x10(5) copies per 400 ng total RNA after reverse transcription). In contrast, the copy number of the class Ia nrdAB transcripts was some 10- to 20-fold less than that of nrdJ in the early-exponential growth phase (2.8x10(4) copies), and decreased markedly at the mid-exponential (4x10(3) copies) and late-exponential phases (1.1x10(3) copies) of growth. A possible role for the involvement of two RNRs during vegetative growth is discussed.

Structure and function of the radical enzyme ribonucleotide reductase

Ribonucleotide reductases (RNRs) catalyze all new production in nature of deoxyribonucleotides for DNA synthesis by reducing the corresponding ribonucleotides. The reaction involves the action of a radical that is produced differently for different classes of the enzyme. Class I enzymes, which are present in eukaryotes and microorganisms, use an iron center to produce a stable tyrosyl radical that is stored in one of the subunits of the enzyme. The other classes are only present in microorganisms. Class II enzymes use cobalamin for radical generation and class III enzymes, which are found only in anaerobic organisms, use a glycyl radical. The reductase activity is in all three classes contained in enzyme subunits that have similar structures containing active site cysteines. The initiation of the reaction by removal of the 3'-hydrogen of the ribose by a transient cysteinyl radical is a common feature of the different classes of RNR. This cysteine is in all RNRs located on the tip of a finger loop inserted into the center of a special barrel structure. A wealth of structural and functional information on the class I and class III enzymes can now give detailed views on how these enzymes perform their task. The class I enzymes demonstrate a sophisticated pattern as to how the free radical is used in the reaction, in that it is only delivered to the active site at exactly the right moment. RNRs are also allosterically regulated, for which the structural molecular background is now starting to be revealed.

The evolution of ribonucleotide reduction revisited

Ribonucleotide reductases (RNRs) catalyze the conversion of both purine and pyrimidine nucleotides to deoxynucleotides in all organisms and provide all the monomeric precursors essential for both DNA replication and repair. RNRs have been divided into three classes on the basis of their unique metallo-cofactors. The exquisitely controlled free radical chemistry used by all RNRs, and the commonality of the structures of the subunits where the nucleotide reduction process occurs, together provide compelling evidence for the importance of chemistry in the divergent evolution of RNRs from a common progenitor.

Methyl-RNA: an evolutionary bridge between RNA and DNA?

Given the apparent limitation of double-stranded RNA (dsRNA) genomes to about 30 kb, together with the complexity of DNA synthesis, it appears difficult for a dsRNA genome to encode all the information required before the transition from an RNA to a DNA genome. Ribonucleotide reductase itself, which synthesizes deoxyribonucleotides from ribonucleotides, requires complex protein radical chemistry, and RNA world genomes may have reached their limits of coding capacity well before such complex enzymes had evolved. The transition from RNA to DNA thus appears to require intermediate steps, and we suggest that the naturally occurring 2'-O-methylated RNA, with chemical properties intermediate between RNA and DNA, is a suitable candidate.

Functional inferences from reconstructed evolutionary biology involving rectified databases--an evolutionarily grounded approach to functional genomics

If bioinformatics tools are constructed to reproduce the natural, evolutionary history of the biosphere, they offer powerful approaches to some of the most difficult tasks in genomics, including the organization and retrieval of sequence data, the updating of massive genomic databases, the detection of database error, the assignment of introns, the prediction of protein conformation from protein sequences, the detection of distant homologs, the assignment of function to open reading frames, the identification of biochemical pathways from genomic data, and the construction of a comprehensive model correlating the history of biomolecules with the history of planet Earth.

Hydrogen abstraction from thiols by adenosyl radicals: chemical precedent for thiyl radical formation, the first catalytic step in ribonucleoside triphosphate reductase from Lactobacillus leichmannii

Aqueous solutions of adenosylcobalamin (AdoCbl) were thermolyzed with excess beta-mercaptoethanol under anaerobic conditions. The product studies reveal that approximately 90% Co-C bond homolysis occurs, to yield Co(II)cobalamin, 5'-deoxyadenosine, and the disulfide product from the combination of two HOCH2CH2S* radicals, 2,2'-dithiodiethanol; there is also approximately 10% Co-C bond heterolysis, yielding Co(III)cobalamin, adenine, and 2,3-dihydroxy-4-pentenal. The kinetic studies show there is a first-order dependence on AdoCbl and zero-order dependence on thiol under the higher [RSH] experimental conditions used, consistent with the rate-determining step at high [RSH] being the generation of Ado*. The kinetic results require that, in enzyme-free AdoCbl solution, adenosyl radical (Ado*) is formed as a discrete intermediate which then abstracts H* from the added thiol. The activation parameters for Co-C bond homolysis in the presence of thiol trap are the same within experimental error as the activation parameters for Co-C bond homolysis without trap, standard delta H(obs) = 29(2) kcal mol(-1) and standard delta S(obs) = -1(5) e.u. The results, in comparison to the rate of Co-C bond homolysis in ribonucleoside triphosphate reductase (RTPR), reveal that RTPR accelerates Co-C bond cleavage in AdoCbl by approximately 10(10+/-1). The recent literature evidence bearing on the exact mechanism of RTPR enzymic cleavage of the Co-C bond of AdoCbl is briefly discussed, notably the fact that this mechanism is presently controversial, but does involve at least coupled (and possibly concerted) Co-C cleavage, -S-H cleavage, and C-H (Ado-H) formation steps.

Glycyl radical enzymes: a conservative structural basis for radicals

The first structures of glycyl radical enzymes, the anaerobic ribonucleotide reductase from bacteriophage T4 and pyruvate formate lyase from Escherichia coli, have been recently determined. This work provides new insights into the structure and chemistry of glycyl radical sites.

Did DNA replication evolve twice independently?

DNA replication is central to all extant cellular organisms. There are substantial functional similarities between the bacterial and the archaeal/eukaryotic replication machineries, including but not limited to defined origins, replication bidirectionality, RNA primers and leading and lagging strand synthesis. However, several core components of the bacterial replication machinery are unrelated or only distantly related to the functionally equivalent components of the archaeal/eukaryotic replication apparatus. This is in sharp contrast to the principal proteins involved in transcription and translation, which are highly conserved in all divisions of life. We performed detailed sequence comparisons of the proteins that fulfill indispensable functions in DNA replication and classified them into four main categories with respect to the conservation in bacteria and archaea/eukaryotes: (i) non-homologous, such as replicative polymerases and primases; (ii) containing homologous domains but apparently non-orthologous and conceivably independently recruited to function in replication, such as the principal replicative helicases or proofreading exonucleases; (iii) apparently orthologous but poorly conserved, such as the sliding clamp proteins or DNA ligases; (iv) orthologous and highly conserved, such as clamp-loader ATPases or 5'-->3' exonucleases (FLAP nucleases). The universal conservation of some components of the DNA replication machinery and enzymes for DNA precursor biosynthesis but not the principal DNA polymerases suggests that the last common ancestor (LCA) of all modern cellular life forms possessed DNA but did not replicate it the way extant cells do. We propose that the LCA had a genetic system that contained both RNA and DNA, with the latter being produced by reverse transcription. Consequently, the modern-type system for double-stranded DNA replication likely evolved independently in the bacterial and archaeal/eukaryotic lineages.

Adenosylcobalamin-dependent ribonucleoside triphosphate reductase from Lactobacillus leichmannii. Rapid, improved purification involving dGTP-based affinity chromatography plus biophysical characterization studies demonstrating enhanced, "crystallographic level" purity

Ribonucleoside triphosphate reductase (RTPR, EC 1.17.4.2) from Lactobacillus leichmannii is a 5'-deoxyadenosylcobalamin-dependent (AdoCbl; Coenzyme B12) enzyme. RTPR is also a prototypical adenosylcobalamin-dependent ribonucleotide reductase, one that, as its name indicates, converts ribonucleoside triphosphates (NTP) to deoxyribonucleoside triphosphates (dNTP). Upon substrate binding to RTPR, AdoCbl's cobalt-carbon bond is cleaved to generate cob(II)alamin, 5'-deoxyadenosine, and the cysteine (C408) derived thiyl radical. Five key cysteines (Cys 119, 408, 419, 731, and 736), from among the ten total cysteines, are involved in RTPR's catalytic mechanism. A critical examination of the RTPR isolation and purification literature suggested that the purification protocol currently used results in RTPR which contains 2040% microheterogeneity, along with minor contamination by other proteins. In addition, no report of crystalline RTPR has ever appeared. The literature indicates that irreversible cysteine oxidation (e.g., to -SO2H or -SO3H) is one highly plausible reason for the microheterogeneity of RTPR. The literature also indicates that improvement in the level of enzyme purity is the most effective next step in coaxing enzymes to crystallize that have previously failed to do so. A shortened, improved purification of RTPR has been developed, one involving a shorter purification time, a lower pH, a higher concentration of the more effective reductant DTT (all designed to help protect the cysteines from oxidation), and a final step utilizing our recently reported, improved dGTP-based affinity chromatography resin. The resultant RTPR is approximately 20-30% higher in both specific activity and in its ability to undergo single turnovers, and is homogeneous by mass spectrometry and dynamic light scattering. Additionally, the revised purification procedure eliminates > 30 proteins present in 2-3% amounts along with damaged RTPR that does not bind properly (i.e. tightly) to the dGTP-affinity resin. Finally, dGTP-based affinity chromatography purified RTPR has yielded the first reported, albeit small, single crystals of RTPR.

Pyruvate formate lyase is structurally homologous to type I ribonucleotide reductase

BACKGROUND

Pyruvate formate lyase (PFL) catalyses a key step in Escherichia coli anaerobic glycolysis by converting pyruvate and CoA to formate and acetylCoA. The PFL mechanism involves an unusual radical cleavage of pyruvate, involving an essential C alpha radical of Gly734 and two cysteine residues, Cys418 and Cys419, which may form thiyl radicals required for catalysis. We undertook this study to understand the structural basis for catalysis.

RESULTS

The first structure of a fragment of PFL (residues 1-624) at 2.8 A resolution shows an unusual barrel-like structure, with a catalytic beta finger carrying Cys418 and Cys419 inserted into the centre of the barrel. Several residues near the active-site cysteines can be ascribed roles in the catalytic mechanism: Arg176 and Arg435 are positioned near Cys419 and may bind pyruvate/formate and Trp333 partially buries Cys418. Both cysteine residues are accessible to each other owing to their cis relationship at the tip of the beta finger. Finally, two clefts that may serve as binding sites for CoA and pyruvate have been identified.

CONCLUSIONS

PFL has striking structural homology to the aerobic ribonucleotide reductase (RNR): the superposition of PFL and RNR includes eight of the ten strands in the unusual RNR alpha/beta barrel as well as the beta finger, which carries key catalytic residues in both enzymes. This provides the first structural proof that RNRs and PFLs are related by divergent evolution from a common ancestor.

Allosteric control of three B12-dependent (class II) ribonucleotide reductases. Implications for the evolution of ribonucleotide reduction

Three separate classes of ribonucleotide reductases are known, each with a distinct protein structure. One common feature of all enzymes is that a single protein generates each of the four deoxyribonucleotides. Class I and III enzymes contain an allosteric substrate specificity site capable of binding effectors (ATP or various deoxyribonucleoside triphosphates) that direct enzyme specificity. Some (but not all) enzymes contain a second allosteric site that binds only ATP or dATP. Binding of dATP to this site inhibits the activity of these enzymes. X-ray crystallography has localized the two sites within the structure of the Escherichia coli class I enzyme and identified effector-binding amino acids. Here, we have studied the regulation of three class II enzymes, one from the archaebacterium Thermoplasma acidophilum and two from eubacteria (Lactobacillus leichmannii and Thermotoga maritima). Each enzyme has an allosteric site that binds ATP or various deoxyribonucleoside triphosphates and that regulates its substrate specificity according to the same rules as for class I and III enzymes. dATP does not inhibit enzyme activity, suggesting the absence of a second active allosteric site. For the L. leichmannii and T. maritima enzymes, binding experiments also indicate the presence of only one allosteric site. Their primary sequences suggest that these enzymes lack the structural requirements for a second site. In contrast, the T. acidophilum enzyme binds dATP at two separate sites, and its sequence contains putative effector-binding amino acids for a second site. The presence of a second site without apparent physiological function leads to the hypothesis that a functional site was present early during the evolution of ribonucleotide reductases, but that its function was lost from the T. acidophilum enzyme. The other two B12 enzymes lost not only the function, but also the structural basis for the site. Also a large subgroup (Ib) of class I enzymes, but none of the investigated class III enzymes, has lost this site. This is further indirect evidence that class II and I enzymes may have arisen by divergent evolution from class III enzymes.

Evolutionary, mechanistic, and predictive analyses of the hydroxymethyldihydropterin pyrophosphokinase family of proteins

A prediction has been prepared ab initio for the secondary structure of the hydroxymethyldihydropterin pyrophosphokinase (HPPK) family of proteins starting from a set of aligned homologous protein sequences. Attempts to identify a fold by threading failed, judging by the inability to find a threading "hit" that had a secondary structure that was plausibly congruent to the predicted secondary structure for the HPPK family. Therefore, a set of tertiary structure models was assembled ab initio, where alternative models were built and used to select between alternative secondary structure models. This prediction report illustrates the importance of non-computational approaches to structure prediction at its present frontier, which is to obtain medium resolution models of tertiary structure.

Ribonucleotide reductases

Ribonucleotide reductases provide the building blocks for DNA replication in all living cells. Three different classes of enzymes use protein free radicals to activate the substrate. Aerobic class I enzymes generate a tyrosyl radical with an iron-oxygen center and dioxygen, class II enzymes employ adenosylcobalamin, and the anaerobic class III enzymes generate a glycyl radical from S-adenosylmethionine and an iron-sulfur cluster. The X-ray structure of the class I Escherichia coli enzyme, including forms that bind substrate and allosteric effectors, confirms previous models of catalytic and allosteric mechanisms. This structure suggests considerable mobility of the protein during catalysis and, together with experiments involving site-directed mutants, suggests a mechanism for radical transfer from one subunit to the other. Despite large differences between the classes, common catalytic and allosteric mechanisms, as well as retention of critical residues in the protein sequence, suggest a similar tertiary structure and a common origin during evolution. One puzzling aspect is that some organisms contain the genes for several different reductases.

Functional analysis of the Escherichia coli genome using the sequence-to-structure-to-function paradigm: identification of proteins exhibiting the glutaredoxin/thioredoxin disulfide oxidoreductase activity

The application of an automated method for the screening of protein activity based on the sequence-to-structure-to-function paradigm is presented for the complete Escherichia coli genome. First, the structure of the protein is identified from its sequence using a threading algorithm, which aligns the sequences to the best matching structure in a structural database and extends sequence analysis well beyond the limits of local sequence identity. Then, the active site is identified in the resulting sequence-to-structure alignment using a "fuzzy functional form" (FFF), a three-dimensional descriptor of the active site of a protein. Here, this sequence-to-structure-to-function concept is applied to analysis of the complete E. coli genome, i.e. all E. coli open reading frames (ORFs) are screened for the thiol-disulfide oxidoreductase activity of the glutaredoxin/thioredoxin protein family. We show that the method can identify the active sites in ten sequences that are known to or proposed to exhibit this activity. Furthermore, oxidoreductase activity is predicted in two other sequences that have not been identified previously. This method distinguishes protein pairs with similar active sites from proteins pairs that are just topological cousins, i.e. those having similar global folds, but not necessarily similar active sites. Thus, this method provides a novel approach for extraction of active site and functional information based on three-dimensional structures, rather than simple sequence analysis. Prediction of protein activity is fully automated and easily extendible to new functions. Finally, it is demonstrated here that the method can be applied to complete genome database analysis.

Mechanistically diverse enzyme superfamilies: the importance of chemistry in the evolution of catalysis

The strategy that nature has used to evolve new catalytic activities from pre-existing enzymes (i.e. retention of substrate binding or of catalytic mechanism) has been controversial. Recent work supports a strategy in which a partial reaction, catalyzed by a progenitor, is retained, and the active-site architecture is modified to allow the intermediate generated to be directed to different products.

Method for prediction of protein function from sequence using the sequence-to-structure-to-function paradigm with application to glutaredoxins/thioredoxins and T1 ribonucleases

The practical exploitation of the vast numbers of sequences in the genome sequence databases is crucially dependent on the ability to identify the function of each sequence. Unfortunately, current methods, including global sequence alignment and local sequence motif identification, are limited by the extent of sequence similarity between sequences of unknown and known function; these methods increasingly fail as the sequence identity diverges into and beyond the twilight zone of sequence identity. To address this problem, a novel method for identification of protein function based directly on the sequence-to-structure-to-function paradigm is described. Descriptors of protein active sites, termed "fuzzy functional forms" or FFFs, are created based on the geometry and conformation of the active site. By way of illustration, the active sites responsible for the disulfide oxidoreductase activity of the glutaredoxin/thioredoxin family and the RNA hydrolytic activity of the T1 ribonuclease family are presented. First, the FFFs are shown to correctly identify their corresponding active sites in a library of exact protein models produced by crystallography or NMR spectroscopy, most of which lack the specified activity. Next, these FFFs are used to screen for active sites in low-to-moderate resolution models produced by ab initio folding or threading prediction algorithms. Again, the FFFs can specifically identify the functional sites of these proteins from their predicted structures. The results demonstrate that low-to-moderate resolution models as produced by state-of-the-art tertiary structure prediction algorithms are sufficient to identify protein active sites. Prediction of a novel function for the gamma subunit of a yeast glycosyl transferase and prediction of the function of two hypothetical yeast proteins whose models were produced via threading are presented. This work suggests a means for the large-scale functional screening of genomic sequence databases based on the prediction of structure from sequence, then on the identification of functional active sites in the predicted structure.

B12-dependent ribonucleotide reductases from deeply rooted eubacteria are structurally related to the aerobic enzyme from Escherichia coli

The ribonucleotide reductases from three ancient eubacteria, the hyperthermophilic Thermotoga maritima (TM), the radioresistant Deinococcus radiodurans (DR), and the thermophilic photosynthetic Chloroflexus aurantiacus, were found to be coenzyme-B12 (class II) enzymes, similar to the earlier described reductases from the archaebacteria Thermoplasma acidophila and Pyrococcus furiosus. Reduction of CDP by the purified TM and DR enzymes requires adenosylcobalamin and DTT. dATP is a positive allosteric effector, but stimulation of the TM enzyme only occurs close to the temperature optimum of 80-90 degrees C. The TM and DR genes were cloned by PCR from peptide sequence information. The TM gene was sequenced completely and expressed in Escherichia coli. The deduced amino acid sequences of the two eubacterial enzymes are homologous to those of the archaebacteria. They can also be aligned to the sequence of the large protein of the aerobic E. coli ribonucleotide reductase that belongs to a different class (class I), which is not dependent on B12. Structure determinations of the E. coli reductase complexed with substrate and allosteric effectors earlier demonstrated a 10-stranded beta/alpha-barrel in the active site. From the conservation of substrate- and effector-binding residues we propose that the B12-dependent class II enzymes contain a similar barrel.

Archaea: what can we learn from their sequences?

Gene-by-gene and traditional biochemical approaches continue to reveal surprising molecular features in the archaeal domain. In addition, the complete sequencing of several archaeal genomes has further confirmed the phenotypic coherence of these micro-organisms at the molecular level. Nevertheless, the phylogeny of Archaea and the nature of the last universal common ancestor are still matters for debate.

Complete genome sequence of Methanobacterium thermoautotrophicum deltaH: functional analysis and comparative genomics

The complete 1,751,377-bp sequence of the genome of the thermophilic archaeon Methanobacterium thermoautotrophicum deltaH has been determined by a whole-genome shotgun sequencing approach. A total of 1,855 open reading frames (ORFs) have been identified that appear to encode polypeptides, 844 (46%) of which have been assigned putative functions based on their similarities to database sequences with assigned functions. A total of 514 (28%) of the ORF-encoded polypeptides are related to sequences with unknown functions, and 496 (27%) have little or no homology to sequences in public databases. Comparisons with Eucarya-, Bacteria-, and Archaea-specific databases reveal that 1,013 of the putative gene products (54%) are most similar to polypeptide sequences described previously for other organisms in the domain Archaea. Comparisons with the Methanococcus jannaschii genome data underline the extensive divergence that has occurred between these two methanogens; only 352 (19%) of M. thermoautotrophicum ORFs encode sequences that are >50% identical to M. jannaschii polypeptides, and there is little conservation in the relative locations of orthologous genes. When the M. thermoautotrophicum ORFs are compared to sequences from only the eucaryal and bacterial domains, 786 (42%) are more similar to bacterial sequences and 241 (13%) are more similar to eucaryal sequences. The bacterial domain-like gene products include the majority of those predicted to be involved in cofactor and small molecule biosyntheses, intermediary metabolism, transport, nitrogen fixation, regulatory functions, and interactions with the environment. Most proteins predicted to be involved in DNA metabolism, transcription, and translation are more similar to eucaryal sequences. Gene structure and organization have features that are typical of the Bacteria, including genes that encode polypeptides closely related to eucaryal proteins. There are 24 polypeptides that could form two-component sensor kinase-response regulator systems and homologs of the bacterial Hsp70-response proteins DnaK and DnaJ, which are notably absent in M. jannaschii. DNA replication initiation and chromosome packaging in M. thermoautotrophicum are predicted to have eucaryal features, based on the presence of two Cdc6 homologs and three histones; however, the presence of an ftsZ gene indicates a bacterial type of cell division initiation. The DNA polymerases include an X-family repair type and an unusual archaeal B type formed by two separate polypeptides. The DNA-dependent RNA polymerase (RNAP) subunits A', A", B', B" and H are encoded in a typical archaeal RNAP operon, although a second A' subunit-encoding gene is present at a remote location. There are two rRNA operons, and 39 tRNA genes are dispersed around the genome, although most of these occur in clusters. Three of the tRNA genes have introns, including the tRNAPro (GGG) gene, which contains a second intron at an unprecedented location. There is no selenocysteinyl-tRNA gene nor evidence for classically organized IS elements, prophages, or plasmids. The genome contains one intein and two extended repeats (3.6 and 8.6 kb) that are members of a family with 18 representatives in the M. jannaschii genome.

Ribonucleotide reductase in the archaeon Pyrococcus furiosus: a critical enzyme in the evolution of DNA genomes?

Ribonucleotide reductase (RNR), the enzyme responsible for deoxyribonucleotide synthesis, has been isolated from Pyrococcus furiosus, a deeply branching hyperthermophilic, strictly anaerobic archaeon. Its gene has been cloned, sequenced, and shown to harbor two insertions encoding inteins. The purified enzyme absolutely requires adenosylcobalamin for activity, a trait that defines it as a member of class II (adenosyl-cobalamin-dependent) prokaryotic RNRs. On the other hand, the archaeal RNR has significant amino acid sequence homology with class I (aerobic non-heme iron-dependent) and class III (anaerobic iron-sulfur-dependent) RNRs present in eukaryotes and bacteria, respectively. It is proposed that this enzyme may be the closest possible relative of the original RNR, which allowed the key "RNA world" to "DNA world" transition, and that the different classes of present-day RNRs are the products of divergent evolution.