Binding of echinomycin to d(GCGC)2 and d(CCGG)2: distinct stacking interactions dictate the sequence-dependent formation of Hoogsteen base pairs

Modulation of Hoogsteen dynamics on DNA recognition

In naked duplex DNA, G–C and A–T Watson-Crick base pairs exist in dynamic equilibrium with their Hoogsteen counterparts. Here, we used nuclear magnetic resonance (NMR) relaxation dispersion and molecular dynamics (MD) simulations to examine how Watson-Crick/Hoogsteen dynamics are modulated upon recognition of duplex DNA by the bisintercalator echinomycin and monointercalator actinomycin D. In both cases, DNA recognition results in the quenching of Hoogsteen dynamics at base pairs involved in intermolecular base-specific hydrogen bonds. In the case of echinomycin, the Hoogsteen population increased 10-fold for base pairs flanking the chromophore most likely due to intermolecular stacking interactions, whereas actinomycin D minimally affected Hoogsteen dynamics at other sites. Modulation of Hoogsteen dynamics at binding interfaces may be a general phenomenon with important implications for DNA–ligand and DNA–protein recognition.

DNA is found in a dynamic equilibrium between standard Watson-Crick (WC) base pairs and non-standard Hoogsteen (HG) base pairs. Here the authors describe the influence of echinomycin and actinomycin D ligands binding on the HG-WC base pair dynamics in DNA.

A historical account of Hoogsteen base-pairs in duplex DNA

In 1957, a unique pattern of hydrogen bonding between N3 and O4 on uracil and N7 and N6 on adenine was proposed to explain how poly(rU) strands can associate with poly(rA)-poly(rU) duplexes to form triplexes. Two years later, Karst Hoogsteen visualized such a noncanonical A-T base-pair through X-ray analysis of co-crystals containing 9-methyladenine and 1-methylthymine. Subsequent X-ray analyses of guanine and cytosine derivatives yielded the expected Watson-Crick base-pairing, but those of adenine and thymine (or uridine) did not yield Watson-Crick base-pairs, instead favoring "Hoogsteen" base-pairing. More than two decades ensued without experimental "proof" for A-T Watson-Crick base-pairs, while Hoogsteen base-pairs continued to surface in AT-rich sequences, closing base-pairs of apical loops, in structures of DNA bound to antibiotics and proteins, damaged and chemically modified DNA, and in polymerases that replicate DNA via Hoogsteen pairing. Recently, NMR studies have shown that base-pairs in duplex DNA exist as a dynamic equilibrium between Watson-Crick and Hoogsteen forms. There is now little doubt that Hoogsteen base-pairs exist in significant abundance in genomic DNA, where they can expand the structural and functional versatility of duplex DNA beyond that which can be achieved based only on Watson-Crick base-pairing. Here, we provide a historical account of the discovery and characterization of Hoogsteen base-pairs, hoping that this will inform future studies exploring the occurrence and functional importance of these alternative base-pairs.

Stability of the I-motif structure is related to the interactions between phosphodiester backbones

The i-motif DNA tetrameric structure is formed of two parallel duplexes intercalated in a head-to-tail orientation, and held together by hemiprotonated cytosine pairs. The four phosphodiester backbones forming the structure define two narrow and wide grooves. The short interphosphate distances across the narrow groove induce a strong repulsion which should destabilize the tetramer. To investigate this point, molecular dynamics simulations were run on the [d(C2)]4 and [d(C4)]4 tetramers in 3'E and 5'E topologies, for which the interaction of the phosphodiester backbones through the narrow groove is different. The analysis of the simulations, using the Molecular Mechanics Generalized Born Solvation Area and Molecular Mechanics Poisson-Boltzmann Solvation Area approaches, shows that it is the van der Waals energy contribution which displays the largest relative difference between the two topologies. The comparison of the solvent-accessible area of each topology reveals that the sugar-sugar interactions account for the greater stability of the 3'E topology. This stresses the importance of the sugar-sugar contacts across the narrow groove which, enforcing the optimal backbone twisting, are essential to the base stacking and the i-motif stability. Tighter interactions between the sugars are observed in the case of N-type sugar puckers.

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.

Solution structure and dynamics of a complex between DNA and the antitumor bisnaphthalimide LU-79553: intercalated ring flipping on the millisecond time scale

