Contiguous Base Pairs Research Articles

Triplex-forming properties and enzymatic incorporation of a base-modified nucleotide capable of duplex DNA recognition at neutral pH.

The sequence-specific recognition of duplex DNA by unmodified parallel triplex-forming oligonucleotides is restricted to low pH conditions due to a necessity for cytosine protonation in the third strand. This has severely restricted their use as gene-targeting agents, as well as for the detection and/or functionalisation of synthetic or genomic DNA. Here I report that the nucleobase 6-amino-5-nitropyridin-2-one (Z) finally overcomes this constraint by acting as an uncharged mimic of protonated cytosine. Synthetic TFOs containing the nucleobase enabled stable and selective triplex formation at oligopurine-oligopyrimidine sequences containing multiple isolated or contiguous GC base pairs at neutral pH and above. Moreover, I demonstrate a universal strategy for the enzymatic assembly of Z-containing TFOs using its commercially available deoxyribonucleotide triphosphate. These findings seek to improve not only the recognition properties of TFOs but also the cost and/or expertise associated with their chemical syntheses.

Mirror-Image Oligonucleotides: History and Emerging Applications.

As chiral molecules, naturally occurring d-oligonucleotides have enantiomers, l-DNA and l-RNA, which are comprised of l-(deoxy)ribose sugars. These mirror-image oligonucleotides have the same physical and chemical properties as that of their native d-counterparts, yet are highly orthogonal to the stereospecific environment of biology. Consequently, l-oligonucleotides are resistant to nuclease degradation and many of the off-target interactions that plague traditional d-oligonucleotide-based technologies; thus making them ideal for biomedical applications. Despite a flurry of interest during the early 1990s, the inability of d- and l-oligonucleotides to form contiguous Watson-Crick base pairs with each other has ultimately led to the perception that l-oligonucleotides have only limited utility. Recently, however, scientists have begun to uncover novel strategies to harness the bio-orthogonality of l-oligonucleotides, while overcoming (and even exploiting) their inability to Watson-Crick base pair with the natural polymer. Herein, a brief history of l-oligonucleotide research is presented and emerging l-oligonucleotide-based technologies, as well as their applications in research and therapy, are presented.

Parallel Guanine Duplex and Cytosine Duplex DNA with Uninterrupted Spines of AgI-Mediated Base Pairs.

Hydrogen bonding between nucleobases produces diverse DNA structural motifs, including canonical duplexes, guanine (G) quadruplexes, and cytosine (C) i-motifs. Incorporating metal-mediated base pairs into nucleic acid structures can introduce new functionalities and enhanced stabilities. Here we demonstrate, using mass spectrometry (MS), ion mobility spectrometry (IMS), and fluorescence resonance energy transfer (FRET), that parallel-stranded structures consisting of up to 20 G-AgI-G contiguous base pairs are formed when natural DNA sequences are mixed with silver cations in aqueous solution. FRET indicates that duplexes formed by poly(cytosine) strands with 20 contiguous C-AgI-C base pairs are also parallel. Silver-mediated G-duplexes form preferentially over G-quadruplexes, and the ability of Ag+ to convert G-quadruplexes into silver-paired duplexes may provide a new route to manipulating these biologically relevant structures. IMS indicates that G-duplexes are linear and more rigid than B-DNA. DFT calculations were used to propose structures compatible with the IMS experiments. Such inexpensive, defect-free, and soluble DNA-based nanowires open new directions in the design of novel metal-mediated DNA nanotechnology.

Molecular Mechanism, Dynamics, and Energetics of Protein-Mediated Dinucleotide Flipping in a Mismatched DNA: A Computational Study of the RAD4-DNA Complex.

