Unusual frequencies of certain alternating purine-pyrimidine runs in natural DNA sequences: relation to Z-DNA

Prokaryotic, eukaryotic and mitochondrial DNA sequences of total Length 300 000 nucleotides have been analyzed to find out whether stretches of alternating purines and pyrimidines are unusual in terms of occurrence, composition and base sequence. Alternating runs longer than 5 nucleotides are significantly under-represented in the natural sequences as compared to random ones. Octanucleotides are the most deficient, occurring at only 60% of the frequency expected in random sequences. An unexpectedly high proportion of these octamers consists of alternating tetramers with the repeat structure (PuPyPuPy)2 or (PyPuPyPu)2. DNA stretches containing such sequences can potentially form a S1 nuclease sensitive slippage (staggered loop) structure, which might serve as a locally unstacked intermediate in the B- to Z-DNA conformational transition.     

The structure of the vimentin gene

The structure of the chromosomal gene encoding the intermediate filament protein vimentin is described. This gene, which is present as a single copy in the hamster genome, comprises about 10 kb of DNA and contains more than 80% of intron sequences. S1 mapping and sequence analysis reveal nine exons with a total length of 1848 nucleotides. For the complete primary structure of hamster vimentin, 464 amino acids are predicted, giving a molecular weight of 53,500 daltons. The intron positions are at codons 186, 206/207, 238/239, 292/293, 334/335, 408, 423, and 451/452. The overall homology with chicken desmin is 60% and is even higher in the central (alpha-helical) regions of both molecules. Cross-hybridization at the DNA level, however, is low. Comparison of the amino acid sequence of vimentin with prekeratin sequences shows that there is lesser homology of primary structure, but both the position and size of alpha-helical regions are strongly conserved. At the 5' end of the gene there is a consensus promoter sequence. The first AUG start codon is found 132 nucleotides downstream of the estimated cap site. The 3' nontranslated sequence shows homologies with the chicken vimentin gene. An interesting feature of the vimentin gene is a stretch of 44 nucleotides of alternating dC and dA within intron 2 that may form left-handed Z-DNA.     

Low molecular weight RNAs with homologies to cellular DNA at sites of adenovirus DNA insertion in hamster or mouse cells.

The adenovirus type 2 (Ad2)-transformed hamster cell line HE5 contains one or very few integrated copies of Ad2 DNA. At the site of insertion of Ad2 DNA, the cellular DNA sequence has been completely preserved and has homologies to small unpolyadenylated, cytoplasmic RNAs of 300 nucleotides in length and to minority populations of smaller RNAs present in HE5 cells and in normal hamster cells. The 300-nucleotide RNA is present on average in approximately 20 copies per cell. This RNA, and shorter RNAs, reveal homologies to the hamster DNA sequence of approximately 400 nucleotides to the right of the site of insertion of Ad2 DNA, which is present in one or very few copies per genome. The nucleotide sequence of the DNA segment homologous to this RNA does not contain open reading frames in excess of a sequence encoding 18 amino acids. Thus, it is unlikely that the small RNAs are actually translated and their function is unknown. The nucleotide sequence does not exhibit similarities to known low mol. wt. RNAs of eukaryotic origin. The low mol. wt. cellular RNA has been found in HE5 cells, in other hamster cell lines and organs, and also in mouse cells. There are differences with respect to size and abundance in the RNAs smaller than 300 nucleotides between HE5 cells and LSH hamster embryo cells. The adenovirus type 12 (Ad12)-induced mouse tumor CBA-12-1-T carries greater than 30 copies of integrated Ad12 DNA. The cellular DNA sequence at the site of Ad12 DNA insertion exhibits homologies to small RNAs (approximately 300 nucleotides long) from mouse cells.(ABSTRACT TRUNCATED AT 250 WORDS)     

Nucleotide discrimination with DNA immobilized in the MspA nanopore

Nanopore sequencing has the potential to become a fast and low-cost DNA sequencing platform. An ionic current passing through a small pore would directly map the sequence of single stranded DNA (ssDNA) driven through the constriction. The pore protein, MspA, derived from Mycobacterium smegmatis, has a short and narrow channel constriction ideally suited for nanopore sequencing. To study MspA's ability to resolve nucleotides, we held ssDNA within the pore using a biotin-NeutrAvidin complex. We show that homopolymers of adenine, cytosine, thymine, and guanine in MspA exhibit much larger current differences than in α-hemolysin. Additionally, methylated cytosine is distinguishable from unmethylated cytosine. We establish that single nucleotide substitutions within homopolymer ssDNA can be detected when held in MspA's constriction. Using genomic single nucleotide polymorphisms, we demonstrate that single nucleotides within random DNA can be identified. Our results indicate that MspA has high signal-to-noise ratio and the single nucleotide sensitivity desired for nanopore sequencing devices.     

