Targeting critical regions in genomic DNA with AT-specific anticancer drugs

Cellular DNA is not a uniform target for DNA-reactive drugs. At the nucleotide level, drugs recognize and bind short motifs of a few base pairs. The location of drug adducts at the genomic level depends on how these short motifs are distributed in larger domains. This aspect, referred to as region specificity, may be critical for the biological outcome of drug action. Recent studies demonstrated that certain minor groove binding (MGB) drugs, such as bizelesin, produce region-specific lesions in cellular DNA. Bizelesin binds mainly T(A/T)(4)A sites, which are on average scarce, but occasionally cluster in distinct minisatellite regions (200-1000 bp of approximately 85-100% AT), herein referred to as AT islands. Bizelesin-targeted AT islands are likely to function as strong matrix attachment regions (MARs), domains that organize DNA loops on the nuclear matrix. Distortion of MAR-like AT islands may be a basis for the observed inhibition of new replicon initiation and the extreme lethality of bizelesin adducts (<10 adducts/cell for cell growth inhibition). Hence, long AT-islands represent a novel class of critical targets for anticancer drugs. The AT island paradigm illustrates the potential of the concept of regional targeting as an essential component of the rational design of new sequence-specific DNA-reactive drugs.     

The genome factor in region-specific DNA damage: the DNA-reactive drug U-78779 prefers mixed A/T-G/C sequences at the nucleotide level but is region-specific for long pure AT islands at the genomic level

Bizelesin is the first anticancer drug capable of damaging specific regions of the genome with clusters of its binding sites T(A/T)(4)A. This study characterized the sequence- and region-specificity of a bizelesin analogue, U-78779, designed to interact with mixed A/T-G/C motifs. At the nucleotide level, U-78779 was found to prefer runs of A/Ts interspersed with 1 or 2 G/C pairs, although 25% of the identified sites corresponded to pure AT motifs similar to bizelesin sites. The in silico computational analysis showed that the preferred mixed A/T-G/C motifs distribute uniformly at the genomic level. In contrast, the secondary, pure AT motifs (A/T)(6)A were found densely clustered in the same long islands of AT-rich DNA that bizelesin targets. Mapping the sites and quantitating the frequencies of U-78779 adducts in model AT island and non-AT island naked DNAs demonstrated that clusters of pure AT motifs outcompete isolated mixed A/T-G/C sites in attracting drug binding. Regional preference of U-78779 for AT island domains was verified also in DNA from drug-treated cells. Thus, while the primary sequence preference gives rise to non-region-specific scattered lesions, the clustering of the minor pure AT binding motifs seems to determine region-specificity of U-78779 in the human genome. The closely correlated cytotoxic activities of U-78779 and bizelesin in several cell lines further imply that both drugs may share common cellular targets. This study underscores the significance of the genome factor in a drug's potential for region-specific DNA damage, by showing that it can take precedence over drug binding preferences at the nucleotide level.     

Region-specific DNA damage by AT-specific DNA-reactive drugs is predicted by drug binding specificity

Bizelesin and adozelesin are DNA-reactive antitumor drugs that alkylate adenines at the 3' ends of their preferred binding sites [5'T(A/T)(4)A3'and 5'(A/T)(3)(-4)A3', respectively]. We used these drugs to examine the determinants for region-specific damage of human genomic DNA. The distribution of bizelesin binding motifs in several regions analyzed "in silico" correlated well with the experimentally determined lesions in these regions assessed by quantitative polymerase chain reaction (QPCR) stop assay. In contrast to the typically low motif density, clusters of potential bizelesin binding sites were found in the matrix-associated regions (MAR domains) of the c-myc and apolipoprotein B (apoB) genes. Accordingly, lesions induced by bizelesin in these domains (2.13 and 7.06 lesions kbp(-1) microM(-1), respectively) markedly exceeded lesions in bulk DNA (0.87 lesions kbp(-1) microM(-1)) or in regions with typically low motif density (e.g., 0.75 and 0.87 lesions kbp(-1) microM(-1) in a beta-globin gene and c-myc origin of replication regions, respectively). Consistent with the more frequent, less localized adozelesin motif, actual lesions induced by adozelesin exceeded by severalfold lesions by bizelesin in four selected regions (within the c-myc and HPRT loci). Whereas adozelesin is likely to affect similar regions as bizelesin, adozelesin's more promiscuous binding probably compromises its relative specificity for such targets. In contrast, findings for bizelesin provide for the first time a proof of principle that a small molecular weight drug can preferentially damage specific regions in cellular DNA. Targeting of critical repetitive sequences, such as AT-rich MAR domains, which allow for clustering of drug binding motif, can be the paradigm for region specificity of small molecular weight agents.     

