1H-NMR studies of a monointercalating drug into a d(CpGpApTpCpG)2 minihelix

The structure of the complex formed between the 7H-pyridocarbazole monomer [[(2-piperidyl)-2,1-ethane-yl] [10-methoxy-7H-pyrido[4,3-c]carbazolium] dimethane sulfonate] and the autocomplementary hexanucleotide d(CpGpApTpCpG)2 in aqueous solution is analyzed by 270- and 400-MHz 1H-nmr. The large upfield shifts observed for both the drug and the self-complementary hexanucleotide protons provide evidence for intercalated complexes. The observation of intermolecular nuclear Overhauser effects between drug and the hexanucleotide protons gives a privileged orientation of the drug in the intercalation site with the quaternarizing ethyl piperidine chain protruding in the major groove. Moreover, the data suggest an intercalation based on the neighbor exclusion site principle in the three alternating sequences.     

Minor groove binding of a bis-quaternary ammonium compound: the crystal structure of SN 7167 bound to d(CGCGAATTCGCG)2.

The X-ray crystal structure of the complex between the synthetic antitumour and antiviral DNA binding ligand SN 7167 and the DNA oligonucleotide d(CGCGAATTCGCG)2 has been determined to an R factor of 18.3% at 2.6 A resolution. The ligand is located within the minor groove and covers almost 6 bp with the 1-methylpyridinium ring extending as far as the C9-G16 base pair and the 1-methylquinolinium ring lying between the G4-C21 and A5-T20 base pairs. The ligand interacts only weakly with the DNA, as evidenced by long range contacts and shallow penetration into the groove. This structure is compared with that of the complex between the parent compound SN 6999 and the alkylated DNA sequence d(CGC[e6G]AATTCGCG)2. There are significant differences between the two structures in the extent of DNA bending, ligand conformation and groove binding.     

Interactions of chromomycin A3 and mithramycin with the sequence d(TAGCTAGCTA)2

Anti-cancer antibiotics, chromomycin A3 (CHR) and mithramycin (MTR) inhibit DNA directed RNA synthesis in vivo by binding reversibly to template DNA in the minor groove with GC base specificity, in the presence of divalent cations like Mg2+. Under physiological conditions, (drug)2Mg2+ complexes formed by the antibiotics are the potential DNA binding ligands. Structures of CHR and MTR differ in their saccharide residues. Scrutiny of the DNA binding properties reveal significant differences in their sequence selectivity, orientation and stoichiometry of binding. Here, we have analyzed binding and thermodynamic parameters for the interaction of the antibiotics with a model oligonucleotide sequence, d(TAGCTAGCTA)2 to understand the role of sugars. The oligomer contains two potential binding sites (GpC) for the ligands. The study illustrates that the drugs bind differently to the sequence. (MTR)2Mg2+ binds to both sites whereas (CHR)2Mg2+ binds to a single site. UV melting profiles for the decanucleotide saturated with the ligands show that MTR bound oligomer is highly stabilized and melts symmetrically. In contrast, with CHR, loss of symmetry in the oligomer following its association with a single (CHR)2Mg2+ complex molecule leads to a biphasic melting curve. Results have been interpreted in the light of saccharide dependent differences in ligand flexibility between the two antibiotics.     

Structural analysis of the binding modes of minor groove ligands comprised of disubstituted benzenes

Two-dimensional homonuclear NMR was used to characterize synthetic DNA minor groove-binding ligands in complexes with oligonucleotides containing three different A-T binding sites. The three ligands studied have a C₂ axis of symmetry and have the same general structural motif of a central para-substituted benzene ring flanked by two meta-substituted rings, giving the molecules a crescent shape. As with other ligands of this shape, specificity seems to arise from a tight fit in the narrow minor groove of the preferred A-T-rich sequences. We found that these ligands slide between binding subsites, behavior attributed to the fact that all of the amide protons in the ligand backbone cannot hydrogen bond to the minor groove simultaneously.     

Interaction of chromium(III) complex of chiral binaphthyl tetradentate ligand with DNA

