Left-handed Z-DNA and in vivo supercoil density in the Escherichia coli chromosome

Conformation-dependent DNA attraction

Understanding how DNA molecules interact with other biomolecules is related to how they utilize their functions and is therefore critical for understanding their structure-function relationships. For a long time, the existence of Z-form DNA (a left-handed double helical version of DNA, instead of the common right-handed B-form) has puzzled the scientists, and the definitive biological significance of Z-DNA has not yet been clarified. In this study, the effects of DNA conformation in DNA-DNA interactions are explored by molecular dynamics simulations. Using umbrella sampling, we find that for both B- and Z-form DNA, surrounding Mg(2+) ions always exert themselves to screen the Coulomb repulsion between DNA phosphates, resulting in very weak attractive force. On the contrary, a tight and stable bound state is discovered for Z-DNA in the presence of Mg(2+) or Na(+), benefiting from their hydrophobic nature. Based on the contact surface and a dewetting process analysis, a two-stage binding process of Z-DNA is outlined: two Z-DNA first attract each other through charge screening and Mg(2+) bridges to phosphate groups in the same way as that of B-DNA, after which hydrophobic contacts of the deoxyribose groups are formed via a dewetting effect, resulting in stable attraction between two Z-DNA molecules. The highlighted hydrophobic nature of Z-DNA interaction from the current study may help to understand the biological functions of Z-DNA in gene transcription.

Gel mobilities of linking-number topoisomers and their dependence on DNA helical repeat and elasticity

Agarose-gel electrophoresis has been used for more than thirty years to characterize the linking-number (Lk) distribution of closed-circular DNA molecules. Although the physical basis of this technique remains poorly understood, the gel-electrophoretic behavior of covalently closed DNAs has been used to determine the local unwinding of DNA by proteins and small-molecule ligands, characterize supercoiling-dependent conformational transitions in duplex DNA, and to measure helical-repeat changes due to shifts in temperature and ionic strength. Those results have been analyzed by assuming that the absolute mobility of a particular topoisomer is mainly a function of the integral number of superhelical turns, and thus a slowly varying function of plasmid molecular weight. In examining the mobilities of Lk topoisomers for a series of plasmids that differ incrementally in size over more than one helical turn, we found that the size-dependent agarose-gel mobility of individual topoisomers with identical values of Lk (but different values of the excess linking number, DeltaLk) vary dramatically over a duplex turn. Our results suggest that a simple semi-empirical relationship holds between the electrophoretic mobility of linking-number topoisomers and their average writhe in solution.

Methods to determine DNA structural alterations and genetic instability

Chromosomal DNA is a dynamic structure that can adopt a variety of non-canonical (i.e., non-B) conformations. In this regard, at least 10 different forms of non-B DNA conformations have been identified; many of them have been found to be mutagenic, and associated with human disease development. Despite the importance of non-B DNA structures in genetic instability and DNA metabolic processes, mechanisms by which instability occurs remain largely undefined. The purpose of this review is to summarize current methodologies that are used to address questions in the field of non-B DNA structure-induced genetic instability. Advantages and disadvantages of each method will be discussed. A focused effort to further elucidate the mechanisms of non-B DNA-induced genetic instability will lead to a better understanding of how these structure-forming sequences contribute to the development of human disease.

Comparative analyses of distributions and functions of Z-DNA in Arabidopsis and rice

Left-handed Z-DNA is an energetically unfavorable DNA structure that could form mostly under certain physiological conditions and was known to be involved in a number of cellular activities such as transcription regulation. We have compared the distributions and functions of Z-DNA in the genomes of Arabidopsis and rice, and observed that Z-DNA occurs in rice at least 9 times more often than in Arabidopsis; similar observations hold for other monocots and dicots. In addition, Z-DNA is significantly enriched in the coding regions of Arabidopsis, and in the high-GC-content regions of rice. Based on our analyses, we speculate that Z-DNA may play a role in regulating the expression of transcription factors, inhibitors, translation repressors, succinate dehydrogenases and glutathione-disulfide reductases in Arabidopsis, and it may affect the expression of vesicle and nucleosome genes and genes involved in alcohol transporter activity, stem cell maintenance, meristem development and reproductive structure development in rice.

Topological insulators inhibit diffusion of transcription-induced positive supercoils in the chromosome of Escherichia coli