DNA damage alters genetic information and adversely affects gene expression pathways leading to various complex genetic disorders and cancers. DNA repair proteins recognize and rectify DNA damage and mismatches with high fidelity. A critical molecular event that occurs during most protein-mediated DNA repair processes is the extrusion of orphaned bases at the damaged site facilitated by specific repairing enzymes. The molecular-level understanding of the mechanism, dynamics, and energetics of base extrusion is necessary to elucidate the molecular basis of protein-mediated DNA damage repair. The present article investigates the molecular mechanism of dinucleotide extrusion in a mismatched DNA (containing a stretch of three contiguous thymidine-thymidine base pairs) facilitated by Radiation sensitive 4 (RAD4), a key DNA repair protein, on an atom-by-atom basis using molecular dynamics (MD) and umbrella-sampling (US) simulations. Using atomistic models of RAD4-free and RAD4-bound mismatched DNA, the free energy profiles associated with extrusion of mismatched partner bases are determined for both systems. The mismatched bases adopted the most stable intrahelical conformation, and their extrusion was unfavorable in RAD4-free mismatched DNA due to the presence of prohibitively high barriers (>12.0 kcal/mol) along the extrusion pathways. Upon binding of RAD4 to the DNA, the global free energy minimum is shifted to the extrahelical state indicating the key role of RAD4-DNA interactions in catalyzing the dinucleotide base extrusion in the DNA-RAD4 complex. The critical residues of RAD4 contributing to the conformational stability of the mismatched bases are identified, and the energetics of insertion of a β-hairpin of RAD4 into the DNA duplex is examined. The conformational energy landscape-based mechanistic insight into RAD4-mediated base extrusion provided here may serve as a useful baseline to understand the molecular basis of xeroderma pigmentosum C (XPC)-mediated DNA damage repair in humans.

Heterochiral DNA Strand-Displacement Circuits.

The absence of a straightforward strategy to interface native d-DNA with its enantiomer l-DNA-oligonucleotides of opposite chirality are incapable of forming contiguous Watson-Crick base pairs with each other-has enforced a "homochiral" paradigm over the field of dynamic DNA nanotechnology. As a result, chirality, a key intrinsic property of nucleic acids, is often overlooked as a design element for engineering of DNA-based devices, potentially limiting the types of behaviors that can be achieved using these systems. Here we introduce a toehold-mediated strand-displacement methodology for transferring information between orthogonal DNA enantiomers via an achiral intermediary, opening the door for "heterochiral" DNA nanotechnology having fully interfaced d-DNA and l-DNA components. Using this approach, we demonstrate several heterochiral DNA circuits having novel capabilities, including autonomous chiral inversion of DNA sequence information and chirality-based computing. In addition, we show that heterochiral circuits can directly interface endogenous RNAs (e.g., microRNAs) with bioorthogonal l-DNA, suggesting applications in bioengineering and nanomedicine. Overall, this work establishes chirality as a design parameter for engineering of dynamic DNA nanotechnology, thereby expanding the types of architectures and behaviors that can be realized using DNA.

The Varied Effects of Tetraalkylammonium Cations on the Properties of DNA

We have been examining the effect of monovalent cations on the melting transition of short model DNA hairpins using capillary electrophoresis. Here we focus on the effects contributed by two disparate monovalent cations, Na and tetrabutylammonium (TBA). Each melting transition yields two parameters, the melting temperature and the amplitude. The dependence of the melting temperature on cation concentration indicates the number of cations released upon melting. Our measurements indicate that more Na is released than TBA upon melting. By contrast, the amplitude of the transition is greater in TBA than in Na. We interpret these observations to indicate that the condensed counter ion cloud about DNA contains fewer TBA than Na, owing to the larger size and apolar surface of TBA. Consequently, the net charge of DNA is greater in TBA than in Na, increasing mobility but lowering the Tm by about 20o due to charge repulsion. This accounts for the ability of TBA to function as a DNA denaturant. A second feature of TBA is that it appears to associate preferentially with AT base pairs and that such association is diminished in sequences containing A-tracts, four or more contiguous AT base pairs

