A molecular dynamics study of the bis-intercalation complexes of echinomycin with d(ACGT)2 and d(TCGA)2: rationale for sequence-specific Hoogsteen base pairing

Hoogsteen base pairs increase the susceptibility of double-stranded DNA to cytotoxic damage

As the Watson-Crick faces of nucleobases are protected in dsDNA, it is commonly assumed that deleterious alkylation damage to the Watson-Crick faces of nucleobases predominantly occurs when DNA becomes single-stranded during replication and transcription. However, damage to the Watson-Crick faces of nucleobases has been reported in dsDNA in vitro through mechanisms that are not understood. In addition, the extent of protection from methylation damage conferred by dsDNA relative to ssDNA has not been quantified. Watson-Crick base pairs in dsDNA exist in dynamic equilibrium with Hoogsteen base pairs that expose the Watson-Crick faces of purine nucleobases to solvent. Whether this can influence the damage susceptibility of dsDNA remains unknown. Using dot-blot and primer extension assays, we measured the susceptibility of adenine-N1 to methylation by dimethyl sulfate (DMS) when in an A-T Watson-Crick versus Hoogsteen conformation. Relative to unpaired adenines in a bulge, Watson-Crick A-T base pairs in dsDNA only conferred ∼130-fold protection against adenine-N1 methylation, and this protection was reduced to ∼40-fold for A(syn)-T Hoogsteen base pairs embedded in a DNA-drug complex. Our results indicate that Watson-Crick faces of nucleobases are accessible to alkylating agents in canonical dsDNA and that Hoogsteen base pairs increase this accessibility. Given the higher abundance of dsDNA relative to ssDNA, these results suggest that dsDNA could be a substantial source of cytotoxic damage. The work establishes DMS probing as a method for characterizing A(syn)-T Hoogsteen base pairs in vitro and also lays the foundation for a sequencing approach to map A(syn)-T Hoogsteen and unpaired adenines genome-wide in vivo.

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.

New insights into Hoogsteen base pairs in DNA duplexes from a structure-based survey

Hoogsteen (HG) base pairs (bps) provide an alternative pairing geometry to Watson-Crick (WC) bps and can play unique functional roles in duplex DNA. Here, we use structural features unique to HG bps (syn purine base, HG hydrogen bonds and constricted C1'-C1' distance across the bp) to search for HG bps in X-ray structures of DNA duplexes in the Protein Data Bank. The survey identifies 106 A•T and 34 G•C HG bps in DNA duplexes, many of which are undocumented in the literature. It also uncovers HG-like bps with syn purines lacking HG hydrogen bonds or constricted C1'-C1' distances that are analogous to conformations that have been proposed to populate the WC-to-HG transition pathway. The survey reveals HG preferences similar to those observed for transient HG bps in solution by nuclear magnetic resonance, including stronger preferences for A•T versus G•C bps, TA versus GG steps, and also suggests enrichment at terminal ends with a preference for 5'-purine. HG bps induce small local perturbations in neighboring bps and, surprisingly, a small but significant degree of DNA bending (∼14°) directed toward the major groove. The survey provides insights into the preferences and structural consequences of HG bps in duplex DNA.

Molecular simulations of drug–receptor complexes in anticancer research

Molecular modeling and computer simulation techniques have matured significantly in recent years and proved their value in the study of drug–DNA, drug–DNA–protein, drug–protein and protein–protein interactions. Evolution in this area has gone hand-in-hand with an increased availability of structural data on biological macromolecules, major advances in molecular mechanics force fields and considerable improvements in computer technologies, most significantly processing speeds, multiprocessor programming and data-storage capacity. The information derived from molecular simulations of drug–receptor complexes can be used to extract structural and energetic information that is usually beyond current experimental possibilities, provide independent accounts of experimentally observed behavior, help in the interpretation of biochemical or pharmacological results, and open new avenues for research by posing novel relevant questions that can guide the design of new experiments. As drug-screening tools, ligand- and fragment-docking platforms stand out as powerful techniques that can provide candidate molecules for hit and lead development. This review provides an overall perspective of the main methods and focuses on some selected applications to both classical and novel anticancer targets.

Molecular modeling of tricyclic compounds with anilino substituents and their intercalation complexes with DNA sequences

