Search

Specific targeting of protein–nucleic acid interactions is an area of current interest, for example, in the regulation of gene-expression. Most transcription factor proteins bind in the DNA major groove; however, we are interested in an approach using small molecules to target the minor groove to control expression by an allosteric mechanism. In an effort to broaden sequence recognition of DNA-targeted-small-molecules to include both A·T and G·C base pairs, we recently discovered that the heterocyclic diamidine, DB2277, forms a strong monomer complex with a DNA sequence containing 5΄-AAAGTTT-3΄. Competition mass spectrometry and surface plasmon resonance identified new monomer complexes, as well as unexpected binding of two DB2277 with certain sequences. Inherent microstructural differences within the experimental DNAs were identified through computational analyses to understand the molecular basis for recognition. These findings emphasize the critical nature of the DNA minor groove microstructure for sequence-specific recognition and offer new avenues to design synthetic small molecules for effective regulation of gene-expression.

Sequence-specific binding to DNA is crucial for targeting transcription factor-DNA complexes to modulate gene expression. The heterocyclic diamidine, DB2277, specifically recognizes a single G•C base pair in the minor groove of mixed base pair sequences of the type AAAGTTT. NMR spectroscopy reveals the presence of major and minor species of the bound compound. To understand the principles that determine the binding affinity and orientation in mixed sequences of DNA, over thirty DNA hairpin substrates were examined by NMR and thermal melting. The NMR exchange dynamics between major and minor species shows that the exchange is much faster than compound dissociation determined from biosensor–surface plasmon resonance. Extensive modifications of DNA sequences resulted in a unique DNA sequence with binding site AAGATA that binds DB2277 in a single orientation. A molecular docking result agrees with the model representing rapid flipping of DB2277 between major and minor species. Imino spectral analysis of a 15N-labeled central G clearly shows the crucial role of the exocyclic amino group of G in sequence-specific recognition. Our results suggest that this approach can be expanded to additional modules for recognition of more sequence-specific DNA complexes. This approach provides substantial information about the sequence-specific, highly efficient, dynamic nature of minor groove binding agents.

Forward genetic screens have successfully identified many genes and continue to be powerful tools for dissecting biological processes in Arabidopsis and other model species. Next-generation sequencing technologies have revolutionized the time-consuming process of identifying the mutations that cause a phenotype of interest. However, due to the cost of such mapping-by-sequencing experiments, special attention should be paid to experimental design and technical decisions so that the read data allows to map the desired mutation. Here, we simulated different mapping-by-sequencing scenarios. We first evaluated which short-read technology was best suited for analyzing gene-rich genomic regions in Arabidopsis and determined the minimum sequencing depth required to confidently call single nucleotide variants. We also designed ways to discriminate mutagenesis-induced mutations from background Single Nucleotide Polymorphisms in mutants isolated in Arabidopsis non-reference lines. In addition, we simulated bulked segregant mapping populations for identifying point mutations and monitored how the size of the mapping population and the sequencing depth affect mapping precision. Finally, we provide the computational basis of a protocol that we already used to map T-DNA insertions with paired-end Illumina-like reads, using very low sequencing depths and pooling several mutants together; this approach can also be used with single-end reads as well as to map any other insertional mutagen. All these simulations proved useful for designing experiments that allowed us to map several mutations in Arabidopsis.

Structural results with minor groove binding agents, such as netropsin, have provided detailed, atomic level views of DNA molecular recognition. Solution studies, however, indicate that there is complexity in the binding of minor groove agents to a single site. Netropsin, for example, has two DNA binding enthalpies in isothermal titration calorimetry (ITC) experiments that indicate the compound simultaneously forms two thermodynamically different complexes at a single AATT site. Two proposals for the origin of this unusual observation have been developed: (i) two different bound species of netropsin at single binding sites and (ii) a netropsin induced DNA hairpin to duplex transition. To develop a better understanding of DNA recognition complexity, the two proposals have been tested with several DNAs and the methods of mass spectrometry (MS), polyacrylamide gel electrophoresis (PAGE) and nuclear magnetic resonance spectroscopy in addition to ITC. All of the methods with all of the DNAs investigated clearly shows that netropsin forms two different complexes at AATT sites, and that the proposal for an induced hairpin to duplex transition in this system is incorrect.