Finding the unexpected--how we identified a second class of introns and the U12-dependent spliceosome.

Finding the unexpected--how we identified a second class of introns and the U12-dependent spliceosome.

Structural Design and Validation of Linkers for Zinc Finger Proteins

The possibility of using engineered DNA binding proteins to either control gene expression (artificial transcription factors) or to alter the sequence or structure of specific genes (Zinc Finger Nucleases) is a highly appealing idea with a wide range of applications in biological research and molecular therapeutics. Zinc finger proteins (ZFPs), belonging to Cys2His2 family constitute the most common DNA binding motifs found in eukaryotes. ZFP normally occur in repeats of two to three zinc finger motifs (ZFMs) to bind 6-9 contiguous DNA base pairs in a sequence specific manner. Several methods of varying complexity are available to engineer ZFPs that can target all the 64 codons in the genome. Although ZFPs are becoming a powerful tool for site specific modification in the genome, several challenges remain before the full potential of ZFPs can be realized. The engineered ZFPs generated using the present design platforms target mostly base triplets with 5' Guanine (GNN, where N is any nucleotide) and the non-GNN or AT rich modules are difficult to target. In the present project we attempt to address this challenge by designing linker regions between the ZFP motifs to target non-contiguous base pairs in the DNA. This will increase the number of targetable DNA sequences by an order of magnitude and will help to realize the full potential of ZFPs. Using structure based methods, we provide an extensive library of possible linker molecules that can be introduced between the individual zinc finger motifs to skip up to 10 base pairs between adjacent zinc finger protein recognition sites in the DNA sequences. We also performed a proof of principle experiment to validate the binding affinity and specificity of one of the computationally designed ZFP to its target DNA sequence.

Akt-Mediated Phosphorylation of Argonaute 2 Downregulates Cleavage and Upregulates Translational Repression of MicroRNA Targets

A high-throughput RNA interference (RNAi) screen targeting 542 genes of the human kinome was usedto discover regulators of RNAi. Here we report that the proto-oncogene Akt-3/PKBγ (Akt3) phosphorylates Argonaute 2 (Ago2) at S387, which downregulates cleavage and upregulates translational repression of endogenous microRNA (miRNA)-targeted messenger RNAs (mRNAs). We further demonstrate that Akt3 coimmunoprecipitates with Ago2 and phosphorylation of Ago2 at S387 facilitates its interaction with GW182 and localization to cytoplasmic processing bodies (P bodies), where miRNA-targeted mRNAs are thought to be stored and degraded. Therefore, Akt3-mediated phosphorylation of Ago2 is a molecular switch between target mRNA cleavage and translational repression activities of Ago2.

A rule of seven in Watson-Crick base-pairing of mismatched sequences

Sequence recognition through base pairing is essential for DNA repair and gene regulation but the basic rules governing this process remain elusive. In particular, the kinetics of annealing between two imperfectly matched strands is not well characterized despite its potential importance in nucleic acids-based biotechnologies and gene silencing. Here we use single molecule fluorescence to visualize the multiple annealing and melting reactions of two untethered strands inside a porous vesicle, allowing us to quantify precisely the annealing and melting rates. The data as a function of mismatch position suggest that seven contiguous base pairs are needed for rapid annealing of DNA and RNA. This phenomenological rule of seven may underlie the requirement of seven nucleotides complementarity to seed gene silencing by small non-coding RNA and may help guide performance improvement in DNA and RNA-based bio- and nano-technologies where off-target effects can be detrimental.

Molecular Characterization of Borrelia persica, the Agent of Tick Borne Relapsing Fever in Israel and the Palestinian Authority

