Strand specificity of ribonucleotide excision repair in Escherichia coli

In Escherichia coli, replication of both strands of genomic DNA is carried out by a single replicase-DNA polymerase III holoenzyme (pol III HE). However, in certain genetic backgrounds, the low-fidelity TLS polymerase, DNA polymerase V (pol V) gains access to undamaged genomic DNA where it promotes elevated levels of spontaneous mutagenesis preferentially on the lagging strand. We employed active site mutants of pol III (pol IIIα_S759N) and pol V (pol V_Y11A) to analyze ribonucleotide incorporation and removal from the E. coli chromosome on a genome-wide scale under conditions of normal replication, as well as SOS induction. Using a variety of methods tuned to the specific properties of these polymerases (analysis of lacI mutational spectra, lacZ reversion assay, HydEn-seq, alkaline gel electrophoresis), we present evidence that repair of ribonucleotides from both DNA strands in E. coli is unequal. While RNase HII plays a primary role in leading-strand Ribonucleotide Excision Repair (RER), the lagging strand is subject to other repair systems (RNase HI and under conditions of SOS activation also Nucleotide Excision Repair). Importantly, we suggest that RNase HI activity can also influence the repair of single ribonucleotides incorporated by the replicase pol III HE into the lagging strand.

5'-End Mapping in Human Mitochondrial DNA

