Deletions of bases in one strand of duplex DNA, in contrast to single-base mismatches, produce highly kinked molecules: possible relevance to the folding of single-stranded nucleic acids

Increasing the length of poly-pyrimidine bulges broadens RNA conformational ensembles with minimal impact on stacking energetics

Helical elements separated by bulges frequently undergo transitions between unstacked and coaxially stacked conformations during the folding and function of noncoding RNAs. Here, we examine the dynamic properties of poly-pyrimidine bulges of varying length (n = 1-4, 7) across a range of Mg²⁺ concentrations using HIV-1 TAR RNA as a model system and solution NMR spectroscopy. In the absence of Mg²⁺, helices linked by bulges with n ≥ 3 residues adopt predominantly unstacked conformations (stacked population <15%), whereas one-bulge and two-bulge motifs adopt predominantly stacked conformations (stacked population >74%). In the presence of 3 mM Mg²⁺, the helices predominantly coaxially stack (stacked population >84%), regardless of bulge length, and the midpoint for the Mg²⁺-dependent stacking transition is within threefold regardless of bulge length. In the absence of Mg²⁺, the difference between free energy of interhelical coaxial stacking across the bulge variants is estimated to be ∼2.9 kcal/mol, based on an NMR chemical shift mapping with stacking being more energetically disfavored for the longer bulges. This difference decreases to ∼0.4 kcal/mol in the presence of Mg²⁺ NMR RDCs and resonance intensity data show increased dynamics in the stacked state with increasing bulge length in the presence of Mg²⁺ We propose that Mg²⁺ helps to neutralize the growing electrostatic repulsion in the stacked state with increasing bulge length thereby increasing the number of coaxial conformations that are sampled. Energetically compensated interhelical stacking dynamics may help to maximize the conformational adaptability of RNA and allow a wide range of conformations to be optimally stabilized by proteins and ligands.

Characterizing the bending and flexibility induced by bulges in DNA duplexes

Advances in DNA nanotechnology have stimulated the search for simple motifs that can be used to control the properties of DNA nanostructures. One such motif, which has been used extensively in structures such as polyhedral cages, two-dimensional arrays, and ribbons, is a bulged duplex, that is, two helical segments that connect at a bulge loop. We use a coarse-grained model of DNA to characterize such bulged duplexes. We find that this motif can adopt structures belonging to two main classes: one where the stacking of the helices at the center of the system is preserved, the geometry is roughly straight, and the bulge is on one side of the duplex and the other where the stacking at the center is broken, thus allowing this junction to act as a hinge and increasing flexibility. Small loops favor states where stacking at the center of the duplex is preserved, with loop bases either flipped out or incorporated into the duplex. Duplexes with longer loops show more of a tendency to unstack at the bulge and adopt an open structure. The unstacking probability, however, is highest for loops of intermediate lengths, when the rigidity of single-stranded DNA is significant and the loop resists compression. The properties of this basic structural motif clearly correlate with the structural behavior of certain nano-scale objects, where the enhanced flexibility associated with larger bulges has been used to tune the self-assembly product as well as the detailed geometry of the resulting nanostructures. We further demonstrate the role of bulges in determining the structure of a "Z-tile," a basic building block for nanostructures.

New insights into the fundamental role of topological constraints as a determinant of two-way junction conformation

Recent studies have shown that topological constraints encoded at the RNA secondary structure level involving basic steric and stereochemical forces can significantly restrict the orientations sampled by helices across two-way RNA junctions. Here, we formulate these topological constraints in greater quantitative detail and use this topological framework to rationalize long-standing but poorly understood observations regarding the basic behavior of RNA two-way junctions. Notably, we show that the asymmetric nature of the A-form helix and the finite length of a bulge provide a physical basis for the experimentally observed directionality and bulge-length amplitude dependence of bulge induced inter-helical bends. We also find that the topologically allowed space can be modulated by variations in sequence, particularly with the addition of non-canonical GU base pairs at the junction, and, surprisingly, by the length of the 5′ and 3′ helices. A survey of two-way RNA junctions in the protein data bank confirms that junction residues have a strong preference to adopt looped-in, non-canonically base-paired conformations, providing a route for extending our bulge-directed framework to internal loop motifs and implying a simplified link between secondary and tertiary structure. Finally, our results uncover a new simple mechanism for coupling junction-induced topological constraints with tertiary interactions.

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.

Topological constraints: using RNA secondary structure to model 3D conformation, folding pathways, and dynamic adaptation

