Computational and in vitro analysis of destabilized DNA regions in the interferon gene cluster: potential of predicting functional gene domains

Recent evidence adds support to a traditional concept according to which the eukaryotic nucleus is organized into functional domains by scaffold or matrix attachment regions (S/MARs). These regions have previously been predicted to have a high potential for stress-induced duplex destabilization (SIDD). Here we report the parallel results of binding (reassociation) and computational SIDD analyses for regions within the human interferon gene cluster on the short arm of chromosome 9 (9p22). To verify and further refine the biomathematical methods, we focus on a 10 kb region in the cluster with the pseudogene IFNWP18 and the interferon alpha genes IFNA10 and IFNA7. In a series of S/MAR binding assays, we investigate the promoter and termination regions and additional attachment sequences that were detected in the SIDD profile. The promoters of the IFNA10 and the IFNA7 genes have a moderate approximately 20% binding affinity to the nuclear matrix; the termination sequences show stronger association (70-80%) under our standardized conditions. No comparable destabilized elements were detected flanking the IFNWP18 pseudogene, suggesting that selective pressure acts on the physicochemical properties detected here. In extended, noncoding regions a striking periodicity is found of rather restricted SIDD minima with scaffold binding potential. By various criteria, the underlying sequences represent a new class of S/MARs, thought to be involved in a higher level organization of the genome. Together, these data emphasize the relevance of SIDD calculations as a valid approach for the localization of structural, regulatory, and coding regions in the eukaryotic genome.     

The analysis of stress-induced duplex destabilization in long genomic DNA sequences.

We present a method for calculating predicted locations and extents of stress-induced DNA duplex destabilization (SIDD) as functions of base sequence and stress level in long DNA molecules. The base pair denaturation energies are assigned individually, so the influences of near neighbors, methylated bases, adducts, or lesions can be included. Sample calculations indicate that copolymeric energetics give results that are close to those derived when full near-neighbor energetics are used; small but potentially informative differences occur only in the calculated SIDD properties of moderately destabilized regions. The method presented here for analyzing long sequences calculates the destabilization properties within windows of fixed length N, with successive windows displaced by an offset distance d(o). The final values of the relevant destabilization parameters for each base pair are calculated as weighted averages of the values computed for each window in which that base pair appears. This approach implicitly assumes that the strength of the direct coupling between remote base pairs that is induced by the imposed stress attenuates with their separation distance. This strategy enables calculations of the destabilization properties of DNA sequences of any length, up to and including complete chromosomes. We illustrate its utility by calculating the destabilization properties of the entire E. coli genomic DNA sequence. A preliminary analysis of the results shows that promoters are associated with SIDD regions in a highly statistically significant manner, suggesting that SIDD attributes may prove useful in the computational prediction of promoter locations in prokaryotes.     

The Functional Role of S/MARs in Episomal Vectors as Defined by the Stress-Induced Destabilization Profile of the Vector Sequences

The scaffold/matrix attachment regions (S/MARs) are chromosomal elements that participate in the formation of chromatin domains and have origin of replication support functions. Because of all these functions, in recent years, they have been used as part of episomal vectors for gene transfer. The S/MAR of the human β-interferon gene has been shown to support efficient episome retention and transgene expression in various mammalian cells. In Jurkat and other cells, DNA plasmid vectors containing Epstein-Barr virus origin of replication (EBV OriP) and the EBV nuclear antigen-1 gene mediate prolonged episome retention in the host cell nucleus, which, however, diminishes over time. In order to enhance retention, we combined this system with an S/MAR element. Unexpectedly, this completely eliminated the capacity of episomes to replicate. Calculation of the stress-induced DNA duplex destabilization profile of the vectors suggested that the S/MAR element had created an increase in molecular stability at the OriP site that may have disturbed replicative potential. In contrast, introduction of an alternative initiation of replication region from the β-globin gene locus, instead of EBV OriP and the EBV nuclear antigen-1 gene, restored replicative capacity and enhanced episome retention mediated by the S/MAR. These effects were associated with a destabilization profile at the initiation of replication region. These data demonstrate a correlation between S/MAR-mediated vector retention and the presence of an unstable duplex at a replication origin, in this particular setting. We consider that the calculation of stress-induced duplex destabilization may be an informative first step in the design of units that replicate extrachromosomally, particularly as the latter present a safer and, therefore, attractive alternative to integrating viral vectors for gene therapy applications.     

