Inverse radiation dose-rate effects on somatic and germ-line mutations and DNA damage rates

The mutagenic effect of low linear energy transfer ionizing radiation is reduced for a given dose as the dose rate (DR) is reduced to a low level, a phenomenon known as the direct DR effect. Our reanalysis of published data shows that for both somatic and germ-line mutations there is an opposite, inverse DR effect, with reduction from low to very low DR, the overall dependence of induced mutations being parabolically related to DR, with a minimum in the range of 0.1 to 1.0 cGy/min (rule 1). This general pattern can be attributed to an optimal induction of error-free DNA repair in a DR region of minimal mutability (MMDR region). The diminished activation of repair at very low DRs may reflect a low ratio of induced (“signal”) to spontaneous background DNA damage (“noise”). Because two common DNA lesions, 8-oxoguanine and thymine glycol, were already known to activate repair in irradiated mammalian cells, we estimated how their rates of production are altered upon radiation exposure in the MMDR region. For these and other abundant lesions (abasic sites and single-strand breaks), the DNA damage rate increment in the MMDR region is in the range of 10% to 100% (rule 2). These estimates suggest a genetically programmed optimatization of response to radiation in the MMDR region.

Distamycin A modulates the sequence specificity of DNA alkylation by duocarmycin A

Duocarmycin A (Duo) normally alkylates adenine N3 at the 3′ end of A+T-rich sequences in DNA. The efficient adenine alkylation by Duo is achieved by its monomeric binding to the DNA minor groove. The addition of another minor groove binder, distamycin A (Dist), dramatically modulates the site of DNA alkylation by Duo, and the alkylation switches preferentially to G residues in G+C-rich sequences. HPLC product analysis using oligonucleotides revealed a highly efficient G–N3 alkylation via the cooperative binding of a heterodimer between Duo and Dist to the minor groove. The three-dimensional structure of the ternary alkylated complex of Duo/Dist/d(CAGGTGGT)·d(ACCACCTG) has been determined by nuclear Overhauser effect (NOE)-restrained refinement using 750 MHz two-dimensional NOE spectroscopy data. The refined NMR structure fully explains the sequence requirement of such modulated alkylations. This is the first demonstration of Duo DNA alkylation through cooperative binding with another structurally different natural product, and it suggests a promising new way to alter or modify the DNA alkylation selectivity in a predictable manner.

Interactions of RecF protein with RecO, RecR, and single-stranded DNA binding proteins reveal roles for the RecF–RecO–RecR complex in DNA repair and recombination

The products of the recF, recO, and recR genes are thought to interact and assist RecA in the utilization of single-stranded DNA precomplexed with single-stranded DNA binding protein (Ssb) during synapsis. Using immunoprecipitation, size-exclusion chromatography, and Ssb protein affinity chromatography in the absence of any nucleotide cofactors, we have obtained the following results: (i) RecF interacts with RecO, (ii) RecF interacts with RecR in the presence of RecO to form a complex consisting of RecF, RecO, and RecR (RecF–RecO–RecR); (iii) RecF interacts with Ssb protein in the presence of RecO. These data suggested that RecO mediates the interactions of RecF protein with RecR and with Ssb proteins. Incubation of RecF, RecO, RecR, and Ssb proteins resulted in the formation of RecF–RecO–Ssb complexes; i.e., RecR was excluded. Preincubation of RecF, RecO, and RecR proteins prior to addition of Ssb protein resulted in the formation of complexes consisting of RecF, RecO, RecR, and Ssb proteins. These data suggest that one role of RecF is to stabilize the interaction of RecR with RecO in the presence of Ssb protein. Finally, we found that interactions of RecF with RecO are lost in the presence of ATP. We discuss these results to explain how the RecF–RecO–RecR complex functions as an anti-Ssb factor.

Multiple independent transpositions of mitochondrial DNA control region sequences to the nucleus

Transpositions of mtDNA sequences to the nuclear genome have been documented in a wide variety of individual taxa, but little is known about their taxonomic frequency or patterns of variation. We provide evidence of nuclear sequences homologous to the mtDNA control region in seven species of diving ducks (tribe Aythyini). Phylogenetic analysis places each nuclear sequence as a close relative of the mtDNA haplotypes of the specie(s) in which it occurs, indicating that they derive from six independent transposition events, all occurring within the last ≈1.5 million years. Relative-rate tests and comparison of intraspecific variation in nuclear and mtDNA sequences confirm the expectation of a greatly reduced rate of evolution in the nuclear copies. By representing mtDNA haplotypes from ancestral populations, nuclear insertions may be valuable in some phylogenetic analyses, but they also confound the accurate determination of mtDNA sequences. In particular, our data suggest that the presumably nonfunctional but more slowly evolving nuclear sequences often will not be identifiable by changes incompatible with function and may be preferentially amplified by PCR primers based on mtDNA sequences from related taxa.