Accompanying recent advances in determining RNA secondary structure is the growing appreciation for the importance of relatively simple topological constraints, encoded at the secondary structure level, in defining the overall architecture, folding pathways, and dynamic adaptability of RNA. A new view is emerging in which tertiary interactions do not define RNA 3D structure, but rather, help select specific conformers from an already narrow, topologically pre-defined conformational distribution. Studies are providing fundamental insights into the nature of these topological constraints, how they are encoded by the RNA secondary structure, and how they interplay with other interactions, breathing new meaning to RNA secondary structure. New approaches have been developed that take advantage of topological constraints in determining RNA backbone conformation based on secondary structure, and a limited set of other, easily accessible constraints. Topological constraints are also providing a much-needed framework for rationalizing and describing RNA dynamics and structural adaptation. Finally, studies suggest that topological constraints may play important roles in steering RNA folding pathways. Here, we review recent advances in our understanding of topological constraints encoded by the RNA secondary structure.

Barcode-like heteroduplex DNA pattern as an aid for rapid identification of anthracnose fungi

We have shown the usefulness of the heteroduplex mobility assay (HMA) for phylogenetic analysis and for the discrimination of closely related Colletotrichum species. Because the heteroduplex mobility of a tested strain shows a unique banding pattern that is the function of the sequence of the referred strain, we further explored the potential use of heteroduplex DNA patterns (HPs) as DNA fingerprints for the identification of these fungi. The 29 Colletotrichum strains previously identified by HMA to be taxonomic members of CG, CA, CM, CC and CL species groups were re-examined with an emphasis on their unique heteroduplex banding patterns. The species attributes of these tested strains were characterized by HMA using ITS fragments amplified from six representative Colletotrichum strains as pairwise compared references. By comparing the unique homoduplex and heteroduplex banding patterns of each tested strain on a polyacrylamide gel with those of the respective reference strain, the species identity of tested strains was determined. The obtained barcode-like HPs classified these 29 Colletotrichum strains into 6 distinctive groups: CG1, CG2, CA, CM, CC and CL. Notably, the HPs differentiated strains CG1 and CG2, which differed in their ITS sequences by only six bases. The presented results revealed that the species-characteristic barcode-like HP classification of ITS regions is a relatively rapid and valuable system for species identification of Colletotrichum species. The potential use of the established barcode-like system for the identification of anthracnose fungi and other fungal pathogens is discussed.

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.

DNA models of trinucleotide frameshift deletions: the formation of loops and bulges at the primer-template junction

Although mechanisms of single-nucleotide residue deletion have been investigated, processes involved in the loss of longer nucleotide sequences during DNA replication are poorly understood. Previous reports have shown that in vitro replication of a 3′-TGC TGC template sequence can result in the deletion of one 3′-TGC. We have used low-energy circular dichroism (CD) and fluorescence spectroscopy to investigate the conformations and stabilities of DNA models of the replication intermediates that may be implicated in this frameshift. Pyrrolocytosine or 2-aminopurine residues, site-specifically substituted for cytosine or adenine in the vicinity of extruded base sequences, were used as spectroscopic probes to examine local DNA conformations. An equilibrium mixture of four hybridization conformations was observed when template bases looped-out as a bulge, i.e. a structure flanked on both sides by duplex DNA. In contrast, a single-loop structure with an unusual unstacked DNA conformation at its downstream edge was observed when the extruded bases were positioned at the primer–template junction, showing that misalignments can be modified by neighboring DNA secondary structure. These results must be taken into account in considering the genetic and biochemical mechanisms of frameshift mutagenesis in polymerase-driven DNA replication.

A loop loop interaction and a K-turn motif located in the lysine aptamer domain are important for the riboswitch gene regulation control

The lysine riboswitch is associated to the lysC gene in Bacillus subtilis, and the binding of lysine modulates the RNA structure to allow the formation of an intrinsic terminator presumably involved in transcription attenuation. The complex secondary structure of the lysine riboswitch aptamer is organized around a five-way junction that undergoes structural changes upon ligand binding. Using single-round transcription assays, we show that a loop-loop interaction is important for lysine-induced termination of transcription. Moreover, upon close inspection of the secondary structure, we find that an unconventional kink-turn motif is present in one of the stems participating in the loop-loop interaction. We show that the K-turn adopts a pronounced kink and that it binds the K-turn-binding protein L7Ae of Archaeoglobus fulgidus in the low nanomolar range. The functional importance of this K-turn motif is revealed from single-round transcription assays, which show its importance for efficient transcription termination. This motif is essential for the loop-loop interaction, and consequently, for lysine binding. Taken together, our results depict for the first time the importance of a K-turn-dependent loop-loop interaction for the transcription regulation of a lysine riboswitch.

The use of unilateral PCR to identify prominent heteroduplexes formed during PCR of the mouse microsatellite locus D17Mit23

Microsatellite markers are useful tools for understanding the evolutionary history of discrete segments of the mammalian genome. We used the microsatellite marker D17Mit23 to study the portion of the mouse genome known as the t complex, a naturally occurring variant of Chromosome 17. We identified an allelic variant of D17Mit23, which is shared by two forms of the t complex, the t haplotypes t(w2) and t(Lub2). Polymerase chain reaction (PCR) analysis of DNA samples from mice that were heterozygous for either t haplotype resulted in gel patterns with prominent bands of higher molecular weight in addition to the bona-fide D17Mit23 alleles. The appearance of these higher molecular weight bands, although consistent with heteroduplex formation, was not diminished through the use of reconditioning PCR. We used a modified form of asymmetric PCR, called "unilateral PCR," to show that the higher molecular weight bands are heteroduplexes and to identify their constituent strands. Certain microsatellite motifs may be especially prone to the production of prominent heteroduplex products, and this may lead to the erroneous genotyping of DNA samples.

