Secondary structure prediction and structure-specific sequence analysis of single-stranded DNA.

DNA sequence analysis by oligonucleotide binding is often affected by interference with the secondary structure of the target DNA. Here we describe an approach that improves DNA secondary structure prediction by combining enzymatic probing of DNA by structure-specific 5'-nucleases with an energy minimization algorithm that utilizes the 5'-nuclease cleavage sites as constraints. The method can identify structural differences between two DNA molecules caused by minor sequence variations such as a single nucleotide mutation. It also demonstrates the existence of long-range interactions between DNA regions separated by >300 nt and the formation of multiple alternative structures by a 244 nt DNA molecule. The differences in the secondary structure of DNA molecules revealed by 5'-nuclease probing were used to design structure-specific probes for mutation discrimination that target the regions of structural, rather than sequence, differences. We also demonstrate the performance of structure-specific 'bridge' probes complementary to non-contiguous regions of the target molecule. The structure-specific probes do not require the high stringency binding conditions necessary for methods based on mismatch formation and permit mutation detection at temperatures from 4 to 37 degrees C. Structure-specific sequence analysis is applied for mutation detection in the Mycobacterium tuberculosis katG gene and for genotyping of the hepatitis C virus.

An invasive cleavage assay for direct quantitation of specific RNAs.

RNA quantitation is becoming increasingly important in basic, pharmaceutical, and clinical research. For example, quantitation of viral RNAs can predict disease progression and therapeutic efficacy. Likewise, gene expression analysis of diseased versus normal, or untreated versus treated, tissue can identify relevant biological responses or assess the effects of pharmacological agents. As the focus of the Human Genome Project moves toward gene expression analysis, the field will require a flexible RNA analysis technology that can quantitatively monitor multiple forms of alternatively transcribed and/or processed RNAs (refs 3,4). We have applied the principles of invasive cleavage and engineered an improved 5'-nuclease to develop an isothermal, fluorescence resonance energy transfer (FRET)-based signal amplification method for detecting RNA in both total RNA and cell lysate samples. This detection format, termed the RNA invasive cleavage assay, obviates the need for target amplification or additional enzymatic signal enhancement. In this report, we describe the assay and present data demonstrating its capabilities for sensitive (<100 copies per reaction), specific (discrimination of 95% homologous sequences, 1 in > or =20,000), and quantitative (1.2-fold changes in RNA levels) detection of unamplified RNA in both single- and biplex-reaction formats.

Mapping of RNA accessible sites by extension of random oligonucleotide libraries with reverse transcriptase.

A rapid and simple method for determining accessible sites in RNA that is independent of the length of target RNA and does not require RNA labeling is described. In this method, target RNA is allowed to hybridize with sequence-randomized libraries of DNA oligonucleotides linked to a common tag sequence at their 5'-end. Annealed oligonucleotides are extended with reverse transcriptase and the extended products are then amplified by using PCR with a primer corresponding to the tag sequence and a second primer specific to the target RNA sequence. We used the combination of both the lengths of the RT-PCR products and the location of the binding site of the RNA-specific primer to determine which regions of the RNA molecules were RNA extendible sites, that is, sites available for oligonucleotide binding and extension. We then employed this reverse transcription with the random oligonucleotide libraries (RT-ROL) method to determine the accessible sites on four mRNA targets, human activated ras (ha-ras), human intercellular adhesion molecule-1 (ICAM-1), rabbit beta-globin, and human interferon-gamma (IFN-gamma). Our results were concordant with those of other researchers who had used RNase H cleavage or hybridization with arrays of oligonucleotides to identify accessible sites on some of these targets. Further, we found good correlation between sites when we compared the location of extendible sites identified by RT-ROL with hybridization sites of effective antisense oligonucleotides on ICAM-1 mRNA in antisense inhibition studies. Finally, we discuss the relationship between RNA extendible sites and RNA accessibility.

Experimental and theoretical analysis of the invasive signal amplification reaction.