The double helical nature of DNA implies that progression of transcription machinery that cannot rotate easily around the DNA axis creates waves of positive supercoils ahead of it and negative supercoils behind it. Using topological reporters that detect local variations in DNA supercoiling, we have characterized the diffusion of transcription-induced (TI) positive supercoils in plasmids or in the chromosome of wild type Escherichia coli cells. Transcription-induced positive supercoils were able to diffuse and affect local supercoiling several kilobases away from the site of origin. By testing the effect of various DNA sequences, these reporters enabled us to identify elements that impede supercoil diffusion, i.e. behave as topological insulators. All the elements tested correspond to DNA gyrase catalytic targets. These results correlate the ability of a DNA sequence to be cleaved by DNA gyrase with topological insulator activity. Implications of the asymmetry in supercoil diffusion for the control of DNA topology are discussed.

Distributions of Z-DNA and nuclear factor I in human chromosome 22: a model for coupled transcriptional regulation

An analysis of the human chromosome 22 genomic sequence shows that both Z-DNA forming regions (ZDRs) and promoter sites for nuclear factor-I (NFI) are correlated with the locations of known and predicted genes across the chromosome and accumulate around the transcriptional start sites of the known genes. Thus, the occurrence of Z-DNA across human genomic sequences mirrors that of a known eukaryotic transcription factor. In addition, 43 of the 383 fully annotated chromosomal genes have ZDRs within 2 nucleosomes upstream of strong NFIs. This suggests a distinct class of human genes that may potentially be transcriptionally regulated by a mechanism that couples Z-DNA with NFI activation, similar to the mechanism previously elucidated for the human colony stimulation factor-I promoter [Liu et al. (2001) Cell, 106, 309-318]. The results from this study will facilitate the design of experimental studies to test the generality of this mechanism for other genes in the cell.

Interaction of the Zalpha domain of human ADAR1 with a negatively supercoiled plasmid visualized by atomic force microscopy

Interest to the left-handed DNA conformation has been recently boosted by the findings that a number of proteins contain the Zalpha domain, which has been shown to specifically recognize Z-DNA. The biological function of Zalpha is presently unknown, but it has been suggested that it may specifically direct protein regions of Z-DNA induced by negative supercoiling in actively transcribing genes. Many studies, including a crystal structure in complex with Z-DNA, have focused on the human ADAR1 Zalpha domain in isolation. We have hypothesized that the recognition of a Z-DNA sequence by the Zalpha(ADAR1) domain is context specific, occurring under energetic conditions, which favor Z-DNA formation. To test this hypothesis, we have applied atomic force microscopy to image Zalpha(ADAR1) complexed with supercoiled plasmid DNAs. We have demonstrated that the Zalpha(ADAR1) binds specifically to Z-DNA and preferentially to d(CG)(n) inserts, which require less energy for Z-DNA induction compared to other sequences. A notable finding is that site-specific Zalpha binding to d(GC)(13) or d(GC)(2)C(GC)(10) inserts is observed when DNA supercoiling is insufficient to induce Z-DNA formation. These results indicate that Zalpha(ADAR1) binding facilities the B-to-Z transition and provides additional support to the model that Z-DNA binding proteins may regulate biological processes through structure-specific recognition.

The intrinsic structure and stability of out-of-alternation base pairs in Z-DNA

Alternating pyrimidine-purine sequences typically form Z-DNA, with the pyrimidines in the anti and purines in the syn conformations. The observation that dC and dT nucleotides can also adopt the syn conformation (i.e. the nucleotides are out-of-alternation) extends the range of sequences that can convert to this left-handed form of DNA. Here, we study the effects of placing two adjacent d(G*C) base pairs as opposed to a single d(G*C) base pair or two d(A*T) base pairs out-of-alternation by comparing the structure of d(m5CGGCm5CG)2with the previously published structures of d(m5CGGGm5CG)*d(m5CGCCm5CG) and d(m5CGATm5CG)2. A high buckle and loss of stacking interactions are observed as intrinsic properties of the out-of-alternation base pairs regardless of sequence and the context of the dinucleotide. From solution titrations, we find that the destabilizing effect of out-of-alternation d(G*C) base pairs are identical whether these base pairs are adjacent or isolated. We can therefore conclude that it is these intrinsic distortions in the structure of the base pairs and not neighboring effects that account for the inability of out-of-alternation base pairs to adopt the left-handed Z conformation.

