Effects of bulge composition and flanking sequence on the kinking of DNA by bulged bases

Increased Flexibility between Stems of Intramolecular Three-Way Junctions by the Insertion of Bulges

Intramolecular junctions are a ubiquitous structure within DNA and RNA; three-way junctions in particular have high strain around the junction because of the lack of flexibility, preventing the junctions from adopting conformations that would allow for optimal folding. In this work, we used a combination of calorimetric and spectroscopic techniques to study the unfolding of four intramolecular three-way junctions. The control three-way junction, 3H, has the sequence d(GAAATTGCGCT5GCGCGTGCT5GCACAATTTC), which has three arms of different sequences. We studied three other three-way junctions in which one (2HS1H), two (HS12HS1), and three (HS1HS1HS1) cytosine bulges were placed at the junction to allow the arms to adopt a wider range of conformations that may potentially relieve strain. Through calorimetric studies, it was concluded that bulges produce only minor effects on the enthalpic and thermal stability at physiological salt concentrations for 2HS1H and HS1HS1HS1. HS12HS1 displays the strongest effect, with the GTGC stem lacking a defined transition. In addition to unfolding thermodynamics, the differential binding of counterions, water, and protons was determined. It was found that with each bulge, there was a large increase in the binding of counterions; this correlated with a decrease in the immobilization of structural water molecules. The increase in counterion uptake upon folding likely displaces binding of structural water, which is measured by the osmotic stress method, in favor of electrostricted waters. The cytosine bulges do not affect the binding of protons; this finding indicates that the bulges are not forming base-triplet stacks. These results indicate that bulges in junctions do not affect the unfolding profile or the enthalpy of oligonucleotides but do affect the number and amount of molecules immobilized by the junction.

Rational design of a quantitative, pH-insensitive, nucleic acid based fluorescent chloride reporter

Chloride plays a major role in cellular homeostasis by regulating the lumenal pH of intracellular organelles. We have described a pH-independent, fluorescent chloride reporter called Clensor that has successfully measured resting chloride in organelles of living cells. Here, we describe the rational design of Clensor. Clensor integrates a chloride sensitive fluorophore called 10,10'-bis[3-carboxypropyl]-9,9'-biacridinium dinitrate (BAC) with the programmability, modularity and targetability available to nucleic acid scaffolds. We show that simple conjugation of BAC to a DNA backbone fails to yield a viable chloride-sensitive reporter. Fluorescence intensity and lifetime investigations on a series of BAC-functionalized structural variants yielded molecular insights that guided the rational design and successful realization of the chloride sensitive fluorescent reporter, Clensor. This study provides some general design principles that would aid the realization of diverse ion-sensitive nucleic acid reporters based on the sensing strategy of Clensor.

Stabilization of hairpins and bulged secondary structures of nucleic acids by single incorporation of α,β-D-CNA featuring a gauche(+) alpha torsional angle

A constrained dinucleotide unit featuring a gauche(+) alpha torsional angle configuration was used to stabilize DNA hairpin or bulged structures.

Conformational trapping of mismatch recognition complex MSH2/MSH3 on repair-resistant DNA loops

Insertion and deletion of small heteroduplex loops are common mutations in DNA, but why some loops are prone to mutation and others are efficiently repaired is unknown. Here we report that the mismatch recognition complex, MSH2/MSH3, discriminates between a repair-competent and a repair-resistant loop by sensing the conformational dynamics of their junctions. MSH2/MSH3 binds, bends, and dissociates from repair-competent loops to signal downstream repair. Repair-resistant Cytosine-Adenine-Guanine (CAG) loops adopt a unique DNA junction that traps nucleotide-bound MSH2/MSH3, and inhibits its dissociation from the DNA. We envision that junction dynamics is an active participant and a conformational regulator of repair signaling, and governs whether a loop is removed by MSH2/MSH3 or escapes to become a precursor for mutation.

RNA helical imperfections regulate activation of the protein kinase PKR: effects of bulge position, size, and geometry

The protein kinase, PKR, is activated by long stretches of double-stranded (ds) RNA. Viruses often make long dsRNA elements with imperfections that still activate PKR. However, due to the complexity of the RNA structure, prediction of whether a given RNA is an activator of PKR is difficult. Herein, we systematically investigated how various RNA secondary structure defects contained within model dsRNA affect PKR activation. We find that bulges increasingly disfavor activation as they are moved toward the center of a duplex and as they are increased in size. Model RNAs designed to conform to cis, trans, or bent global geometries through strategic positioning of one or more bulges decreased activation of PKR relative to perfect dsRNA, although cis-bulged RNAs activated PKR much more potently than trans-bulged RNAs. Activation studies on bulge-containing chimeric duplexes support a model wherein PKR monomers interact adjacently, rather than through-space, for activation on bulged substrates. Last, unusually low ionic strength induced substantial increases in PKR activation in the presence of bulged RNAs suggesting that discrimination against bulges is higher under biological ionic strength conditions. Overall, this study provides a set of rules for understanding how secondary structural defects affect PKR activity.

Binding of a Flexibly-linked Dinuclear Ruthenium(II) Complex to Adenine-bulged DNA Duplexes

Using 1H NMR spectroscopy and molecular modelling, the DNA binding of a chiral dinuclear ruthenium(ii) complex {Δ,Δ-[{Ru(phen)2}2(μ-bb7)]4+; phen = 1,10-phenanthroline, bb7 = 1,7-bis[4(4′-methyl-2,2′-bipyridyl)]-heptane} involving a bridging ligand containing a flexible aliphatic chain has been studied. The binding of the ruthenium(ii) complex was examined with the non-self-complementary duplexes d(CCGAGAATCGGCC):d(GGCCGATTCCGG) (containing a single adenine bulge: designated SB) and d(CCGAGCCGTGCC):d(GGCACGAGCCGG) (containing two adenine bulge sites separated by two base-pairs: designated DB). The NMR data indicated that the ruthenium(ii) complex bound at the bulge site of SB, with one ruthenium centre located at the bulge site with the second metal centre binding with lower affinity and selectivity in the duplex region adjacent to the bulge site. Less specific binding is inferred from chemical shift changes of nucleotide protons two to five base pairs from the single adenine bulge. The ruthenium(ii) complex selectively bound the DB duplex with one metal centre located at each bulge site. The NMR results also suggested that the metal complex binding induced greater changes to the structure of the SB duplex, compared with the DB duplex. Modelling indicates the bridging ligand allowed each ruthenium(ii) metal centre to bind one adenine bulge of the doubly-bulged duplex without disrupting the DNA structure, using the additional torsional flexibility conferred by the aliphatic bridging ligand. However, the second ruthenium(ii) metal centre is not able to bind in the minor groove of the singly-bulged duplex without disrupting the structure, as the metal centre is too bulky. The results of this study suggest dinuclear ruthenium(ii) complexes have considerable potential as probes for DNA and RNA sequences that contain two bulge sites separated by a small number of base-pairs.