Heterocyclic diamidines are compounds with antiparasitic properties that target the minor groove of kinetoplast DNA. The mechanism of action of these compounds is unknown, but topological changes to DNA structures are likely to be involved. In this study, we have developed a polyacrylamide gel electrophoresis-based screening method to determine topological effects of heterocyclic diamidines on four minor groove target sequences: AAAAA, TTTAA, AAATT and ATATA. The AAAAA and AAATT sequences have the largest intrinsic bend, whereas the TTTAA and ATATA sequences are relatively straight. The changes caused by binding of the compounds are sequence dependent, but generally the topological effects on AAAAA and AAATT are similar as are the effects on TTTAA and ATATA. A total of 13 compounds with a variety of structural differences were evaluated for topological changes to DNA. All compounds decrease the mobility of the ATATA sequence that is consistent with decreased minor groove width and bending of the relatively straight DNA into the minor groove. Similar, but generally smaller, effects are seen with TTTAA. The intrinsically bent AAAAA and AAATT sequences, which have more narrow minor grooves, have smaller mobility changes on binding that are consistent with increased or decreased bending depending on compound structure.

To determine what topological changes antiparasitic heterocyclic dications can have on kinetoplast DNA, we have constructed ligation ladders, with phased A5 and ATATA sequences in the same flanking sequence context, as models. Bending by the A5 tract is observed, as expected, while the ATATA sequence bends DNA very little. Complexes of these DNAs with three diamidines containing either furan, thiophene or selenophene groups flanked by phenylamidines were investigated along with netropsin. With the bent A5 ladder the compounds caused either a slight increase or decrease in the bending angle. Surprisingly, however, with ATATA all of the compounds caused significant bending, to values close to or even greater than the A5 bend angle. Results with a mixed cis sequence, which has one A5 and one ATATA, show that the compounds bend ATATA in the same direction as a reference A5 tract, that is, into the minor groove. These results are interpreted in terms of a groove structure for A5 which is largely pre-organized for a fit to the heterocyclic amidines. With ATATA the groove is intrinsically wider and must close to bind the compounds tightly. The conformational change at the binding site then leads to significant bending of the alternating DNA sequence.

Alternative models have been presented to provide explanations for the sequence-dependent variation of the DNA minor groove width. In a structural model groove narrowing in A-tracts results from direct, short-range interactions among DNA bases. In an electrostatic model, the narrow minor groove of A-tracts is proposed to respond to sequence-dependent localization of water and cations. Molecular dynamics simulations on partially methylphosphonate substituted helical chains of d(TATAGGCCTATA) and d(CGCGAATTCGCG) duplexes have been carried out to help evaluate the effects of neutralizing DNA phosphate groups on the minor groove width. The results show that the time-average minor groove width of the GGCC duplex becomes significantly more narrow on neutralizing the phosphate backbone with methylphosphonates. The minor groove of the AATT sequence is normally narrow and the methylphosphonate substitutions have a smaller but measurable affect on this sequence. These results and models provide a system that can be tested by experiment and they support the hypothesis that the electrostatic environment around the minor groove affects the groove width in a sequence-dependent dynamic and time-average manner.