The invasive signal amplification reaction is a sensitive method for single nucleotide polymorphism detection and quantitative determination of viral load and gene expression. The method requires the adjacent binding of upstream and downstream oligonucleotides to a target nucleic acid (either DNA or RNA) to form a specific substrate for the structure-specific 5' nucleases that cleave the downstream oligonucleotide to generate signal. By running the reaction at an elevated temperature, the downstream oligonucleotide cycles on and off the target leading to multiple cleavage events per target molecule without temperature cycling. We have examined the performance of the FEN1 enzymes from Archaeoglobus fulgidus and Methanococcus jannaschii and the DNA polymerase I homologues from Thermus aquaticus and Thermus thermophilus in the invasive signal amplification reaction. We find that the reaction has a distinct temperature optimum which increases with increasing length of the downstream oligonucleotide. Raising the concentration of either the downstream oligonucleotide or the enzyme increases the reaction rate. When the reaction is configured to cycle the upstream instead of the downstream oligonucleotide, only the FEN1 enzymes can support a high level of cleavage. To investigate the origin of the background signal generated during the invasive reaction, the cleavage rates for several nonspecific substrates that arise during the course of a reaction were measured and compared with the rate of the specific reaction. We find that the different 5' nuclease enzymes display a much greater variability in cleavage rates on the nonspecific substrates than on the specific substrate. The experimental data are compared with a theoretical model of the invasive signal amplification reaction.

Sensitive detection of DNA polymorphisms by the serial invasive signal amplification reaction.

The invasive signal amplification reaction has been previously developed for quantitative detection of nucleic acids and discrimination of single-nucleotide polymorphisms. Here we describe a method that couples two invasive reactions into a serial isothermal homogeneous assay using fluorescence resonance energy transfer detection. The serial version of the assay generates more than 10(7) reporter molecules for each molecule of target DNA in a 4-h reaction; this sensitivity, coupled with the exquisite specificity of the reaction, is sufficient for direct detection of less than 1,000 target molecules with no prior target amplification. Here we present a kinetic analysis of the parameters affecting signal and background generation in the serial invasive signal amplification reaction and describe a simple kinetic model of the assay. We demonstrate the ability of the assay to detect as few as 600 copies of the methylene tetrahydrofolate reductase gene in samples of human genomic DNA. We also demonstrate the ability of the assay to discriminate single base differences in this gene by using 20 ng of human genomic DNA.

RNA template-dependent 5' nuclease activity of Thermus aquaticus and Thermus thermophilus DNA polymerases.

DNA replication and repair require a specific mechanism to join the 3'- and 5'-ends of two strands to maintain DNA continuity. In order to understand the details of this process, we studied the activity of the 5' nucleases with substrates containing an RNA template strand. By comparing the eubacterial and archaeal 5' nucleases, we show that the polymerase domain of the eubacterial enzymes is critical for the activity of the 5' nuclease domain on RNA containing substrates. Analysis of the activity of chimeric enzymes between the DNA polymerases from Thermus aquaticus (TaqPol) and Thermus thermophilus (TthPol) reveals two regions, in the "thumb" and in the "palm" subdomains, critical for RNA-dependent 5' nuclease activity. There are two critical amino acids in those regions that are responsible for the high activity of TthPol on RNA containing substrates. Mutating glycine 418 and glutamic acid 507 of TaqPol to lysine and glutamine, respectively, increases its RNA-dependent 5' nuclease activity 4-10-fold. Furthermore, the RNA-dependent DNA polymerase activity is controlled by a completely different region of TaqPol and TthPol, and mutations in this region do not affect the 5' nuclease activity. The results presented here suggest a novel substrate binding mode of the eubacterial DNA polymerase enzymes, called a 5' nuclease mode, that is distinct from the polymerizing and editing modes described previously. The application of the enzymes with improved RNA-dependent 5' nuclease activity for RNA detection using the invasive signal amplification assay is discussed.

The kinetics of oligonucleotide replacements.

The formation of a duplex between two nucleic acid strands is restricted if one of the strands forms an intra- or intermolecular secondary structure. The formation of the new duplex requires the dissociation and replacement of the initial structure. To understand the mechanism of this type of kinetics we studied the replacement of a labeled DNA oligonucleotide probe bound to a complementary DNA target with an unlabeled probe of the same sequence. The replacement kinetics were measured using a gel-shift assay for 12, 14 and 16-nucleotide probes as a function of temperature and concentration of the unlabeled probe. The results demonstrate that the overall replacement rate is a combination of two kinetic pathways: dissociative and sequential displacement. The dissociative pathway occurs by the spontaneous dissociation of the initial duplex followed by association of the target and unlabeled probe. The sequential displacement pathway requires only the partial melting of the initial duplex to allow for the formation of a branched nucleation complex with the unlabeled probe, followed by the complete displacement of the labeled probe by migration of the branch point. The contribution from the dissociative pathway is predominant at temperatures close to the melting point of the labeled probe, whereas the contribution from the displacement pathway prevails at lower temperatures and when the concentration of the replacing unlabeled probe is high. The results show that at physiological conditions, duplex formation between a single-stranded oligonucleotide probe and a structured region of a target molecule occurs mainly by the sequential-displacement mechanism.