Intrinsic curvature in the VP1 gene of SV40: comparison of solution and gel results

DNA restriction fragments that are stably curved are usually identified by polyacrylamide gel electrophoresis because curved fragments migrate more slowly than normal fragments containing the same number of basepairs. In free solution, curved DNA molecules can be identified by transient electric birefringence (TEB) because they exhibit rotational relaxation times that are faster than those of normal fragments of the same size. In this article, the results observed in free solution and in polyacrylamide gels are compared for a highly curved 199-basepair (bp) restriction fragment taken from the VP1 gene in Simian Virus 40 (SV40) and various sequence mutants and insertion derivatives. The TEB method of overlapping fragments was used to show that the 199-bp fragment has an apparent bend angle of 46 +/- 2 degrees centered at sequence position 1922 +/- 2 bp. Four unphased A- and T-tracts and a mixed A3T4-tract occur within a span of approximately 60 bp surrounding the apparent bend center; for brevity, this 60-bp sequence element is called a curvature module. Modifying any of the A- or T-tracts in the curvature module by site-directed mutagenesis decreases the curvature of the fragment; replacing all five A- and T-tracts by random-sequence DNA causes the 199-bp mutant to adopt a normal conformation, with normal electrophoretic mobilities and birefringence relaxation times. Hence, stable curvature in this region of the VP1 gene is due to the five unphased A- and T- tracts surrounding the apparent bend center. Discordant solution and gel results are observed when long inverted repeats are inserted within the curvature module. These insertion derivatives migrate anomalously slowly in polyacrylamide gels but have normal, highly flexible conformations in free solution. Discordant solution and gel results are not observed if the insert does not contain a long inverted repeat or if the long inverted repeat is added to the 199-bp fragment outside the curvature module. The results suggest that long inverted repeats can form hairpins or cruciforms when they are located within a region of the helix backbone that is intrinsically curved, leading to large mobility anomalies in polyacrylamide gels. Hairpin/cruciform formation is not observed in free solution, presumably because of rapid conformational exchange. Hence, DNA restriction fragments that migrate anomalously slowly in polyacrylamide gels are not necessarily stably curved in free solution.

Sequencing-based detection of low-frequency human immunodeficiency virus type 1 drug-resistant mutants by an RNA/DNA heteroduplex generator-tracking assay

Drug-resistant viruses may be present as minority variants during early treatment failures or following discontinuation of failed antiretroviral regimens. A limitation of the traditional direct PCR population sequencing method is its inability to detect human immunodeficiency virus type 1 (HIV-1) variants present at frequencies lower than 20%. A drug resistance genotyping assay based on the isolation and DNA sequencing of minority HIV protease variants is presented here. A multiple-codon-specific heteroduplex generator probe was constructed to improve the separation of HIV protease genes varying in sequence at 12 codons associated with resistance to protease inhibitors. Using an RNA molecule as probe allowed the simple sequencing of protease variants isolated as RNA/DNA heteroduplexes with different electrophoretic mobilities. The protease gene RNA heteroduplex generator-tracking assay (RNA-HTA) was tested on plasma quasispecies from 21 HIV-1-infected persons in whom one or more protease resistance mutations emerged during therapy or following initiation of salvage regimens. In 11 of 21 cases, RNA-HTA testing of virus from the first episode of virologic failure identified protease resistance mutations not seen by population-based PCR sequencing. In 8 of these 11 cases, all of the low-frequency drug resistance mutations detected exclusively by RNA-HTA during the first episode became detectable by population-based PCR sequencing at the later time point. Distinct sets of protease mutations could be linked on different genomes in patients with high-frequency protease gene lineages. The enhanced detection of minority drug resistance variants using a sequencing-based assay may improve the efficacy of genotype-assisted salvage therapies.

Mechanics of DNA flexibility visualized by selective 2'-amine acylation at nucleotide bulges

We used selective acylation of 2'-amine-substituted nucleotides to visualize local backbone conformations that occur preferentially at bulged sites in DNA duplexes. 2'-Amine acylation reports local nucleotide flexibility because unconstrained 2'-amino nucleotides more readily reach a reactive conformation in which the amide-forming transition state is stabilized by interactions between the amine nucleophile and the adjacent 3'-phosphodiester group. Bulged 2'-amine-substituted cytidine nucleotides react approximately 20-fold more rapidly than nucleotides constrained by base-pairing at 35 degrees C. In contrast, base-paired 2'-amine-substituted nucleotides flanked by a 5' or 3' bulge react two- or six-fold more rapidly, respectively, than the perfectly paired duplex. The relative lack of 2'-amine reactivity for nucleotides adjacent to a DNA bulge emphasizes, first, that structural perturbations do not extend significantly into the flanking duplex structure. Second, the exquisite sensitivity towards very local perturbations in nucleic acid structure suggests that 2'-amine acylation can be used to chemically interrogate deletion mutations in DNA. Finally, these data support the mechanical interpretation that the reactive ribose conformation for 2'-amine acylation requires that the base lies out of the helix and in the major groove, a mechanistic insight useful for designing 2'-amine-based sensors.