Energetic signatures of single base bulges: thermodynamic consequences and biological implications

DNA bulges are biologically consequential defects that can arise from template-primer misalignments during replication and pose challenges to the cellular DNA repair machinery. Calorimetric and spectroscopic characterizations of defect-containing duplexes reveal systematic patterns of sequence-context dependent bulge-induced destabilizations. These distinguishing energetic signatures are manifest in three coupled characteristics, namely: the magnitude of the bulge-induced duplex destabilization (DeltaDeltaG(Bulge)); the thermodynamic origins of DeltaDeltaG(Bulge) (i.e. enthalpic versus entropic); and, the cooperativity of the duplex melting transition (i.e. two-state versus non-two state). We find moderately destabilized duplexes undergo two-state dissociation and exhibit DeltaDeltaG(Bulge) values consistent with localized, nearest neighbor perturbations arising from unfavorable entropic contributions. Conversely, strongly destabilized duplexes melt in a non-two-state manner and exhibit DeltaDeltaG(Bulge) values consistent with perturbations exceeding nearest-neighbor expectations that are enthalpic in origin. Significantly, our data reveal an intriguing correlation in which the energetic impact of a single bulge base centered in one strand portends the impact of the corresponding complementary bulge base embedded in the opposite strand. We discuss potential correlations between these bulge-specific differential energetic profiles and their overall biological implications in terms of DNA recognition, repair and replication.

[Bent dsDNA with defined geometric characteristics in terms of complexes of bridged oligonucleotides]

An opportunity of designing nontypical double-stranded DNA structures containing nonnatural inserts in a regular nucleotide DNA sequence has been investigated. The looped nucleotide inserts on the basis of adenylates and thymidilates, and nonnucleotide inserts on the basis of phosphodiesters of diethyleneglycol, 1,10-decanediol, and 3-hydroxy-2-hudroxymethyltetrahydrofuran were introduced into the backbone of a 32-mer native DNA duplex. These inserts formed the internal loops in the modified double-stranded DNA fragments which were shown to lead to bending of the linear duplex structure by 16 to 83 degrees. The dependencies of the bend angle of dsDNA on the composition and the length of the looped regions were determined. It was established that the bend of the irregular region of dsDNA depended on the electrostatic interaction of the phosphate residues. The tension in the complex structure could be reduced by the introduction of additional nucleotide units opposite the loop, which led to some relaxation of the bent helix. The resulting parameters of the bend values were shown to be in a good agreement with the published data obtained by NMR spectroscopy. It was demonstrated that the variation of the nature or the length of the insert allowed one to regulate the level of the local perturbation of the duplex structure and, thereby, influence both the bend level of the double helix and the destabilization of the modified complex.

Analysis of branched nucleic acid structure using comparative gel electrophoresis

Electrophoresis in polyacrylamide gels provides a simple yet powerful means of analyzing the relative disposition of helical arms in branched nucleic acids. The electrophoretic mobility of DNA or RNA with a central discontinuity is determined by the angle subtended between the arms radiating from the branchpoint. In a multi-helical branchpoint, comparative gel electrophoresis can provide a relative measure of all the inter-helical angles and thus the shape and symmetry of the molecule. Using the long-short arm approach, the electrophoretic mobility of all the species with two helical arms that are longer than all others is compared. This can be done as a function of conditions, allowing the analysis of ion-dependent folding of branched DNA and RNA species. Notable successes for the technique include the four-way (Holliday) junction in DNA and helical junctions in functionally significant RNA species such as ribozymes. Many of these structures have subsequently been proved correct by crystallography or other methods, up to 10 years later in the case of the Holliday junction. Just as important, the technique has not failed to date. Comparative gel electrophoresis can provide a window on both fast and slow conformational equilibria such as conformer exchange in four-way DNA junctions. But perhaps the biggest test of the approach has been to deduce the structures of complexes of four-way DNA junctions with proteins. Two recent crystallographic structures show that the global structures were correctly deduced by electrophoresis, proving the worth of the method even in these rather complex systems. Comparative gel electrophoresis is a robust method for the analysis of branched nucleic acids and their complexes.

A new approach to varietal identification in plants by microsatellite high resolution melting analysis: application to the verification of grapevine and olive cultivars

BACKGROUND

Microsatellites are popular molecular markers in many plant species due to their stable and highly polymorphic nature. A number of analysis methods have been described but analyses of these markers are typically performed on cumbersome polyacrylamide gels or more conveniently by capillary electrophoresis on automated sequencers. However post-PCR handling steps are still required. High resolution melting can now combine detailed sequence analysis with the closed-tube benefits of real-time PCR and is described here as a novel way to verify the identity of plant varieties such as grapevine and olive.

RESULTS

DNA melting profiles for various plant variety and rootstock samples were compared to profiles for certified reference samples. Two closely related grapevine rootstocks differing by as little as a single di-nucleotide repeat could be rapidly differentiated while there was high reproducibility of melting profiles for identical cultivars.

CONCLUSION

This novel microsatellite analysis method allows high sample throughput with greatly reduced time to results for varietal certification and is amenable to other microsatellite analyses.

Nucleosome Positioning Determinants