Here, we describe an assay that enables mapping of 5'-ends across the genome using next-generation sequencing on an Illumina platform, 5'-End-sequencing (5'-End-seq). We use this method to map free 5'-ends in mtDNA isolated from fibroblasts. This method can be used to answer key questions regarding DNA integrity, DNA replication mechanisms and to identify priming events, primer processing, nick processing, and double strand break processing on the entire genome.

Tracking Escherichia coli DNA polymerase V to the entire genome during the SOS response

Ribonucleotides are frequently incorporated into DNA and can be used as a marker of DNA replication enzymology. To investigate on a genome-wide scale, how E. coli pol V accesses undamaged chromosomal DNA during the SOS response, we mapped the location of ribonucleotides incorporated by steric gate variants of pol V across the entire E. coli genome. To do so, we used strains that are deficient in ribonucleotide excision repair (ΔrnhB), deficient in pol IV DNA polymerase, constitutively express all SOS-regulated genes [lexA(Def)] and constitutively “activated” RecA* (recA730). The strains also harbor two steric gate variants of E. coli pol V (Y11A, or F10L), or a homolog of pol V, (pol V_R391-Y13A). Ribonucleotides are frequently incorporated by the pol V-Y11A and pol V_R391-Y13A variants, with a preference to the lagging strand. In contrast, the pol V-F10L variant incorporates less ribonucleotides and no strand preference is observed. Sharp transitions in strand specificity are observed at the replication origin (oriC), while a gradient is observed at the termination region. To activate RecA* in a recA⁺ strain, we treated the strains with ciprofloxacin and genome-wide mapped the location of the incorporated ribonucleotides. Again, the pol V-Y11A steric gate variant exhibited a lagging strand preference. Our data are consistent with a specific role for pol V in lagging strand DNA synthesis across the entire E. coli genome during the SOS response.

DNA polymerase η contributes to genome-wide lagging strand synthesis

DNA polymerase η (pol η) is best known for its ability to bypass UV-induced thymine-thymine (T-T) dimers and other bulky DNA lesions, but pol η also has other cellular roles. Here, we present evidence that pol η competes with DNA polymerases α and δ for the synthesis of the lagging strand genome-wide, where it also shows a preference for T-T in the DNA template. Moreover, we found that the C-terminus of pol η, which contains a PCNA-Interacting Protein motif is required for pol η to function in lagging strand synthesis. Finally, we provide evidence that a pol η dependent signature is also found to be lagging strand specific in patients with skin cancer. Taken together, these findings provide insight into the physiological role of DNA synthesis by pol η and have implications for our understanding of how our genome is replicated to avoid mutagenesis, genome instability and cancer.

RNase H1 directs origin-specific initiation of DNA replication in human mitochondria

Human mitochondrial DNA (mtDNA) replication is first initiated at the origin of H-strand replication. The initiation depends on RNA primers generated by transcription from an upstream promoter (LSP). Here we reconstitute this process in vitro using purified transcription and replication factors. The majority of all transcription events from LSP are prematurely terminated after ~120 nucleotides, forming stable R-loops. These nascent R-loops cannot directly prime mtDNA synthesis, but must first be processed by RNase H1 to generate 3′-ends that can be used by DNA polymerase γ to initiate DNA synthesis. Our findings are consistent with recent studies of a knockout mouse model, which demonstrated that RNase H1 is required for R-loop processing and mtDNA maintenance in vivo. Both R-loop formation and DNA replication initiation are stimulated by the mitochondrial single-stranded DNA binding protein. In an RNase H1 deficient patient cell line, the precise initiation of mtDNA replication is lost and DNA synthesis is initiated from multiple sites throughout the mitochondrial control region. In combination with previously published in vivo data, the findings presented here suggest a model, in which R-loop processing by RNase H1 directs origin-specific initiation of DNA replication in human mitochondria.

Human mitochondria contain a double-stranded DNA genome that codes for key components of the oxidative phosphorylation system. The mitochondrial DNA (mtDNA) is replicated by a replication machinery distinct from that operating in the nucleus and mutations affecting individual replication factors have been associated with an array of rare, human diseases. In the present work, we demonstrate that RNase H1 directs origin-specific initiation of DNA replication in human mitochondria and that disease-causing mutations may impair this process. A unique feature of mtDNA replication is that primers required for initiation of leading-strand DNA replication are produced by the mitochondrial transcription machinery. A substantial fraction of all transcription events is prematurely terminated about 120 nucleotides downstream of the promoter and the RNA remains firmly associated with the genome, forming R-loops. Interestingly, the free 3′-end of these R-loops cannot directly prime initiation of DNA synthesis, but must first be processed by RNase H1. The process is stimulated by the mitochondrial single-stranded DNA binding protein and faithfully reconstitutes replication events mapped in vivo. In combination with mapping of replication events in fibroblasts derived from patients with mutations in RNASEH1, our findings point to a possible model for replication initiation in human mitochondria similar to that previously described in the E. coli plasmid, ColE1.

Elevated pyrimidine dimer formation at distinct genomic bases underlies promoter mutation hotspots in UV-exposed cancers

Sequencing of whole cancer genomes has revealed an abundance of recurrent mutations in gene-regulatory promoter regions, in particular in melanoma where strong mutation hotspots are observed adjacent to ETS-family transcription factor (TF) binding sites. While sometimes interpreted as functional driver events, these mutations are commonly believed to be due to locally inhibited DNA repair. Here, we first show that low-dose UV light induces mutations preferably at a known ETS promoter hotspot in cultured cells even in the absence of global or transcription-coupled nucleotide excision repair (NER). Further, by genome-wide mapping of cyclobutane pyrimidine dimers (CPDs) shortly after UV exposure and thus before DNA repair, we find that ETS-related mutation hotspots exhibit strong increases in CPD formation efficacy in a manner consistent with tumor mutation data at the single-base level. Analysis of a large whole genome cohort illustrates the widespread contribution of this effect to recurrent mutations in melanoma. While inhibited NER underlies a general increase in somatic mutation burden in regulatory elements including ETS sites, our data supports that elevated DNA damage formation at specific genomic bases is at the core of the prominent promoter mutation hotspots seen in skin cancers, thus explaining a key phenomenon in whole-genome cancer analyses.

Mapping Ribonucleotides Incorporated into DNA by Hydrolytic End-Sequencing

Ribonucleotides embedded within DNA render the DNA sensitive to the formation of single-stranded breaks under alkali conditions. Here, we describe a next-generation sequencing method called hydrolytic end sequencing (HydEn-seq) to map ribonucleotides inserted into the genome of Saccharomyce cerevisiae strains deficient in ribonucleotide excision repair. We use this method to map several genomic features in wild-type and replicase variant yeast strains.

Simultaneous Mapping and Quantitation of Ribonucleotides in Human Mitochondrial DNA

Established approaches to estimate the number of ribonucleotides present in a genome are limited to the quantitation of incorporated ribonucleotides using short synthetic DNA fragments or plasmids as templates and then extrapolating the results to the whole genome. Alternatively, the number of ribonucleotides present in a genome may be estimated using alkaline gels or Southern blots. More recent in vivo approaches employ Next-generation sequencing allowing genome-wide mapping of ribonucleotides, providing the position and identity of embedded ribonucleotides. However, they do not allow quantitation of the number of ribonucleotides which are incorporated into a genome. Here we describe how to simultaneously map and quantitate the number of ribonucleotides which are incorporated into human mitochondrial DNA in vivo by Next-generation sequencing. We use highly intact DNA and introduce sequence specific double strand breaks by digesting it with an endonuclease, subsequently hydrolyzing incorporated ribonucleotides with alkali. The generated ends are ligated with adapters and these ends are sequenced on a Next-generation sequencing machine. The absolute number of ribonucleotides can be calculated as the number of reads outside the recognition site per average number of reads at the recognition site for the sequence specific endonuclease. This protocol may also be utilized to map and quantitate free nicks in DNA and allows adaption to map other DNA lesions that can be processed to 5´-OH ends or 5´-phosphate ends. Furthermore, this method can be applied to any organism, given that a suitable reference genome is available. This protocol therefore provides an important tool to study DNA replication, 5´-end processing, DNA damage, and DNA repair.

Gene duplications and losses among vertebrate deoxyribonucleoside kinases of the non-TK1 Family

Deoxyribonucleoside kinases (dNKs) salvage deoxyribonucleosides (dNs) and catalyze the rate limiting step of this salvage pathway by converting dNs into corresponding monophosphate forms. These enzymes serve as an excellent model to study duplicated genes and their evolutionary history. So far, among vertebrates only four mammalian dNKs have been studied for their substrate specificity and kinetic properties. However, some vertebrates, such as fish, frogs, and birds, apparently possess a duplicated homolog of deoxycytidine kinase (dCK). In this study, we characterized a family of dCK/deoxyguanosine kinase (dGK)-like enzymes from a frog Xenopus laevis and a bird Gallus gallus. We showed that X. laevis has a duplicated dCK gene and a dGK gene, whereas G. gallus has a duplicated dCK gene but has lost the dGK gene. We cloned, expressed, purified, and subsequently determined the kinetic parameters of the dCK/dGK enzymes encoded by these genes. The two dCK enzymes in G. gallus have broader substrate specificity than their human or X. laevis counterparts. Additionally, the duplicated dCK enzyme in G. gallus might have become mitochondria. Based on our study we postulate that changing and adapting substrate specificities and subcellular localization are likely the drivers behind the evolution of vertebrate dNKs.

Nucleotide pools dictate the identity and frequency of ribonucleotide incorporation in mitochondrial DNA

Previous work has demonstrated the presence of ribonucleotides in human mitochondrial DNA (mtDNA) and in the present study we use a genome-wide approach to precisely map the location of these. We find that ribonucleotides are distributed evenly between the heavy- and light-strand of mtDNA. The relative levels of incorporated ribonucleotides reflect that DNA polymerase γ discriminates the four ribonucleotides differentially during DNA synthesis. The observed pattern is also dependent on the mitochondrial deoxyribonucleotide (dNTP) pools and disease-causing mutations that change these pools alter both the absolute and relative levels of incorporated ribonucleotides. Our analyses strongly suggest that DNA polymerase γ-dependent incorporation is the main source of ribonucleotides in mtDNA and argues against the existence of a mitochondrial ribonucleotide excision repair pathway in human cells. Furthermore, we clearly demonstrate that when dNTP pools are limiting, ribonucleotides serve as a source of building blocks to maintain DNA replication. Increased levels of embedded ribonucleotides in patient cells with disturbed nucleotide pools may contribute to a pathogenic mechanism that affects mtDNA stability and impair new rounds of mtDNA replication.

Human mitochondria contain a small double-stranded DNA genome (mtDNA) of only 16,569 base pairs (bp) that encodes 13 essential subunits of the oxidative phosphorylation system. Depletion of mtDNA and different types of mtDNA mutations cause mitochondrial disease, and are also implicated in biological ageing. For almost half a century it has been known that mtDNA contains ribonucleotides, but their identity and precise location are not known. The source of these ribonucleotides and their relevance for mitochondrial genome stability in healthy individuals and in patients with mitochondrial defects has not been addressed. We have used a combination of next-generation sequencing, and in vivo and in vitro biochemistry to address some of these questions. Our findings demonstrate that DNA polymerase γ-dependent incorporation is the main source of ribonucleotides in mtDNA and argues against the existence of ribonucleotide excision repair pathways in human mitochondria. Our data also reveal that when dNTP pools are limiting, ribonucleotides serves as a second line of building blocks for DNA synthesis. We also demonstrate increased levels of embedded ribonucleotides in patient cells with disturbed nucleotide pools, which may constitute a new pathogenic mechanism that affects mtDNA stability and impairs later rounds of mtDNA replication.

Measuring ribonucleotide incorporation into DNA in vitro and in vivo

Ribonucleotides are incorporated into genomes by DNA polymerases, they can be removed, and if not removed, they can have deleterious and beneficial consequences. Here, we describe an assay to quantify stable ribonucleotide incorporation by DNA polymerases in vitro, and an assay to probe for ribonucleotides in each of the two DNA strands of the yeast nuclear genome.

Plants salvage deoxyribonucleosides in mitochondria

Deoxyribonucleoside kinases phosphorylate deoxyribonucleosides into the corresponding 5'-monophosphate deoxyribonucleosides to supply the cell with nucleic acid precursors. In mitochondrial fractions of the model plant Arabidopsis thaliana, we detected deoxyadenosine and thymidine kinase activities, while the cytosol fraction contained six-fold lower activity and chloroplasts contained no measurable activities. In addition, a mitochondrial fraction isolated from the potato Solanum tuberosum contained thymidine kinase and deoxyadenosine kinase activities. We conclude that an active salvage of deoxyribonucleosides in plants takes place in their mitochondria. In general, the observed localization of the plant dNK activities in the mitochondrion suggests that plants have a different organization of the deoxyribonucleoside salvage compared to mammals.

Ribonucleotides are signals for mismatch repair of leading-strand replication errors

To maintain genome stability, mismatch repair of nuclear DNA replication errors must be directed to the nascent strand, likely by DNA ends and PCNA. Here we show that the efficiency of mismatch repair in Saccharomyces cerevisiae is reduced by inactivating RNase H2, which nicks DNA containing ribonucleotides incorporated during replication. In strains encoding mutator polymerases, this reduction is preferential for repair of mismatches made by leading-strand DNA polymerase ε as compared to lagging-strand DNA polymerase δ. The results suggest that RNase-H2-dependent processing of ribonucleotides transiently present in DNA after replication may direct mismatch repair to the continuously replicated nascent leading strand.

Thymidine kinase 1 regulatory fine-tuning through tetramer formation

Thymidine kinase 1 (TK1) provides a crucial precursor, deoxythymidine monophosphate, for nucleic acid synthesis, and the activity of TK1 increases by up to 200-fold during the S-phase of cell division in humans. An important part of the regulatory checkpoints is the ATP and enzyme concentration-dependent transition of TK1 from a dimer with low catalytic efficiency to a tetramer with high catalytic efficiency. This regulatory fine-tuning serves as an additional control to provide a balanced pool of nucleic acid precursors in the cell. We subcloned and over-expressed 10 different TK1s, originating from widely different organisms, and characterized their kinetic and oligomerization properties. Whilst bacteria, plants and Dictyostelium only exhibited dimeric TK1, we found that all animals had a tetrameric TK1. However, a clear ATP-dependent switch between dimer and tetramer was found only in higher vertebrates and was especially pronounced in mammalian and bird TK1s. We suggest that the dimer form is the original form and that the tetramer originated in the animal lineage after the split of Dictyostelium and the lineages leading to invertebrates and vertebrates. The efficient switching mechanism was probably first established in warm-blooded animals when they separated from the rest of the vertebrates.

Ribonucleotide incorporation, proofreading and bypass by human DNA polymerase δ

In both budding and fission yeast, a large number of ribonucleotides are incorporated into DNA during replication by the major replicative polymerases (Pols α, δ and ɛ). They are subsequently removed by RNase H2-dependent repair, which if defective leads to replication stress and genome instability. To extend these studies to humans, where an RNase H2 defect results in an autoimmune disease, here we compare the ability of human and yeast Pol δ to incorporate, proofread, and bypass ribonucleotides during DNA synthesis. In reactions containing nucleotide concentrations estimated to be present in mammalian cells, human Pol δ stably incorporates one rNTP for approximately 2000 dNTPs, a ratio similar to that for yeast Pol δ. This result predicts that human Pol δ may introduce more than a million ribonucleotides into the nuclear genome per replication cycle, an amount recently reported to be present in the genome of RNase H2-defective mouse cells. Consistent with such abundant stable incorporation, we show that the 3'-exonuclease activity of yeast and human Pol δ largely fails to edit ribonucleotides during polymerization. We also show that, like yeast Pol δ, human Pol δ pauses as it bypasses ribonucleotides in DNA templates, with four consecutive ribonucleotides in a DNA template being more problematic than single ribonucleotides. In conjunction with recent studies in yeast and mice, this ribonucleotide incorporation may be relevant to impaired development and disease when RNase H2 is defective in mammals. As one tool to investigate ribonucleotide incorporation by Pol δ in human cells, we show that human Pol δ containing a Leu606Met substitution in the polymerase active site incorporates 7-fold more ribonucleotides into DNA than does wild type Pol δ.

Two thymidine kinases and one multisubstrate deoxyribonucleoside kinase salvage DNA precursors in Arabidopsis thaliana

Deoxyribonucleotides are the building blocks of DNA and can be synthesized via de novo and salvage pathways. Deoxyribonucleoside kinases (EC 2.7.1.145) salvage deoxyribonucleosides by transfer of a phosphate group to the 5' of a deoxyribonucleoside. This salvage pathway is well characterized in mammals, but in contrast, little is known about how plants salvage deoxyribonucleosides. We show that during salvage, deoxyribonucleosides can be phosphorylated by extracts of Arabidopsis thaliana into corresponding monophosphate compounds with an unexpected preference for purines over pyrimidines. Deoxyribonucleoside kinase activities were present in all tissues during all growth stages. In the A. thaliana genome, we identified two types of genes that could encode enzymes which are involved in the salvage of deoxyribonucleosides. Thymidine kinase activity was encoded by two thymidine kinase 1 (EC 2.7.1.21)-like genes (AtTK1a and AtTK1b). Deoxyadenosine, deoxyguanosine and deoxycytidine kinase activities were encoded by a single AtdNK gene. T-DNA insertion lines of AtTK1a and AtTK1b mutant genes had normal growth, although AtTK1a AtTK1b double mutants died at an early stage, which indicates that AtTK1a and AtTK1b catalyze redundant reactions. The results obtained in the present study suggest a crucial role for the salvage of thymidine during early plant development.

Pasteurella multocida thymidine kinase 1 efficiently activates pyrimidine nucleoside analogs

In the Pasteurella multocida genome only one putative deoxyribonucleoside kinase encoding gene, for thymidine kinase 1 (PmTK1), was identified. The PmTK1 gene was sub-cloned into Escherichia coli KY895 and it sensitized the host towards 2',2'-difluoro-deoxycytidine (gemcitabine, dFdC), 3'-azido-thymidine (AZT) and 5-fluoro-deoxyuridine (5F-dU). PmTK1 was over-expressed and purified with two different tags. Apparently, deoxyuridine (dU), and not thymidine (dT), is the preferred substrate. We suggest that PmTK1s could be employed as a species-specific activator of uracil-based nucleoside antibiotics.

Nucleoside analogues are activated by bacterial deoxyribonucleoside kinases in a species-specific manner

OBJECTIVES

To investigate the bactericidal activity of antiviral and anticancer nucleoside analogues against a variety of pathogenic bacteria and characterize the activating enzymes, deoxyribonucleoside kinases (dNKs).

METHODS

Several FDA-approved nucleoside analogue drugs were screened for their potential bactericidal activity against several clinical bacterial isolates and type strains. We identified and subcloned the genes coding for putative deoxyribonucleoside kinases in Escherichia coli, Pasteurella multocida, Salmonella enterica, Yersinia enterocolitica, Bacillus cereus, Clostridium perfringens and Listeria monocytogenes. These genes were tested for their ability to increase the susceptibility of a dNK-deficient E. coli strain to various analogues. We overexpressed, purified and characterized the substrate specificity and kinetic properties of the recombinant enzymes from S. enterica and B. cereus.

RESULTS

The tested Gram-negative bacteria were susceptible to 3'-azido-3'-deoxythymidine (AZT) in the concentration range 0.032-31.6 microM except for a single E. coli isolate and two Pseudomonas aeruginosa isolates which were resistant to the tested AZT concentrations. Purified recombinant S. enterica thymidine kinase phosphorylated AZT efficiently with a Km of 73.3 microM and k(cat)/Km of 6.6 x 10(4) s(-1) M(-1) and is the activator of this drug in vivo. 2',2'-Difluoro-2'-deoxycytidine (gemcitabine) was a potent antibiotic against Gram-positive bacteria in the concentration range between 0.001 and 1.0 microM. The B. cereus deoxyadenosine kinase had a Km for gemcitabine of 33.5 microM and k(cat)/Km of 5.1 x 10(3) s(-1) M(-1) and activates gemcitabine in vivo. S. enterica and B. cereus are now amongst the first bacteria with a completely characterized set of dNK enzymes.

CONCLUSIONS

Bacterial dNKs efficiently activate nucleoside analogues in a species-specific manner. Therefore, nucleoside analogues have a potential to be employed as antibiotics in the fight against emerging multiresistant bacteria.

Deoxyribonucleoside kinases activate nucleoside antibiotics in severely pathogenic bacteria

Common bacterial pathogens are becoming progressively more resistant to traditional antibiotics, representing a major public-health crisis. Therefore, there is a need for a variety of antibiotics with alternative modes of action. In our study, several nucleoside analogs were tested against pathogenic staphylococci and streptococci. We show that pyrimidine-based nucleoside analogs, like 3'-azido-3'-deoxythymidine (AZT) and 2',2'-difluoro-2'deoxycytidine (gemcitabine), are specifically activated by the endogenous bacterial deoxyribonucleoside kinases, leading to cell death. Deoxyribonucleoside kinase-deficient Escherichia coli strains become highly susceptible to nucleoside analogs when they express recombinant kinases from Staphylococcus aureus or Streptococcus pyogenes. We further demonstrate that recombinant S. aureus deoxyadenosine kinase efficiently phosphorylates the anticancer drug gemcitabine in vitro and is therefore the key enzyme in the activation pathway. When adult mice were infected intraperitoneally with a fatal dose of S. pyogenes strain AP1 and afterwards received gemcitabine, they failed to develop a systemic infection. Nucleoside analogs may therefore represent a promising alternative for combating pathogenic bacteria.

Thymidine kinase diversity in bacteria

Thymidine kinases (TKs) appear to be almost ubiquitous and are found in nearly all prokaryotes, eukaryotes, and several viruses. They are the key enzymes in thymidine salvage and activation of several anti-cancer and antiviral drugs. We show that bacterial TKs can be subdivided into 2 groups. The TKs from Gram-positive bacteria are more closely related to the eukaryotic TK1 enzymes than are TKs from Gram-negative bacteria.

Thymidine kinases in archaea

Twenty-six fully sequenced archaeal genomes were searched for genes coding for putative deoxyribonucleoside kinases (dNKs). We identified only 5 human-like thymidine kinase 1 genes (TK1s) and none for non-TK1 kinases. Four TK1s were identified in the Euryarchaea and one was found in the Crenarchaea, while none was found in Nanoarchaeum. The identified TK1s have high identity to Gram-positive bacteria TK1s. The TK1s from archaea, Gram-positive bacteria and eukaryotes share the same common ancestor, while the TK1s from Gram-negative bacteria belong to a less-related subgroup. It seems that a functional deoxyribonucleoside salvage pathway is not crucial for the archaeal cell.

Structural studies of thymidine kinases from Bacillus anthracis and Bacillus cereus provide insights into quaternary structure and conformational changes upon substrate binding

Thymidine kinase (TK) is the key enzyme in salvaging thymidine to produce thymidine monophosphate. Owing to its ability to phosphorylate nucleoside analogue prodrugs, TK has gained attention as a rate-limiting drug activator. We describe the structures of two bacterial TKs, one from the pathogen Bacillus anthracis in complex with the substrate dT, and the second from the food-poison-associated Bacillus cereus in complex with the feedback inhibitor dTTP. Interestingly, in contrast with previous structures of TK in complex with dTTP, in this study dTTP occupies the phosphate donor site and not the phosphate acceptor site. This results in several conformational changes compared with TK structures described previously. One of the differences is the way tetramers are formed. Unlike B. anthracis TK, B. cereus TK shows a loose tetramer. Moreover, the lasso-domain is in open conformation in B. cereus TK without any substrate in the active site, whereas in B. anthracis TK the loop conformation is closed and thymidine occupies the active site. Another conformational difference lies within a region of 20 residues that we refer to as phosphate-binding beta-hairpin. The phosphate-binding beta-hairpin seems to be a flexible region of the enzyme which becomes ordered upon formation of hydrogen bonds to the alpha-phosphate of the phosphate donor, dTTP. In addition to descriptions of the different conformations that TK may adopt during the course of reaction, the oligomeric state of the enzyme is investigated.