The catalytic subunit of DNA-dependent protein kinase selectively regulates p53-dependent apoptosis but not cell-cycle arrest

DNA damage induced by ionizing radiation (IR) activates p53, leading to the regulation of downstream pathways that control cell-cycle progression and apoptosis. However, the mechanisms for the IR-induced p53 activation and the differential activation of pathways downstream of p53 are unclear. Here we provide evidence that the catalytic subunit of DNA-dependent protein kinase (DNA-PKcs) serves as an upstream effector for p53 activation in response to IR, linking DNA damage to apoptosis. DNA-PKcs knockout (DNA-PKcs−/−) mice were exposed to whole-body IR, and the cell-cycle and apoptotic responses were examined in their thymuses. Our data show that IR induction of apoptosis and Bax expression, both mediated via p53, was significantly suppressed in the thymocytes of DNA-PKcs−/− mice. In contrast, IR-induced cell-cycle arrest and p21 expression were normal. Thus, DNA-PKcs deficiency selectively disrupts p53-dependent apoptosis but not cell-cycle arrest. We also confirmed previous findings that p21 induction was attenuated and cell-cycle arrest was defective in the thymoctyes of whole body-irradiated Atm−/− mice, but the apoptotic response was unperturbed. Taken together, our results support a model in which the upstream effectors DNA-PKcs and Atm selectively activate p53 to differentially regulate cell-cycle and apoptotic responses. Whereas Atm selects for cell-cycle arrest but not apoptosis, DNA-PKcs selects for apoptosis but not cell-cycle arrest.

Knowledge-based analysis of microarray gene expression data by using support vector machines

We introduce a method of functionally classifying genes by using gene expression data from DNA microarray hybridization experiments. The method is based on the theory of support vector machines (SVMs). SVMs are considered a supervised computer learning method because they exploit prior knowledge of gene function to identify unknown genes of similar function from expression data. SVMs avoid several problems associated with unsupervised clustering methods, such as hierarchical clustering and self-organizing maps. SVMs have many mathematical features that make them attractive for gene expression analysis, including their flexibility in choosing a similarity function, sparseness of solution when dealing with large data sets, the ability to handle large feature spaces, and the ability to identify outliers. We test several SVMs that use different similarity metrics, as well as some other supervised learning methods, and find that the SVMs best identify sets of genes with a common function using expression data. Finally, we use SVMs to predict functional roles for uncharacterized yeast ORFs based on their expression data.

p53-induced DNA bending and twisting: p53 tetramer binds on the outer side of a DNA loop and increases DNA twisting

DNA binding activity of p53 is crucial for its tumor suppressor function. Our recent studies have shown that four molecules of the DNA binding domain of human p53 (p53DBD) bind the response elements with high cooperativity and bend the DNA. By using A-tract phasing experiments, we find significant differences between the bending and twisting of DNA by p53DBD and by full-length human wild-type (wt) p53. Our data show that four subunits of p53DBD bend the DNA by 32–36°, whereas wt p53 bends it by 51–57°. The directionality of bending is consistent with major groove bends at the two pentamer junctions in the consensus DNA response element. More sophisticated phasing analyses also demonstrate that p53DBD and wt p53 overtwist the DNA response element by ≈35° and ≈70°, respectively. These results are in accord with molecular modeling studies of the tetrameric complex. Within the constraints imposed by the protein subunits, the DNA can assume a range of conformations resulting from correlated changes in bend and twist angles such that the p53–DNA tetrameric complex is stabilized by DNA overtwisting and bending toward the major groove at the CATG tetramers. This bending is consistent with the inherent sequence-dependent anisotropy of the duplex. Overall, the four p53 moieties are placed laterally in a staggered array on the external side of the DNA loop and have numerous interprotein interactions that increase the stability and cooperativity of binding. The novel architecture of the p53 tetrameric complex has important functional implications including possible p53 interactions with chromatin.

Regulation of the pAD1 sex pheromone response of Enterococcus faecalis by direct interaction between the cAD1 peptide mating signal and the negatively regulating, DNA-binding TraA protein

The Enterococcus faecalis conjugative plasmid pAD1 (60 kb) encodes a mating response to the recipient-produced peptide sex pheromone cAD1. The response involves two key plasmid-encoded regulatory proteins: TraE1, which positively regulates all or most structural genes relating to conjugation, and TraA, which binds DNA and negatively regulates expression of traE1. In vitro studies that included development of a DNA-associated protein-tag affinity chromatography technique showed that TraA (37.9 kDa) binds directly to cAD1 near its carboxyl-terminal end and, as a consequence, loses its affinity for DNA. Analyses of genetically modified TraA proteins indicated that truncations within the carboxyl-terminal 9 residues significantly affected the specificity of peptide-directed association/dissociation of DNA. The data support earlier observations that transposon insertions near the 3′ end of traA eliminated the ability of cells to respond to cAD1.