Since conformation of the molecule plays a vital role in the activity of drug, we have investigated the DNA interaction of a chromium(III) complex with ligands in two conformations. Chromium(III) complexes derived from chiral binaphthyl Schiff base ligands, viz. R- and S-2,2'-bis(salicylideneamino) 1,1'-binaphthyl, have been synthesized and characterized by mass, IR, and electronic spectra. The interaction of these R- and S-binaphthyl Schiff base chromium(III) complexes with CT-DNA was investigated with the goal of examining whether the chirality has an influence on the chromium(III)-DNA binding properties. The difference in chirality of the ligand did not show any striking difference in binding properties. The binding constants for R and S conformers were estimated to be 18 (+/-0.4) x 10(3) and 9.4 (+/-0.3) x 10(3) M(-1), respectively, through spectroscopic titrations. All the experimental results are suggestive that both the isomers are DNA groove binders. The results of steady-state as well as time-resolved fluorescence experiments, however, suggest that the R conformer has restricted mobility when bound to DNA because it is more deeply buried in the groove of DNA compared to the S isomer.     

(+)-CC-1065 as a structural probe of Mu transposase-induced bending of DNA: overcoming limitations of hydroxyl-radical footprinting.

Phage Mu transposase (A-protein) is primarily responsible for transposition of the Mu genome. The protein binds to six att sites, three at each end of Mu DNA. At most att sites interaction of a protein monomer with DNA is seen to occur over three minor and two consecutive major grooves and to result in bending up to about 90 degrees. To probe the directionality and locus of these A-protein-induced bends, we have used the antitumor antibiotic (+)-CC-1065 as a structural probe. As a consequence of binding within the minor groove, (+)-CC-1065 is able to alkylate N3 of adenine in a sequence selective manner. This selectivity is partially determined by conformational flexibility of the DNA sequence, and the covalent adduct has a bent DNA structure in which narrowing of the minor groove has occurred. Using this drug in experiments in which either gel retardation or DNA strand breakage are used to monitor the stability of the A-protein--DNA complex or the (+)-CC-1065 alkylation sites on DNA (att site L3), we have demonstrated that of the three minor grooves implicated in the interaction with A-protein, the peripheral two are 'open' or accessible to drug bonding following protein binding. These drug-bonding sites very likely represent binding at at least two A-protein-induced bending sites. Significantly, the locus of bending at these sites is spaced approximately two helical turns apart, and the bending is proposed to occur by narrowing of the minor groove of DNA. The intervening minor groove between these two peripheral sites is protected from (+)-CC-1065 alkylation. The results are discussed in reference to a proposed model for overall DNA bending in the A-protein att L3 site complex. This study illustrates the utility of (+)-CC-1065 as a probe for protein-induced bending of DNA, as well as for interactions of minor groove DNA bending proteins with DNA which may be masked in hydroxyl radical footprinting experiments.     

Probing site specificity of DNA binding metallointercalators by NMR spectroscopy and molecular modeling

The molecular recognition of oligonucleotides by chiral ruthenium complexes has been probed by NMR spectroscopy using the template Delta-cis-alpha- and Delta-cis-beta-[Ru(RR-picchxnMe(2)) (bidentate)](2+), where the bidentate ligand is one of phen (1,10-phenanthroline), dpq (dipyrido[3,2-f:2',3'-h]quinoxaline), or phi (9,10-phenanthrenequinone diimine) and picchxnMe(2)() is N,N'-dimethyl-N,N'-di(2-picolyl)-1,2-diaminocyclohexane. By varying only the bidentate ligand in a series of complexes, it was shown that the bidentate alone can alter binding modes. DNA binding studies of the Delta-cis-alpha-[Ru(RR-picchxnMe(2))(phen)](2+) complex indicate fast exchange kinetics on the chemical shift time scale and a "partial intercalation" mode of binding. This complex binds to [d(CGCGATCGCG)](2) and [d(ATATCGATAT)](2) at AT, TA, and GA sites from the minor groove, as well as to the ends of the oligonucleotide at low temperature. Studies of the Delta-cis-beta-[Ru(RR-picchxnMe(2))(phen)](2+) complex with [d(CGCGATCGCG)](2) showed that the complex binds only weakly to the ends of the oligonucleotide. The interaction of Delta-cis-alpha-[Ru(RR-picchxnMe(2))(dpq)](2+) with [d(CGCGATCGCG)](2) showed intermediate exchange kinetics and evidence of minor groove intercalation at the GA base step. In contrast to the phen and dpq complexes, Delta-cis-alpha- and Delta-cis-beta-[Ru(RR-picchxnMe(2))(phi)](2+) showed evidence of major groove binding independent of the metal ion configuration. DNA stabilization induced by complex binding to [d(CGCGATCGCG)](2) (measured as DeltaT(m)) increases in the order phen < dpq and DNA affinity in the order phen < dpq < phi. The groove binding preferences exhibited by the different bidentate ligands is explained with the aid of molecular modeling experiments.     

