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.

Detection of alkali metal ions in DNA crystals using state-of-the-art X-ray diffraction experiments

The observation of light metal ions in nucleic acids crystals is generally a fortuitous event. Sodium ions in particular are notoriously difficult to detect because their X-ray scattering contributions are virtually identical to those of water and Na+…O distances are only slightly shorter than strong hydrogen bonds between well-ordered water molecules. We demonstrate here that replacement of Na+ by K+, Rb+ or Cs+ and precise measurements of anomalous differences in intensities provide a particularly sensitive method for detecting alkali metal ion-binding sites in nucleic acid crystals. Not only can alkali metal ions be readily located in such structures, but the presence of Rb+ or Cs+ also allows structure determination by the single wavelength anomalous diffraction technique. Besides allowing identification of high occupancy binding sites, the combination of high resolution and anomalous diffraction data established here can also pinpoint binding sites that feature only partial occupancy. Conversely, high resolution of the data alone does not necessarily allow differentiation between water and partially ordered metal ions, as demonstrated with the crystal structure of a DNA duplex determined to a resolution of 0.6 Å.

Mechanism of stimulation of the DNA glycosylase activity of hOGG1 by the major human AP endonuclease: bypass of the AP lyase activity step

The generation of reactive oxygen species in the cell provokes, among other lesions, the formation of 8-oxo-7,8-dihydroguanine (8-oxoG) in DNA. Due to mispairing with adenine during replication, 8-oxoG is highly mutagenic. To minimise the mutagenic potential of this oxidised purine, human cells have a specific 8-oxoG DNA glycosylase/AP lyase (hOGG1) that initiates the base excision repair (BER) of 8-oxoG. We show here that in vitro this first enzyme of the BER pathway is relatively inefficient because of a high affinity for the product of the reaction it catalyses (half-life of the complex is >2 h), leading to a lack of hOGG1 turnover. However, the glycosylase activity of hOGG1 is stimulated by the major human AP endonuclease, HAP1 (APE1), the enzyme that performs the subsequent step in BER, as well as by a catalytically inactive mutant (HAP1-D210N). In the presence of HAP1, the AP sites generated by the hOGG1 DNA glycosylase can be occupied by the endonuclease, avoiding the re-association of hOGG1. Moreover, the glycosylase has a higher affinity for a non-cleaved AP site than for the cleaved DNA product generated by HAP1. This would shift the equilibrium towards the free glycosylase, making it available to initiate new catalytic cycles. In contrast, HAP1 does not affect the AP lyase activity of hOGG1. This stimulation of only the hOGG1 glycosylase reaction accentuates the uncoupling of its glycosylase and AP lyase activities. These data indicate that, in the presence of HAP1, the BER of 8-oxoG residues can be highly efficient by bypassing the AP lyase activity of hOGG1 and thus excluding a potentially rate limiting step.

ALFRED: an allele frequency database for diverse populations and DNA polymorphisms—an update