A previous report demonstrated that one site in a nucleosome assembled onto a synthetic positioning sequence known as Fragment 67 is hypersensitive to permanganate. The site is required for positioning activity and is located 1.5 turns from the dyad, which is a region of high DNA curvature in the nucleosome. Here, the permanganate sensitivity of the nucleosome positioning Fragment 601 was examined in order to expand the dataset of nucleosome sequences containing KMnO(4) hypersensitive sites. The hyperreactive T residue in the six sites detected as well as the one in Fragment 67 and three in the 5 S rDNA positioning sequence were contained within a TA step. Seven of the ten sequences were of the form CTAGPuG or the related sequence TTAAPu. These motifs were also found in the binding sites of several transcriptional regulatory proteins that kink DNA. In order to assess the significance of these sites, the 10 bp positioning determinant in Fragment 67 was removed and replaced by the nine sequences from the 5 S rDNA and Fragment 601. The results demonstrated that these derivative fragments promoted high nucleosome stability and positioning as compared to a control sequence that contained an AT step in place of the TA step. The importance of the TA step was further tested by making single base-pair substitutions in Fragment 67 and the results revealed that stability and positioning activity followed the order: TA>TG>TT>/=TC approximately GG approximately GA approximately AT. Sequences flanking the TA step were also shown to be critical for nucleosome stability and positioning. Nucleosome positioning was restored to near wild-type levels with (CTG)(3), which can form slipped stranded structures and with one base bulges that kink DNA. The results of this study suggest that local DNA structures are important for positioning and that single base-pair changes at these sites could have profound effects on those genomic functions that depend on ordered nucleosomes.

Assessment of the sequence dependency for the binding of 2-aminonaphthyridine to the guanine bulge

We have studied the sequence dependent binding of 2-amino-1,8-naphthyridine derivative 1 to a single guanine bulge. The free energy changes for the binding to a guanine bulge with different sequence contexts (5'X_Y3'/3'X'GY'5') were determined by a curve fitting of the thermal denaturation profile of DNA in the presence and absence of 1. The data showed that (i) the binding of 1 to a guanine bulge is stronger for those flanking the G-C base pair than A-T base pair, (ii) the guanine 3' side to 1 in the complex is especially effective for the complex stabilization, and (iii) the increase of T(m) in the presence of 1 is not a good estimate for the sequence dependent binding. The most efficient 1-binding was observed for the sequence of G_G/CGC. Molecular modeling simulations suggested that stacking interaction between the 3' side guanine and 1 is the molecular basis for the strong binding to G_G/CGC.

Heteroduplex mobility assay for detection of new avian influenza virus variants

Highly pathogenic avian influenza (HPAI) in poultry causes high morbidity and mortality, and it is a List A disease of the Office International des Epizooties. An outbreak of HPAI in commercial poultry not only causes direct disease losses but often results in trade restrictions for the affected country. Because HPAI viruses can mutate from H5 and H7 low pathogenic avian influenza viruses, it is necessary to monitor and control even the low pathogenic form of the virus. We report a practical approach for screening large numbers of isolates that uses amplification by reverse transcriptase-polymerase chain reaction of a segment of the hemagglutinin (HA) gene (536-560 bp) of H7 avian influenza viruses followed by the heteroduplex mobility assay (HMA). The HMA test compares the amplified polymerase chain reaction product from unknown samples with reference isolates, which allows the identification of new variants. The HMA test results were compared with sequence analysis of the isolates used in the study. On the basis of the HMA, we could identify several new variant viruses present in the live bird markets in the northeastern United States. New strains gave a distinct pattern of bands in the gels in accordance with the different heteroduplexes formed when their HA region amplification products were incubated together with the same amplification product of a reference strain. These differences correlate with phylogenetic analysis from sequence data.

Distinguishing "looped-out" and "stacked-in" DNA bulge conformation using fluorescent 2-aminopurine replacing a purine base

The conformation of a bulged DNA base, whether looped-out of the DNA helix or stacked-in between the flanking bases, can be distinguished using fluorescence spectroscopy of an inserted fluorescent base. If 2-aminopurine, a structural analog of adenine and guanine, is placed in duplex DNA as the bulged base replacing an adenine or guanine, it loops out of the DNA helix into solution. This is determined by the decrease or increase of 2-aminopurine fluorescence during DNA thermomelting: if the 2-aminopurine base stacks into the helix, its fluorescence increases or remains about the same during DNA duplex melting, but if the 2-aminopurine base loops out of the helix, its fluorescence decreases upon melting of the DNA duplex.

Conformations of an adenine bulge in a DNA octamer and its influence on DNA structure from molecular dynamics simulations

Molecular dynamics simulations have been applied to the DNA octamer d(GCGCA-GAAC). d(GTTCGCGC), which has an adenine bulge at the center to determine the pathway for interconversion between the stacked and extended forms. These forms are known to be important in the molecular recognition of bulges. From a total of ~35 ns of simulation time with the most recent CHARMM27 force field a variety of distinct conformations and subconformations are found. Stacked and fully looped-out forms are in excellent agreement with experimental data from NMR and x-ray crystallography. Furthermore, in a number of conformations the bulge base associates with the minor groove to varying degrees. Transitions between many of the conformations are observed in the simulations and used to propose a complete transition pathway between the stacked and fully extended conformations. The effect on the surrounding DNA sequence is investigated and biological implications of the accessible conformational space and the suggested transition pathway are discussed, in particular for the interaction of the MS2 replicase operator RNA with its coat protein.

Differentiation of cassava-infecting begomoviruses using heteroduplex mobility assays

Cassava mosaic disease is an enormous problem affecting the viability and productivity of cassava in all the developing regions in Africa. The disease is now known to be caused by a number of cassava-infecting begomoviruses. New viruses and viral strains continue to be discovered and due to the lack of cloning and sequencing facilities in many African laboratories, a simple, rapid and sensitive technique is needed for screening of cassava plantations. Here we report on the development of a heteroduplex mobility assay (HMA) which could be used for the testing of viral-infected cassava. The assay involves amplifying the highly conserved core region of the coat protein gene of field isolates followed by denaturing and annealing with a number of reference strains. The HMA profiles in this study were able to differentiate four different viral species and 11 different virus strains, and showed a good correlation with sequencing results and phylogenetic comparisons with other sequenced cassava viruses. This technique is sensitive and rapid and has the added advantage of being able to detect mixtures of viruses in field-grown cassava.

Viral genome characterisation by the heteroduplex mobility and heteroduplex tracking assays