Molecular dynamics of an in vacuo model of duplex d(CGCGAATTCGCG) in the B-form based on the amber 3.0 force field.

The characteristics of 100 ps of molecular dynamics (MD) on the DNA dodecamer d(CGCGAATTCGCG) at 300 K are described and investigated. The simulation is based on an in vacuo model of the oligomer and the AMBER 3.0 force field configured in the manner of Singh, U. C., S. J. Weiner, and P. A. Kollman, (1985, Proc. Natl. Acad. Sci. USA. 82:755-759). The analysis of the results was carried out using the "curves, dials, and windows" procedure (Ravishanker, G., S. Swaminathan, D. L. Beveridge, R. Lavery, and H. Sklenar. 1989. J. Biomol. Struct. Dyn. 6:669-699). The results indicate this dynamical model to be a provisionally stable double helix which lies at approximately 3.2 A rms deviation from the canonical B-form. There is, however, a persistent nonplanarity in the base pair orientations which resemble that observed in canonical A-DNA. The major groove width is seen to narrow during the course of the simulation and the minor groove expands, contravariant to the alterations in groove width seen in the crystal structure of the native dodecamer (Drew, H. R., R. M. Wing, T. Takano, C. Broka, S. Tanaka, I. Itakura, and R. E. Dickerson, 1981. Proc. Natl. Acad. Sci. USA. 78:2179-2183). The propeller twist in the bases, the sequence dependence of the base pair roll and aspects of bending in the helix axis are in some degree of agreement with the crystal structure. The patterns in DNA bending are observed to follow Zhurkin theory (Zhurkin, V. B. 1985. J. Biomol. Struct. Dyn. 2:785-804.). The relationship between the dynamical model and structure in solution is discussed.     

Manipulative interplay of two adozelesin molecules with d(ATTAAT)₂achieving ligand-stacked Watson-Crick and Hoogsteen base-paired duplex adducts

Previous structural studies of the cyclopropapyrroloindole (CPI) antitumor antibiotics have shown that these ligands bind covalently edge-on into the minor groove of double-stranded DNA. Reversible covalent modification of the DNA via N3 of adenine occurs in a sequence-specific fashion. Early nuclear magnetic resonance and molecular modeling studies with both mono- and bis-alkylating ligands indicated that the ligands fit tightly within the minor groove, causing little distortion of the helix. In this study, we propose a new binding model for several of the CPI-based analogues, in which the aromatic secondary rings form π-stacked complexes within the minor groove. One of the adducts, formed with adozelesin and the d(ATTAAT)(2) sequence, also demonstrates the ability of these ligands to manipulate the DNA of the binding site, resulting in a Hoogsteen base-paired adduct. Although this type of base pairing has been previously observed with the bisfunctional CPI analogue bizelesin, this is the first time that such an observation has been made with a monoalkylating nondimeric analogue. Together, these results provide a new model for the design of CPI-based antitumor antibiotics, which also has a significant bearing on other structurally related and structurally unrelated minor groove-binding ligands. They indicate the dynamic nature of ligand-DNA interactions, demonstrating both DNA conformational flexibility and the ability of two DNA-bound ligands to interact to form stable covalent modified complexes.     

Mixed mode of ligand-DNA binding results in S-shaped binding curves

S-shaped binding curves often characterize interactions of ligands with nucleic acid molecules as analyzed by different physico-chemical and biophysical techniques. S-shaped experimental binding curves are usually interpreted as indicative of the positive cooperative interactions between the bound ligand molecules. This paper demonstrates that S-shaped binding curves may occur as a result of the "mixed mode" of DNA binding by the same ligand molecule. Mixed mode of the ligand-DNA binding can occur, for example, due to 1) isomerization or dimerization of the ligands in solution or on the DNA lattice, 2) their ability to intercalate the DNA and to bind it within the minor groove in different orientations. DNA-ligand complexes are characterized by the length of the ligand binding site on the DNA lattice (so-called "multiple-contact" model). We show here that if two or more complexes with different lengths of the ligand binding sites could be produced by the same ligand, the dependence of the concentration of the complex with the shorter length of binding site on the total concentration of ligand should be S-shaped. Our theoretical model is confirmed by comparison of the calculated and experimental CD binding curves for bis-netropsin binding to poly(dA-dT) poly(dA-dT). Bis-netropsin forms two types of DNA complexes due to its ability to interact with the DNA as monomers and trimers. Experimental S-shaped bis-netropsin-DNA binding curve is shown to be in good correlation with those calculated on the basis of our theoretical model. The present work provides new insight into the analysis of ligand-DNA binding curves.     