The antibiotic AT2433-B1 belongs to a therapeutically important class of antitumor agents. This natural product contains an indolocarbazole aglycone connected to a unique disaccharide consisting of a methoxyglucose and an amino sugar subunit, 2,4-dideoxy-4-methylamino-L-xylose. The configuration of the amino sugar distinguishes AT2433-B1 from its diastereoisomer iso-AT2433-B1. Here we have investigated the interaction of these two disaccharide indolocarbazole derivatives with different DNA sequences by means of DNase I footprinting and surface plasmon resonance (SPR). Accurate binding measurements performed at 4 and 25 degrees C using the BIAcore SPR method revealed that AT2433-B1 binds considerably more tightly to a hairpin oligomer containing a [CG](4) block than to an oligomer with a central [AT](4) tract. The kinetic analysis shows that the antibiotic dissociates much more slowly from the GC sequence compared to the AT one. Preferential binding of AT2433-B1 to GC-rich sequences in DNA was independently confirmed by DNase I footprinting experiments performed with a 117 bp DNA restriction fragment. The specific binding sequence 5'-AACGCCAG identified from the footprints was then converted into a biotin-labeled DNA hairpin duplex and compound interactions with this specific sequence were characterized by high resolution BIAcore SPR experiments. Such a combined approach provided a detailed understanding of the molecular basis of DNA recognition. The discovery that the glycosyl antibiotic AT2433-B1 preferentially recognizes defined sequences offers novel opportunities for the future design of sequence-specific DNA-reading small molecules.

An imidazole-containing polyamide trimer, f-ImImIm, where f is a formamido group, was recently found using NMR methods to recognize T*G mismatched base pairs. In order to characterize in detail the T*G recognition affinity and specificity of imidazole-containing polyamides, f-ImIm, f-ImImIm and f-PyImIm were synthesized. The kinetics and thermodynamics for the polyamides binding to Watson-Crick and mismatched (containing one or two T*G, A*G or G*G mismatched base pairs) hairpin oligonucleotides were determined by surface plasmon resonance and circular dichroism (CD) methods. f-ImImIm binds significantly more strongly to the T*G mismatch-containing oligonucleotides than to the sequences with other mismatched or with Watson-Crick base pairs. Compared with the Watson-Crick CCGG sequence, f-ImImIm associates more slowly with DNAs containing T*G mismatches in place of one or two C*G base pairs and, more importantly, the dissociation rate from the T*G oligonucleotides is very slow (small k(d)). These results clearly demonstrate the binding selectivity and enhanced affinity of side-by-side imidazole/imidazole pairings for T*G mismatches and show that the affinity and specificity increase arise from much lower k(d) values with the T*G mismatched duplexes. CD titration studies of f-ImImIm complexes with T*G mismatched sequences produce strong induced bands at approximately 330 nm with clear isodichroic points, in support of a single minor groove complex. CD DNA bands suggest that the complexes remain in the B conformation.

ETS transcription factors mediate a wide array of cellular functions and are attractive targets for pharmacological control of gene regulation. We report the inhibition of the ETS-family member PU.1 with a panel of novel heterocyclic diamidines. These diamidines are derivatives of furamidine (DB75) in which the central furan has been replaced with selenophene and/or one or both of the bridging phenyl has been replaced with benzimidazole. Like all ETS proteins, PU.1 binds sequence specifically to 10-bp sites by inserting a recognition helix into the major groove of a 5'-GGAA-3' consensus, accompanied by contacts with the flanking minor groove. We showed that diamidines target the minor groove of AT-rich sequences on one or both sides of the consensus and disrupt PU.1 binding. Although all of the diamidines bind to one or both of the expected sequences within the binding site, considerable heterogeneity exists in terms of stoichiometry, site-site interactions and induced DNA conformation. We also showed that these compounds accumulate in live cell nuclei and inhibit PU.1-dependent gene transactivation. This study demonstrates that heterocyclic diamidines are capable of inhibiting PU.1 by targeting the flanking sequences and supports future efforts to develop agents for inhibiting specific members of the ETS family.