Structure of FokI has implications for DNA cleavage

FokI is a member an unusual class of restriction enzymes that recognize a specific DNA sequence and cleave nonspecifically a short distance away from that sequence. FokI consists of an N-terminal DNA recognition domain and a C-terminal cleavage domain. The bipartite nature of FokI has led to the development of artificial enzymes with novel specificities. We have solved the structure of FokI to 2.3 Å resolution. The structure reveals a dimer, in which the dimerization interface is mediated by the cleavage domain. Each monomer has an overall conformation similar to that found in the FokI–DNA complex, with the cleavage domain packing alongside the DNA recognition domain. In corroboration with the cleavage data presented in the accompanying paper in this issue of Proceedings, we propose a model for FokI DNA cleavage that requires the dimerization of FokI on DNA to cleave both DNA strands.

Determination of gene organization in individual haplotypes by analyzing single DNA fragments from single spermatozoa

To determine human Ig heavy chain variable region (VH) gene segment organization on individual homologous chromosomes, an efficient approach has been developed. Single spermatozoa were used as subjects for the study. Upon sperm lysis, VH regions in each sperm were randomly sheared into fragments by the random Brownian force. The fragments were separated from each other by aliquoting the lysate into a certain number of tubes. The gene segments in the VH1 and VH4 families in each tube were identified by denaturing gradient gel electrophoresis after PCR amplification. The polymorphic VH sequences were used to determine the parental origins of the analyzed sperm. VH segment organization in the parental haplotypes was determined by aligning the overlapping fragments from the spermatozoa with the corresponding haplotypes. Based on this comparison between the resulting haplotype maps and the composite map reported previously, the VH region on chromosome 14 could be subdivided into four portions. The numbers and compositions of the VH gene segments differ considerably among the maps in two portions, but are highly conserved in the other two. The data also indicate that the VH region on chromosome 15 may contain a large duplicated block with copy number varying among haplotypes. The approach used in the present study may be used to construct high-resolution haplotype maps without molecular cloning.

Analysis of variable (diversity) joining recombination in DNAdependent protein kinase (DNA-PK)-deficient mice reveals DNA-PK-independent pathways for both signal and coding joint formation

Previous studies have suggested that ionizing radiation causes irreparable DNA double-strand breaks in mice and cell lines harboring mutations in any of the three subunits of DNA-dependent protein kinase (DNA-PK) (the catalytic subunit, DNA-PKcs, or one of the DNA-binding subunits, Ku70 or Ku86). In actuality, these mutants vary in their ability to resolve double-strand breaks generated during variable (diversity) joining [V(D)J] recombination. Mutant cell lines and mice with targeted deletions in Ku70 or Ku86 are severely compromised in their ability to form coding and signal joints, the products of V(D)J recombination. It is noteworthy, however, that severe combined immunodeficient (SCID) mice, which bear a nonnull mutation in DNA-PKcs, are substantially less impaired in forming signal joints than coding joints. The current view holds that the defective protein encoded by the murine SCID allele retains enough residual function to support signal joint formation. An alternative hypothesis proposes that DNA-PKcs and Ku perform different roles in V(D)J recombination, with DNA-PKcs required only for coding joint formation. To resolve this issue, we examined V(D)J recombination in DNA-PKcs-deficient (SLIP) mice. We found that the effects of this mutation on coding and signal joint formation are identical to the effects of the SCID mutation. Signal joints are formed at levels 10-fold lower than in wild type, and one-half of these joints are aberrant. These data are incompatible with the notion that signal joint formation in SCID mice results from residual DNA-PKcs function, and suggest a third possibility: that DNA-PKcs normally plays an important but nonessential role in signal joint formation.

A terbenzimidazole that preferentially binds and conformationally alters structurally distinct DNA duplex domains: A potential mechanism for topoisomerase I poisoning