Flexible docking of DNA fragments and actinocin derivatives

We have applied molecular docking methods to systems containing nucleic acids as targets and biologically active substances as ligands. The complexes of DNA fragments and actinocin derivatives with different lengths of aminoalkyl side chains were obtained by molecular docking. It was observed that actinocin derivatives could form energetically favourable complexes with DNA both as intercalators and minor groove binders. It was shown that small changes in the binding energy (～1?kcal/mol) could result in complexes with substantially different structure. The complexes of actinocin derivatives and DNA fragments were stabilized by hydrogen bonding upon intercalation and minor groove binding. It was found that the change of solvent-accessible surface area upon binding of the actinocin derivative to DNA linear increased with the growth of methylene groups' number in ligand side chains. The solvation energy change upon binding of actinocin derivatives to DNA calculated by the WSAS method was favourable in the case of small uncharged ligands and unfavourable for positively charged ligands.     

Biophysical analysis of DNA modified by 1,2-diaminocyclohexane platinum(II) complexes.

Modification of DNA and double-stranded deoxyoligonucleotides with antitumour 1,2-diamino-cyclohexanedinitroplatinum(II) (Pt-dach) complexes was investigated with the aid of physico-chemical methods and chemical probes of nucleic acid conformation. The three Pt-dach complexes were used which differed in isomeric forms of the dach nonleaving ligand-Pt(1R,2R-dach), Pt(1S,2S-dach) and Pt(1R,2S-dach) complexes. The latter complex has lower antitumour activity than the other two Pt-dach complexes. Pt(1R,2S-dach) complex exhibits the slowest kinetics of its binding to DNA and of the conversion of monofunctional binding to bifunctional lesions. The anomalously slow electrophoretic mobility of multimers of the platinated and ligated oligomers suggests that bifunctional binding of Pt-dach complexes to a d(GG) site within double-stranded oligonucleotides induces bending of the oligomer. In addition, chemical probing of double-helical deoxyoligonucleotides modified by the Pt-dach complexes at the d(GG) sites reveals that Pt(1R,2S-dach) complex induces more extensive conformational changes in the oligomer than Pt(1R,2R-dach) and Pt(1S,2S-dach) complexes. It is proposed that different effects of the Pt-dach complexes on DNA observed in this work arise mainly from a steric crowding of the axially oriented cyclohexane ring in the DNA adduct of Pt(1R,2S-dach) complex.     

Stereoelectronic factors in the interaction with DNA of small aromatic molecules substituted with a short cationic chain: importance of the polarity of the aromatic system of the molecule.

We have performed a quantitative analysis of the interaction with DNA of several unfused aromatic compounds synthesized in our laboratory and substituted with one or two short cationic chains. These and similar literature compounds, for which DNA binding data are available, bind with DNA by partial intercalation of the aromatic system, groove interaction of the linker chain, and groove electrostatic interactions of the terminal cationic group. Several independent quantitative and qualitative approaches show consistently that the strength of the interaction of the aromatic unit of the molecule with DNA binding sites depends on the direction and magnitude of polarity of the aromatic system. The phenomenon is explained in terms of the greatest negative potential in the DNA grooves, a concept extensively elaborated by Pullman and Pullman [cf. Lavery, R. and Pullman, B. [(1985) J. Biomol. Struct. Dyn. 2, 1021-1032] and references therein]. Classical, fused-ring planar intercalators do not follow the polarity-DNA affinity correlation, presumably because the intercalative forces depend more strongly on polarizability than on polarity of the aromatic system.     

Mixed Mode of Ligand-DNA Binding Results in S-Shaped Binding Curves

Abstract

S-shaped binding curves often characterize interactions of ligands with nucleic acid molecules as analyzed by different physicochemical and biophysical techniques. S-shaped experimental binding curves are usually interpreted as indicative of the positive cooperative interactions between the bound ligand molecules. This paper demonstrates that S-shaped binding curves may occur as a result of the “mixed mode” of DNA binding by the same ligand molecule. Mixed mode of the ligand-DNA binding can occur, for example, due to 1) isomerization or dimerization of the ligands in solution or on the DNA lattice, 2) their ability to intercalate the DNA and to bind it within the minor groove in different orientations. DNA- ligand complexes are characterized by the length of the ligand binding site on the DNA lattice (so-called “multiple-contact” model). We show here that if two or more complexes with different lengths of the ligand binding sites could be produced by the same ligand, the dependence of the concentration of the complex with the shorter length of binding site on the total concentration of ligand should be S-shaped. Our theoretical model is confirmed by comparison of the calculated and experimental CD binding curves for bis-netropsin binding to poly(dA-dT) poly(dA-dT). Bis-netropsin forms two types of DNA complexes due to its ability to interact with the DNA as monomers and trimers. Experimental S-shaped bis-netropsin-DNA binding curve is shown to be in good correlation with those calculated on the basis of our theoretical model. The present work provides new insight into the analysis of ligand-DNA binding curves.     