A comparison of eubacterial and archaeal structure-specific 5'-exonucleases.

The 5'-exonuclease domains of the DNA polymerase I proteins of Eubacteria and the FEN1 proteins of Eukarya and Archaea are members of a family of structure-specific 5'-exonucleases with similar function but limited sequence similarity. Their physiological role is to remove the displaced 5' strands created by DNA polymerase during displacement synthesis, thereby creating a substrate for DNA ligase. In this paper, we define the substrate requirements for the 5'-exonuclease enzymes from Thermus aquaticus, Thermus thermophilus, Archaeoglobus fulgidus, Pyrococcus furiosus, Methanococcus jannaschii, and Methanobacterium thermoautotrophicum. The optimal substrate of these enzymes resembles DNA undergoing strand displacement synthesis and consists of a bifurcated downstream duplex with a directly abutted upstream duplex that overlaps the downstream duplex by one base pair. That single base of overlap causes the enzymes to leave a nick after cleavage and to cleave several orders of magnitude faster than a substrate that lacks overlap. The downstream duplex needs to be 10 base pairs long or greater for most of the enzymes to cut efficiently. The upstream duplex needs to be only 2 or 3 base pairs long for most enzymes, and there appears to be interaction with the last base of the primer strand. Overall, the enzymes display very similar substrate specificities, despite their limited level of sequence similarity.

Electron microscopic studies of sequence-specific recognition of duplex DNA by different ligands.

The following ligands were used to study sequence specific recognition of duplex DNA by electron microscopic techniques: methyltransferases BspR1 and EcoR124 (recognition sequences GGCC and GAAN7RTCG, respectively), a biotinylated deoxyoligonucleotide 5'-CTCTCTCTCTCTCT-3' capable of forming triplex DNA, and PNA oligomer H-T10-LysNH2. For each ligand the best conditions for electron microscopic (EM) detection of stable specific complex formation were determined. It was demonstrated that EM allowed us to determine the position of the individual target site with an error of 15-20 bp, the relative affinities for individual target sites and kinetic parameters of the binding. These results open new possibilities for EM investigations of sequence-specific interactions with a wide range of other ligands of a similar nature. They also imply that a wide range of different sequences can be unambiguously and precisely mapped by EM and greatly extend the scope of EM applications for physical mapping of genomic DNA.

Photofootprinting of DNA triplexes.

We have used a photofootprinting assay to study intermolecular and intramolecular DNA triplexes. The assay is based on the fact that the DNA duplex is protected against photodamage (specifically, against the formation of the (6-4) pyrimidine photoproducts) within a triplex structure. We have shown that this is the case for PyPuPu (YRR) as well as PyPuPy (YRY) triplexes. Using the photofootprinting assay, we have studied the triplex formation under a variety of experimentally defined conditions. At acid pH, d(C)n.d(G)n.d(C)n and d(CT)n.d(GA)n.d(CT)n triplexes are detected by this method. The d(CT)n.d(GA)n.d(CT)n triplexes are additionally stabilized by divalent cations and spermidine. PyPuPu triplexes are pH-independent and are stabilized by divalent cations, such as Mg++ and Zn++. The effect depends on the type of cation and on the DNA sequence. The d(CT)n.d(GA)n.d(GA)n triplex is stabilized by Zn++, but not by Mg++, whereas the d(C)n.d(G)n.d(G)n triplex is stabilized by Mg++. In H-DNA, virtually the entire pyrimidine chain is protected against photodimerization, whereas only half of the pyrimidine chain participating in a triplex is protected in the CGG intramolecular triplex.

Protection against UV-induced pyrimidine dimerization in DNA by triplex formation.