Tallimustine lesions in cellular DNA are AT sequence-specific but not region-specific.

Tallimustine (FCE 24517) is an AT-specific alkylating antitumor derivative of distamycin. This study examined levels of tallimustine lesions in intracellular DNA, their sequence- and region-specificity, and the long-range distribution of the drug binding motif. Tallimustine adducts in DNA converted to strand breaks by heating allowed the quantitation of drug lesions. In bulk DNA of intact human leukemia CEM cells, tallimustine formed 0.15 +/- 0.04 and 0.64 +/- 0.18 lesions/kbp at 5 and 50 microM, respectively. These lesions represent monoadducts as no interstrand cross-links or DNA-protein cross-links were detected. Tallimustine adducts in intracellularly treated DNA showed a general preference for sequences with T-tracts, suggesting a propensity for intrinsically bent motifs. Major drug-adducted sites identified by repetitive primer extension, included 5'-TTTTGPu-3' and 5'-TTTTGC-3' motif. Despite the high specificity at the nucleotide level, tallimustine did not differentiate among bulk DNA and three discrete AT-rich regions of genomic DNA examined by quantitative PCR stop assay with lesion frequencies ranging from 0.23 to 0.39 lesions/kbp at 25 microM drug. In comparisons of lesion frequencies and cytotoxicity, tallimustine adducts are approximately 50 times more lethal than relatively nonsequence specific cisplatin adducts but are >100 times less lethal than lesions by an unrelated AT-specific drug, bizelesin. However, the 5'-TTTTGPu-3' motifs targeted by tallimustine are relatively infrequent and scattered throughout the genome. In contrast, the motifs 5'-T(A/T)(4)A-3' motifs targeted by bizelesin, while also infrequent, cluster in defined AT-rich islands. The lack of region-specificity may be the reason tallimustine adducts, despite high AT-specificity at the nucleotide level, are less lethal than region-specific bizelesin adducts.     

Matrix attachment region (MAR) properties and abnormal expansion of AT island minisatellites in FRA16B fragile sites in leukemic CEM cells

AT-rich minisatellites (AT islands) are sites of genomic instability in cancer cells and targets for extremely lethal AT-specific drugs, such as bizelesin. Here we investigated the AT islands in the FRA16B fragile site region for their possible roles in the organization of DNA on the nuclear matrix. The FRA16B AT island nominally spans ~3 kb of mostly >90% A/T DNA. In silico analysis indicates that this domain exhibits characteristics of nuclear matrix attachment regions (MARs): an exceptionally intense computed ‘MAR potential’ and profound duplex destabilization and flexibility. FRA16B repeats specifically bind to isolated nuclear matrices, which indicates their in vitro MAR function. This binding is several-fold greater than that of a known MAR in the c-myc gene. AT islands in fragile sites FRA16B and FRA16D are significantly more abundant in CEM cells that are hypersensitive to bizelesin compared to normal WI-38 cells. FRA16B overabundance in CEM is due to an ~10-fold expansion of FRA16B repeats. The expanded FRA16B minisatellites in CEM cells preferentially localize to the nuclear matrix-associated DNA indicating their in vivo MAR function. The unexpanded repeats in WI-38 cells localize to the loop DNA. The c-myc MAR is also matrix-associated in CEM cells while localizing to loop DNA in WI-38 cells. These results are the first to demonstrate that AT islands in fragile sites can function as MARs both in vitro and in vivo. The ability of FRA16B-mediated MAR sites to rearrange depending on the repeat expansion status could be relevant to both genomic instability of cancer cells and their sensitivity to AT-island targeting drugs.     

Induction of AT-specific DNA-interstrand crosslinks by bizelesin in genomic and simian virus 40 DNA

Bizelesin is a bifunctional AT-specific DNA alkylating drug. Our study characterized the ability of bizelesin to induce interstrand crosslinks, a potential lethal lesion. In genomic DNA of BSC-1 cells, bizelesin formed from approx. 0.3 to 6.03+/-0.85 interstrand crosslinks per 106 base pairs, at 5-100 nM drug concentration, respectively, comparable to the number of total adducts previously determined in the same system (J.M. Woynarowski, M.M. McHugh, L.S. Gawron, T.A. Beerman, Biochemistry 34 (1995) 13042-13050). Bizelesin did not induce DNA-protein crosslinks or strand breaks. A model defined target, intracellular simian virus 40 (SV40) DNA, was employed to map at the nucleotide level sites of bizelesin adducts, including potential interstrand crosslinks. Preferential adduct formation was observed at AT tracts which are abundant in the SV40 matrix associated region and the origin of replication. Many sites, including each occurrence of 5'-T(A/T)4A-3', co-mapped on both DNA strands suggesting interstrand crosslinks, although monoadducts were also formed. Bizelesin adducts in naked SV40 DNA were found at similar sites. The localization of bizelesin-induced crosslinks in AT-rich tracts of replication-related regions is consistent with the potent anti-replicative properties of bizelesin. Given the apparent lack of other types of lesions in genomic DNA, interstrand crosslinks localized in AT-rich tracts, and to some extent perhaps also monoadducts, are likely to be lethal effects of bizelesin.     