Crystal structure of d(GCGCGCG) with 5'-overhang G residues

The crystal structure of the DNA heptamer d(GCGCGCG) has been solved at 1.65 A resolution by the molecular replacement method and refined to an R-value of 0.184 for 3598 reflections. The heptamer forms a Z-DNA d(CGCGCG)2 with 5'-overhang G residues instead of an A-DNA d(GCGCGC)2 with 3'-overhang G residues. The overhang G residues from parallel strands of two adjacent duplexes form a trans reverse Hoogsteen G x G basepair that stacks on the six Z-DNA basepairs to produce a pseudocontinuous helix. The reverse Hoogsteen G x G basepair is unusual in that the displacement of one G base relative to the other allows them to participate in a bifurcated (G1)N2 . . . N7(G8) and an enhanced (G8)C8-H . . . O6(G1) hydrogen bond, in addition to the two usual hydrogen bonds. The 5'-overhang G residues are anti and C2'-endo while the 3'-terminal G residues are syn and C2'-endo. The conformations of both G residues are different from the syn/C3'-endo for the guanosine in a standard Z-DNA. The two cobalt hexammine ions bind to the phosphate groups in both GpC and CpG steps in Z(I) and Z(II) conformations. The water structure motif is similar to the other Z-DNA structures.

Flexible DNA: genetically unstable CTG.CAG and CGG.CCG from human hereditary neuromuscular disease genes

The properties of duplex CTG.CAG and CGG.CCG, which are involved in the etiology of several hereditary neurodegenerative diseases, were investigated by a variety of methods, including circularization kinetics, apparent helical repeat determination, and polyacrylamide gel electrophoresis. The bending moduli were 1.13 x 10(-19) erg.cm for CTG and 1.27 x 10(-19) erg.cm for CGG, approximately 40% less than for random B-DNA. Also, the persistence lengths of the triplet repeat sequences were approximately 60% the value for random B-DNA. However, the torsional moduli and the helical repeats were 2.3 x 10(-19) erg.cm and 10.4 base pairs (bp)/turn for CTG and 2.4 x 10(-19) erg.cm and 10.3 bp/turn for CGG, respectively, all within the range for random B-DNA. Determination of the apparent helical repeat by the band shift assay indicated that the writhe of the repeats was different from that of random B-DNA. In addition, molecules of 224-245 bp in length (64-71 triplet repeats) were able to form topological isomers upon cyclization. The low bending moduli are consistent with predictions from crystallographic variations in slide, roll, and tilt. No unpaired bases or non-B-DNA structures could be detected by chemical and enzymatic probe analyses, two-dimensional agarose gel electrophoresis, and immunological studies. Hence, CTG and CGG are more flexible and highly writhed than random B-DNA and thus would be expected to act as sinks for the accumulation of superhelical density.

Base-base and deoxyribose-base stacking interactions in B-DNA and Z-DNA: a quantum-chemical study

Base-stacking interactions in canonical and crystal B-DNA and in Z-DNA steps are studied using the ab initio quantum-chemical method with inclusion of electron correlation. The stacking energies in canonical B-DNA base-pair steps vary from -9.5 kcal/mol (GG) to -13.2 kcal/mol (GC). The many-body nonadditivity term, although rather small in absolute value, influences the sequence dependence of stacking energy. The base-stacking energies calculated for CGC and a hypothetical TAT sequence in Z-configuration are similar to those in B-DNA. Comparison with older quantum-chemical studies shows that they do not provide even a qualitatively correct description of base stacking. We also evaluate the base-(deoxy)ribose stacking geometry that occurs in Z-DNA and in nucleotides linked by 2',5'-phosphodiester bonds. Although the molecular orbital analysis does not rule out the charge-transfer n-pi* interaction of the sugar 04' with the aromatic base, the base-sugar contact is stabilized by dispersion energy similar to that of stacked bases. The stabilization amounts to almost 4 kcal/mol and is thus comparable to that afforded by normal base-base stacking. This enhancement of the total stacking interaction could contribute to the propensity of short d(CG)n sequences to adopt the Z-conformation.

Cloning, characterization, and properties of seven triplet repeat DNA sequences