The kink-turn motif in RNA is dimorphic, and metal ion-dependent

The kink-turn (K-turn) is a new motif in RNA structure that was identified by examination of the crystal structures of the ribosome. We examined the structural and dynamic properties of this element in free solution. The K-turn RNA exists in a dynamic equilibrium between a tightly kinked conformation and a more open structure similar to a simple bulge bend. The highly kinked form is stabilized by the noncooperative binding of metal ions, but a significant population of the less-kinked form is present even in the presence of relatively high concentrations of divalent metal ions. The conformation of the tightly kinked population is in excellent agreement with that of the K-turn structures observed in the ribosome by crystallography. The end-to-end FRET efficiency of this species agrees closely with that of the ribosomal K-turn, and the direction of the bend measured by comparative gel electrophoresis also corresponds very well. These results show that the tightly kinked conformation of the K-turn requires stabilization by other factors, possibly by protein binding, for example. The K-turn is therefore unlikely to be of itself a primary organizing feature in RNA.

Turnover of env variable region 1 and 2 genotypes in subjects with late-stage human immunodeficiency virus type 1 infection

The env gene of human immunodeficiency virus type 1 (HIV-1) includes some of the most genetically diverse regions of the viral genome, which are called variable regions 1 through 5 (V1 through V5). We have developed a heteroduplex tracking assay to detect changes in variable regions 1 and 2 of env (V1/V2-HTA). Using sequences from two molecular clones as probes, we have studied the nature of longitudinal virus population changes in a cohort of HIV-1-infected subjects. Viral sequences present in 21 subjects with late-stage HIV-1 infection were initially screened for stability of the virus population by V1/V2-HTA. The virus populations at entry comprised an average of five coexisting V1/V2 genotypic variants (as identified by HTA). Eight of the 21 subjects were examined in detail because of the dynamic behavior of their env variants over an approximately 9-month period. In each of these cases we detected a single discrete transition of V1/V2 genotypes based on monthly sampling. The major V1/V2 genotypes (those present at >10% abundance) from the eight subjects were cloned and sequenced to define the nature of V1/V2 variability associated with a discrete transition. Based on a comparison of V1/V2 genotypic variants present at entry with the newly emerged variants we categorized the newly emerged variants into two groups: variants without length differences and variants with length differences. Variants without length differences had fewer nucleotide substitutions, with the changes biased to either V1 or V2, suggestive of recent evolutionary events. Variants with length differences included ones with larger numbers of changes that were distributed, suggestive of recall of older genotypes. Most length differences were located in domains where the codon motif AVT (V = A, G, C) had become enriched and fixed. Finally, recombination events were detected in two subjects, one of which resulted in the reassortment of V1 and V2 regions. We suggest that turnover in V1/V2 populations was largely driven by selection on either V1 or V2 and that escape was accomplished either through changes focused in the region under selection or by the appearance of a highly divergent variant.

Recognition of bulged DNA by a neocarzinostatin product via an induced fit mechanism

The binding of the wedge-shaped isostructural analogue of the biradical species of the chromphore of antitumor antibiotic neocarzinostatin to sequence-specific bulged DNAs results in alterations in ellipticity of the DNAs. Circular dichroism (CD) spectroscopic results suggest that the drug specifically recognizes bulges of DNA via a combination of conformational selection and induced fit, not by binding to a preorganized site. Analysis of circular dichroism spectra indicates that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the bulge cavity. The effective recognition of the bulge site on duplex DNA appears to depend to a significant extent on the bent groove space being flexible enough to be able to adopt the geometrically optimal conformation compatible with the wedge-shaped drug molecule, rather than involving 'lock and key' recognition. The spectroscopic results indicate a change of DNA conformation, consistent with an allosteric binding model. Spectroscopic studies with various bulged DNAs also reveal that the binding strength directly correlates with the stability of the bulge structures.

Heteroduplex formation in SMN gene dosage analysis