Although 9-anilinoacridines are among the best studied antitumoral intercalators, there are few studies about the effect of isosteric substitution of a benzene moiety for a heterocycle ring in the acridine framework. According to these studies, this approach may lead to effective cytotoxic agents, but good cytotoxic activity depends on structural requirements in the aniline ring which differ from those in 9-anilinoacridines. The present paper deals with molecular modeling studies of some 9-anilino substituted tricyclic compounds and their intercalation complexes (in various DNA sequences) resulting from docking the compounds into various DNA sequences. As expected, the isosteric substitution in 9-anilinoacridines influences the LUMO energy values and orbital distribution, the dipole moment, electrostatic charges and the conformation of the anilino ring. Other important differences are observed during the docking studies, for example, changes in the spatial arrangement of the tricyclic nucleus and the anilino ring at the intercalation site. Semiempirical calculations of the intercalation complexes show that the isosteric replacement of a benzene ring in the acridine nucleus affects not only DNA affinity but also base pair selectivity. These findings explain, at least partially, the different structural requirements observed in several 9-anilino substituted tricyclic compounds for cytotoxic activity. Thus, the data presented here may guide the rational design of new agents with different DNA binding properties and/or a cytotoxic profile by isosteric substitution of known intercalators.

Antitumor Activity, X-ray Crystal Structure, and DNA Binding Properties of Thiocoraline A, a Natural Bisintercalating Thiodepsipeptide

The marine natural product thiocoraline A displayed approximately equal cytotoxic activity at nanomolar concentrations in a panel of 12 human cancer cell lines. X-ray diffraction analyses of orthorhombic crystals of this DNA-binding drug revealed arrays of docked pairs of staple-shaped molecules in which one pendent hydroxyquinoline chromophore from each cysteine-rich molecule appears intercalated between the two chromophores of a facing molecule. This arrangement is in contrast to the proposed mode of binding to DNA that shows the two drug chromophores clamping two stacked base pairs, in agreement with the nearest-neighbor exclusion principle. Proof of DNA sequence recognition was obtained from both classical DNase I footprinting experiments and determination of the melting temperatures of several custom-designed fluorescently labeled oligonucleotides. A rationale for the DNA-binding behavior was gained when models of thiocoraline clamping a central step embedded in several octanucleotides were built and studied by means of unrestrained molecular dynamics simulations in aqueous solution.

RNA unrestrained molecular dynamics ensemble improves agreement with experimental NMR data compared to single static structure: a test case

Nuclear magnetic resonance (NMR) provides structural and dynamic information reflecting an average, often non-linear, of multiple solution-state conformations. Therefore, a single optimized structure derived from NMR refinement may be misleading if the NMR data actually result from averaging of distinct conformers. It is hypothesized that a conformational ensemble generated by a valid molecular dynamics (MD) simulation should be able to improve agreement with the NMR data set compared with the single optimized starting structure. Using a model system consisting of two sequence-related self-complementary ribonucleotide octamers for which NMR data was available, 0.3 ns particle mesh Ewald MD simulations were performed in the AMBER force field in the presence of explicit water and counterions. Agreement of the averaged properties of the molecular dynamics ensembles with NMR data such as homonuclear proton nuclear Overhauser effect (NOE)-based distance constraints, homonuclear proton and heteronuclear (1)H-(31)P coupling constant (J) data, and qualitative NMR information on hydrogen bond occupancy, was systematically assessed. Despite the short length of the simulation, the ensemble generated from it agreed with the NMR experimental constraints more completely than the single optimized NMR structure. This suggests that short unrestrained MD simulations may be of utility in interpreting NMR results. As expected, a 0.5 ns simulation utilizing a distance dependent dielectric did not improve agreement with the NMR data, consistent with its inferior exploration of conformational space as assessed by 2-D RMSD plots. Thus, ability to rapidly improve agreement with NMR constraints may be a sensitive diagnostic of the MD methods themselves.

Structure-activity relationships for the anti-HIV activity of flavonoids

Quantitative structure-activity relationship (QSAR) models are useful in providing a biochemical understanding of the biological activity of natural and synthetic chemicals based solely on molecular structure. One- to three-parameter multiregression equations were generated for three different groups of flavonoids to model experimental flavonoid-induced cytotoxicity and inhibition of human immunodeficiency virus I replication in lymphocyte-infected cells, using quantum chemical and geometrical parameters derived from optimized molecular structures. Both biological properties were basically dependent on electronic parameters describing charge distribution on the two fused rings of the flavonoid molecule. Atomic charges in C3 and the carbonyl carbon as well as the dipolar moment were important electronic descriptors to define the studied biological properties of flavonoids.