The identification of the Tick Borne Relapsing Fever (TBRF) agent in Israel and the Palestinian Authority relies on the morphology and the association of Borrelia persica with its vector Ornithodoros tholozani. Molecular based data on B. persica are very scarce as the organism is still non-cultivable. In this study, we were able to sequence three complete 16S rRNA genes, 12 partial flaB genes, 18 partial glpQ genes, 16 rrs-ileT intergenic spacers (IGS) from nine ticks and ten human blood samples originating from the West Bank and Israel. In one sample we sequenced 7231 contiguous base pairs that covered completely the region from the 5′end of the 16S rRNA gene to the 5′end of the 23S rRNA gene comprising the whole 16S rRNA (rrs), and the following genes: Ala tRNA (alaT), Ile tRNA (ileT), adenylosuccinate lyase (purB), adenylosuccinate synthetase (purA), methylpurine-DNA glycosylase (mag), hypoxanthine-guanine phosphoribosyltransferase (hpt), an hydrolase (HAD superfamily) and a 135 bp 5′ fragment of the 23S rRNA (rrlA) genes. Phylogenic sequence analysis defined all the Borrelia isolates from O. tholozani and from human TBRF cases in Israel and the West Bank as B. persica that clustered between the African and the New World TBRF species. Gene organization of the intergenic spacer between the 16S rRNA and the 23S rRNA was similar to that of other TBRF Borrelia species and different from the Lyme disease Borrelia species. Variants of B. persica were found among the different genes of the different isolates even in the same sampling area.

ChemInform Abstract: Design of a Nonnatural Deoxyribonucleoside for Recognition of GC Base Pairs by Oligonucleotide-Directed Triple Helix Formation.

The deoxyribonucleoside 1-(2-deoxy-beta-D-ribofuranosyl)-3-methyl-5-amino-1H-pyrazolo4,3-d]py rimidin-7-one (P1) was designed such that two specific hydrogen bonds would form with guanine (G) of a Watson-Crick guanine-cytosine (GC) base pair in the major groove of double-helical DNA. One edge of the P1 heterocycle mimics N3-protonated cytosine, which would circumvent the pH dependence observed for the formation of triple helices containing C + GC base triplets. P1 was synthesized in five steps and incorporated by automated methods in pyrimidine oligodeoxyribonucleotides. From affinity cleaving analyses, the stabilities of base triplets decrease in the order P1.GC >> P1.CG >> P1.AT approximately P1.TA (pH 7.4, 35-degrees-C). P1 binds GC base pairs within a pyrimidinte triple-helix motif as selectively and strongly as C but over an extended pH range. Oligodeoxyribonucleotides containing P1 residues were shown to bind within plasmid DNA a single 15 base pair site containing five GC base pairs at pH 7.8 and a single 16 base pair site containing six contiguous GC base pairs at pH 7.4.

Meso-tetrakis(N-methylpyridinium-4-yl)porphyrin at the Minor Groove of Contiguous Adenine-Thymine Base Pairs

Three possible binding modes of cationic meso-tetrakis(N-methylpyridinium-4-yl)porphyrin (TMPyP) to d[(GCATATATGC)2] duplex were investigated by the molecular dynamics (MD) simulation. Among the three binding modes namely, “along the groove”, “across the groove” and “face on the groove”, the “across the groove” model exhibited the largest negative binding free energy and the DNA backbone remained as the B form. In this model, the molecular plain of the TMPyP tilts 45o with respect to the DNA helix axis and is largely exposed to the solvent. TMPyP was stabilized mainly by the interaction between the positively charged neighboring pyridinium moieties of TMPyP and negatively charged phosphate groups of DNA. The result obtained in this work by MD and the report (Jin, B. et al., J. Am. Chem. Soc. 2005, 127, 2417.) that the spectral properties of poly[d(A-T)2] bound TMPyP in the presence and absence of the minor groove binding drug 4',6- diamidino-2-phenylindole are similar, we propose that TMPyP bind across the minor groove of the AT rich- DNA.

Synthesis and Evaluation of an Intercalator–Polyamide Hairpin Designed to Target the Inverted CCAAT Box 2 in the Topoisomerase IIα Promoter

