Sequence organization in Xenopus DNA studied by the electron microscope.

Xenopus laevis DNA was extracted from red blood cells and sheared to a mean length of 2780 nucleotides. The DNA was stripped of foldback-containing fragments and incubated to C₀t 10 (mol · s · l⁻¹), allowing most repetitive sequences to form duplex structures. Duplex-containing fragments were eluted from an hydroxylapatite column and visualized for electron microscopy by spreading from 57% formamide according to the modified Kleinschmidt technique of Davis et al. (1971). The mean length of the fragments observed was 2445 nucleotides. A total of 1700 DNA strands were photographed and studied. Less than 5% of the total strand length was in uninterpretable structures. Every molecule falling within the confines of the plates was included in the sample. Over 50% of the total strand length in the sample was found in structures bearing at least one interspersed repetitive sequence duplex terminated by four single-strand regions. The fraction of DNA present in duplex regions was almost exactly that predicted if the duplex regions represent all the interspersed middle repetitive sequence in the Xenopus genome. Direct measurement of visualized duplexes shows that the mean length of interspersed repetitive sequence elements in this genome is 345 nucleotides. Duplex length was shown to be independent of the length of the strands bearing the duplexes. These observations provide direct confirmation of the length of approximately 300 nucleotides indicated for interspersed repetitive sequences by earlier physical-chemical studies 011 Xenopus DNA. In strands carrying two duplexes terminated by single-strand regions the interduplex, or single-copy sequence element length could be measured. Sequence interspersion curves generated from these data are roughly consistent with those derived earlier from measurements of hydroxylapatite binding as a function of fragment length.     

Force-induced DNA slippage

DNA containing repetitive sequences displays richer dynamics than heterogeneous sequences. In the genome the number of repeat units of repetitive sequences, known as microsatellites, may vary during replication by DNA slippage and their expansion gives rise to serious disorders. We studied the mechanical properties of repetitive DNA using dynamic force spectroscopy and found striking differences compared with ordinary heterogeneous sequences. Repetitive sequences dissociate at lower forces and elongate above a certain threshold force. This yield force was found to be rate dependent. Following the rapid stretching of the DNA duplex, the applied force relaxes by stepwise elongation of this duplex. Conversely, contraction of the DNA duplex can be observed at low forces. The stepwise elongation and shortening is initiated by single slippage events, and single-molecule experiments might help to explain the molecular mechanisms of microsatellites formation. In addition to the biological importance, the remarkable properties of repetitive DNA can be useful for different nanomechanical applications.     

The paperclip triplex: understanding the role of apex residues in tight turns

In this study, we investigate the role of the apex nucleotides of the two turns found in the intramolecular "paperclip" type triplex DNA formed by 5'-TCTCTCCTCTCTAGAGAG-3'. Our previously published structure calculations show that residues C7-A18 form a hairpin turn via Watson-Crick basepairing and residues T1-C6 bind into the major groove of the hairpin via Hoogsteen basepairing resulting in a broad turn of the T1-T12 5'-pyrimidine section of the DNA. We find that only the C6C7/G18 apex triad (and not the T12A13/T1 apex triad) is required for intramolecular triplex formation, is base independent, and occurs whether the purine section is located at the 5' or 3' end of the sequence. NMR spectroscopy and molecular dynamics simulations are used to investigate a bimolecular complex (which retains only the C6C7/G18 apex) in which a pyrimidine strand 5'- TCTCTCCTCTCT-3' makes a broad fold stabilized by the purine strand 5'-AGAGAG-3' via Watson Crick pairing to the T8-T12 and Hoogsteen basepairing to T1-T5 of the pyrimidine strand. Interestingly, this investigation shows that this 5'-AGAGAG-3' oligo acts as a new kind of triplex forming oligonucleotide, and adds to the growing number of triplex forming oligonucleotides that may prove useful as therapeutic agents.     

Structural properties of hybrid triplex of polycation deoxyribonucleic S-methylthiourea (DNmt) strands with a complementary DNA strand, probed by nanosecond molecular dynamics