Rebek imides and their adenine complexes: preferences for Hoogsteen binding in the solid state and in solution

Rebek imides (3), formed from Kemp's triacid, were developed in the mid-1980's as model receptors for adenine derivatives. We report here the first account of their hydrogen-bonding preferences upon binding 9-ethyladenine (1a) in the solid state. Structural analysis begins with simple imides 3a-e that form discrete dimers, while bis-imide 4 forms ribbon-like structures in the crystalline phase. The hydrogen-bonding interface within each of the representative assemblies features short intermolecular N(3)imide...O(8*)imide* distances (ca. 2.95 A), indicative of two-point hydrogen bonding. Imides 3f-h could be co-crystallized with 1a; single-crystal X-ray analysis of the resulting complexes reveals hydrogen-bonding geometries nearly identical to those observed in nucleobase complexes of adenine and pyrimidine derivatives. Imides 3f and 3g form 2:1 ternary assemblies with 1a; the complex of the former, (3f)2 x 1a, displays both Watson-Crick- and Hoogsteen-type hydrogen bonding, whereas the complex of the latter, (3g)2 x 1a, shows the Hoogsteen motif and imide hydrogen bonding to N(3) of the purine base (N(3)adenine...N(3'')imide = 3.07(1) A). Imide 3h forms a 1:1 complex with 1a (3h x 1a x CHCl3) and displays Hoogsteen binding exclusively. All of the 3 x 1a assemblies show C(adenine)...O(imide) distances (3.38-3.75 A) that are consistent with C-H...O hydrogen bonding. Base-pairing preferences for the Rebek imides are further explored in solution by 1H NMR one-dimensional NOE experiments and by computational means; in all cases the Hoogsteen motif is modestly favored relative to its Watson-Crick counterpart.

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.

DNA recognition by quinoxaline antibiotics: use of base-modified DNA molecules to investigate determinants of sequence-specific binding of triostin A and TANDEM

The methodology of DNAase I footprinting has been adapted to investigate the sequence-specific binding of two quinoxaline drugs to DNA fragments containing natural and modified bases. In order to help comprehend the molecular origin of selectivity in the bis-intercalation of triostin A and TANDEM at CpG and TpA sites respectively, we have specifically examined the effect of the 2-amino group of guanine on their sequence specificity by using DNA in which that group has been either removed from guanine, added to adenine or both. Previous studies suggested that the recognition of particular nucleotide sequences by these drugs might be dependent upon the placement of the purine 2-amino group, serving as a positive or a negative effector for triostin A and TANDEM respectively. However, the footprinting data reported here indicate that this is not entirely correct, since they show that the 2-amino group of guanine is important for the binding of triostin A to DNA but has absolutely no influence on the interaction of TANDEM with TpA steps. Apparently the binding of triostin A to CpG sites is primarily due to hydrogen bonding interaction between the cyclic peptide of the antibiotic and the 2-amino group of guanine residues, whereas the selective binding of TANDEM to TpA sites is not hydrogen-bond driven and probably originates mainly from steric and/or hydrophobic interactions, perhaps involving indirect recognition of a suitable minor groove structure.

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.

Dissociation kinetics of echinomycin from CpG binding sites in different sequence environments

We have examined the kinetics of dissociation of echinomycin from CpG sites in several DNA fragments containing synthetic DNA inserts by a variation of the footprinting technique. Complexes of the ligand with radiolabeled DNA fragments were dissociated by adding an excess of unlabeled calf thymus DNA. Samples were removed from this mixture at subsequent time intervals and subjected to DNase I footprinting. The rate of disappearance of the footprints varied considerably between the various CpG sites. At 20 degrees C, echinomycin dissociates more slowly from CpG sites flanked by (AT)n (t1/2 approximately 40 min) and (CA)n.(TG)n (t1/2 approximately 11 min) than by An.Tn (t1/2 < 3 min). In each sequence context the dissociation from ACGT is slower than that from TCGA. (TAA)4CG(TTA)4 also represents a very good binding site (t1/2 approximately 35 min), which is less sensitive to changes in temperature than most other sites. Within sequences (AT)10(G/C)4(AT)10, the dissociation from CGGC is slower than that from CCCG or CCGC.

Site selective bis-intercalation of a homodimeric thiazole orange dye in DNA oligonucleotides