Most spinal muscular atrophy patients lack both copies of SMN1 exon 7 and most carriers have only one copy of SMN1 exon 7. We investigated the effect of SMN1/SMN2 heteroduplex formation on SMN gene dosage analysis, which is an assay to determine copy number of SMN1 exon 7 that utilizes multiplex quantitative polymerase chain reaction (PCR) with DraI digestion to differentiate SMN1 from SMN2. Heteroduplex formation in PCR is a well-described phenomenon. In addition to demonstrating the presence of heteroduplexes by sequence analysis of purified SMN1 bands, we compared the SMN1 signals in various genotype groups (total n = 260) to those in a group lacking SMN2 (n = 13), and we estimated the relative amounts of SMN1/SMN2 heteroduplexes. The SMN1 signal increased as SMN2 copy number increased despite a constant SMN1 copy number, although not all pairwise comparisons showed a statistically significant difference in the SMN1 signal. In conclusion, SMN1/SMN2 heteroduplexes form in SMN gene dosage analysis, falsely increasing the SMN1 signal. External controls for SMN gene dosage analysis should be chosen carefully with regard to SMN2 copy number. The effect of heteroduplex formation should be considered when performing quantitative multiplex PCR.

Heteroduplex mobility assay (HMA) pre-screening: an improved strategy for the rapid identification of inserts selected from phage-displayed peptide libraries

Phage-displayed peptide libraries represent an efficient tool to isolate peptides that bind a given target molecule. After several selection rounds, generally a large pool of target binding phages is obtained. Conventional analysis of the selected phage population involves extensive sequencing of many clones, most of which can be identical. We have adapted the Heteroduplex Mobility Assay (HMA) for pre-screening of phage inserts that were amplified by direct colony PCR of ELISA-positive clones. This strategy allowed for the rapid and reproducible assignment of insert sequences to different 'heteroduplex migration groups'. Sequence analysis of only one representative of each HMA migration group then completes the characterisation of the binding phage population. In our model experiments, only 16% of HMA pre-screened clones required further sequence analysis.

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.

Fluorescence energy transfer analysis of DNA structures containing several bulges and their interaction with CAP

DNA molecules with three bulges separated by double-stranded helical sections of B-DNA were constructed to be used as substrates for DNA-protein binding assays. Fluorescence resonance energy transfer (FRET) between dye molecules attached to the 5'-ends of the DNA molecules is used to monitor the protein binding. The A5 bulge, which consists of five unpaired adenine nucleotides, alters the direction of the helical axis by approximately 80 to 90 at every bulge site. Computer molecular modeling facilitated a pre-selection of suitable helix lengths that bring the labeled ends of the three-bulge DNA molecules (60 to 70 base-pairs long) into close proximity. The FRET experiments verified that the labeled ends of the helices of these long molecules were indeed close. A series of FRET experiments was carried out with two A5 and two A7 bulge molecules. The relative positions of the bulges were varied along the central helical DNA sequence (between the bulges) in order to determine the relative angular juxtapositions of the outlying helical arms flanking the central helical region. The global structural features of the DNA molecules are manifested in the FRET data. The FRET experiments, especially those of the two-bulge series, could be interpreted remarkably well with molecular models based on the NMR structure of the A5 bulge. These models assume that the DNA molecules do not undergo large torsional conformational fluctuations at the bulge sites. The magnitude of the FRET efficiency attests to a relatively rigid structure for many of the long 5'-end-labeled molecules. The changes in the FRET efficiency of three-bulge structures containing the specific binding sequence of the catabolite activator protein (CAP) demonstrated significant deformation of the DNA upon binding of CAP. No direct interaction of CAP with the dyes was observed.

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%.

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.

Rapid determination of transgene copy number in stably-transfected mammalian cells by competitive PCR

We describe here an application of the competitive PCR technique to the analysis of copy number of recombinant rat parathyroid hormone-related protein (rPTHrP) gene in stably-transfected murine erythroleukemia (MEL) cell lines. A single-copy reference gene (endogenous mouse PTHrP gene or mPTHrP) is used as an internal control. This control gene, present in the genome of MEL cells, shares the same primer binding sites as the rPTHrP cDNA but contains an internal PvuII site, which allows resolution of the amplified products after restriction enzyme digestion by polyacrylamide gel electrophoresis (PAGE). The transgene copy number is determined by the ratio of band intensity of the rPTHrP product to that of the mPTHrP product. Using this method, we have determined the copy number of the rPTHrP transgene from isolated genomic DNA, and compared the results with those obtained from Southern blot analysis. In addition, we have demonstrated that the procedure can be applied very simply to whole MEL cells without DNA extractions and that as few as 10(4) cells are required for the analysis.

A rapid method for detecting and distinguishing clinically important yeasts by heteroduplex mobility assays (HMAs)

Heteroduplex Mobility Assays (HMAs) of 700 bp amplified products of a 17S rDNA region were used to identify and differentiate seven yeast species of clinical importance Candida albicans, Torulopsis (Candida) glabrata, Candida tropicalis, Candida parapsilosis, Candida (Pichia) guilliermondii and Hansenula (Pichia) anomala. Distance of heteroduplex migration (dHE) was found to be negatively correlated to the number of nucleotide differences between amplified DNA sequences.

Thermostable flap endonuclease from the archaeon, Pyrococcus horikoshii, cleaves the replication fork-like structure endo/exonucleolytically