Using a combination of nuclear magnetic resonance (NMR) spectroscopy experiments and molecular dynamics, we have analyzed the structure and dynamics of a complex between the bisnaphthalimide drug LU-79553 and the DNA duplex d(ATGCAT)(2). LU-79553 is a DNA-binding topoisomerase II inhibitor that is particularly effective against human solid tumors that are refractory to other drugs. We have found that the two naphthalimide chromophores of the drug bisintercalate at the TpG and CpA steps of the DNA hexanucleotide, stacking mainly with the purine G and A bases from opposite strands. The 3, 7-diazanonylene linker lies in the major groove of the DNA molecule, with its two amino groups hydrogen-bonded to the symmetry-related guanine bases. Unexpectedly, we have detected an unprecedented exchange process between two equivalent and intercalated states of the naphthalimide rings in the drug-DNA complex. The interconversion process takes place by rotational ring flipping, has an activation energy of 22 kcal mol(-)(1) for the two rings, and does not affect the aminoalkyl linker region of the drug. The exchange rate is intermediate to fast on the chemical shift time scale at 36 degrees C (1800 s(-)(1)) but slow at 2 degrees C (20 s(-)(1)). We have also observed limited flexibility for the drug linker on the picosecond time scale on the basis of NMR data and a time-averaged restrained molecular dynamics simulation. The implications of the structural and dynamic features of the DNA-LU-79553 complex on the binding specificity and on the antitumor activity of bisnaphthalimide agents are discussed.

The folding of centromeric DNA strands into intercalated structures: a physicochemical and computational study

We have carried out a physicochemical and computational analysis on the stability of the intercalated structures formed by cytosine-rich DNA strands. In the computational study, the electrostatic energy components have been calculated using a Poisson-Boltzmann model, and the non-polar energy components have been computed with a van der Waals function and/or a term dependent on the solvent-accessible surface area of the molecules. The results have been compared with those obtained for Watson-Crick duplexes and with thermodynamic data derived from UV experiments. We have found that intercalated DNA is mainly stabilized by very favorable electrostatic interactions between hydrogen-bonded protonated and neutral cytosines, and by non-polar forces including the hydrophobic effect and enhanced van der Waals contacts. Cytosine protonation electrostatically promotes the association of DNA strands into a tetrameric structure. The electrostatic interactions between stacked C.C+ pairs are strongly attenuated by the reaction field of the solvent, and are modulated by a complex interplay of geometric and protonation factors. The forces stabilizing intercalated DNA must offset an entropic penalty due to the uptake of protons for cytosine protonation, at neutral pH, and also the electrostatic contribution to the solvation free energy. The latter energy component is less favorable for protonated DNA due to the partial neutralization of the negative charge of the molecule, and probably affects other protonated DNA and RNA structures such as C+-containing triplexes.

Stacking interactions and intercalative DNA binding

The DNA-binding properties of many ligands can be rationalized on the basis of their structural and electronic complementarity with the functional groups present in the minor and major grooves of particular DNA sequences. Specific hydrogen bonding patterns are particularly useful for the purpose of sequence recognition. Less obvious, however, is the influence of base composition on the conformational preferences of individual base steps and on the binding of intercalating moieties which become sandwiched between contiguous base pairs. Improved knowledge of stacking interactions may lead to a better understanding of the architecture and inherent flexibility of particular DNA sequences and may provide insight into the principles that dictate the structural changes and specificity patterns observed in the binding of some intercalating ligands to DNA.

Centromeric pyrimidine strands fold into an intercalated motif by forming a double hairpin with a novel T:G:G:T tetrad: solution structure of the d(TCCCGTTTCCA) dimer

The solution structures of the oligodeoxynucleotides d(CCCGTTTCC) and d(TCCCGTTTCCA) have been determined by two-dimensional NMR spectroscopy. These oligomers are part of a DNA box in human centromeric alpha satellite targeted by the centromere protein B (CENP-B). Both CENP-B and its recognition box in alphoid DNA are conserved in mammals, suggesting an important biological role. At acidic pH, d(CCCGTTTCC), d(TCCCGTTTCCA) and the full d(TCCCGTTTCCAACGAAG) CENP-B box strand all fold and dimerize in solution forming a stable bimolecular structure containing two GTTT hairpin loops that interact through a novel T : G : G : T tetrad. The stem region of the dimer is a four-stranded intercalated motif in which the hairpin monomers are parallel and held together by C : C+ hydrogen-bonding and intercalation. The loops are at the same end of the dimer and lie across the narrow grooves of the tetraplex. They are remarkably structured and stabilized by base-base cross-stacking, sugar-base stacking, and parallel G:G and antiparallel G:T pairing. In the d(TCCCGTTTCCA)2 structure, the intercalated motif is continued at the other end of the dimer with unpaired but stacked adenine and thymine bases. The possible biological implications of these structures are discussed.