The replacement of phosphodiester linkages of the polyanion DNA with S-methylthiourea linkers provides the polycation deoxyribonucleic S-methylthiourea (DNmt). Molecular dynamics studies to 1,220 ps of the hybrid triplex formed from octameric DNmt strands d(Tmt)8 with a complementary DNA oligomer strand d(Ap)8 have been carried out with explicit water solvent and Na+Cl- counterions under periodic boundary conditions using the CHARMM force field and the Ewald summation method. The Watson-Crick and Hoogsteen hydrogen-bonding patterns of the A/T tracts remained intact without any structural restraints for triplex structures throughout the simulation. The duplex portion of the triplex structure equilibrated at a B-DNA conformation in terms of the helical rise and other helical parameters. The dynamic structures of the DNmt x DNA x DNmt triplex were determined by examining histograms from the last 800 ps of the dynamics run. These included the hydrogen-bonding pattern (sequence recognition), three-centered bifurcating occurrences, minor groove width variations, and bending of tracts for the hybrid triplex structures. Together with the Watson-Crick hydrogen-bondings, the strong Hoogsteen hydrogen-bondings, the partially maintained three-centered bifurcatings in the Watson-Crick pair, and the medium-strength three-centered bifurcatings in the Hoogsteen pair suggest that the hybrid triplex is energetically favorable as compared to a duplex with similar base stacking, van der Waals interactions, and helical parameters. This is in agreement with our previously reported thermodynamic study, in which only triplex structures were observed in solution. The bending angle measured between the local axis vectors of the first and last helical axis segments is about 20 degrees for the Watson-Crick portion of the averaged structure. Propeller twist (associated with three-centered hydrogen-bonding) up to -30 degrees, native to DNA AT base pairing, was also observed for the triplex structure. The sugar pseudorotation phase angles and the ring rotation angles for the DNA strand are within the C3'-endo domain and C2'-endo domain for the DNmt strand. Water spines are observed in both minor and major grooves throughout the dynamics run. The molecular dynamics simulations of the structural properties of DNmt x DNA x DNmt hybrid triplex is compared to the DNG x DNA x DNG hybrid triplex (In DNG the -O-(PO2-)-O- linkers in DNA is replaced by -NH-C(=N+H2)-NH-).     

Affinity generation of single-stranded DNA for dideoxy sequencing following the polymerase chain reaction

A method to rapidly generate single stranded DNA for dideoxy sequencing following the polymerase chain reaction is described. By incorporating biotin in one of the amplification primers, we are able to physically separate the two DNA strands produced in the polymerase chain reaction. After amplification, the mixture is passed through a column containing streptavidin agarose. The strand produced by the biotinylated primer is bound in this matrix. The unbiotinylated strand is eluted with 0.2 N NaOH and sequenced by the dideoxy method. This method was utilized to sequence mitochondrial DNA from crude genomic DNA and to determine the sequences of four clones containing human mitochondrial DNA as a test of its accuracy. The use of biotin-facilitated separation permitted us to amplify and sequence DNA samples in a single day.     

Molecular dynamics simulation study of DNA triplex formed by mixed sequences in solution

The unrestrained molecular dynamics simulation of the triple helical DNA with mix sequences d(GACTGGTGAC).d(CTGACCACTG)*d (GACTGGTGAC), using the particle mesh Ewald sum, is presented here. The Ewald summation method effectively eliminates the usualcut-of of the long range interactions and allowed us to evaluate the full effect of the electrostatic forces. The AMBER5.0 force field has been used during the simulation in solvent. The MD results support a dynamically stable model of DNA triplex over the entire length of the trajectory. The duplex structure assumes the conformation, which is very close to B-DNA. In mixed sequences the purine bases occurs in both strand of DNA duplex. The bases of third strand do not favor the Hoogsteen or/and reverse Hoogsteen type of Hydrogen bonding but they form hydrogen bonds with the bases of both the strand of DNA duplex. The orientation of the third strand is parallel to one of the strand of duplex and all nucleotides (C, A, G & T) show isomorphic behavior with respect to the DNA duplex. The conformation of all the three strands is almost same except few exceptions. Due to interaction of third strand the conformational change in the duplex structure and a finite amount of displacement in the W-C base pairs have been observed. The conformational variation of the back bone torsion angles and helicoidal parameters, groove widths have been discussed. The sequence dependent effects on local conformation, helicoidal and morphological structure, width of the grooves of DNA helix may have important implication for understanding the functional energetics and specificity of interactions of DNA and its triplexes with proteins, pharmaceutical agents and other ligands.     

三核苷酸双链重复序列扩展合成特性及其机理

简单重复序列在各种生物基因组中广泛存在,同分子进化、遗传多样性、分子标记和某些遗传性疾病等密切相关．本文以全部32种18 bp三核苷酸双链重复序列作为研究对象,对它们在聚合酶作用下的扩展合成进行了系统研究．探讨了反应温度、序列本身等对扩展效率和产物长度的影响．结果显示,几乎所有的序列都能扩展变长．但序列对其扩展效率有很大影响：GC含量较多的短链,尤其是一条链中同时有G和C的短链扩展效率较高;全由AC和GT组成的双链也较易扩展．短链的最适扩展温度与短链GC组成呈一定的正相关关系．琼脂糖凝胶电泳和聚丙烯酰胺凝胶电泳分析发现,大部分产物单一性良好,且产物分子质量的大小随反应时间线性增加;随着反应温度升高,产物差异性增大．最后分析了双链重复序列的“复制滑移”扩展机理,有望为进一步研究重复序列的分子进化和基因检测中的非特异性扩增等奠定基础．     