Interaction of echinomycin with An.Tn. and (AT)n regions flanking its CG binding site.

We have prepared DNA fragments containing the sequences A15CGT15, T15CGA15 and T(AT)8CG(AT)15 cloned within the SmaI site of the pUC19 polylinker. These have been used as substrates in footprinting experiments with DNase I and diethylpyrocarbonate probing the effects of echinomycin, binding to the central CG, on the structure of the surrounding sequences. No clear DNase I footprints are seen with T15CGA15 though alterations in the nuclease susceptibility of surrounding regions suggest that the ligand is binding, albeit weakly at this site. All the other fragments show the expected footprints around the CG site. Regions of An and Tn are rendered much more reactive to DNase I and adenines on the 3'-side of the CG become hyperreactive to diethylpyrocarbonate. Regions of alternating AT show unusual changes in the presence of the ligand. At low concentrations (5 microM) cleavage of TpA is enhanced, whereas at higher concentrations a cleavage pattern with a four base pair repeat is evident. A similar pattern is seen with micrococcal nuclease. Modification by diethylpyrocarbonate is strongest at alternate adenines which are staggered in the 5'-direction across the two strands. We interpret these changes by suggesting secondary drug binding within regions of alternating AT, possibly to the dinucleotide ApT. DNase I footprinting experiments performed at 4 degrees C revealed neither enhancements nor footprints for flanking regions of homopolymeric A and T suggesting that the conformational changes are necessary consequence of drug binding.     

Footprinting studies of DNA-sequence recognition by nogalamycin.

We have studied the DNA sequence binding preference of the antitumour antibiotic nogalamycin by DNase-I footprinting using a variety of DNA fragments. The DNA fragments were obtained by cloning synthetic oligonucleotides into longer DNA fragments and were designed to contain isolated ligand-binding sites surrounded by repetitive sequences such as (A)n.(T)n and (AT)n. Within regions of (A)n.(T)n, clear footprints are observed with low concentrations of nogalamycin (< 5 microM), with apparent binding affinities for tetranucleotide sequences which decrease in the order TGCA > AGCT = ACGT > TCGA. In contrast, within regions of (AT)n, the ligand binds best to AGCT; binding to TCGA and TGCA is no stronger than to alternating AT. Within (ATT)n, the preference is for ACGT > TCGA. Although each of these binding sites contains all four base pairs, there is no apparent consensus sequence, suggesting that the selectivity is affected by local DNA dynamic and structural effects. At higher drug concentrations (> 25 microM), nogalamycin prevents DNAse-I cleavage of (AT)n but shows no interaction with regions of (AC)n.(GT)n. Regions of (A)n.(T)n, which are poorly cut by DNase I, show enhanced rates of cleavage in the presence of low concentrations of nogalamycin, but are protected from cleavage at higher concentrations. We suggest that this arises because drug binding to adjacent regions distorts the DNA to a structure which is more readily cut by the enzyme and which is better able to bind further ligand molecules.     

Characterization of randomly-obtained matrix attachment regions (MARs) from higher plants.

Matrix attachment regions (MARs) can be operationally defined as DNA fragments that bind to the nuclear matrix. We have created a library of randomly obtained MARs from tobacco (Nicotiana tobacum) by cloning DNA fragments that co-isolate with nuclear matrixes prepared by a method involving lithium diiodosalicylate. The interactions of several of the cloned MARs with nuclear matrixes were tested by an in vitro binding assay in which genomic DNA was used as competitor. Based on this assay, the MARs were classified as strong, medium, and weak binders. Examples of each of the binding classes were further studied by in vitro binding using self- and cross-competition. Estimates of dissociation constants for several MARs ranged from 6 to 11 nM and correlated inversely with binding strength. The number of binding sites per matrix for several MARs ranged from 4 x 10(5) to 9 x 10(5) and correlated directly with binding strength. We conclude that binding strength, as we have measured it, is a function of both numbers of binding sites and affinity for the sites. The tobacco MARs were sequenced and analyzed for overall AT content, for distribution of AT-rich regions, and for the abundance of several MAR-related motifs. Previously identified MAR motifs correlate to various degrees with binding strength. Notably, the Drosophila topoisomerase II motif does not correlate with binding strength of the tobacco MARs. A newly identified motif, the "90%AT Box," correlates better with binding strength than any of the previously identified motifs we investigated.     