Several neuromuscular and neurodegenerative diseases are caused by genetically unstable triplet repeat sequences (CTG.CAG, CGG.CCG, or AAG.CTT) in or near the responsible genes. We implemented novel cloning strategies with chemically synthesized oligonucleotides to clone seven of the triplet repeat sequences (GTA.TAC, GAT.ATC, GTT.AAC, CAC.GTG, AGG.CCT, TCG.CGA, and AAG.CTT), and the adjoining paper (Ohshima, K., Kang, S., Larson, J. E., and Wells, R. D.(1996) J. Biol. Chem. 271, 16784-16791) describes studies on TTA.TAA. This approach in conjunction with in vivo expansion studies in Escherichia coli enabled the preparation of at least 81 plasmids containing the repeat sequences with lengths of approximately 16 up to 158 triplets in both orientations with varying extents of polymorphisms. The inserts were characterized by DNA sequencing as well as DNA polymerase pausings, two-dimensional agarose gel electrophoresis, and chemical probe analyses to evaluate the capacity to adopt negative supercoil induced non-B DNA conformations. AAG.CTT and AGG.CCT form intramolecular triplexes, and the other five repeat sequences do not form any previously characterized non-B structures. However, long tracts of TCG.CGA showed strong inhibition of DNA synthesis at specific loci in the repeats as seen in the cases of CTG.CAG and CGG.CCG (Kang, S., Ohshima, K., Shimizu, M., Amirhaeri, S., and Wells, R. D.(1995) J. Biol. Chem. 270, 27014-27021). This work along with other studies (Wells, R. D.(1996) J. Biol. Chem. 271, 2875-2878) on CTG.CAG, CGG.CCG, and TTA.TAA makes available long inserts of all 10 triplet repeat sequences for a variety of physical, molecular biological, genetic, and medical investigations. A model to explain the reduction in mRNA abundance in Friedreich's ataxia based on intermolecular triplex formation is proposed.

Analysis of left-handed Z-DNA formation in short d(CG)n sequences in Escherichia coli and Halobacterium halobium plasmids. Stabilization by increasing repeat length and DNA supercoiling but not salinity

To evaluate the relative importance of alternating d(CG) sequence length, DNA supercoiling, and salt in left-handed Z-DNA formation, plasmids containing short d(CG)n sequences (n = 3-17) with the capability of replicating in either Escherichia coli or the halophilic archaeum Halobacterium halobium were constructed. Z-DNA conformation in the d(CG)n sequences was assayed by (i) a band shift assay using the Z-DNA-specific Z22 monoclonal antibody (ZIBS assay); (ii) an S1 nuclease cleavage-primer extension assay to map B-Z junctions; and (iii) a BssHII restriction inhibition assay. Using the ZIBS assay on plasmids purified from E. coli, the transition from B-DNA to Z-DNA occurred from d(CG)4, to d(CG)5, with 20% of d(CG)4, and 90% of d(CG)5 in Z-DNA conformation. These findings were consistent with the results of S1 nuclease cleavage observed at B-Z junctions flanking d(CG)4 and d(CG)5 sequences. Resistance to BssHII restriction endonuclease digestion was observed only in supercoiled plasmids containing d(CG)8 or longer sequences, indicating that shorter d(CG)n sequences are in dynamic equilibrium between B- and Z-DNA conformations. When a plasmid containing d(CG)4, was isolated from a topA mutant of E. coli, it contained 25% greater linking deficiency and 40% greater Z-DNA conformation in the alternating d(CG) region. In plasmids purified from H. halobium, which showed 30% greater linking deficiency than from E. coli, 20-40% greater Z-DNA formation was found in d(CG)4-6 sequences. Surprisingly, no significant difference in Z-DNA formation could be detected in d(CG)3-17 sequences in plasmids from either E. coli or H. halobium in the NaCl concentration range of 0.1-4 M.

Cooperative relaxation of supercoils and periodic transcriptional initiation within polymerase batteries