The synthesis and DNA-binding properties of a novel naphthalimide-polyamide hairpin (3) designed to target the inverted CCAAT box 2 (ICB2) site on the topoisomerase IIalpha (topoIIalpha) promoter are described. The polyamide component of 3 was derived from the minor-groove binder, 2, and tailored to bind to the 5'-TTGGT sequence found in and flanking ICB2. The propensity of mitonafide 4 to intercalate between G-C base pairs was exploited by the incorporation of a naphthalimide moiety at the N terminus of 2. Hybrid 3 targeted 5'-CGATTGGT and covered eight contiguous base pairs, which included the underlined ICB2 site. DNase I footprinting analysis with the topoIIalpha promoter sequence demonstrated that 3 bound selectively to the ICB2 and ICB3 sites. Thermal-denaturation studies confirmed these results, and the highest degree of stabilization was found for ICB2 and -3 in preference to ICB1 (4.1, 4.6, and 0.6 degrees C, respectively). CD studies confirmed minor-groove binding and suggested a 1:1 binding stoichiometry. Emission-titration experiments established intercalative binding. Surface plasmon resonance results showed strong binding to ICB2 (2.5x10(7) M(-1)) with no observable binding to ICB1. Furthermore, the binding constant of 3 to ICB2 was larger than that of the parent polyamide 2. The increased binding affinity was primarily due to a reduction in the dissociation-rate constant of the polyamide-DNA complex, which can be attributed to the N-terminal naphthalimide moiety. In addition, the binding site of 3 was larger than that of 2, which innately improved sequence selectivity. We conclude that the polyamide-naphthalimide 3 selectively binds to the ICB2 site by simultaneous intercalation and minor-groove binding, and warrants further investigation as a model compound for the regulation of topoIIalpha gene expression.

Programmable Oligomers for Minor Groove DNA Recognition

The four Watson-Crick base pairs of DNA can be distinguished in the minor groove by pairing side-by-side three five-membered aromatic carboxamides, imidazole (Im), pyrrole (Py), and hydroxypyrrole (Hp), four different ways. On the basis of the paradigm of unsymmetrical paired edges of aromatic rings for minor groove recognition, a second generation set of heterocycle pairs, imidazopyridine/pyrrole (Ip/Py) and hydroxybenzimidazole/pyrrole (Hz/Py), revealed that recognition elements not based on analogues of distamycin could be realized. A new set of end-cap heterocycle dimers, oxazole-hydroxybenzimidazole (No-Hz) and chlorothiophene-hydroxybenzimidazole (Ct-Hz), paired with Py-Py are shown to bind contiguous base pairs of DNA in the minor groove, specifically 5'-GT-3' and 5'-TT-3', with high affinity and selectivity. Utilizing this technology, we have developed a new class of oligomers for sequence-specific DNA minor groove recognition no longer based on the N-methyl pyrrole carboxamides of distamycin.

DNA-binding affinities and sequence selectivity of quaternary benzophenanthridine alkaloids sanguinarine, chelerythrine, and nitidine

A comparative study on the intercalating binding of sanguinarine, chelerythrine, and nitidine with CT DNA, poly(dG–dC)·poly(dG–dC), poly(dA–dT)·poly(dA–dT), and seven sequence-designed double-stranded oligodeoxynucleotides has been performed using fluorometric and spectrophotometric techniques, aiming at providing insights into their sequence selectivity for DNA-binding. The results show that both sanguinarine and nitidine bind preferentially to DNA containing alternating GC base pairs [d(TGCGCA) 2], while chelerythrine exhibits quite distinct sequence selectivity from sanguinarine, which shows a high specificity for DNA containing contiguous GC base pairs [5′-TGGGGA-3′/3′-ACCCCT-5′].

A multivariate prediction model for microarray cross-hybridization.