Introns of the chicken ovalbumin gene promote nucleosome alignment in vitro.

A defined in vitro chromatin assembly system was used to examine the nucleosome alignment induced by histone H5 throughout a 12 kilobase pair chicken genomic DNA fragment containing the ovalbumin gene. In contrast with total fragmented chicken DNA and several anonymous cloned fragments, much of the gene permitted histone H5 to space nucleosomes at physiological intervals in an extended array. Nucleosomes at the 3'-end of the gene and on approximately 4 kilobase pairs of 5'-flanking ovalbumin sequence did not become aligned to appreciable extents. Analysis of cloned 2-3 kilobase pair subfragments suggested that a strong nucleosome alignment signal, specifying a 196 +/- 5 base pair repeat exists in intron E. A second discrete region of the gene, which mapped approximately to intron A, exhibited nucleosome alignment with a spacing periodicity of about 200 base pairs. The ovalbumin cDNA did not permit nucleosome alignment. These findings suggest that some of the introns contain signals that direct nucleosome alignment over the ovalbumin gene in a way conducive to its regulation.     

Characterization of frequencies and distribution of single nucleotide insertions/deletions in the human genome

Most of the studies on single nucleotide variations are on substitutions rather than insertions/deletions. In this study, we examined the distribution and characteristics of single nucleotide insertions/deletions (SNindels), using data available from dbSNP for all the human chromosomes. There are almost 300,000 SNindels in the database, of which only 0.8% are validated. They occur at the frequency of 0.887 per 10 kb on average for the whole genome, or approximately 1 for every 11,274 bp. More than half occur in regions with mononucleotide repeats the longest of which is 47 bases. Overall the mononucleotide repeats involving C and G are much shorter than those for A and T. About 12% are surrounded by palindromes. There is general correlation between chromosome size and total number for each chromosome. Inter-chromosomal variation in density ranges from 0.6 to 21.7 per kilobase. The overall spectrum shows very high proportion of SNindel of types -/A and -/T at over 81%. The proportion of -/A and -/T SNindels for each chromosome is correlated to its AT content. Less than half of the SNindels are within or near known genes and even fewer (<0.183%) in coding regions, and more than 1.4% of -/C and -/G are in coding compared to 0.2% for -/A and -/T types. SNindels of -/A and -/T types make up 80% of those found within untranslated regions but less than 40% of those within coding regions. A separate analysis using the subset of 2324 validated SNindels showed slightly less AT bias of 74%, SNindels not within mononucleotide repeats showed even less AT bias at 58%. Density of validated SNindels is 0.007/10 kb overall and 90% are found within or near genes. Among all chromosomes, Y has the lowest numbers and densities for all SNindels, validated SNindels, and SNindels not within repeats.     

AHM1, a novel type of nuclear matrix-localized, MAR binding protein with a single AT hook and a J domain-homologous region

Interactions between the nuclear matrix and special regions of chromosomal DNA called matrix attachment regions (MARs) have been implicated in various nuclear functions. We have identified a novel protein from wheat, AT hook-containing MAR binding protein1 (AHM1), that binds preferentially to MARs. A multidomain protein, AHM1 has the special combination of a J domain-homologous region and a Zn finger-like motif (a J-Z array) and an AT hook. For MAR binding, the AT hook at the C terminus was essential, and an internal portion containing the Zn finger-like motif was additionally required in vivo. AHM1 was found in the nuclear matrix fraction and was localized in the nucleoplasm. AHM1 fused to green fluorescent protein had a speckled distribution pattern inside the nucleus. AHM1 is most likely a nuclear matrix component that functions between intranuclear framework and MARs. J-Z arrays can be found in a group of (hypothetical) proteins in plants, which may share some functions, presumably to recruit specific Hsp70 partners as co-chaperones.     

A map of nuclear matrix attachment regions within the breast cancer loss-of-heterozygosity region on human chromosome 16q22.1