The heteroduplex mobility assay (HMA) is a means of comparing two PCR amplicons or, in the variation known as the heteroduplex tracking assay (HTA), a means of estimating the quasispecies diversity of a viral genome. Heteroduplex assays have many applications including subtyping viral genomes, screening for low frequency variants in a population, scanning the relative genetic diversity across a genome and screening for recombinant clones. They can be used to detect dual infections, superinfections, contaminated blood products and laboratory contaminations. PCR amplicons of about 65% sequence similarity or greater will form heteroduplexes under appropriate conditions, and phylogenetic trees can be drawn from heteroduplex mobility data. While homoduplexes indicate more than 98% similarity between two DNA sequences, heteroduplexes indicate at least seven mismatches in a 500-bp amplicon, or a three-base pair gap in 1000-bp. Minority variants comprising 1% to 5% of the genome population can be detected and quantified by HTA. Thus far, heteroduplex assays have been described for HIV and other lentiviruses, hepatitis C and G viruses, Norwalk-like viruses, influenza, measles and poliovirus. They could be applied to a wide range of other viral species.

Position and degree of mismatches and the mobility of DNA heteroduplexes

Heteroduplex mobility assay (HMA) is a fast and inexpensive method for determining relatedness between DNA sequences. Rapidly evolving viruses such as HIV-1 develop marked sequence differences in their genomes over the course of the epidemic and infection in a single individual. HMA can be used to monitor both processes. Here, we systematically evaluated the influence of single base mismatches on heteroduplex mobility. The impact of mismatches at nine different positions in 559 bp double-stranded DNA molecules, within a background of overall sequence divergence ranging from 1.97 to 9.65%, was evaluated in both non-denaturing and partially-denaturing acrylamide gels. We found that the electrophoretic mobility of heteroduplexes was proportional to the level of mismatch when that level exceeded 4.5%. Overall, mismatches near the center of the fragment and clustered mismatches tended to have an exaggerated influence on the mobility of heteroduplexes. Thus, the use of HMA for quantitative inference of genetic distances under the conditions we describe is of greatest utility at levels of mismatch >5%.

Heteroduplex mobility assay detects DNA mutations for differentiation of closely related phytoplasma strains

We present the first use of DNA heteroduplex mobility assay (HMA) to detect the point mutations including substitutions and deletions/insertions in 16S rDNA of aster yellows phytoplasma (AY27) and to differentiate phytoplasmas collected from field samples of clover proliferation (CP) and alfalfa witches'-broom (AWB). The phytoplasmal 16S rDNA fragment was amplified from AY27 by polymerase chain reaction (PCR) and cloned into a plasmid vector. The cloned DNA fragment was subjected to in vitro mutation to produce 1- to 4-base substitutions and 1- to 3-base deletions. The mutated 16S rDNA fragments were analyzed by HMA. The results showed that a single two-base substitution or a single-base deletion/insertion in the 529 bp DNA fragment was directly detected and that a DNA divergence at a level of as low as 0.2% was detectable by HMA. Heteroduplex mobilities were affected by the number and composition of the phytoplasma DNA bases in mismatches or gaps and were proportional to the degree of DNA divergences. Gaps caused greater retardation in heteroduplex mobility than mismatches did. HMA was highly sensitive in detecting the mixed infections of phytoplasmas. In analyses of CP and AWB field samples collected in Alberta, two CP and one AWB phytoplasma isolates were differentiated from others by HMA but not by restriction fragment length polymorphism (RFLP). Therefore, HMA provides a simple, rapid, highly sensitive and analytical method to detect and estimate the genetic divergence of phytoplasmas when other methods such as RFLP are not readily applicable.

Induced heteroduplex genotyping of TNF-alpha, IL-1beta, IL-6 and IL-10 polymorphisms associated with transcriptional regulation

We describe the construction and use of 7 induced heteroduplex generators, reagents for the rapid and unequivocal genotyping of nucleotide sequence polymorphism in TNF-alpha, IL-1beta, IL-6 and IL-10. Polymorphisms detected are those previously associated with regulation of gene transcription: TNF-alpha positions -308 and -238; IL-1beta position +3953; IL-6 position -174; and IL-10 positions -1082, -819 and -592. The reagents were used for analysis of allele and haplotype frequencies in a population of healthy Caucasian volunteer blood donors.

Methods for subtyping and molecular comparison of human viral genomes

The development over the past two decades of molecular methods for manipulation of RNA and DNA has afforded molecular virologists the ability to study viral genomes in detail that has heretofore not been possible. There are many molecular techniques now available for typing and subtyping of viruses. The available methods include restriction fragment length polymorphism analysis, Southern blot analysis, oligonucleotide fingerprint analysis, reverse hybridization, DNA enzyme immunoassay, RNase protection analysis, single-strand conformation polymorphism analysis, heteroduplex mobility assay, nucleotide sequencing, and genome segment length polymorphism analysis. The methods have certain advantages and disadvantages which should be considered in their application to specific viruses or for specific purposes. These techniques are likely to become more widely used in the future for epidemiologic studies and for investigations into the pathophysiology of virus infections.

Structure of branched DNA molecules: gel retardation and atomic force microscopy studies

DNA heteroduplexes as models for slipped strand DNA have been analyzed by polyacrylamide gel migration and atomic force microscopy (AFM). All heteroduplexes containing one hairpin or loop have reduced electrophoretic mobilities compared with that expected for their molecular weights. The retarded gel mobility correlates with the formation of a sharp kink detected by AFM. Increasing the hairpin length from 7 bp to 50 bp results in a monotonous decrease in gel mobility of heteroduplexes. This secondary retardation effect appears to depend only on the hairpin size since the AFM data show no dependence of the kink angle on the hairpin length. Heteroduplex isomers with a loop or hairpin in opposite strands migrate with distinct mobilities. Analysis of gel migration of heteroduplexes with altered hairpin orientations as well as of truncated heteroduplexes indicates that the difference in mobility is due to an inherent curvature in one of the long arms. This is confirmed by the end-to-end distance measurements from AFM images. In addition, significant variation of the end-to-end distances is consistent with a dynamic structure of heteroduplexes at the three-way junction. Double heteroduplexes containing one hairpin in each of the complementary strands also separate in a gel as two isomers. Their appearance in AFM showed a complicated pattern of flat representations of the three-dimensional structure and may indicate a certain degree of interaction between complementary parts of the hairpins that are several helical turns apart.