Sequence organization of the rat genome by electron microscopy.

The size and arrangement of repetitive and inverted repeat (foldback) sequences in rat DNA were studied by visualization of hybrid and heteroduplex structures in the electron microscope. The self-reassociation of repetitive sequence-bearing DNA strands often results in the formation of four-ended "H" structures, whose duplex regions equal the repetitive sequence length and can be measured in the electron microscope. In this way, it was determined that the average size of the class of numerous short repetitive sequences is 0.40 +/- 0.15 kbp. Heteroduplex structures were prepared between long whole DNA single strands and short repeat-sequence-bearing strands. The analysis of these structures confirms that the size of the repetitive sequences in 0.4 kbp on average. Length measurements between adjacent duplexes show that the average spacing between two interspersed repeats is at least 1.5-1.8 kbp. By examining 29.4-kbp single strands after brief renaturation, the size and distribution of foldback sequences were determined. There are 1.9 X 10(5) foldback apirs per rat genome, spaced an average of 9.7 kbp apart according to our measurement. Repetitive, inverted repeat and unique sequences are interspersed with each other in at least half the genome.     

Targeting individual subunits of the FokI restriction endonuclease to specific DNA strands

Many restriction endonucleases are dimers that act symmetrically at palindromic DNA sequences, with each active site cutting one strand. In contrast, FokI acts asymmetrically at a non-palindromic sequence, cutting ‘top’ and ‘bottom’ strands 9 and 13 nucleotides downstream of the site. FokI is a monomeric protein with one active site and a single monomer covers the entire recognition sequence. To cut both strands, the monomer at the site recruits a second monomer from solution, but it is not yet known which DNA strand is cut by the monomer bound to the site and which by the recruited monomer. In this work, mutants of FokI were used to show that the monomer bound to the site made the distal cut in the bottom strand, whilst the recruited monomer made in parallel the proximal cut in the top strand. Procedures were also established to direct FokI activity, either preferentially to the bottom strand or exclusively to the top strand. The latter extends the range of enzymes for nicking specified strands at specific sequences, and may facilitate further applications of FokI in gene targeting.     

Coding capacity of complementary DNA strands.

A Fortran computer algorithm has been used to analyze the nucleotide sequence of several structural genes. The analysis performed on both coding and complementary DNA strands shows that whereas open reading frames shorter than 100 codons are randomly distributed on both DNA strands, open reading frames longer than 100 codons ("virtual genes") are significantly more frequent on the complementary DNA strand than on the coding one. These "virtual genes" were further investigated by looking at intron sequences, splicing points, signal sequences and by analyzing gene mutations. On the basis of this analysis coding and complementary DNA strands of several eukaryotic structural genes cannot be distinguished. In particular we suggest that the complementary DNA strand of the human epsilon-globin gene might indeed code for a protein.     

Molecular Dynamics Simulation Study of DNA Triplex Formed by Mixed Sequences in Solution

Abstract

The unrestrained molecular dynamics simulation of the triple helical DNA with mix sequences d(GACTGGTGAC)?d(CTGACCACTG)*d (GACTGGTGAC), using the particle mesh Ewald sum, is presented here. The Ewald summation method effectively eliminates the usual “cut-off” of the long—range interactions and allowed us to evaluate the full effect of the electrostatic forces. The AMBER5.0 force field has been used during the simulation in solvent. The MD results support a dynamically stable model of DNA triplex over the entire length of the trajectory. The duplex structure assumes the conformation, which is very close to B-DNA. In mixed sequences the purine bases occurs in both strand of DNA duplex. The bases of third strand do not favor the Hoogsteen or/and reverse Hoogsteen type of Hydrogen bonding but they form hydrogen bonds with the bases of both the strand of DNA duplex. The orientation of the third strand is parallel to one of the strand of duplex and all nucleotides (C, A, G & T) show isomorphic behavior with respect to the DNA duplex. The conformation of all the three strands is almost same except few exceptions. Due to interaction of third strand the conformational change in the duplex structure and a finite amount of displacement in the W-C base pairs have been observed. The conformational variation of the back bone torsion angles and helicoidal parameters, groove widths have been discussed. The sequence—dependent effects on local conformation, helicoidal and morphological structure, width of the grooves of DNA helix may have important implication for understanding the functional energetics and specificity of interactions of DNA and its triplexes with proteins, pharmaceutical agents and other legends.     