Structure-affinity relationships for the binding of actinomycin D to DNA

Molecular models of the complexes between actinomycin D and 14 different DNA hexamers were built based on the X-ray crystal structure of the actinomycin-d(GAAGCTTC)2 complex. The DNA sequences included the canonical GpC binding step flanked by different base pairs, nonclassical binding sites such as GpG and GpT, and sites containing 2,6-diamino-purine. A good correlation was found between the intermolecular interaction energies calculated for the refined complexes and the relative preferences of actinomycin binding to standard and modified DNA. A detailed energy decomposition into van der Waals and electrostatic components for the interactions between the DNA base pairs and either the chromophore or the peptidic part of the antibiotic was performed for each complex. The resulting energy matrix was then subjected to principal component analysis, which showed that actinomycin D discriminates among different DNA sequences by an interplay of hydrogen bonding and stacking interactions. The structure-affinity relationships for this important antitumor drug are thus rationalized and may be used to advantage in design of novel sequence-specific DNA-binding agents.

Molecular dynamics simulations of the bis-intercalated complexes of ditercalinium and Flexi-Di with the hexanucleotide d(GCGCGC)2: theoretical analysis of the interaction and rationale for the sequence binding specificity

The X-ray crystal structures of the complexes of ditercalinium and Flexi-Di with d(CGCG)2 have been studied by computational chemistry methods in an attempt to rationalize their distinct structural features. In addition, the complexes of these two bisintercalating drugs with d(GCGCGC)2 have been modeled and subjected to 0.5 ns of molecular dynamics simulations in explicit solvent with the aim of evaluating the relative importance of hydrogen bonding and stacking interactions in the sequence binding specificity of these compounds. According to our calculations, the electrostatic term is attractive for the stacking interactions between the pyridocarbazole chromophores of these drugs and the base pairs that make up the sandwiched GpC step. On the contrary, this energy term is repulsive for the base pairs that make up the boundaries of the bisintercalation site. This differential electrostatic binding energy component, which is shown to have a strong orientational dependence, could lie at the origin of the observed binding preferences of these drugs. In addition, both the Lennard-Jones and the electrostatic energy terms contribute to stabilizing the underwound central GpC step. The attractive electrostatic interactions between the linkers and the major groove are in concert with the stacking specificities for the sandwiched GpC step, which is thus very effectively stapled by the drugs. The hydrogen-bonding potential of the linkers, however, appears to be reduced in an aqueous medium due to competing interactions with water. Binding of either ditercalinium or Flexi-Di to d(GCGCGC)2 appears to favor the A-type conformation that this DNA molecule most likely adopts in the free state. The possible relevance of these findings to the process of bis-intercalation and to the pharmacological action of these compounds is discussed.

UV light-induced crosslinking of short DNA duplex strands: nucleotide sequence preferences and a prominent role of the duplex ends

By various doses of UV light, we irradiated 32 short DNA duplexes having between 12 and 40 nucleotide residues in length, and monitored the induced crosslinks between the complementary DNA strands by denaturing polyacrylamide gel electrophoresis. The experiments revealed that the crosslinking was strongest with the alternating sequence of T and A and weaker with the alternating sequence of T and G (C and A in the complementary strand). On the other hand, GC blocks of any sequence provided undetectable amounts of interstrand crosslinks even at the highest doses of UV irradiation. The amount of crosslinked strands logarithmically increased with the UV dose but it did not depend on the oligonucleotide concentration, ionic strength, divalent magnesium or manganese cations and pH at least within the examined regions of the experimental conditions, unless the oligonucleotide denatured or isomerized into a unimolecular foldback. The extent of crosslinking also did not depend on the (dT-dA)n duplex length to indicate that the crosslink was predominantly localized at a specific duplex locus. Experiments with (dT-dA)8 "mutants" in which AT pairs were systematically replaced by GC pairs at various molecule positions, revealed that the crosslinking predominantly occurred at the oligo(dT-dA) duplex ends. The crosslinking is a direct method to detect duplexes of DNA, which is here, for example, demonstrated with the heteroduplex of (dT-dA)12 and (dT-dA)16.