Structural analysis of slipped-strand DNA (S-DNA) formed in (CTG)n. (CAG)n repeats from the myotonic dystrophy locus

The mechanism of disease-associated trinucleotide repeat length variation may involve slippage of the triplet-containing strand at the replication fork, generating a slipped-strand DNA structure. We recently reported formation in vitro of slipped-strand DNA (S-DNA) structures when DNAs containing triplet repeat blocks of myotonic dystrophy or fragile X diseases were melted and allowed to reanneal to form duplexes. Here additional evidence is presented that is consistent with the existence of S-DNA structures. We demonstrate that S-DNA structures can form between two complementary strands containing equal numbers of repeats. In addition, we show that both the propensity for S-DNA formation and the structural complexity of S-DNAs formed increase with increasing repeat length. S-DNA structures were also analyzed by electron microscopy, confirming that the two strands are slipped out of register with respect to each other and confirming the structural polymorphism expected within long tracts of trinucleotide repeats. For (CTG)50.(CAG)50 two distinct populations of slipped structures have been identified: those involving </=10 repeats per slippage, which appear as bent/kinked DNA molecules, and those involving >10 repeats, which have multiple loops or hairpins indicative of complex alternative DNA secondary structures.

Analysis of the stability of looped-out and stacked-in conformations of an adenine bulge in DNA using a continuum model for solvent and ions

A combination of conformational search, energy minimization, and energetic evaluation using a continuum solvent treatment has been employed to study the stability of various conformations of the DNA fragment d(CGCAGAA)/d(TTCGCG) containing a single adenine bulge. The extra-helical (looped-out) bulge conformation derived from a published x-ray structure and intra-helical (stacked bulge base) model structures partially based on nuclear magnetic resonance (NMR) data were used as start structures for the conformational search. Solvent-dependent contributions to the stability of the conformations were calculated from the solvent exposed molecular surface area and by using the finite difference Poisson-Boltzmann approach. Three classes (I-III) of bulge conformations with calculated low energies can be distinguished. The lowest-energy conformations were found in class I, corresponding to structures with the bulge base stacked between flanking helices, and class II, composed of structures forming a triplet of the bulge base and a flanking base pair. All extra-helical bulge structures, forming class III, were found to be less stable compared with the lowest energy structures of class I and II. The results are consistent with NMR data on an adenine bulge in the same sequence context indicating an intra-helical or triplet bulge conformation in solution. Although the total energies and total electrostatic energies of the low-energy conformations show only relatively modest variations, the energetic contributions to the stability were found to vary significantly among the classes of bulge structures. All intra-helical bulge structures are stabilized by a more favorable Coulomb charge-charge interaction but destabilized by a larger electrostatic reaction field contribution compared with all extra-helical and most triplet bulge structures. Van der Waals packing interactions and nonpolar surface-area-dependent contributions appear to favor triplet class II structures and to a lesser degree also the intra-helical stacked bulge conformations. The large conformational variation found for class III conformers might add a favorable entropic contribution to the stability of the extra-helical bulge form.

Evolutionary variants of the human immunodeficiency virus type 1 V3 region characterized by using a heteroduplex tracking assay

Syncytium-inducing (SI) variants of human immunodeficiency virus type 1 (HIV-1) are evolutionary variants that are associated with rapid CD4+ cell loss and rapid disease progression. The heteroduplex tracking assay (HTA) was used to detect evolutionary V3 variants by amplifying the V3 sequences from viral RNA derived from 50 samples of patient plasma. For this V3-specific HTA (V3-HTA), heteroduplexes were formed between the patient V3 sequences and a probe with the subtype B consensus V3 sequence. Evolution was then measured by divergence from the consensus. The presence of evolutionary variants was correlated with SI detection data on the same samples from the MT-2 cell culture assay. Evolutionary variants were correlated with the SI phenotype in 88% of the samples, and 96% of the SI samples contained evolutionary variants. In most cases the evolutionary V3 variants represented discrete clonal outgrowths of virus. Sequence analysis of the six discordant samples that did not show this correlation indicated that three non-syncytium-inducing (NSI) samples had V3 sequences that had evolved away from the consensus sequence but not toward an SI genotype. A fourth sample showed little evolution away from the consensus but was SI, which indicates that not all SI variants require basic substitutions in V3. The other two samples had SI-like genotypes and NSI phenotypes, suggesting that V3-HTA was able to detect SI emergence in these samples in the absence of their detection in vitro. V3-HTA was also used to confirm SI variant selection in MT-2 cells and to examine the possibility of variant selection during virus culture in peripheral blood cells.

A practical approach to genetic screening for influenza virus variants

This report describes a quick genetic approach to the screening of influenza virus variants. Multiplex reverse transcription (MRT) and multiplex PCR (MPCR) were used to amplify and differentiate the hemagglutinin (HA) genes of different types and subtypes of influenza viruses. Heteroduplex mobility assay (HMA) was then used to differentiate strains within the same type and subtype. Three primers complementary to the consensus 3' termini of viral genomic RNA segments of human influenza virus types A, B, and C were used in a single MRT reaction to prime the synthesis of cDNA of all the viral genome segments. The cDNA was then amplified in an MPCR containing primers for the HA genes of the H1 and H3 subtypes of type A, the HA gene of type B, and the counterpart of type C virus. Amplicons of the different types and subtypes differ in size, thus allowing typing and subtyping. The regions amplified cover most of the HA1 portion of the HA genes and therefore amplicons of variants identified by the described HMA can be sequenced directly. In the HMA, the amplicon of an individual strain was mixed with that of a reference strain and heteroduplexes derived from mismatches migrated more slowly than homoduplexes of the same size in electrophoresis, with the mobility shift pattern indicating the divergence of amplicons. The whole process from viral RNA extraction to HMA can be completed within 2 days and thus provides a quick screening before further analysis by hemagglutination inhibition testing and sequencing. In addition, all segments of the viral genome can be amplified from a single MRT reaction, which can yield valuable sources of products for future genetic analyses. This method should facilitate genetic screening and characterization of influenza virus variants.