S/MAR与基因表达

在真核生物的细胞核内,基因组是通过ＤＮＡ的核骨架附着（ＳＡＲ）或称核基质附着区（ＭＡＲ）（简记为Ｓ／ＭＡＲ）锚定在核骨架网状系统上的．Ｓ／ＭＡＲ既有一定的特征,又有多样性,研究认为它参与了ＤＮＡ复制调控和转录调控等多种核内生化过程,通过重组,在目的基因一侧或两侧带上Ｓ／ＭＡＲ后作基因转染或基因动植物,发现整合后的基因表达有时可增强几倍,甚至上万倍和／或显示位置独立效应,有些研究还报道,Ｓ／ＭＡＲ能     

Susceptibility to superhelically driven DNA duplex destabilization: a highly conserved property of yeast replication origins

Strand separation is obligatory for several DNA functions, including replication. However, local DNA properties such as A+T content or thermodynamic stability alone do not determine the susceptibility to this transition in vivo. Rather, superhelical stresses provide long-range coupling among the transition behaviors of all base pairs within a topologically constrained domain. We have developed methods to analyze superhelically induced duplex destabilization (SIDD) in genomic DNA that take into account both this long-range stress-induced coupling and sequence-dependent local thermodynamic stability. Here we apply this approach to examine the SIDD properties of 39 experimentally well-characterized autonomously replicating DNA sequences (ARS elements), which function as replication origins in the yeast Saccharomyces cerevisiae. We find that these ARS elements have a strikingly increased susceptibility to SIDD relative to their surrounding sequences. On average, these ARS elements require 4.78 kcal/mol less free energy to separate than do their immediately surrounding sequences, making them more than 2,000 times easier to open. Statistical analysis shows that the probability of this strong an association between SIDD sites and ARS elements arising by chance is approximately 4 × 10⁻¹⁰. This local enhancement of the propensity to separate to single strands under superhelical stress has obvious implications for origin function. SIDD properties also could be used, in conjunction with other known origin attributes, to identify putative replication origins in yeast, and possibly in other metazoan genomes.     

AT-rich islands in genomic DNA as a novel target for AT-specific DNA-reactive antitumor drugs

Interstrand cross-links at T(A/T)4A sites in cellular DNA are associated with hypercytotoxicity of an anticancer drug, bizelesin. Here we evaluated whether these lethal effects reflect targeting critical genomic regions. An in silico analysis of human sequences showed that T(A/T)4A motifs are on average scarce and scattered. However, significantly higher local motif densities were identified in distinct minisatellite regions (200-1000 base pairs of approximately 85-100% AT), herein referred to as "AT islands." Experimentally detected bizelesin lesions agree with these in silico predictions. Actual bizelesin adducts clustered within the model AT island naked DNA, whereas motif-poor sequences were only sparsely adducted. In cancer cells, bizelesin produced high levels of lesions (approximately 4.7-7.1 lesions/kilobase pair/microM drug) in several prominent AT islands, compared with markedly lower lesion levels in several motif-poor loci and in bulk cellular DNA (approximately 0.8-1.3 and approximately 0.9 lesions/kilobase pair/microM drug, respectively). The identified AT islands exhibit sequence attributes of matrix attachment regions (MARs), domains that organize DNA loops on the nuclear matrix. The computed "MAR potential" and propensity for supercoiling-induced duplex destabilization (both predictive of strong MARs) correlate with the total number of bizelesin binding sites. Hence, MAR-like AT-rich non-coding domains can be regarded as a novel class of critical targets for anticancer drugs.     

The positive aspects of stress: strain initiates domain decondensation (SIDD).