Cyclobutane and [6-4]-pyrimidine dimers are major photoproducts of ultraviolet-irradiated DNA. The yield of these photoproducts is dependent on the sequence and structure of the DNA. By analysing the photofootprints of fragments produced by cleavage of the DNA chain near [6-4]-pyrimidine dimers, we show here that a homopurine-homopyrimidine insert (with either d(TC)x or d(C)n) in plasmid pUC19 is, as expected, a good target for UV-induced pyrimidine-dimer formation. But we find that dimerization is virtually completely suppressed when the pyrimidine oligonucleotides d(TC)y or d(C)m are added to DNA carrying d(TC)x- or d(C)n-containing inserts, respectively. This effect is dependent on the type of oligonucleotide used and is site-specific. The protection occurs under acidic conditions that favour the formation of intermolecular triplexes between the homopurine-homopyrimidine inserts and homologous oligopyrimidines. We therefore conclude that triplex formation effectively protects the DNA duplex from UV-induced damage (pyrimidine dimerization). This observation makes the photofootprinting assay a very promising method for studying intermolecular and intramolecular triplexes (H-form DNA) both in vitro and in vivo.

Energetics of the B-H transition in supercoiled DNA carrying d(CT)x.d(AG)x and d(C)n.d(G)n inserts.

We have studied the B-H transition in the d(AG)x inserts of varying length under superhelical stress. The new data and previously published results for the d(G)31 insert are treated within a phenomenological model of the B-H transition, making it possible to obtain, for the first time, the energy parameters of the B-H transition in the d(AG)x and d(G)n sequences.

An unusual DNA structure detected in a telomeric sequence under superhelical stress and at low pH.

Telomeric sequences of DNA, which are found at the ends of linear chromosomes, have been attracting attention as potential sites for the formation of unusual DNA structures. They consist of (GnTm) or (GnATm) motifs (n greater than or equal to m) and, in the single-stranded state, form hairpins stabilized by non-canonical G.G pairs. In the duplex state and under superhelical stress they exhibit hypersensitivity to SI nuclease which by analogy with homopurine-homopyrimidine sequences may reflect the formation of an unusual structure. To determine whether this is the case we have inserted into a plasmid the Tetrahymena telomeric motif (G4T2).(A2C4) and probed it by two-dimensional gel electrophoresis, chemical modification and oligonucleotide binding. Our data demonstrate that, under superhelical stress and at low pH, the insert does indeed adopt a novel DNA conformation. We have concluded that in this structure the C-rich strand forms a hairpin stabilized by non-Watson-Crick base pairs C.C+ and A.A+, whereas the G-rich strand remains unstructured. We term this new DNA structure the (C,A)-hairpin.

Chemical probing of homopurine-homopyrimidine mirror repeats in supercoiled DNA.

We have recently shown that under superhelical stress and/or acid pH the homopurine-homopyrimidine tracts conforming to the mirror symmetry (H palindromes) form a novel DNA structure, the H form. According to our model, the H form includes (1) a triplex formed by half of the purine strand and by the homopyrimidine hairpin and (2) the unstructured other half of the purine strand. We used four specially designed sequences to demonstrate that, in accordance with our model, the mirror symmetry is essential for facile formation of the H form detected by two-dimensional gel electrophoresis. Here we report that, under conditions favouring the H-form extrusion, guanines of the 3' half of the purine strand are protected against alkylation by dimethylsulphate, whereas adenines of the 5' half of the purine strand react with diethyl pyrocarbonate. These data indicate that the 3' half of the homopurine strand is within the triplex whereas the 5' half is unstructured, in full agreement with our model.

A stable complex between homopyrimidine oligomers and the homologous regions of duplex DNAs.

When plasmid DNA duplexes carrying the regular homopurine-homopyrimidine inserts (dGA)n, (dTC)n and (dG)n, (dC)n are preincubated with homologous labeled oligo(dPy) ((dTC)n and (dC)n respectively) at acid pH, the label co-electrophoreses with the duplex DNA. Thus, a very strong complex is formed. Complementary oligo(dPu) does not form a complex under these conditions. No binding is observed for oligo(dPy) with non-homologous inserts as well as with vector plasmids without inserts. The complex is formed equally well with linear, nicked or superhelical DNA. The complex is not detected at pH greater than 6. Complex formation leads to very little, if any, unwinding of the duplex. The observed complex appears to be the Py.Pu.Py triplex consisting of TAT and CGC base-triads with protonated cytosines. Two-dimensional gel electrophoresis patterns show that the presence of homologous oligo(dPy) destabilizes the formation of the H form.

DNA H form requires a homopurine-homopyrimidine mirror repeat.