Methyl-directed repair of mismatched small heterologous sequences in cell extracts from Escherichia coli

The methyl-directed DNA repair efficiency of a set of M13mp18 heteroduplexes containing 1-8 or 22 unpaired bases was determined by using an in vitro DNA mismatch repair assay. The unpaired bases of each heteroduplex residing at overlapping recognition sites of two restriction endonucleases allow independent assay of repair on either DNA strand. Our results showed that the repair of small nucleotide heterologies in Escherichia coli extracts was very similar to base-base mismatch repair, being strand-specific and highly biased to the unmethylated strand. The in vitro activity was also dependent on products of mutH, mutL, mutS, and uvrD loci and was equally efficient on nucleotide insertions and deletions. The repair levels of small heterologies were affected by base composition of the heterologies. However, the extent of repair of heteroduplexes containing small heterologous sequences was found to decrease with an increase in the number of unpaired bases. Heteroduplexes containing an extra nucleotide of 22 bases provoked very low level of methyl-directed repair.

Tracking changes in HIV-1 envelope quasispecies using DNA heteroduplex analysis

A DNA heteroduplex tracking assay (HTA) using single-stranded probes is described. This assay provides a rapid means of resolving genetic variants coamplified by PCR and of measuring the level of particular variants in complex populations. To confidently detect minor quasispecies changes, the importance of maximizing template input into nested PCR (nPCR) and of duplicating nPCR and HTA to ensure correct population sampling is highlighted. The sensitivity of detection of rare variants within a genetically mixed population using single-stranded DNA probes is shown to be 1:500. The effects of nucleotide substitution at different locations on heteroduplex electrophoretic mobility are used to illustrate the limits of HTA for mutation detection. This simple assay may be used to track the evolution of HIV as well as to address issues of contamination and transmission.

Interaction of MutS protein with the major and minor grooves of a heteroduplex DNA

Thermus aquaticus MutS protein is a DNA mismatch repair protein that recognizes and binds to heteroduplex DNAs containing mispaired or unpaired bases. Using enzymatic and chemical probe methods, we have examined the binding of Taq MutS protein to a heteroduplex DNA having a single unpaired thymidine residue. DNase I footprinting identifies a symmetrical region of protection 24-28 nucleotides long centered on the unpaired base. Methylation protection and interference studies establish that Taq MutS protein makes contacts with the major groove of the heteroduplex in the immediate vicinity of the unpaired base. Hydroxyl radical and 1, 10-phenanthroline-copper footprinting experiments indicate that MutS also interacts with the minor groove near the unpaired base. Together with the identification of key phosphate groups detected by ethylation interference, these data reveal critical contact points residing in the major and minor grooves of the heteroduplex DNA.

DNA mobility anomalies are determined primarily by polyacrylamide gel concentration, not gel pore size

The dependence of DNA mobility anomalies on gel pore size has been studied in polyacrylamide gels with a wide variety of compositions, using molecular weight ladders containing multiple copies of normal (12B) and anomalously slowly migrating (12A) 147-base pair restriction fragments from plasmid pBR322 as the migrating probe molecules. If the gel pore size is increased by decreasing the total acrylamide concentration (%T) at constant cross-linker ratio (%C), the usual method of increasing gel pore size, the mobility anomalies decrease with increasing gel pore radius as though the 12A multimers were retarded by a sieving mechanism. However, the decrease in the mobility anomalies is independent of whether the apparent gel pore radius is larger or smaller than the DNA radius of gyration, suggesting that gel pore size is not the controlling variable. If the acrylamide concentration is held constant and the gel pore size is increased by decreasing %C at constant %T, the mobility anomalies of the largest 12A multimers (6 mers and higher) decrease with increasing gel pore radius, because of sieving effects, until the effective gel pore radius becomes approximately equal to the DNA radius of gyration, after which the mobility anomalies level off and become independent of gel pore size. The mobility anomalies exhibited by 5-mers and smaller multimers of fragment 12A are independent of gel pore radius in all gels with constant %T. Similar results are observed with a molecular weight ladder containing phased A-tracts from the kinetoplast bending locus. Since the anomalous electrophoretic mobilities depend primarily on the total acrylamide concentration in the gel, and not on the apparent gel pore radius, increases in the magnitude of the mobility anomalies with increasing gel concentration (and decreasing gel pore radius) cannot be taken as evidence for DNA curvature.

DNA structural forms

In the years that have passed since the publication of Wolfram Saenger's classic book on nucleic acid structure (Saenger, 1984), a considerable amount of new data has been accumulated on the range of conformations which can be adopted by DNA. Many unusual species have joined the DNA zoo, including new varieties of two, three and four stranded helices. Much has been learnt about intrinsic DNA curvature, dynamics and conformational transitions and many types of damaged or deformed DNA have been investigated. In this article, we will try to summarise this progress, pointing out the scope of the various experimental techniques used to study DNA structure, and, where possible, trying to discern the rules which govern the behaviour of this subtle macromolecule. The article is divided into six major sections which begin with a general discussion of DNA structure and then present successively, B-DNA, DNA deformations, A-DNA, Z-DNA and DNARNA hybrids. An extensive set of references is included and should serve the reader who wishes to delve into greater detai.

Differentiation between vaccine-related and wild-type polioviruses using a heteroduplex mobility assay

A heteroduplex mobility assay (HMA) was developed for intratypic differentiation between poliovirus isolates. The assay is based on polymerase chain reaction (PCR) amplification of a 480 base pair fragment which encodes a variable segment of VP1, followed by denaturation and reannealing of the resulting single strands with those from reference Sabin targets. Mismatches between wild-type and Sabin vaccine templates result in the formation of detectable heteroduplexes of reduced electrophoretic mobility. Poliovirus strains confirmed previously as wild-type or vaccine-like by PCR and sequencing were all correctly identified using the HMA. Mixtures of both wild-type and vaccine-like strains in a single isolate could also be detected using this technique. The results of this study demonstrate that heteroduplex analysis is a simple, rapid, and sensitive means for differentiating between vaccine-like and wild-type poliovirus isolates.

Recognition of DNA insertion/deletion mismatches by an activity in Saccharomyces cerevisiae