The terbenzimidazoles are a class of synthetic ligands that poison the human topoisomerase I (TOP1) enzyme and promote cancer cell death. It has been proposed that drugs of this class act as TOP1 poisons by binding to the minor groove of the DNA substrate of TOP1 and altering its structure in a manner that results in enzyme-mediated DNA cleavage. To test this hypothesis, we characterize and compare the binding properties of a 5-phenylterbenzimidazole derivative (5PTB) to the d(GA4T4C)2 and d(GT4A4C)2 duplexes. The d(GA4T4C)2 duplex contains an uninterrupted 8-bp A⋅T domain, which, on the basis of x-ray crystallographic data, should induce a highly hydrated “A-tract” conformation. This duplex also exhibits anomalously slow migration in a polyacrylamide gel, a feature characteristic of a noncanonical global conformational state frequently described as “bent.” By contrast, the d(GT4A4C)2 duplex contains two 4-bp A⋅T tracts separated by a TpA dinucleotide step, which should induce a less hydrated “B-like” conformation. This duplex also migrates normally in a polyacrylamide gel, a feature further characteristic of a global, canonical B-form duplex. Our data reveal that, at 20°C, 5PTB exhibits an ≈2.3 kcal/mol greater affinity for the d(GA4T4C)2 duplex than for the d(GT4A4C)2 duplex. Significantly, we find this sequence/conformational binding specificity of 5PTB to be entropic in origin, an observation consistent with a greater degree of drug binding-induced dehydration of the more solvated d(GA4T4C)2 duplex. By contrast with the differential duplex affinity exhibited by 5PTB, netropsin and 4′,6-diamidino-2-phenylindole (DAPI), two AT-specific minor groove binding ligands that are inactive as human TOP1 poisons, bind to both duplexes with similar affinities. The electrophoretic behaviors of the ligand-free and ligand-bound duplexes are consistent with 5PTB-induced bending and/or unwinding of both duplexes, which, for the d(GA4T4C)2 duplex, is synergistic with the endogenous sequence-directed electrophoretic properties of the ligand-free duplex state. By contrast, the binding to either duplex of netropsin or DAPI induces little or no change in the electrophoretic mobilities of the duplexes. Our results demonstrate that the TOP1 poison 5PTB binds differentially to and alters the structures of the two duplexes, in contrast to netropsin and DAPI, which bind with similar affinities to the two duplexes and do not significantly alter their structures. These results are consistent with a mechanism for TOP1 poisoning in which drugs such as 5PTB differentially target conformationally distinct DNA sites and induce structural changes that promote enzyme-mediated DNA cleavage.

Dynamic modeling of gene expression data

We describe the time evolution of gene expression levels by using a time translational matrix to predict future expression levels of genes based on their expression levels at some initial time. We deduce the time translational matrix for previously published DNA microarray gene expression data sets by modeling them within a linear framework by using the characteristic modes obtained by singular value decomposition. The resulting time translation matrix provides a measure of the relationships among the modes and governs their time evolution. We show that a truncated matrix linking just a few modes is a good approximation of the full time translation matrix. This finding suggests that the number of essential connections among the genes is small.

Stimulation of human 8-oxoguanine-DNA glycosylase by AP-endonuclease: potential coordination of the initial steps in base excision repair

8-Oxoguanine-DNA glycosylase 1 (OGG1), with intrinsic AP lyase activity, is the major enzyme for repairing 7,8-dihydro-8-oxoguanine (8-oxoG), a critical mutagenic DNA lesion induced by reactive oxygen species. Human OGG1 excised the damaged base from an 8-oxoG·C-containing duplex oligo with a very low apparent kcat of 0.1 min–1 at 37°C and cleaved abasic (AP) sites at half the rate, thus leaving abasic sites as the major product. Excision of 8-oxoG by OGG1 alone did not follow Michaelis–Menten kinetics. However, in the presence of a comparable amount of human AP endonuclease (APE1) the specific activity of OGG1 was increased ∼5-fold and Michaelis–Menten kinetics were observed. Inactive APE1, at a higher molar ratio, and a bacterial APE (Nfo) similarly enhanced OGG1 activity. The affinity of OGG1 for its product AP·C pair (Kd ∼ 2.8 nM) was substantially higher than for its substrate 8-oxoG·C pair (Kd ∼ 23.4 nM) and the affinity for its final β-elimination product was much lower (Kd ∼ 233 nM). These data, as well as single burst kinetics studies, indicate that the enzyme remains tightly bound to its AP product following base excision and that APE1 prevents its reassociation with its product, thus enhancing OGG1 turnover. These results suggest coordinated functions of OGG1 and APE1, and possibly other enzymes, in the DNA base excision repair pathway.

Polarity of human replication protein A binding to DNA

Replication protein A (RPA), the nuclear single-stranded DNA binding protein is involved in DNA replication, nucleotide excision repair (NER) and homologous recombination. It is a stable heterotrimer consisting of subunits with molecular masses of 70, 32 and 14 kDa (p70, p32 and p14, respectively). Gapped DNA structures are common intermediates during DNA replication and NER. To analyze the interaction of RPA and its subunits with gapped DNA we designed structures containing 9 and 30 nucleotide gaps with a photoreactive arylazido group at the 3′-end of the upstream oligonucleotide or at the 5′-end of the downstream oligonucleotide. UV crosslinking and subsequent analysis showed that the p70 subunit mainly interacts with the 5′-end of DNA irrespective of DNA structure, while the subunit orientation towards the 3′-end of DNA in the gap structures strongly depends on the gap size. The results are compared with the data obtained previously with the primer–template systems containing 5′- or 3′-protruding DNA strands. Our results suggest a model of polar RPA binding to the gapped DNA.