Direct modulation of gene expression by targeting oncogenic transcription factors is a new area of research for cancer treatment. ERG, an ETS-family transcription factor, is commonly over-expressed or translocated in leukaemia and prostate carcinoma. In this work, we selected the di-(thiophene-phenyl-amidine) compound DB1255 as an ERG/DNA binding inhibitor using a screening test of synthetic inhibitors of the ERG/DNA interaction followed by electrophoretic mobility shift assays (EMSA) validation. Spectrometry, footprint and biosensor-surface plasmon resonance analyses of the DB1255/DNA interaction evidenced sequence selectivity and groove binding as dimer. Additional EMSA evidenced the precise DNA-binding sequence required for optimal DB1255/DNA binding and thus for an efficient ERG/DNA complex inhibition. We further highlighted the structure activity relationships from comparison with derivatives. In cellulo luciferase assay confirmed this modulation both with the constructed optimal sequences and the Osteopontin promoter known to be regulated by ERG and which ERG-binding site was protected from DNaseI digestion on binding of DB1255. These data showed for the first time the ERG/DNA complex modulation, both in vitro and in cells, by a heterocyclic diamidine that specifically targets a portion of the ERG DNA recognition site.

The development of small molecules to control gene expression could be the spearhead of future-targeted therapeutic approaches in multiple pathologies. Among heterocyclic dications developed with this aim, a phenyl-furan-benzimidazole dication DB293 binds AT-rich sites as a monomer and 5′-ATGA sequence as a stacked dimer, both in the minor groove. Here, we used a protein/DNA array approach to evaluate the ability of DB293 to specifically inhibit transcription factors DNA-binding in a single-step, competitive mode. DB293 inhibits two POU-domain transcription factors Pit-1 and Brn-3 but not IRF-1, despite the presence of an ATGA and AT-rich sites within all three consensus sequences. EMSA, DNase I footprinting and surface-plasmon-resonance experiments determined the precise binding site, affinity and stoichiometry of DB293 interaction to the consensus targets. Binding of DB293 occurred as a cooperative dimer on the ATGA part of Brn-3 site but as two monomers on AT-rich sites of IRF-1 sequence. For Pit-1 site, ATGA or AT-rich mutated sequences identified the contribution of both sites for DB293 recognition. In conclusion, DB293 is a strong inhibitor of two POU-domain transcription factors through a cooperative binding to ATGA. These findings are the first to show that heterocyclic dications can inhibit major groove transcription factors and they open the door to the control of transcription factors activity by those compounds.

Deoxyribooligonucleotides containing 19 repeating bases of A, T or U were prepared with normal phosphodiester (dA19, dT19, dU19) or methylphosphonate (dA*19, dT*19, dU*19) linkages. Complexes of these strands have been investigated at 1:1 and 1:2 molar ratios (purine:pyrimidine) by thermal melting and gel electrophoresis. There are dramatic sequence dependent differences in stabilities of complexes containing methylphosphonate strands. Duplexes of dA*19 with dT19 or dU19 have sharp melting curves, increased Tm values, and slopes of Tm versus log (sodium ion activity) plots reduced by about one half relative to their unmodified 'parent' duplexes. Duplexes of dA19 with either dT*19 or dU*19, however, have broader melting curves, reduced Tm values at most salt concentrations and slopes of less than one tenth the values for the unmodified duplexes. Duplex stabilization due to reduced phosphate charge repulsion is offset in the pyrimidine methylphosphonate complexes by steric and other substituent effects. Triple helical complexes with dA19 + 2dT19 and dA19 + 2dU19, which can be detected by biphasic melting curves and gel electrophoresis, are stable at increased Na+ or Mg+2 concentrations. Surprisingly, however, no triple helix forms, even at very high salt concentrations, when any normal strand(s) is replaced by a methylphosphonate strand. Since triple helical complexes with methylphosphonates have been reported for shorter oligomers, inhibition with larger oligomers may vary due to their length and extent of substitution.