An activity in nuclear extracts of S.cerevisiae binds specifically to heteroduplexes containing four to nine extra bases in one strand. The specificity of this activity (IMR, for insertion mismatch recognition) in band shift assays was confirmed by competition experiments. IMR is biochemically and genetically distinct from the MSH2 dependent, single base mismatch binding activity. The two activities migrate differently during electrophoresis, they are differentially competable and their spectra of mispair binding are distinct. Furthermore, IMR activity is observed in extracts from an msh2- msh3- msh4- strain. IMR exhibits specificity for insertion mispairs in two different sequence contexts. Binding is influenced by the structure of the mismatch since an insertion with a hairpin configuration is not recognized by this activity. IMR does not result from single-strand binding because single-stranded probes to not yield IMR complex and single-stranded competitors are unable to displace insertion heteroduplexes from the complex. Similar results with intrinsically bent duplexes make it unlikely that recognition is conferred by a bend alone. Heteroduplexes bound by IMR do not contain any obvious damage. These findings are consistent with the idea that yeast contains a distinct recognition factor, IMR that is specific for insertion/deletion mismatches.

Rapid molecular epidemiology of human immunodeficiency virus transmission

Close sequence homology between strains of HIV-1 have been used to corroborate cases of epidemiologically identified transmission. As an alternative to extensive DNA sequence analysis, genetic relateness between pairs of HIV quasispecies was estimated using the reduced electrophoretic mobilities of HIV-1 envelope DNA heteroduplexes through polyacrylamide gels. All six infections acquired in a dental practice in the late 1980s and four of six infections acquired through blood product transfusions and sexual contact in 1984-1985 could be rapidly identified. A rising level of genetic diversity within HIV-1 subtype B facilitated the detection of later transmission events. Transmission linkages could be detected up to 4 years following infection. The simple and rapid technique of DNA heteroduplex tracking can therefore assist epidemiological investigations of HIV transmission and potentially of other genetically variable infectious agents.

Kinking of DNA and RNA by base bulges

Images

Kinking of DNA and RNA helices by bulged nucleotides observed by fluorescence resonance energy transfer

Fluorescence resonance energy transfer (FRET) has been used to demonstrate the bending of DNA and RNA helices for three series of double-stranded molecules containing bulge loops of unopposed adenosine nucleotides (An, n = 0-9). Fluorescein and rhodamine were covalently attached to the 5' termini of the two component strands. Three different methods were applied to measure the FRET efficiencies. The extent of energy transfer within each series increases as the number of bulged nucleotides varies from 1 to 7, indicating a shortening of the end-to-end distance. This is consistent with a bending of DNA and RNA helices that is greater for larger bulges. The FRET efficiency for DNA molecules with A9 bulges is lower than the efficiency for the corresponding A7 bulged molecules, although the A9 molecules exhibit increased electrophoretic retardation. Ranges of bending angles can be estimated from the FRET results.

Human immunodeficiency virus type 1 evolution in vivo tracked by DNA heteroduplex mobility assays

High mutation rates and strong selective pressures imposed on human immunodeficiency viruses in vivo result in the formation of pools of genetic variants known as quasispecies. DNA heteroduplex mobility and tracking analyses were used to monitor the generation of HIV sequence diversity, to estimate quasispecies complexity, and to assess the turnover of genetic variants to approach an understanding of the relationship between viral quasispecies evolution in vivo and disease progression. Proviral DNA pools were nearly homogeneous soon after sexual transmission. The emergence and clearance of individual variants then occurred at different rates in different individuals. High quasispecies complexity was found in long-term-infected, asymptomatic individuals, while rapid CD4+ cell decline and AIDS were often, but not always, associated with lower quasispecies complexity. Proviral genetic variation was often low following in vitro culture, because of the outgrowth of one or a few variants that often became more abundant only later as proviruses in peripheral blood mononuclear cells. These studies provide insight into the dynamics of human immunodeficiency virus sequence changes in vivo and illustrate the utility of heteroduplex analysis for the study of phenomena associated with rapid genetic changes.

Conformation-sensitive gel electrophoresis for rapid detection of single-base differences in double-stranded PCR products and DNA fragments: evidence for solvent-induced bends in DNA heteroduplexes

Several techniques have recently been developed to detect single-base mismatches in DNA heteroduplexes that contain one strand of wild-type and one strand of mutated DNA. Here we tested the hypothesis that an appropriate system of mildly denaturing solvents can amplify the tendency of single-base mismatches to produce conformational changes, such as bends in the double helix, and thereby increase the differential migration of DNA heteroduplexes and homoduplexes during gel electrophoresis. The best separations of heteroduplexes and homoduplexes were obtained with a standard 6% polyacrylamide gel polymerized in 10% ethylene glycol/15% formamide/Tris-taurine buffer. As predicted by the hypothesis of solvent-induced bends, when the concentration of either ethylene glycol or formamide was increased, the differential migration decreased. Also, single-base mismatches within 50 bp of one end of a heteroduplex did not produce differential migration. Sixty of 68 single-base mismatches in a series of PCR products were detected in some 59 different sequence contexts. The eight mismatches not detected were either within 50 bp of the nearest end of the PCR product or in isolated high-melting-temperature domains. Therefore, it was possible to predict in advance the end regions and sequence contexts in which mismatches may be difficult to detect. The procedure can be applied to any PCR products of 200-800 bp and requires no special equipment or preparation of samples.

The structure of branched DNA species

Branched DNA molecules provide a challenging set of structural problems. Operationally we define branched DNA species as molecules in which double helical segments are interrupted by abrupt discontinuities, and we draw together a number of different kinds of structure in the class, including helical junctions of different orders, and base bulges (Fig. 1).

Solution structure of the exocyclic 1,N2-propanodeoxyguanosine adduct opposite deoxyadenosine in a DNA nonamer duplex at pH 8.9. model of pH-dependent conformational transition