Transcription and DNA supercoiling are known to be linked by a cause-effect relationship that operates in both directions. It is proposed here that this two-way relationship may be exploited by the E. coli genome to facilitate constitutive transcription of supercoil-sensitive genes by polymerase batteries made up of uniformly spaces RNA polymerase elongation complexes. Specifically, it is argued that (1) polymerases transcribing DNA in tandem cooperate to relax each other's transcription-driven positive supercoils; and (2) negative supercoils driven upstream by elongation complexes tend to be 'harnessed' and used to cooperatively (and periodically) initiate fresh transcription from promoters. Harnessing of transcription-driven negative supercoils is thought to be achieved through the erection of protein barriers to the rotational upstream propagation of supercoils from transcription events. The possible relevance of such cooperation amongst polymerases to the activation of transcription by DNA-binding protein factors is emphasized. Some testable predictions are made and implications are discussed.

An intramolecular triplex in the human gamma-globin 5'-flanking region is altered by point mutations associated with hereditary persistence of fetal hemoglobin

The properties of an intramolecular triplex formed in vitro at the 5'-flanking region of the human gamma-globin genes were studied by chemical and physical probes. Chemical modifications performed with osmium tetroxide, chloroacetaldehyde, and diethyl pyrocarbonate revealed the presence of non-paired nucleotides on the "coding strand" at positions -209 through -217. These reactivities were induced by negative supercoiling, low pH, and magnesium ions. Downstream point mutations associated with hereditary persistence of fetal hemoglobin (HPFH) altered the extent of the modifications and some of the patterns. Specifically, C-202-->G and C-202-->T significantly decreased the reactivities, whereas the patterns were increased and altered in the T-198-->C. C-196-->T and C-195-->G caused local decreases in reactivity. Modifications at the upstream flanking duplex were modulated by the composition of the vector sequence. In summary, our data indicates the formation of an intramolecular triplex between nucleotides -209 to -217 of the "non-coding strand" and the downstream sequence containing the HPFH mutations. All of the HPFH point mutations altered the structure. More than one sequence alignment is possible for each of the triplexes. In addition, a consequence of some of the point mutations may be to facilitate slippage of the third strand relative to the Watson-Crick duplex.

Accelerated hybridization of oligonucleotides to duplex DNA

We report two strategies for accelerating the hybridization of oligonucleotides to DNA. We demonstrate that oligodeoxyribonucleotides and peptide nucleic acid oligomers hybridize to inverted repeats within duplex DNA by D-loop formation. Oligonucleotides and duplex template form an active complex, which can be recognized by T7 DNA polymerase to prime polymerization. Quantitation of polymerization products allowed the rate of hybridization to be estimated, and peptide nucleic acid oligomers and oligonucleotide-protein adducts anneal with association constants 500- and 12,000-fold greater, respectively, than the analogous unmodified oligonucleotides. Together, these results indicate that sequences within duplex DNA can be targeted by Watson-Crick base pairing and that chemical modifications can dramatically enhance the rate of strand association. These findings should facilitate targeting of oligomers for priming DNA polymerization, the detection of diagnostic sequences, and the disruption of gene expression. The observed acceleration of hybridization may offer a new perspective on the ability of RecA or other proteins to accelerate strand invasion.

Differential effects of simple repeating DNA sequences on gene expression from the SV40 early promoter

The influence of simple repeat sequences, cloned into different positions relative to the SV40 early promoter/enhancer, on the transient expression of the chloramphenicol acetyltransferase (CAT) gene was investigated. Insertion of (G)29.(C)29 in either orientation into the 5'-untranslated region of the CAT gene reduced expression in CV-1 cells 50-100 fold when compared with controls with random sequence inserts. Analysis of CAT-specific mRNA levels demonstrated that the effect was due to a reduction of CAT mRNA production rather than to posttranscriptional events. In contrast, insertion of the same insert in either orientation upstream of the promoter-enhancer or downstream of the gene stimulated gene expression 2-3-fold. These effects could be reversed by cotransfection of a competitor plasmid carrying (G)25.(C)25 sequences. The results suggest that a G.C-binding transcription factor modulates gene expression in this system and that promoter strength can be regulated by providing protein-binding sites in trans. Although constructs containing longer tracts of alternating (C-G), (T-G), or (A-T) sequences inhibited CAT expression when inserted in the 5'-untranslated region of the CAT gene, the amount of CAT mRNA was unaffected. Hence, these inhibitions must be due to posttranscriptional events, presumably at the level of translation. These effects of microsatellite sequences on gene expression are discussed with respect to recent data on related simple repeat sequences which cause several human genetic diseases.