BackgroundExpression microarray analysis is one of the most popular molecular diagnostic techniques in the post-genomic era. However, this technique faces the fundamental problem of potential cross-hybridization. This is a pervasive problem for both oligonucleotide and cDNA microarrays; it is considered particularly problematic for the latter. No comprehensive multivariate predictive modeling has been performed to understand how multiple variables contribute to (cross-) hybridization.ResultsWe propose a systematic search strategy using multiple multivariate models [multiple linear regressions, regression trees, and artificial neural network analyses (ANNs)] to select an effective set of predictors for hybridization. We validate this approach on a set of DNA microarrays with cytochrome p450 family genes. The performance of our multiple multivariate models is compared with that of a recently proposed third-order polynomial regression method that uses percent identity as the sole predictor. All multivariate models agree that the 'most contiguous base pairs between probe and target sequences,' rather than percent identity, is the best univariate predictor. The predictive power is improved by inclusion of additional nonlinear effects, in particular target GC content, when regression trees or ANNs are used.ConclusionA systematic multivariate approach is provided to assess the importance of multiple sequence features for hybridization and of relationships among these features. This approach can easily be applied to larger datasets. This will allow future developments of generalized hybridization models that will be able to correct for false-positive cross-hybridization signals in expression experiments.

Cenozoic biogeography and evolution in direct-developing frogs of Central America (Leptodactylidae: Eleutherodactylus) as inferred from a phylogenetic analysis of nuclear and mitochondrial genes

We report the first phylogenetic analysis of DNA sequence data for the Central American component of the genus Eleutherodactylus (Anura: Leptodactylidae: Eleutherodactylinae), one of the most ubiquitous, diverse, and abundant components of the Neotropical amphibian fauna. We obtained DNA sequence data from 55 specimens representing 45 species. Sampling was focused on Central America, but also included Bolivia, Brazil, Jamaica, and the USA. We sequenced 1460 contiguous base pairs (bp) of the mitochondrial genome containing ND2 and five neighboring tRNA genes, plus 1300 bp of the c- myc nuclear gene. The resulting phylogenetic inferences were broadly concordant between data sets and among analytical methods. The subgenus Craugastor is monophyletic and its initial radiation was potentially rapid and adaptive. Within Craugastor, the earliest splits separate three northern Central American species groups, milesi, augusti, and alfredi, from a clade comprising the rest of Craugastor. Within the latter clade, the rhodopis group as formerly recognized comprises three deeply divergent clades that do not form a monophyletic group; we therefore restrict the content of the rhodopis group to one of two northern clades, and use new names for the other northern ( mexicanus group) and one southern clade ( bransfordii group). The new rhodopis and bransfordii groups together form the sister taxon to a clade comprising the biporcatus, fitzingeri, mexicanus, and rugulosus groups. We used a Bayesian MCMC approach together with geological and biogeographic assumptions to estimate divergence times from the combined DNA sequence data. Our results corroborated three independent dispersal events for the origins of Central American Eleutherodactylus: (1) an ancestor of Craugastor entered northern Central America from South American in the early Paleocene, (2) an ancestor of the subgenus Syrrhophus entered northern Central America from the Caribbean at the end of the Eocene, and (3) a wave of independent dispersal events from South America coincided with formation of the Isthmus of Panama during the Pliocene. We elevate the subgenus Craugastor to the genus rank.

DNA Sequence Recognition by the Antitumor Drug Ditercalinium