The flap endonuclease gene homologue from the hyperthermophilic archaeon, Pyrococcus horikoshii, was overexpressed in Escherichia coli and purified. The results of gel filtration indicated that this protein was a 41-kDa monomer. P. horikoshii flap endonuclease (phFEN) cleaves replication fork-like substrates (RF) and 5' double-strand flap structures (DF) using both flap endonuclease and 5'-3'-exonuclease activities. The mammalian flap endonuclease (mFEN) is a single-strand flap-specific endonuclease (Harrington, J. J., and Lieber, M. R. (1994) EMBO J. 13, 1235-1246), but the action patterns of phFEN appear to be quite different from those of mFEN at this point. The DF-specific flap endonuclease and 5'-exonuclease activities have not yet been reported. Therefore, this is the first report of the specific endo/exonuclease activities of phFEN. The DF-specific 5'-exonuclease activity degraded the downstream primer of 3' single-flap structure and was 15 times higher than the activities against nicked substrates without 3' flap strand. DF-specific flap endonuclease cleaved the 5' double-flap strand in DF and the lagging strand in RF at the junction portion. Because the RF appears to be the intermediate structure, due to the arrest of the replication fork, the double strand breaks after the arrests of the replication forks are probably caused by phFEN.

Conformational analysis of single-base bulges in A-form DNA and RNA using a hierarchical approach and energetic evaluation with a continuum solvent model

The analysis and prediction of non-canonical structural motifs in RNA is of great importance for an understanding of the function and design of RNA structures. A hierarchical method has been employed to generate a large variety of sterically possible conformations for a single-base adenine bulge structure in A -form DNA and RNA. A systematic conformational search was performed on the isolated bulge motif and neighboring nucleotides under the constraint to fit into a continuous helical structure. These substructures were recombined with double-stranded DNA or RNA. Energy minimization resulted in more than 300 distinct bulge conformations. Energetic evaluation using a solvation model based on the finite-difference Poisson-Boltzmann method identified three basic classes of low-energy structures. The three classes correspond to conformations with the bulge base stacked between flanking nucleotides (I), location of the bulge base in the minor groove (II) and conformations with a continuous stacking of the flanking helices and a looped out bulge base (III). For the looped out class, two subtypes (IIIa and IIIb) with different backbone geometries at the bulge site could be distinguished. The conformation of lowest calculated energy was a class I structure with backbone torsion angles close to those in standard A -form RNA. Conformations very close to the extra-helical looped out bulge structure determined by X-ray crystallography were also among the low-energy structures. In addition, topologies observed in other experimentally determined bulge structures have been found among low-energy conformers. The implicit solvent model was further tested by comparing an uridine and adenine bulge flanked by guanine:cytosine base-pairs, respectively. In agreement with the experimental observation, a looped out form was found as the energetically most favorable form for the uridine bulge and a stacked conformation in case of the adenine bulge. The inclusion of solvation effects especially electrostatic reaction field contributions turned out to be critically important in order to select realistic low-energy bulge structures from a large number of sterically possible conformations. The results indicate that the approach might be useful to model the three-dimensional structure of non-canonical motifs embedded in double-stranded RNA, in particular, to restrict the number of possible conformations to a manageable number of conformers with energies below a certain threshold.

Structural flexibility of a DNA hairpin located in the long terminal repeat of the Drosophila 1731 retrotransposon

The structure of the DNA binding site of the Nuclear single-stranded Binding Factor (NssBF), located in the long terminal repeat of the Drosophila 1731 retrotransposon, was investigated by melting temperature experiments, chemical probing and fluorescence measurements using a macrocyclic bis-acridine. The most probable structure of this element, named Bc, mainly involves two hairpins in equilibrium at pH 6.0 at low concentration. The hairpins differ in their apical loop size; 4 and 8 nt. The structural flexibility of Bc probably derives from the three consecutive CATA repeats complementary to the GTAT nucleotides of the palindrome. In contrast, the Bc complementary strand adopts a single hairpin. Since Bc is implicated in repression of transcription via binding of two specific factors, its structural flexibility could be associated with this process.

High-resolution analysis of T-cell receptor beta-chain repertoires using DNA heteroduplex tracking: generally stable, clonal CD8+ expansions in all healthy young adults