Comparison of analyses of DNA curvature

Popular programs for characterizing DNA structure include Curves 5.1 (Lavery, R. and Sklenar, H., J. Biomol. Struct. Dyn. 6, 63-91, 1988; Lavery, R. and Sklenar, H., J. Biomol. Struct. Dyn. 6, 655-67, 1989) and Freehelix98 (Dickerson, R. E., Nucleic Acids Res. 26, 1906-1926, 1998), along with the more recent 3DNA (X. J. Lu, Z. Shakked and W. K. Olson., J. Mol. Biol. 300, 819-840 (2000). Given input of structural coordinates, all of these programs return values of the local helical parameters, such as roll, tilt, twist, etc. The first two programs also provide characterization of global curvature. Madbend (Strahs, D. and Schlick, T., J. Mol. Biol. 301, 643-663, 2000), a program that computes global curvature from local roll, tilt, and twist parameters, can be applied to the output of all three structural programs. We have compared the curvature predicted by the three programs with and without the use of Madbend. Global bend magnitudes and directions as well as values of helical kinks were calculated for four high-resolution DNA structures and four model DNA helices. Global curvature determined by Curves 5.1 without Madbend was found to differ from values obtained using Freehelix98 with or without Madbend or 3DNA and Curves 5.1 with Madbend. Using model helices, this difference was attributed the fact that Curves 5.1 is the only program sensitive to changes in axial displacement, such as shift and slide. Madbend produced robust values of bend magnitude and direction, and displayed little sensitivity to axis displacement or the source of local helical parameters. Madbend also appears to be the method of choice for bending comparisons of high-resolution structures with results from cyclization kinetics, a method that measures DNA curvature as a vectorial sum of local roll and tilt angles.     

The use of diaminopurine to investigate structural properties of nucleic acids and molecular recognition between ligands and DNA.

2,6-Diaminopurine (DAP) is an analogue of adenine which can be converted to nucleotides that serve as substrates for incorporation into nucleic acids by polymerases in place of (d)AMP. It pairs with thymidine (or uracil), engaging in three hydrogen bonds of the Watson-Crick type. The result of DAP incorporation is to add considerable stability to the double helix and to impart other structural features, such as an altered groove width and disruption of the normal spine of hydration. DNA containing DAP may or may not be recognized by restriction endonucleases; RNA containing DAP may not engage in normal splicing. The DAP.T pair affects the local flexibility of DNA and impedes the interaction with helix bending proteins. By providing a non-canonical hydrogen bond donor in the minor groove and/or blocking access to the floor of that groove it strongly affects interactions with small molecules such as antibiotics and anticancer drugs. Examples which illustrate altered recognition of nucleotide sequences in DAP-containing DNA are presented: changed sites of cutting by bleomycin, photocleavage by uranyl nitrate and footprinting with mithramycin. Using DNA in which both A-->DAP and G-->Inosine substitutions have been made it is possible to assess precisely the role of the purine 2-amino group in ligand-DNA recognition.     

Aryl substituted ruthenium bis-terpyridine complexes: intercalation and groove binding with DNA

The non-covalent interaction of five novel ruthenium(II) bis-terpyridine complexes with calf thymus DNA and, where appropriate, with poly[d(G-C)](2) and poly[d(A-T)](2) is described. Each complex is functionalised with aryl tail groups in the 4' position of the terpyridine ligands ((i) 9-anthracenyl, (ii) 4,4'-biphenyl, (iii) beta-naphthyl, (iv) 9-phenanthrenyl, and (v) 1-pyrenyl). Circular dichroism and linear dichroism show that the binding of three of the complexes (phenanthrenyl, anthracenyl and pyrenyl) at low metal complex concentration is dominated by intercalation of the aryl tail groups between the DNA bases. The complex with the biphenyl tail predominantly exhibits groove binding with no significant tail intercalation. The naphthyl derivative binds both by intercalation and a non-intercalative mode even at low metal complex concentrations. At high metal complex concentrations, aggregation of the complexes on the DNA is observed. Resonance light scattering indicates that the aggregates are of low nuclearity along the groove.