Structural Properties of Hybrid Triplex of Polycation Deoxyribonucleic S-methylthiourea (DNmt) Strands with a Complementary DNA Strand,Probed by Nanosecond Molecular Dynamics

Abstract

The replacement of phosphodiester linkages of the polyanion DNA with S-methylthiourea linkers provides the polycation deoxyribonucleic S-methylthiourea (DNmt). Molecular dynamics studies to 1,220 ps of the hybrid triplex formed from octameric DNmt strands d(Tmt)₈ with a complementary DNA oligomer strand d(Ap)₈ have been carried out with explicit water solvent and Na counterions under periodic boundary conditions using the CHARMM force field and the Ewald summation method. The Watson-Crick and Hoogsteen hydrogen-bonding patterns of the A/T tracts remained intact without any structural restraints for triplex structures throughout the simulation. The duplex portion of the triplex structure equilibrated at a B-DNA conformation in terms of the helical rise and other helical parameters. The dynamic structures of the DNmt·DNA·DNmt triplex were determined by examining histograms from the last 800 ps of the dynamics run. These included the hydrogen-bonding pattern (sequence recognition), three-centered bifurcating occurrences, minor groove width variations, and bending of tracts for the hybrid triplex structures. Together with the Watson-Crick hydrogen-bondings, the strong Hoogsteen hydrogen-bondings, the partially maintained three-centered bifurcatings in the Watson-Crick pair, and the medium-strength three-centered bifurcatings in the Hoogsteen pair suggest that the hybrid triplex is energetically favorable as compared to a duplex with similar base stacking, van der Waals interactions, and helical parameters. This is in agreement with our previously reported thermody- namic study, in which only triplex structures were observed in solution. The bending angle measured between the local axis vectors of the first and last helical axis segments is about 20° for the Watson-Crick portion of the averaged structure. Propeller twist (associated with three-centered hydrogen-bonding) up to ?30°, native to DNA AT base pairing, was also observed for the triplex structure. The sugar pseudorotation phase angles and the ring rotation angles for the DNA strand are within the C3′-endo domain and C2′-endo domain for the DNmt strand. Water spines are observed in both minor and major grooves throughout the dynamics run. The molecular dynamics simulations of the structural properties of DNmt·DNA·DNmt hybrid triplex is compared to the DNG·DNA·DNG hybrid triplex (In DNG the -O-(PO₂-)-O- linkers in DNA is replaced by -NH-C(=N₂)-NH-).     

A nonlinear dynamic model of DNA with a sequence-dependent stacking term

No simple model exists that accurately describes the melting behavior and breathing dynamics of double-stranded DNA as a function of nucleotide sequence. This is especially true for homogenous and periodic DNA sequences, which exhibit large deviations in melting temperature from predictions made by additive thermodynamic contributions. Currently, no method exists for analysis of the DNA breathing dynamics of repeats and of highly G/C- or A/T-rich regions, even though such sequences are widespread in vertebrate genomes. Here, we extend the nonlinear Peyrard–Bishop–Dauxois (PBD) model of DNA to include a sequence-dependent stacking term, resulting in a model that can accurately describe the melting behavior of homogenous and periodic sequences. We collect melting data for several DNA oligos, and apply Monte Carlo simulations to establish force constants for the 10 dinucleotide steps (CG, CA, GC, AT, AG, AA, AC, TA, GG, TC). The experiments and numerical simulations confirm that the GG/CC dinucleotide stacking is remarkably unstable, compared with the stacking in GC/CG and CG/GC dinucleotide steps. The extended PBD model will facilitate thermodynamic and dynamic simulations of important genomic regions such as CpG islands and disease-related repeats.     

Sequence and properties of the human KB cell and mouse L cell D-loop regions of mitochondrial DNA.

The sequences of the displacement-loop (D-loop) regions of mitochondrial DNA (mtDNA) from mouse L cells and human KB cells have been determined and provide physical maps to aid in the identification of sequences involved in the regulation of replication and expression of mammalian mtDNA. Both D-loop regions are bounded by the genes for tRNAPhe and tRNAPro. This region contains the most highly divergent sequences in mtDNA with the exceptions of three small conserved sequence blocks near the 5' ends of D-loop strands, a 225 nucleotide conserved sequence block in the center of the D-loop strand template region, and a short sequence associated with the 3' ends of D-loop strands. A sequence similar to that associated with the 3' termini of D-loop strands overlaps one of the conserved sequence blocks near the 5' ends of D-loop strands. The large, central conserved sequence probably does not code for a protein since no open reading frames are discretely conserved. Numerous symmetric sequences and potential secondary structures exist in these sequences, but none appear to be clearly conserved between species.     