Regular homopurine-homopyrimidine tracts, (dG-dA)n(dT-dC)n and (dG)n(dC)n, undergo a superhelix-induced, strongly pH-dependent, structural transition into a novel DNA conformation, the H form. We have suggested that the H form can arise in any homopurine-homopyrimidine mirror repeat (H palindrome). We have now tested this prediction using a tailored series of plasmids carrying the inserts AAGGGAGAAXGGGGTATAGGGGYAAGAGGGAA, where X and Y may be either A or G, and subject them to two-dimensional gel electrophoresis. In support of our hypothesis, the inserts exhibited facile transitions into the H form for X = Y = G, or X = Y = A, whereas the transition was much more difficult or impossible for the two non-palindromes (X = A, Y = G or X = G, Y = A). We present evidence that the H form is the structural basis for S1-nuclease hypersensitivity.

Structure of (dG)n.(dC)n under superhelical stress and acid pH.

We have recently shown that a (GA)n.(TC)n tract undergoes a sharp structural transition under superhelical stress (V.I. Lyamichev, S.M. Mirkin and M.D. Frank-Kamenetskii, J. Biomol. Struct. Dyn. 2,327 (1985]. Unlike the well studied transitions to the cruciform and to the Z form, this novel transition was strongly pH-dependent. We have found the (dG)n.(dC)n insert to undergo a pH-dependent structural transition similar to that of the (GA)n.(TC)n tract. These new data meet our earlier expectations and disagree with the data of D.E. Pulleyblank, D.B. Haniford and A.R. Morgan, Cell 42, 271 (1985). We conclude that a novel DNA structure (the H-form) is typical of homopurine-homopyrimidine mirror repeats (the H palindromes) under superhelical stress and/or acid pH. In the H-form the homopyrimidine strand forms a hairpin while half of the homopurine strand interacts with the hairpin forming a triplex, the other half of the homopurine strand being unstructured (V.I. Lyamichev, S.M. Mirkin and M.D. Frank-Kamenetskii, J. Biomol. Struct. Dyn. 2, 3, 667 (1986].

The energetics of the B-Z transition in DNA.

The paper deals with the energetics of the transition to left-handed Z form in DNA with an arbitrary base sequence. There is a brief outline of the statistical-mechanical model of the B-Z transition allowing for three possible states of each base pair. The parameters of the model can be determined by comparing the theory with experimental data for the B-Z transition in inserts with given sequences in circular DNA. The model contains six energy parameters, most of which have been determined before. In order to find the remaining parameters of the model and test its adequacy, a number of oligonucleotide sequences were synthesized and inserted into the pUC 19 plasmid. Two-dimensional gel electrophoresis was used to determine the superhelical density at which the inserts adopt the Z form. A statistical-mechanical treatment of these data yielded a complete set of six energy parameters for the B-Z transition. The theoretical assumption that the free energy of Z-form pairs does not depend on the type of adjacent pairs proved to be in agreement with the experimental data.

Structures of homopurine-homopyrimidine tract in superhelical DNA.

For homopurine-homopyrimidine tracts in superhelical DNA, we propose a structure involving Watson-Crick and Hoogsteen paired triple helixes, hairpin loops and unstructured domains. Topologically, the whole structure is equivalent to an open region. The proposed structure is consistent with available S1 cleavage, pH and alkylation data and energetics under superhelical stress; this new structure is a much more probable candidate than the one proposed by us recently (V.I. Lyamichev, S.M. Mirkin & M.D. Frank-Kamenetskii, J. Biomole. Str. Dyns 3, 327-338, 1985).

A pH-dependent structural transition in the homopurine-homopyrimidine tract in superhelical DNA.

We have inserted the 509-bp-long fragment of sea urchin P. miliaris histone gene spacer region into plasmid pUC19. The fragment contains the 60-bp-long homopurine-homopyrimidine tract that is known to be hypersensitive to the S1 endonuclease. Using two-dimensional gel electrophoresis we have observed a sharp structural transition in the insert with increasing DNA superhelicity. As in the cases of cruciform and Z form formation, the observed transition partly relaxes the superhelical stress. In contrast with the other two well documented transitions, the observed transition strongly depends on pH. At pH7 and above the transition occurs at negative superhelicities exceeding the physiological range (- sigma greater than 0.08). For pH6 the transition occurs at -sigma = 0.055, whereas for pH4.3 it takes place at -sigma = 0.001. A comprehensive analysis of the obtained data has made it possible to define the nature of the observed transition. We conclude that under superhelical stress or/and at low pH homopurinehomopyrimidine tracts adopt a novel spatial structure called the H form.