The conventional string-based bioinformatic methods of genomic sequence analysis are often insufficient to identify DNA regulatory elements, since many of these do not have a recognizable motif. Even in case a sequence pattern is known to be associated with an element it may only partially mediate its function. This suggests that properties not correlated with the details of base sequence contribute to regulation. One of these attributes is the DNA strand-separation potential, known as SIDD (stress-induced duplex destabilization) which facilitates the access of tracking proteins and the formation of local secondary structures. Using the type 1 interferon gene cluster as a paradigm, we demonstrate that the imprints in a SIDD profile coincide with chromatin domain borders and with DNAse I hypersensitive sites to which regulatory potential could be assigned. The approach permits the computer-guided identification of yet unknown, mostly remote sites and the design of artificial elements with predictable properties for multiple applications.     

Structure and Functions of Nuclear Matrix Associated Regions (S/MARs)

Modern concepts on the chromatin loop–domain organization and the role of the DNA regions specifically binding the nuclear matrix or nuclear scaffold (S/MARs) during its formation, maintenance, and regulation are discussed. Some S/MAR structural features, properties of binding the nuclear matrix, and probable mechanisms of their involvement in the gene regulation of activity are considered.     

Genomic imprinting controls matrix attachment regions in the Igf2 gene

Genomic imprinting at the Igf2/H19 locus originates from allele-specific DNA methylation, which modifies the affinity of some proteins for their target sequences. Here, we show that AT-rich DNA sequences located in the vicinity of previously characterized differentially methylated regions (DMRs) of the imprinted Igf2 gene are conserved between mouse and human. These sequences have all the characteristics of matrix attachment regions (MARs), which are known as versatile regulatory elements involved in chromatin structure and gene expression. Combining allele-specific nuclear matrix binding assays and real-time PCR quantification, we show that retention of two of these Igf2 MARs (MAR0 and MAR2) in the nuclear matrix fraction depends on the tissue and is specific to the paternal allele. Furthermore, on this allele, the Igf2 MAR2 is functionally linked to the neighboring DMR2 while, on the maternal allele, it is controlled by the imprinting-control region. Our work clearly demonstrates that genomic imprinting controls matrix attachment regions in the Igf2 gene.     

A matrix/scaffold attachment region binding protein: identification, purification, and mode of binding.

Matrix/scaffold attachment regions (MARs/SARs) partition chromatin into functional loop domains. Here we have identified a chicken protein that selectively binds to MARs from the chicken lysozyme locus and to MARs from Drosophila, mouse, and human genes. This protein, named ARBP (for attachment region binding protein), was purified to homogeneity and shown to bind to MARs in a cooperative fashion. ARBP is an abundant nuclear protein and a component of the internal nuclear network. Deletion mutants indicate that multiple AT-rich sequences, if contained in a minimal approximately 350 bp MAR fragment, can lead to efficient binding of ARBP. Furthermore, dimerization mutants show that, to bind ARBP efficiently, MAR sequences can act synergistically over large distances, apparently with the intervening DNA looping out. The binding characteristics of ARBP to MARs reproduce those of unfractionated matrix preparations, suggesting that ARBP is an important nuclear element for the generation of functional chromatin loops.     

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.     

Comparative study and prediction of DNA fragments associated with various elements of the nuclear matrix

Scaffold/matrix-associated region (S/MAR) sequences are DNA regions that are attached to the nuclear matrix, and participate in many cellular processes. The nuclear matrix is a complex structure consisting of various elements. In this paper we compared frequencies of simple nucleotide motifs in S/MAR sequences and in sequences extracted directly from various nuclear matrix elements, such as nuclear lamina, cores of rosette-like structures, synaptonemal complex. Multivariate linear discriminant analysis revealed significant differences between these sequences. Based on this result we have developed a program, ChrClass (Win/NT version, ftp.bionet.nsc.ru/pub/biology/chrclass/chrclass.zip), for the prediction of the regions associated with various elements of the nuclear matrix in a query sequence. Subsequently, several test samples were analyzed by using two S/MAR prediction programs (a ChrClass and MAR-Finder) and a simple MRS criterion (S/MAR recognition signature) indicating the presence of S/MARs. Some overlap between the predictions of all MAR prediction tools has been found. Simultaneous use of the ChrClass, MRS criterion and MAR-Finder programs may help to obtain a more clearcut picture of S/MAR distribution in a query sequence. In general, our results suggest that the proportion of missed S/MARs is lower for ChrClass, whereas the proportion of wrong S/MARs is lower for MAR-Finder and MRS.