Polyamides containing an N-terminal formamido (f) group bind to the minor groove of DNA as staggered, antiparallel dimers in a sequence-specific manner. The formamido group increases the affinity and binding site size, and it promotes the molecules to stack in a staggered fashion thereby pairing itself with either a pyrrole (Py) or an imidazole (Im). There has not been a systematic study on the DNA recognition properties of the f/Py and f/Im terminal pairings. These pairings were analyzed here in the context of f-ImPyPy, f-ImPyIm, f-PyPyPy and f-PyPyIm, which contain the central pairing modes, -ImPy- and -PyPy-. The specificity of these triamides towards symmetrical recognition sites allowed for the f/Py and f/Im terminal pairings to be directly compared by SPR, CD and DeltaT (M) experiments. The f/Py pairing, when placed next to the -ImPy- or -PyPy- central pairings, prefers A/T and T/A base pairs to G/C base pairs, suggesting that f/Py has similar DNA recognition specificity to Py/Py. With -ImPy- central pairings, f/Im prefers C/G base pairs (>10 times) to the other Watson-Crick base pairs; therefore, f/Im behaves like the Py/Im pair. However, the f/Im pairing is not selective for the C/G base pair when placed next to the -PyPy- central pairings.

To complement available structure and binding results and to develop a detailed understanding of the basis for selective molecular recognition of T.G mismatches in DNA by imidazole containing polyamides, a full thermodynamic profile for formation of the T.G-polyamide complex has been determined. The amide-linked heterocycles f-ImImIm and f-PyImIm (where f is formamido group, Im is imidazole and Py is pyrrole) were studied by using biosensor-surface plasmon resonance (SPR) and isothermal titration calorimetry (ITC) with a T.G mismatch containing DNA hairpin duplex and a similar DNA with only Watson-Crick base pairs. Large negative binding enthalpies for all of the polyamide-DNA complexes indicate that the interactions are enthalpically driven. SPR results show slower complex formation and stronger binding of f-ImImIm to the T.G than to the match site. The thermodynamic analysis indicates that the enhanced binding to the T.G site is the result of better entropic contributions. Negative heat capacity changes for the complex are correlated with calculated solvent accessible surface area changes and indicate hydrophobic contributions to complex formation. DNase I footprinting analysis in a long DNA sequence provided supporting evidence that f-ImImIm binds selectively to T.G mismatch sites.

The binding of two symmetric bis-benzimidazole compounds, 2,2-bis-[4′-(3″-dimethylamino-1″-propyloxy)phenyl]-5,5-bi-1H-benzimidazole and its piperidinpropylphenyl analog, to the minor groove of DNA, have been studied by DNA footprinting, surface plasmon resonance (SPR) methods and molecular dynamics simulations in explicit solvent. The footprinting and SPR methods find that the former compound has enhanced affinity and selectivity for AT sequences in DNA. The molecular modeling studies have suggested that, due to the presence of the oxygen atom in each side chain of the former compound, a water molecule is immobilized and effectively bridges between side chain and DNA base edges via hydrogen bonding interactions. This additional contribution to ligand–DNA interactions would be expected to result in enhanced DNA affinity, as is observed.

The interaction of berberinium chloride with DNA has been investigated using spectrophotometry, viscometric titrations with sonicated and closed circular superhelical DNA, and flow polarized fluorescence. The binding results for berberinium were found to fit the neighbor exclusion model. The two viscometric titrations and flow polarized fluorescence results exclusion model. The two viscometric titrations and flow polarized fluorescence results also indicated that berberinium binds to DNA by intercalation. Titration of sonicated DNA with berberinium produced viscosity increases which were less than those obtained with quinacrine and the titration of superhelical DNA indicated a significantly smaller unwinding angle for intercalation of berberinium than for quinacrine. These results can be interpreted in terms of a model in which (i) berberinium is partially intercalated into the double helix, or (ii) the alkaloid is more completely intercalated into the double helix, but causes bending of the helix due to the slight nonplanarity of the berberinium ring system, or (iii) a combination of (i) and (22).