The distribution of restriction fragment lengths for non-overlapping restriction sites in a random DNA sequence model.

Overlapping subsequences in a DNA sequence are not independent even if independence is supposed for the single nucleotides. Therefore the often used geometric distribution for the length of restriction fragments is not exact. The exact distribution of this random variable is derived for non-overlapping restriction sites in a DNA sequence with an infinite (or very large) number of nucleotides. Correction to the finite case is easy. It is shown that the simple geometric distribution is a good approximation as long as the basic probability for the occurrence of the recognition sequence at a given site is small.     

De novo telomere addition by Tetrahymena telomerase in vitro.

Previous molecular genetic studies have shown that during programmed chromosomal healing, telomerase adds telomeric repeats directly to non-telomeric sequences in Tetrahymena, forming de novo telomeres. However, the biochemical mechanism underlying this process is not well understood. Here, we show for the first time that telomerase activity is capable in vitro of efficiently elongating completely non-telomeric DNA oligonucleotide primers, consisting of natural telomere-adjacent or random sequences, at low primer concentrations. Telomerase activity isolated from mated or vegetative cells had indistinguishable specificities for nontelomeric and telomeric primers. Consistent with in vivo results, the sequence GGGGT... was the predominant initial DNA sequence added by telomerase in vitro onto the 3' end of the non-telomeric primers. The 3' and 5' sequences of the primer both influenced the efficiency and pattern of de novo telomeric DNA addition. Priming of telomerase by double-stranded primers with overhangs of various lengths showed a requirement for a minimal 3' overhang of 20 nucleotides. With fully single-stranded non-telomeric primers, primer length up to approximately 30 nucleotides strongly affected the efficiency of telomeric DNA addition. We propose a model for the primer binding site of telomerase for non-telomeric primers to account for these length and structural requirements. We also propose that programmed de novo telomere addition in vivo is achieved through a hitherto undetected intrinsic ability of telomerase to elongate completely non-telomeric sequences.     

The inhibition of restriction endonucleases due to Z-DNA in negatively supercoiled plasmid.

Plasmid pGC20 containing the (dGC)9 insert in SmaI recognition site has been used to study the inhibition of cleavage by different restriction endonuclease due to Z-DNA formation in (dCG)10 sequence of the negatively supercoiled plasmid. Data obtained indicate the different sensitivity of restriction endonucleases to DNA conformational perturbations resulted from the Z-DNA formation. Therefore, the inhibition of DNA cleavage by a particular restriction endonuclease cannot serve as a criterion for the estimation of the length of B-Z junctions in circular supercoiled DNAs.     

Concordance of experimentally mapped or predicted Z-DNA sites with positions of selected alternating purine-pyrimidine tracts.

The recent electronmicroscopic and biochemical mapping of Z-DNA sites in phi X174, SV40, pBR322 and PM2 DNAs has been used to determine two sets of criteria for identification of potential Z-DNA sequences in natural DNA genomes. The prediction of potential Z-DNA tracts and corresponding statistical analysis of their occurrence have been made on a sample of 14 DNA genomes. Alternating purine and pyrimidine tracts longer than 5 base pairs in length and their clusters (quasi alternating fragments) in the 14 genomes studied are under-represented compared to the expectation from corresponding random sequences. The fragments [d(G X C)]n and [d(C X G)]n (n greater than or equal to 3) in general do not occur in circular DNA genomes and are under-represented in the linear DNAs of phages lambda and T7, whereas in linear genomes of adenoviruses they are strongly over-represented. With minor exceptions, potential Z-DNA sites are also under-represented compared to random sequences. In the 14 genomes studied, predicted Z-DNA tracts occur in non-coding as well as in protein coding regions. The predicted Z-DNA sites in phi X174, SV40, pBR322 and PM2 correspond well with those mapped experimentally. A complete listing together with a compact graphical representation of alternating purine-pyrimidine fragments and their Z-forming potential are presented.     

Crystal structure of a Z-DNA fragment containing thymine/2-aminoadenine base pairs

The Z-DNA structure has been shown to form in two crystals made from self-complementary DNA hexamers d(CGTDCG) and d(CDCGTG) which contain thymine/2-aminoadenine (TD) base pairs. The latter structure has been solved and refined to 1.3 A resolution and it shows only small conformational changes due to the introduction of the TD base pairs in comparison with the structure of d(CG)3. Spectroscopic studies with these compounds demonstrate that DNA molecules containing 2-aminoadenine residues form Z-DNA slightly more easily than do those containing adenine nucleotides, but not as readily as the parent sequence containing only guanine-cytosine base pairs.     