ALFRED (the ALelle FREquency Database) is designed to store and disseminate frequencies of alleles at human polymorphic sites for multiple populations, primarily for the population genetics and molecular anthropology communities. Currently ALFRED has information on over 180 polymorphic sites for more than 70 populations. Since our initial release of the database we have focussed on increasing the quantity and quality of data, making reciprocal links between ALFRED and other related databases, and providing useful tools to make the data more comprehensible to the end user. ALFRED is accessible from the Kidd Lab home page (http://info.med.yale.edu/genetics/kkidd/) or from ALFRED directly (http://alfred.med.yale.edu/alfred/index.asp).

The Homeodomain Resource: sequences, structures, DNA binding sites and genomic information

The Homeodomain Resource is an annotated collection of non-redundant protein sequences, three-dimensional structures and genomic information for the homeodomain protein family. Release 3.0 contains 795 full-length homeodomain-containing sequences, 32 experimentally-derived structures and 143 homeo­box loci implicated in human genetic disorders. Entries are fully hyperlinked to facilitate easy retrieval of the original records from source databases. A simple search engine with a graphical user interface is provided to query the component databases and assemble customized data sets. A new feature for this release is the addition of DNA recognition sites for all human homeodomain proteins described in the literature. The Homeodomain Resource is freely available through the World Wide Web at http://genome.nhgri.nih.gov/homeodomain.

Mechanism of topology simplification by type II DNA topoisomerases

Type II DNA topoisomerases actively reduce the fractions of knotted and catenated circular DNA below thermodynamic equilibrium values. To explain this surprising finding, we designed a model in which topoisomerases introduce a sharp bend in DNA. Because the enzymes have a specific orientation relative to the bend, they act like Maxwell's demon, providing unidirectional strand passage. Quantitative analysis of the model by computer simulations proved that it can explain much of the experimental data. The required sharp DNA bend was demonstrated by a greatly increased cyclization of short DNA fragments from topoisomerase binding and by direct visualization with electron microscopy.

Mapping the intrinsic curvature and flexibility along the DNA chain

The energy of DNA deformation plays a crucial and active role in its packaging and its function in the cell. Considerable effort has gone into developing methodologies capable of evaluating the local sequence-directed curvature and flexibility of a DNA chain. These studies thus far have focused on DNA constructs expressly tailored either with anomalous flexibility or curvature tracts. Here we demonstrate that these two structural properties can be mapped also along the chain of a “natural” DNA with any sequence on the basis of its scanning force microscope (SFM) images. To know the orientation of the sequence of the investigated DNA molecules in their SFM images, we prepared a palindromic dimer of the long DNA molecule under study. The palindromic symmetry also acted as an internal gauge of the statistical significance of the analysis carried out on the SFM images of the dimer molecules. It was found that although the curvature modulus is not efficient in separating static and dynamic contributions to the curvature of the population of molecules, the curvature taken with its direction (its sign in two dimensions) permits the direct separation of the intrinsic curvature from the flexibility contributions. The sequence-dependent flexibility seems to vary monotonically with the chain's intrinsic curvature; the chain rigidity was found to modulate as its local thermodynamic stability and does not correlate with the dinucleotide chain rigidities evaluation made from x-ray data by other authors.

Herpes simplex virus DNA packaging sequences adopt novel structures that are specifically recognized by a component of the cleavage and packaging machinery

The product of the herpes simplex virus type 1 UL28 gene is essential for cleavage of concatemeric viral DNA into genome-length units and packaging of this DNA into viral procapsids. To address the role of UL28 in this process, purified UL28 protein was assayed for the ability to recognize conserved herpesvirus DNA packaging sequences. We report that DNA fragments containing the pac1 DNA packaging motif can be induced by heat treatment to adopt novel DNA conformations that migrate faster than the corresponding duplex in nondenaturing gels. Surprisingly, these novel DNA structures are high-affinity substrates for UL28 protein binding, whereas double-stranded DNA of identical sequence composition is not recognized by UL28 protein. We demonstrate that only one strand of the pac1 motif is responsible for the formation of novel DNA structures that are bound tightly and specifically by UL28 protein. To determine the relevance of the observed UL28 protein–pac1 interaction to the cleavage and packaging process, we have analyzed the binding affinity of UL28 protein for pac1 mutants previously shown to be deficient in cleavage and packaging in vivo. Each of the pac1 mutants exhibited a decrease in DNA binding by UL28 protein that correlated directly with the reported reduction in cleavage and packaging efficiency, thereby supporting a role for the UL28 protein–pac1 interaction in vivo. These data therefore suggest that the formation of novel DNA structures by the pac1 motif confers added specificity on recognition of DNA packaging sequences by the UL28-encoded component of the herpesvirus cleavage and packaging machinery.

Predicting crossover generation in DNA shuffling

We introduce a quantitative framework for assessing the generation of crossovers in DNA shuffling experiments. The approach uses free energy calculations and complete sequence information to model the annealing process. Statistics obtained for the annealing events then are combined with a reassembly algorithm to infer crossover allocation in the reassembled sequences. The fraction of reassembled sequences containing zero, one, two, or more crossovers and the probability that a given nucleotide position in a reassembled sequence is the site of a crossover event are estimated. Comparisons of the predictions against experimental data for five example systems demonstrate good agreement despite the fact that no adjustable parameters are used. An in silico case study of a set of 12 subtilases examines the effect of fragmentation length, annealing temperature, sequence identity and number of shuffled sequences on the number, type, and distribution of crossovers. A computational verification of crossover aggregation in regions of near-perfect sequence identity and the presence of synergistic reassembly in family DNA shuffling is obtained.

Increased in vivo apoptosis in cells lacking mitochondrial DNA gene expression

We have attempted to determine whether loss of mtDNA and respiratory chain function result in apoptosis in vivo. Apoptosis was studied in embryos with homozygous disruption of the mitochondrial transcription factor A gene (Tfam) and tissue-specific Tfam knockout animals with severe respiratory chain deficiency in the heart. We found massive apoptosis in Tfam knockout embryos at embryonic day (E) 9.5 and increased apoptosis in the heart of the tissue-specific Tfam knockouts. Furthermore, mtDNA-less (ρ0) cell lines were susceptible to apoptosis induced by different stimuli in vitro. The data presented here provide in vivo evidence that respiratory chain deficiency predisposes cells to apoptosis, contrary to previous assumptions based on in vitro studies of cultured cells. These results suggest that increased apoptosis is a pathogenic event in human mtDNA mutation disorders. The finding that respiratory chain deficiency is associated with increased in vivo apoptosis may have important therapeutic implications for human disease. Respiratory chain deficiency and cell loss and/or apoptosis have been associated with neurodegeneration, heart failure, diabetes mellitus, and aging. Furthermore, chemotherapy and radiation treatment of cancer are intended to induce apoptosis in tumor cells. It would therefore be of interest to determine whether manipulation of respiratory chain function can be used to inhibit or enhance apoptosis in these conditions.

Methylation by a mutant T2 DNA [N6-adenine] methyltransferase expands the usage of RecA-assisted endonuclease (RARE) cleavage

Properties of a mutant bacteriophage T2 DNA [N6-adenine] methyltransferase (T2 Dam MTase) have been investigated for its potential utilization in RecA-assisted restriction endonuclease (RARE) cleavage. Steady-state kinetic analyses with oligonucleotide duplexes revealed that, compared to wild-type T4 Dam, both wild-type T2 Dam and mutant T2 Dam P126S had a 1.5-fold higher kcat in methylating canonical GATC sites. Additionally, T2 Dam P126S showed increased efficiencies in methylation of non-canonical GAY sites relative to the wild-type enzymes. In agreement with these steady-state kinetic data, when bacteriophage λ DNA was used as a substrate, maximal protection from restriction nuclease cleavage in vitro was achieved on the sequences GATC, GATN and GACY, while protection of GACR sequences was less efficient. Collectively, our data suggest that T2 Dam P126S can modify 28 recognition sequences. The feasibility of using the mutant enzyme in RARE cleavage with BclI and EcoRV endonucleases has been shown on phage λ DNA and with BclI and DpnII endonucleases on yeast chromosomal DNA embedded in agarose.

GenEST, a powerful bidirectional link between cDNA sequence data and gene expression profiles generated by cDNA-AFLP

The release of vast quantities of DNA sequence data by large-scale genome and expressed sequence tag (EST) projects underlines the necessity for the development of efficient and inexpensive ways to link sequence databases with temporal and spatial expression profiles. Here we demonstrate the power of linking cDNA sequence data (including EST sequences) with transcript profiles revealed by cDNA-AFLP, a highly reproducible differential display method based on restriction enzyme digests and selective amplification under high stringency conditions. We have developed a computer program (GenEST) that predicts the sizes of virtual transcript-derived fragments (TDFs) of in silico-digested cDNA sequences retrieved from databases. The vast majority of the resulting virtual TDFs could be traced back among the thousands of TDFs displayed on cDNA-AFLP gels. Sequencing of the corresponding bands excised from cDNA-AFLP gels revealed no inconsistencies. As a consequence, cDNA sequence databases can be screened very efficiently to identify genes with relevant expression profiles. The other way round, it is possible to switch from cDNA-AFLP gels to sequences in the databases. Using the restriction enzyme recognition sites, the primer extensions and the estimated TDF size as identifiers, the DNA sequence(s) corresponding to a TDF with an interesting expression pattern can be identified. In this paper we show examples in both directions by analyzing the plant parasitic nematode Globodera rostochiensis. Various novel pathogenicity factors were identified by combining ESTs from the infective stage juveniles with expression profiles of ∼4000 genes in five developmental stages produced by cDNA-AFLP.