There is abundant evidence that the DNA in eukaryotic cells is organized into loop domains that represent basic structural and functional units of chromatin packaging. To explore the DNA domain organization of the breast cancer loss-of-heterozygosity region on human chromosome 16q22.1, we have identified a significant portion of the scaffold/matrix attachment regions (S/MARs) within this region. Forty independent putative S/MAR elements were assigned within the 16q22.1 locus. More than 90% of these S/MARs are AT rich, with GC contents as low as 27% in 2 cases. Thirty-nine (98%) of the S/MARs are located within genes and 36 (90%) in gene introns, of which 15 are in first introns of different genes. The clear tendency of S/MARs from this region to be located within the introns suggests their regulatory role. The S/MAR resource constructed may contribute to an understanding of how the genes in the region are regulated and of how the structural architecture and functional organization of the DNA are related.     

核基质附着区与转基因表达

核基质附着区 (matrixattachmentregions,MARs)是与核基质 (或核骨架 )特异结合的DNA序列 ,属于非编码序列 ,富含AT。通过与核基质的结合 ,它能使染色质形成独立的环状结构 ,调控基因的转录和表达 ,减少由于位置效应引起的转基因沉默。MARs在提高转基因表达水平、消除转基因个体间表达水平的差异、抑制转基因沉默等方面起着重要的作用。就MAR的一些结构功能特征及其在基因工程中的应用等方面进行了阐述。     

Molecular cloning of Snyder-Theilen feline leukemia and sarcoma viruses: comparative studies of feline sarcoma virus with its natural helper virus and with Moloney murine sarcoma virus

Extrachromosomal DNA obtained from mink cells acutely infected with the Snyder-Theilen (ST) strain of feline sarcoma virus (feline leukemia virus) [FeSV(FeLV)] was fractionated electrophoretically, and samples enriched for FeLV and FeSV linear intermediates were digested with EcoRI and cloned in lambda phage. Hybrid phages were isolated containing either FeSV or FeLV DNA "inserts" and were characterized by restriction enzyme analysis, R-looping with purified 26 to 32S viral RNA, and heteroduplex formation. The recombinant phages (designated lambda FeSV and lambda FeLV) contain all of the genetic information represented in FeSV and FeLV RNA genomes but lack one extended terminally redundant sequence of 750 bases which appears once at each end of parental linear DNA intermediates. Restriction enzyme and heteroduplex analyses confirmed that sequences unique to FeSV (src sequences) are located at the center of the FeSV genome and are approximately 1.5 kilobase pairs in length. With respect to the 5'-3' orientation of genes in viral RNA, the order of genes in the FeSV genome is 5'-gag-src-env-c region-3'; only 0.9 kilobase pairs of gag and 0.6 kilobase pairs of env-derived FeLV sequences are represented in ST FeSV. Heteroduplex analyses between lambda FeSV or lambda FeLV DNA and Moloney murine sarcoma virus DNA (strain m1) were performed under conditions of reduced stringency to demonstrate limited regions of base pair homology. Two such regions were identified: the first occurs at the extreme 5' end of the leukemia and both sarcoma viral genomes, whereas the second corresponds to a 5' segment of leukemia virus "env" sequences conserved in both sarcoma viruses. The latter sequences are localized at the 3' end of FeSV src and at the 5' end of murine sarcoma virus src and could possibly correspond to regions of helper virus genomes that are required for retroviral transforming functions.     

DNA binding of methyl-CpG-binding protein MeCP2 in human MCF7 cells

MeCP2 has been identified as a chromatin-associated protein that recognizes MAR elements as well as methyl-CpGs. To characterize target sequences of MeCP2 in human cells, we employed two complementary methods. First, by use of a preparative chromatin immunoprecipitation protocol, we created from MCF7 cells a library enriched with sequences bound to MeCP2. A total of 154 representative clones were sequenced and analyzed. A large fraction of clones was found to be associated with retrotransposons, mostly with Alu repeats. A subgroup of four clones is derived from putative MARs; one clone is associated with a CpG island, and four clones contain alphoid repeats. Classical satellite DNAs II and III are not represented among clones, although they are heavily methylated. Second, using indirect immunofluorescence microscopy, we show that MeCP2 staining of human metaphase chromosomes has a dotted to knobby appearance with a reduced level of staining of centromeric regions of some chromosomes. On the other hand, an anti-5-methylcytosine antibody preferentially stained the juxtacentromeric regions of chromosomes 1, 9, and 16, which habor highly methylated, classical satellite DNAs, and methylated alphoid sequences in centromeric regions of several other chromosomes with reduced intensity. In interphase MCF7 cells, the distribution of MeCP2 exhibits a granular appearance throughout the nucleus. This distribution does not parallel that of methylated cytosine and heterochromatin. The selective binding behavior of MeCP2 revealed by these results (preference for murine major satellite DNA, Alu sequences, MARs, and CpG islands) is explained by its ability to recognize the sequence information (guanine bases) adjacent to CpG (TpG) as demonstrated in previous footprinting experiments.