Exploration of long and short repetitive sequence relationships in the sea urchin genome.

Long and short repetitive sequences of sea urchin DNA were prepared by reassociation of 2000 nucleotide long fragments to Cot 4 and digestion with the single strand specific nuclease S1. The S1 resistant duplexes were separated into long repetitive and short repetitive fractions on Agarose A50. The extent of shared sequences was studied by reassociating a labeled preparation of short repetitive DNA with an excess of unlabeled long repetitive DNA. Less than 10% of the long repetitive DNA preparation was able to reassociate with the short repetitive DNA. Thus the long and short repetitive elements appear to be principally independent sequence classes in sea urchin DNA. Precisely reassociating repetitive DNA was prepared by four successive steps of reassociation and thermal chromatography on hydroxyapatite. This fraction (3% of the genome) was reassociated by itself or with a great excess of total sea urchin DNA. The thermal stability of the products was identical in both cases (Tm=81 degrees C), indicating that precisely repeated sequences do not have many imprecise copies in sea urchin DNA.     

Association of double-stranded DNA fragments into multistranded DNA structures.

We have previously observed that double-stranded DNA fragments containing a tract of the tandemly repeated sequence poly(CA). poly(TG) can associate in vitro to form stable complexes of low electrophoretic mobility, which are recognized with high specificity by proteins HMG1 and HMG2. The formation of such complexes has since been observed to depend on interactions of DNA with polypropylene surfaces, with the suggestion that the formation of low mobility complexes might be the result of strand dissociation followed by misaligned reassociation of the repetitive sequences. The data presented here show that at high ionic strength the interactions of DNA with polypropylene are sufficiently strong for DNA to remain bound to the polypropylene surface, which suggests that DNA might also be involved in interactions with hydrophobic molecules in vivo. Under such conditions, low-mobility complexes are found only in the material adsorbed to the polypropylene surface, and all DNA fragments are able to form low-mobility structures, whether or not they contain repetitive sequences. Preventing the separation of strands by ligating hairpin loop oligonucleotides at both ends of the fragments does not prevent the formation of low-mobility complexes. Our results suggest two different pathways for the formation of complexes. In the first, dissociation is followed by misaligned reassociation of repetitive sequences, yielding duplexes with single-stranded end regions that associate to form multimeric complexes. In the second, repetitive as well as nonrepetitive DNA molecules bound to polypropylene adopt a conformation with locally unwound regions, which allows interactions between neighboring duplexes adsorbed on the surface, resulting in the formation of low-mobility complexes.     

Combinatorics of peptide sextets encoded by a single microgene

Genetic information stored in DNA sequences is translated into protein by linking a triplet nucleotide sequence and an amino acid. Because the frames of the triplets can be configured in three ways, a total of six polypeptides, each with a different sequence, can be produced from a single double-stranded DNA molecule. We recently developed the MolCraft system [reviewed in K. Shiba, J. Mol. Catal. B 18 (2004) xxx], which enables us to make combinatorial polymers of three peptides translated from one strand of a double-stranded DNA molecule. To explore all the information that a single double-stranded DNA molecule encodes, we have now developed a new system, La-MolCraft, in which all six reading frames encoded by both strands are combinatorially polymerized using loop-mediated isothermal amplification of DNA (LAMP) [Nucl. Acids Res. 28 (2000) E63].     

Analysis of rat repetitive DNA sequences.

Parameters of repetitive sequence organization have been measured in the rat genome. Experiments using melting, hydroxylapatite binding, and single strand specific nuclease digestion have been used to measure the number, length, and arrangement of repeated DNA sequences. Renaturation and melting or S1 nuclease digestion of 1.0 kbp DNA fragment show about 20% of rat DNA sequences are 3000-fold repeated. Renatured duplexes from 4.0 kbp DNA fragments display two repetitive size fractions after nuclease digestion. About 60% of the repeated sequences are 0.2-0.4 kbp long while the remainder are longer than 1.5 kbp. The arrangement of the repeated sequences has been measured by hydroxylapatite fractionation of DNA fragments of varying lengths bearing a repeated sequence. Repeated DNA sequences are interspersed among 2.5 kbp long nonrepeated sequences throughout more than 70% of the rat genome. There are approximately 350 different 3000-fold short repeated sequences in the rat interspersed among 600,000 nonrepeated DNA sequences.