We have used one and two dimensional 1H NMR spectroscopy to characterize the binding of a homodimeric thiazole orange dye, 1,1'-(4,4,8,8-tetramethyl-4,8-diaza-undecamethylene)-bis-4- (3-methyl-2,3-dihydro-(benzo-1,3-thiazole)-2-methylidene)-quinolin ium tetraiodide (TOTO), to various double stranded DNA oligonucleotides. TOTO binds strongly to all the oligonucleotides used, but usually more than one complex is observed and exchange between different binding sites broadens the lines in the NMR spectra. Complete precipitation occurs when TOTO is bound to small oligonucleotides. Binding to larger oligonucleotides occurs by bis-intercalation. The 1:1 complex of TOTO with the oligonucleotide d(CCGACTGATGC):d (GCATCAGTCGG) gave only one complex that was shown to be a bis-intercalation in the CTGA:TCAG binding site. The binding to this site was also characterized by studying the TOTO complex with the d(CCGCTGAGC):d(GCTCAGCGG) oligonucleotide. NOE connectivities and molecular modelling were used to characterize the complex. The 1:1 complex of TOTO with the oligonucleotide d(CCGCTAGCG):d(CGCTAGCGG) containing a CTAG:CTAG binding site was similarly characterized by NMR. It was concluded that the binding of TOTO to larger oligonucleotides is site selective with CTAG:CTAG as the preferred binding site.

NMR investigation of Hoogsteen base pairing in quinoxaline antibiotic--DNA complexes: comparison of 2:1 echinomycin, triostin A and [N-MeCys3,N-MeCys7] TANDEM complexes with DNA oligonucleotides

Hoogsteen base pairs have been demonstrated to occur in base pairs adjacent to the CpG binding sites in complexes of triostin A and echinomycin with a variety of DNA oligonucleotides. To understand the relationship of these unusual base pairs to the sequence specificity of these quinoxaline antibiotics, the conformation of the base pairs flanking the YpR binding sites of the 2:1 drug-DNA complexes of triostin A with [d(ACGTACGT)]2 and of the TpA specific [N-MeCys3, N-MeCys7] TANDEM with [d(ATACGTAT)]2 have been studied by 1H NMR spectroscopy. In both the 2:1 triostin A-DNA complex and the 2:1 [N-MeCys3, N-MeCys7] TANDEM-DNA complex, the terminal A.T base pairs are Hoogsteen base paired with the 5' adenine in the syn conformation. This indicates that both TpA specific and CpG specific quinoxaline antibiotics are capable of inducing Hoogsteen base pairs in DNA. However, in both 2:1 complexes, Hoogsteen base pairing is limited to the terminal base pairs. In the 2:1 triostin A complex, the internal adenines are anti and in the 2:1 [N-MeCys3, N-MeCys7] TANDEM-DNA complex, the internal guanines are anti regardless of pH, which indicates that the central base pairs of both complexes form Watson-Crick base pairs. This indicates that the sequence dependent nature of Hoogsteen base pairing is the same in TpA specific and CpG specific quinoxaline antibiotic-DNA complexes. We have calculated a low resolution three-dimensional structure of the 2triostin A-[d(ACGTACGT)]2 complex and compared it with other CpG specific quinoxaline antibiotic-DNA complexes. The role of stacking in the formation of Hoogsteen base pairs in these complexes is discussed.

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

Molecular dynamics simulations have been used to explore the behavior of the complexes of echinomycin with the DNA tetramers d(GCGC)2 and d(CCGG)2 in which the terminal bases have been paired according to either a Hoogsteen or a Watson-Crick hydrogen bonding scheme. The energy of the four resulting complexes has been monitored along the dynamics trajectories and the interaction energy between echinomycin and DNA has been decomposed into contributions arising from the planar aromatic systems and the depsipeptide part of the antibiotic. Our calculations predict a large increase in overall stabilization upon protonation of the terminal cytosines and subsequent Hoogsteen pair formation in the complex of echinomycin with d(GCGC)2 but not with d(CCGG)2, in agreement with the experimental evidence [Gao and Patel, Quart. Rev. Biophys. 22, 93-138 (1989)]. The conformational preferences appear to arise mainly from differential stacking interactions in which the electrostatic component is shown to play a dominant role. Differences in hydrogen bonding patterns are also found among the complexes and these are compared in relation to available crystal structures. The binding of echinomycin to DNA appears as a complex process involving many interrelated variables.