The accurate measurement of T-cell receptor (TCR) repertoire changes requires the analysis of a representative sampling of complex T-cell populations. The number and frequency of clonally expanded TCR beta-chain transcripts bearing distinct CDR3 sequences were accurately determined using a simple DNA heteroduplex tracking assay. This method allowed major and minor clonal expansions (> or = 1% of a Vbeta subfamily's transcripts) to be rapidly and reproducibly quantified. Oligoclonal CD8 + cell expansions were detected in all young adults tested, while CD4 + cells generally expressed more polyclonal beta-chain repertoires. The same pattern of CD8 + cells oligoclonality and CD4 + cells polyclonality was observed in asymptomatic HIV-1 infected individuals with high CD4 + cell counts. CD8 + CD45RA + and CD8 + CD45RO + cell fractions both displayed oligoclonal, although distinct, TCR beta chain repertoires while CD8 + cells from umbilical cord blood were generally polyclonal. Oligoclonal CD8 + cell repertoires from young adults were generally stable over a period of weeks, although minor, transient, clonal expansions could also be detected in the absence of symptomatic infections. DNA heteroduplex tracking analysis provided a higher level of sensitivity for the detection of TCR beta chain transcript expansions than CDR3 length (spectrotyping/immunoscope) analysis. DNA heteroduplex tracking of TCR beta-chain transcripts is therefore a simple and sensitive method for assessing the level of clonality and for measuring changes in the TCR beta chain repertoire of different T-cell populations.

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.

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.

Recognition and manipulation of branched DNA structure by junction-resolving enzymes

The junction-resolving enzymes are a class of nucleases that introduce paired cleavages into four-way DNA junctions. They are important in DNA recombination and repair, and are found throughout nature, from eubacteria and their bacteriophages through to higher eukaryotes and their viruses. These enzymes exhibit structure-selective binding to DNA junctions; although cleavage may be more or less sequence-dependent, binding affinity is purely related to the branched structure of the DNA. Binding and cleavage events can be separated for a number of the enzymes by mutagenesis, and mutant proteins that are defective in cleavage while retaining normal junction-selective binding have been isolated. Critical acidic residues have been identified in several resolving enzymes, suggesting a role in the coordination of metal ions that probably deliver the hydrolytic water molecule. The resolving enzymes all bind to junctions in dimeric form, and the subunits introduce independent cleavages within the lifetime of the enzyme-junction complex to ensure resolution of the four-way junction. In addition to recognising the structure of the junction, recent data from four different junction-resolving enzymes indicate that they also manipulate the global structure. In some cases this results in severe distortion of the folded structure of the junction. Understanding the recognition and manipulation of DNA structure by these enzymes is a fascinating challenge in molecular recognition.

The acidic tail of the high mobility group protein HMG-D modulates the structural selectivity of DNA binding

HMG-D is one of the Drosophila counterparts of the vertebrate HMG1/2 class of abundant chromosomal proteins and contains three domains: an HMG domain followed by a basic region and a short acidic carboxyterminal tail. We show that the HMG domain of HMG-D does not bind to deformed DNA structures such as DNA bulges, cis-platinated DNA or four-way junctions but does bind tightly to DNA microcircles, suggesting that in vivo the natural ligands of this domain are tightly bent DNA loops. The flanking basic region substantially increases the DNA-binding activity of the HMG domain to DNA ligands other than microcircles. We demonstrate that the acidic tail alters the structural selectivity of DNA binding by increasing the affinity for deformed DNA and decreasing the affinity for linear B-DNA. Finally, we show that the acidic tail increases the efficiency of constraining preformed negative supercoils but conversely decreases the efficiency of supercoiling relaxed DNA in the presence of topoisomerase I.

Structure of HIV-1 TAR RNA in the absence of ligands reveals a novel conformation of the trinucleotide bulge

Efficient transcription from the human immunodeficiency virus (HIV) promoter depends on binding of the viral regulatory protein Tat to a cis-acting RNA regulatory element, TAR. Tat binds at a trinucleotide bulge located near the apex of the TAR stem-loop structure. An essential feature of Tat-TAR interaction is that the protein induces a conformational change in TAR that repositions the functional groups on the bases and the phosphate backbone that are critical for specific intermolecular recognition of TAR RNA. We have previously determined a high resolution structure for the bound form of TAR RNA using heteronuclear NMR. Here, we describe a high resolution structure of the free TAR RNA based on 871 experimentally determined restraints. In the free TAR RNA, bulged residues U23 and C24 are stacked within the helix, while U25 is looped out. This creates a major distortion of the phosphate backbone between C24 and G26. In contrast, in the bound TAR RNA, each of the three residues from the bulge are looped out of the helix and U23 is drawn into proximity with G26 through contacts with an arginine residue that is inserted between the two bases. Thus, TAR RNA undergoes a transition from a structure with an open and accessible major groove to a much more tightly packed structure that is folded around basic side chains emanating from the Tat protein.

Scanning tunneling microscopy of mercapto-hexyl-oligonucleotides attached to gold

6-mercapto hexyl-oligonucleotides bind to a gold surface strongly enough to permit imaging by a scanning tunneling microscope (STM). STM images showed worm-like chains that were approximately 12-(A-wide for single-stranded DNA and 20-(A-wide for double-stranded DNA. The chain lengths corresponded to 3.4 +/- 0.4 A per basepair for double-stranded DNA and 2.2 +/- 0.4 A per base for single-stranded DNA. This unexpectedly short length for single-stranded DNA was confirmed using oligomers with both single- and double-stranded regions. When the attachment of the samples was weakened (by imaging in water or scraping with the STM tip) the images changed to pairs of "blobs," apparently reflecting the attachment points of the molecules to the gold surface. Given this interpretation, images of DNA containing a five-base bulge imply that the bulge bends the oligomer by 90 degrees +/- 20 degrees.

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.

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.

Molecular recognition of DNA structure by proteins that mediate genetic recombination

The latter half of genetic recombination is mediated by proteins that recognise the structure of the four-way DNA junction, and manipulate this structure. In solution the four-way junction adopts a stacked X-structure in the presence of metal ions. The folding is brought about by the pairwise coaxial stacking of helices in a right-handed antiparallel X-shaped structure. The four-way junction is cleaved by structure-selective resolving enzymes that have been isolated from a wide variety of sources, from eubacteria and their phages through to mammals. In addition, another class of proteins accelerate the branch migration of the junction. These proteins all appear to be divisible into a component that recognises structure and another that carries out a reaction on the junction. Thus the ability of structure-selective binding to the four-way DNA junction is a key feature of enzymes important in genetic recombination.

Modeling large RNAs and ribonucleoprotein particles using molecular mechanics techniques

There is a growing body of low-resolution structural data that can be utilized to devise structural models for large RNAs and ribonucleoproteins. These models are routinely built manually. We introduce an automated refinement protocol to utilize such data for building low-resolution three-dimensional models using the tools of molecular mechanics. In addition to specifying the positions of each nucleotide, the protocol provides quantitative estimates of the uncertainties in those positions, i.e., the resolution of the model. In typical applications, the resolution of the models is about 10-20 A. Our method uses reduced representations and allows us to refine three-dimensional structures of systems as big as the 16S and 23S ribosomal RNAs, which are about one to two orders of magnitude larger than nucleic acids that can be examined by traditional all-atom modeling methods. Nonatomic resolution structural data--secondary structure, chemical cross-links, chemical and enzymatic footprinting patterns, protein positions, solvent accessibility, and so on--are combined with known motifs in RNA structure to predict low-resolution models of large RNAs. These structural constraints are imposed on the RNA chain using molecular mechanics-type potential functions with parameters based on the quality of experimental data. Surface potential functions are used to incorporate shape and positional data from electron microscopy image reconstruction experiments into our models. The structures are optimized using techniques of energy refinement to get RNA folding patterns. In addition to providing a consensus model, the method finds the range of models consistent with the data, which allows quantitative evaluation of the resolution of the model. The method also identifies conflicts in the experimental data. Although our protocol is aimed at much larger RNAs, we illustrate these techniques using the tRNA structure as an example and test-bed.

Strand-invasion of duplex DNA by peptide nucleic acid oligomers

Polyamide oligomers, termed peptide nucleic acids (PNAs), bind with high affinity to both DNA and RNA and offer both antisense and antigene approaches for regulating gene expression. When a PNA binds to a complementary sequence in a double-stranded DNA, one strand of the duplex is displaced, and a stable D-loop is formed. Unlike oligodeoxynucleotides for which binding polarity is determined by the deoxyribose sugar, the unrestrained polyamide backbone of the PNA could permit binding to a DNA target in an orientation-independent manner. We now provide evidence that PNAs can, in fact, bind to their complementary sequence in DNA independent of the DNA-strand polarity--that is, a PNA binds to DNA in both "parallel" and "antiparallel" fashion. With a mixed-sequence 15-mer PNA, kinetic studies of PNA.DNA interactions revealed that D-loop formation was rapid and the complex was stable for several hours. However, when measured either by gel-mobility-shift analysis or RNA polymerase II-elongation termination, D-loop formation was salt dependent, but PNA-strand dissociation was not salt dependent. We observed that D-loop-containing DNA fragments had anomalous gel mobilities that varied as a function of the position of the D-loop relative to the DNA termini. On the basis of permutation analysis, the decreased mobility of the PNA.DNA complex was attributed to a bend in the DNA at or near the D-loop.

Structures of bulged three-way DNA junctions

We have studied a series of three-way DNA junctions containing unpaired bases on one strand at the branch-point of the junctions. The global conformation of the arms of the junctions has been analysed by means of polyacrylamide gel electrophoresis, as a function of conditions. We find that in the absence of added metal ions, all the results for all the junctions can be accounted for by extended structures, with the largest angle being that between the arms defined by the strand containing the extra bases. Upon addition of magnesium (II) or hexamine cobalt (III) ions, the electrophoretic patterns change markedly, indicative of ion-dependent folding transitions for some of the junctions. For the junction lacking the unpaired bases, the three inter-arm angles appear to be quite similar, suggesting an extended structure. However, the addition of unpaired bases permits the three-way junction to adopt a significantly different structure, in which one angle becomes smaller than the other two. These species also exhibit marked protection against osmium addition to thymine bases at the point of strand exchange. These results are consistent with a model in which two of the helical arms undergo coaxial stacking in the presence of magnesium ions, with the third arm defining an angle that depends upon the number of unpaired bases.