Inhibition of Mycobacterium smegmatis topoisomerase I by specific oligonucleotides

DNA topoisomerase I from Mycobacterium smegmatis unlike many other type I topoisomerases is a site specific DNA binding protein. We have investigated the sequence specific DNA binding characteristics of the enzyme using specific oligonucleotides of varied length. DNA binding, oligonucleotide competition and covalent complex assays show that the substrate length requirement for interaction is much longer ( approximately 20 nucleotides) in contrast to short length substrates (eight nucleotides) reported for Escherichia coli topoisomerase I and III. P1 nuclease and KMnO(4) footprinting experiments indicate a large protected region spanning about 20 nucleotides upstream and 2-3 nucleotides downstream of the cleavage site. Binding characteristics indicate that the enzyme interacts efficiently with both single-stranded and double-stranded substrates containing strong topoisomerase I sites (STS), a unique property not shared by any other type I topoisomerase. The oligonucleotides containing STS effectively inhibit the M. smegmatis topoisomerase I DNA relaxation activity.     

Topoisomer gel retardation: detection of anti-Z-DNA antibodies bound to Z-DNA within supercoiled DNA minicircles.

Small DNA fragments of approximately 350 bp in length, either with or without d(CG)n tracts, are ligated into underwound DNA minicircles to generate topoisomeric rings with different topological linking numbers, Lk. These minicircles, differing by an Lk of one, can be separated by acrylamide gel electrophoresis. Furthermore, electrophoresis can be used to reveal DNA double helix conformational changes that are induced by supercoiling, such as left-handed Z-DNA. When anti-Z-DNA antibodies are added to such minicircles, their binding leads to a selective retardation of the electrophoretic migration of the Z-DNA containing circles. This effect is not seen with relaxed minicircles and those with insufficient torsional stress to induce a conformational transition. Thus the technique of 'topoisomer gel retardation' presents a very sensitive assay for the identification of proteins that selectively bind to DNA conformations stabilized by negative DNA supercoiling.     

Patterns in spontaneous mutation revealed by human-baboon sequence comparison

We have analyzed the alignment of a long homologous region of the human and baboon genomes (approximately 1.5 Mb). We show that the frequency of gaps between aligned segments decreases slowly with gap length, indicating that several successive nucleotides are often deleted or inserted in one event. By contrast, runs of consecutive mismatches decrease rapidly in frequency with increasing length, following an exponential distribution, indicating that nucleotides are mostly substituted one at a time. Nucleotide substitutions are clumped at the scales of <10 and 1000-10,000 nucleotides, but show almost no aggregation at the scales of <10-100 and over approximately 50,000 nucleotides. Apparently, two rather different factors make the substitution rate not exactly uniform along the DNA sequence. Comparison of regions of very similar genomes that are approximately selectively neutral makes it possible to study spontaneous mutation at a new level of resolution.     

Crystal Structure of A Z-DNA Fragment Containing Thymine/2-Aminoadenine Base Pairs

Abstract

The Z-DNA structure has been shown to form in two crystals made from self-complementary DNA hexamers d(CGTDCG) and d(CDCGTG) which contain thymine/2-ammoadenine (TD) base pairs. The latter structure has been solved and refined to 1.3 Å resolution and it shows only small conformational changes due to the introduction of the TD base pairs in comparison with the structure of d(CG)3. Spectroscopic studies with these compounds demonstrate that DNA molecules containing 2-aminoadenine residues form Z-DNA slightly more easily than do those containing adenine nucleotides, but not as readily as the parent sequence containing only guanine-cytosine base pairs.     

A computer aided thermodynamic approach for predicting the formation of Z-DNA in naturally occurring sequences.

The ease with which a particular DNA segment adopts the left-handed Z-conformation depends largely on the sequence and on the degree of negative supercoiling to which it is subjected. We describe a computer program (Z-hunt) that is designed to search long sequences of naturally occurring DNA and retrieve those nucleotide combinations of up to 24 bp in length which show a strong propensity for Z-DNA formation. Incorporated into Z-hunt is a statistical mechanical model based on empirically determined energetic parameters for the B to Z transition accumulated to date. The Z-forming potential of a sequence is assessed by ranking its behavior as a function of negative superhelicity relative to the behavior of similar sized randomly generated nucleotide sequences assembled from over 80,000 combinations. The program makes it possible to compare directly the Z-forming potential of sequences with different base compositions and different sequence lengths. Using Z-hunt, we have analyzed the DNA sequences of the bacteriophage phi X174, plasmid pBR322, the animal virus SV40 and the replicative form of the eukaryotic adenovirus-2. The results are compared with those previously obtained by others from experiments designed to locate Z-DNA forming regions in these sequences using probes which show specificity for the left-handed DNA conformation.