The antitumor drug ditercalinium is a rare example of a noncovalent DNA-binding ligand that forms bisintercalation complexes via the major groove of the double helix. Previous structural studies have revealed that the two connected pyridocarbazolium chromophores intercalate into DNA with the positively charged bis(ethylpiperidinium) linking chain oriented to the wide groove side of the helix. Although the interaction of ditercalinium with short oligonucleotides containing 4-6 contiguous GC base pairs has been examined in detail by biophysical and theoretical approaches, the sequence preference for ditercalinium binding to long DNA fragments that offer a wide variety of binding sites has been investigated only superficially. Here we have investigated both sequence preferences and possible molecular determinants of selectivity in the binding of ditercalinium to DNA, primarily using methods based upon DNase I footprinting. A range of multisite DNA substrates, including several natural restriction fragments and different PCR-generated fragments containing unconventional bases (2,6-diaminopurine, inosine, uridine, 5-fluoro- and 5-methylcytosine, 7-deazaguanine, 7-deazaadenine, and N(7)-cyanoboranoguanine), have been employed to show that ditercalinium selectively recognizes certain GC-rich sequences in DNA and to identify some of the factors which affect its DNA-binding sequence selectivity. Specifically, the footprinting data have revealed that the 2-amino group on the purines or the 5-methyl group on the pyrimidines is not essential for the formation of ditercalinium-DNA complexes whereas the major groove-oriented N(7) of guanine does appear as a key element in the molecular recognition process. The loss of N(7) at guanines but not adenines is sufficient to practically abolish sequence-selective binding of ditercalinium to DNA. Thus, as expected for a major groove binding drug, the N(7) of guanine is normally required for effective complex formation with GC base pairs, but interestingly the substitution of the N(7) with a relatively bulky cyanoborane group does not markedly affect the sequence recognition process. Therefore, the hydrogen bond accepting capability at N(7) of guanines is not sufficient to explain the GC-selective drug-DNA association, and the implications of these findings are considered.

Influence of compound structure on affinity, sequence selectivity, and mode of binding to DNA for unfused aromatic dications related to furamidine

In the course of a program aimed at developing sequence-specific gene-regulatory small organic molecules, we have investigated the DNA interactions of a new series of nine diphenylfuran dications related to the antiparasitic drug furamidine (DB75). Two types of structural modifications were tested: the terminal amidine groups of DB75 were shifted from the para to the meta position, and the amidines were replaced with imidazoline or dimethyl-imidazoline groups, to test the importance of both the position and nature of positively charged groups on DNA interactions. The interactions of these compounds with DNA and oligonucleotides were studied by a combination of biochemical and biophysical techniques. Absorption and CD measurements suggested that the drugs bind differently to AT and GC sequences in DNA. The para-para dications, like DB75, bind into the minor groove of poly(dAT)(2) and intercalate between the base pairs of poly(dGC)(2), as revealed by electric linear dichroism experiments. In contrast, the meta-meta compounds exhibit a high tendency to intercalate into DNA whatever the target sequence. The lack of sequence selectivity of the meta-meta compounds containing amidines or dimethyl-imidazoline groups was also evident from DNase I footprinting and surface plasmon resonance (SPR) experiments. Accurate binding measurements using the BIAcore SPR method revealed that all nine compounds bind with similar affinity to an immobilized GC sequence DNA hairpin but exhibit very distinct affinities for the corresponding AT hairpin oligonucleotide. The minor groove-binding para-para compounds have a high specificity for AT sequences. The biophysical data clearly indicate that shifting the cationic substituents from the para to the meta position results in a loss of specificity and change in binding mode. The strong AT selectivity of the para-para compounds was independently confirmed by DNase I footprinting experiments performed with a range of DNA restrictions fragments. In terms of AT selectivity, the compounds rank in the order para-para > para-meta > meta-meta. The para dications bind preferentially to sequences containing four contiguous AT base pairs. Additional footprinting experiments with substrates containing the 16 possible [A.T](4) blocks indicated that the presence of a TpA step within an [A.T] (4) block generally reduces the extent of binding. The diverse methods, from footprinting to SPR to dichroism, provide a consistent model for the interactions of the diphenylfuran dications with DNA of different sequences. Altogether, the results attest unequivocally that the binding mode for unfused aromatic cations can change completely depending on substituent position and DNA sequence. These data provide a rationale to explain the relationships between sequence selectivity and mode of binding to DNA for unfused aromatic dications related to furamidine.