The solution structure of the complementary d(C1-A2-T3-G4-X5-G6-T7-A8-C9).d(G10-T11-A12-C13-A14-C15-A16-T17-G18) DNA duplex (designated X.A 9-mer), which contains a 1,N2-propanodeoxyguanosine exocyclic adduct X5 opposite deoxyadenosine A14 at the center, is pH dependent [Kouchakdjian, M., Eisenberg, M., Live, D., Marinelli, E., Grollman, A., & Patel, D.J. (1990) Biochemistry 29, 4456-4465]. In our previous paper [Huang, P., & Eisenberg, M. (1992) Biochemistry 31, 6518-6532] we established the three-dimensional structure of this X.A 9-mer duplex at pH 5.8 by use of restrained molecular dynamics followed by NOE-based back-calculation refinement. The present paper discusses the structure at pH 8.9 and the pH-dependent conformational transition between the structures at pH 5.8 and at pH 8.9. The structure at pH 8.9 is calculated starting from five different conformations. The final structures converge and agree well with the experimental NOE intensities. These structures are essentially B-type DNA (with X5 and A14 in the BII conformation while the other residues are in the most commonly described BI conformation) and display an approximate 27 degrees kink at the center of the helix. At the kink site, X5 is positioned in the major groove with the exocyclic ring directed toward the G6.C13 base pair, unstacked from the flanking base G6 and exposed to the solvent. A14, opposite the lesion, remains stacked with its neighbor C15, but not with C13. The kinked helix can accommodate the rotation of the bulky X5 about its glycosidic bond. We propose here a model for the pH-dependent transition. Our model explains the conformational change, which includes the anti and syn rotation of the bulky adduct around its glycosidic bond, with a minimal energy barrier and with an overall kink of the DNA helix. These new findings, fully consistent with the NMR experimental data, were revealed only after restrained dynamics refinement. Distance-restrained energy minimization by itself was insufficient, as shown by the previous NMR study.

Simple, rapid, and accurate determination of deletion mutations by automated DNA sequencing of heteroduplex fragments of the adenomatous polyposis coli (APC) gene generated by PCR amplification

Germline mutations in patients with familial adenomatous polyposis were analyzed by polymerase chain reaction (PCR) amplification of the adenomatous polyposis coli gene. PCR products from heterozygous patients for deletions of this gene formed four distinct bands on polyacrylamide gel electrophoresis. The four fragments were subsequently purified and both strands of each fragment were directly sequenced, using an automated DNA sequencer and the same primers as those for PCR amplification. It was found that the two slower migrating fragments were "bulge" heteroduplexes, while the other two were homoduplexes made up of two wild-type strands and two deletion-mutant strands, respectively. The sites of deletions in the adenomatous polyposis coli gene could be exactly determined in four of the five patients. In an attempt to identify deletion-carriers of familial adenomatous polyposis at the presymptomatic stage, a family study was also carried out, and two children were found to have the same mutations as those of their affected parents. The direct sequencing of heteroduplex fragments generated during PCR amplification is a potentially useful method for detecting mutations of not only the adenomatous polyposis coli gene but also many other genes of genetic diseases.

Three-dimensional structures of bulge-containing DNA fragments

The three-dimensional structure of a DNA tridecamer d(CGCAGAATTCGCG)2 containing bulged adenine bases was determined by single crystal X-ray diffraction methods, at 120 K, to 2.6 A resolution. The structure is a B-DNA type double helix with a single duplex in the asymmetric unit. One of the bulged adenine bases loops out from the double helix, while the other stacks in to it. This is in contrast to our preliminary finding, which indicated that both adenine bases were looped out. This revised model was confirmed by the use of a covalently bound heavy-atom derivative. The conformation of the looped-out bulge hardly disrupts base stacking interactions of the bases flanking it. This is achieved by the backbone making a "loop-the-loop" curve with the extra adenine flipping over with respect to the other nucleotides in the strand. The looped-out base intercalates into the stacked-in bulge site of a symmetrically related duplex. The looped-out and stacked-in bases form an A.A reversed Hoogsteen base-pair that stacks between the surrounding base-pairs, thus stabilizing both bulges. The double helix is frayed at one end with the two "melted" bases participating in intermolecular interactions. A related structure, of the same tridecamer, after soaking the crystals with proflavin, was determined to 3.2 A resolution. The main features of this B-DNA duplex are basically similar to the native tridecamer but differ in detail especially in the conformation of the bulged-out base. Accommodation of a large perturbation such as that described here with minimal disruption of the double helix shows both the flexibility and resiliency of the DNA molecule.

Comparative NMR study of A(n)-bulge loops in DNA duplexes: intrahelical stacking of A, A-A, and A-A-A bulge loops

We have prepared a series of deoxyoligonucleotide duplexes of the sequence d(G-C-A-T-C-G-X-G-C-T-A-C-G).d(C-G-T-A-G-C-C-G-A-T-G-C), in which X represents either one (A), two (A-A), or three (A-A-A) unpaired adenine basis. Using two-dimensional proton and phosphorus NMR spectroscopy, we have characterized conformational features of these bulge-loop duplexes in solution. We find that Watson-Crick hydrogen bonding is intact for all 12 base pairs, including the GC bases that flank the bulge loop. Observation of NOE connectivities in both H2O and D2O allows us to unambiguously localize all of the bulged adenine residues to intrahelical positions within the duplex. This is in contrast to an earlier model for multiple-base bulge loops in DNA [Bhattacharyya, A., & Lilley, D. M. J. (1989) Nucleic Acids Res. 17, 6821-6840], in which all but the most 5' bulged base are looped out into solution. We find that insertion of two or three bases into the duplex results in the disruption of specific sequential NOEs for the base step across from the bulge loop site on the opposite strand. This disruption is characterized by a partial shearing apart of these bases, such that certain sequential NOEs for this base step are preserved. We observe a downfield-shifted phosphorus resonance, which we assign in the A-A-A bulge duplex to the 3' side of the last bulged adenine residue. Proton and phosphorus chemical shift trends within the An-bulge duplex series indicate that there is an additive effect on the structural perturbations caused by additional unpaired bases within the bulge loop. This finding parallels previous observations [Bhattacharyya, A., & Lilley, D. M. J. (1989) Nucleic Acids Res. 17, 6821-6840; Hsieh, C.-H., & Griffith, J. D. (1989) Proc. Natl. Acad. Sci. U.S.A. 86, 4833-4837] on the magnitude of the induced bending of DNA duplexes by multiple-base bulge loops.