1H NMR study of [d(GCGATCGC)]2 and its interaction with minor groove binding 4',6-diamidino-2-phenylindole

Two-dimensional NMR spectroscopy has been applied to study the solution binding of 4',6-diamidino-2-phenylindole (DAPI) to synthetic DNA duplex [d(GCGATCGC)]2. The structure of the complex at a molar ratio of 1:1 drug:duplex has been investigated. NMR results indicate that DAPI binds selectively in the minor groove of the DNA region containing only two A:T base pairs. The results disagree with conclusions drawn from footprinting experiments and show that the presence of the G3NH2 group in the minor groove does not prevent the binding. A molecular model is proposed that closely resembles the crystal structure previously published for the interaction of DAPI with the dodecamer [d(CGCGAATTCGCG)]2, containing four A:T base pairs in the binding site. In this model, DAPI lies in the minor groove, nearly isohelical, with its aromatic rings adjacent to H4' protons of T5 and C6 deoxyribose and the NH indole group oriented toward the DNA axis. The binding does not perturb the B-type conformation of the duplex, and the DNA oligomer conserves its 2-fold symmetry, indicating that fast exchange dynamics exist between the two stereochemically equivalent binding sites of the palindromic sequence. The binding constant and the exchange rate between free and bound species were also measured by NMR spectroscopy.     

Improved targeting of an anti-epidermal growth factor receptor variant III monoclonal antibody in tumor xenografts after labeling using N-succinimidyl 5-iodo-3-pyridinecarboxylate

The crystal structure is reported of a complex between the dodecanucleotide sequence d(CGCGAATTCGCG)2and an analogue of the DNA binding drug Hoechst 33258, in which the piperazine ring has been replaced by an amidinium group and the phenol ring by a phenylamidinium group. The structure has been refined to an R factor of 19.5% at 2.2 A resolution. The drug is held in the minor groove by five strong hydrogen bonds, together with bridging water molecules at both ends. There are few other contacts with the floor of the groove, indicating a lack of isohelicity with the groove and suggesting (i) that the observed high DNA affinity of this drug is primarily due to the array of hydrogen bonds and (ii) that these more than compensate for its poor isohelicity.     

Sequence specific modulation of DNA restriction enzyme cleavage by minor groove binders

The inhibition of restriction endonuclease cleavage by a series of bisquaternary ammonium derivatives (BQA-derivatives) which bind to the minor groove of DNA has been studied. The derivatives considered included six sequence-selective binders (SN 6570, SN 6999, SN 6050, SN 6132, SN 6131 and SN 18071) and four non-specific binders (SN 6113, SN 5754, SN 6324 and SN 4094) and can be distinguished by their activity on restriction endonucleases. Digestion experiments with pUC19 DNA were monitored electrophoretically using the transition of the covalently closed circular (ccc) DNA into the linear double stranded (lds) one. Only the sequence-specific binders inhibit the cleavage activity of restriction endonucleases EcoRI, SspI and DraI with four and six dAdT-base pairs within their restriction sites, while the activity of SalI and BamHI with less than four dAdT-sequences was unaffected. In contrast, the non-specific binding ligands were incapable of suppressing enzyme digestion. The inhibition of the restriction endonuclease PvuII indicates that ligand binding in close vicinity to the cleavage sites is also involved in the enzyme inhibition. The dAdT-content in proximity to the palindromic sequences of three DraI cutting sites in pUC19 DNA explains why the derivative SN 6053 protects these sequences in different manners. Gel shift experiments indicated that BQA-derivatives inhibit the DNA-enzyme complex formation if the ligand was added to the DNA before the enzyme. In contrast, complex formation between DNA and enzyme remained unchanged when the enzyme was added first.     

Differential DNA recognition by the enantiomers of 1-Rh(MGP)2 phi: a combination of shape selection and direct readout

The enantiomers of the symmetric metallointercalator complex 1-Rh(MGP)2phi5+ [MGP = 4-(guanidylmethyl)-1,10-phenanthroline; phi = phenanthrenequinone diimine] bound to DNA decamer duplexes containing their respective 6 bp recognition sequences have been investigated using 1H NMR. Shape selection due to the chirality of the metal center and hydrogen-bonding contacts of ancillary guanidinium groups to 3'-G N7 atoms define the recognition by complexes which bind by intercalation to duplex DNA. The titration of Lambda-Rh into the self-complementary decamer containing the recognition sequence (5'-GACATATGTC-3', L1) resulted in one symmetric bound conformation observed in the 1H NMR spectrum, indicating that the DNA duplex retains its symmetry in the presence of the metal complex. Upfield chemical shifts of duplex imino protons and the disruption of the NOE base-sugar contacts defined the central T5-A6 intercalation site. The downfield shift of the G8 imino proton supports the conclusion that the pendant guanidinium arms make simultaneous H-bonding contacts to the N7 atoms of 3'-G8 bases on either side of the site. A variable-temperature study of a partially titrated sample (2:3 Lambda-Rh/L1) showed the exchange rate (kobs) at 298 K to be 68 s-1 and the activation barrier to exchange (DeltaG of association) to be 2.7 kcal/mol, a value comparable to the stacking energy of one base step. The results presented coupled with biochemical data are therefore consistent with binding models in which Lambda-1-Rh(MGP)2phi5+ (Lambda-Rh) traps the recognition site 5'-CATATG-3' in an unwound state, permitting intercalation centrally and hydrogen bonding to guanines at the first and sixth base pair positions. The data suggest a different model of binding and recognition by Delta-Rh. The titration of Delta-Rh into a DNA decamer containing the 6 bp recognition site (D1, 5'-CGCATCTGAC-3'; D2, 5'-GTCAGATGCG-3') resulted in two, distinct conformers, in slow exchange on the NMR time scale. The rate of exchange between the two conformers (kobs) at 298 K is 37 s-1, most likely due to partial dissociation between binding modes. The slower rate relative to Lambda-Rh association reflects the relative rigidity of the D1 and/or D2 sequence in comparison to L1. NOE cross-peaks between the intercalating phi ligand and protons of T5-C6, as well as the upfield shifts observed for imino protons at this step, serve to define the central T5-C6 step as the single site of intercalation. The downfield shift of the 3'-G imino protons indicates the complex makes hydrogen bond contacts with these bases. The complex, which is too small to span a 6 bp B-form DNA sequence, nonetheless makes major groove contacts with 3'-G bases to either side of the site. Notably, both 3'-guanine bases are necessary to impart site specificity and slow dissociation kinetics with the 5'-CATCTG-3' site, as evidenced by the extremely exchange-broadened two-dimensional NOESY spectra of Delta-Rh bound to modified duplexes containing N7-deazaguanine at either G8 or G18; the loss of one major groove contact completely abolishes specificity for 5'-CATCTG-3'. DNA chemical shifts upon binding and intermolecular NOE contacts therefore support a model in which Delta-Rh intercalates in one of two canted binding conformations. Within this model, each intercalation mode allows one guanidinium-guanine hydrogen bond at a time, while bringing the other arm close to the phosphate backbone.     

Structural basis of DNA folding and recognition in an AMP-DNA aptamer complex: distinct architectures but common recognition motifs for DNA and RNA aptamers complexed to AMP

BACKGROUND: Structural studies by nuclear magnetic resonance (NMR) of RNA and DNA aptamer complexes identified through in vitro selection and amplification have provided a wealth of information on RNA and DNA tertiary structure and molecular recognition in solution. The RNA and DNA aptamers that target ATP (and AMP) with micromolar affinity exhibit distinct binding site sequences and secondary structures. We report below on the tertiary structure of the AMP-DNA aptamer complex in solution and compare it with the previously reported tertiary structure of the AMP-RNA aptamer complex in solution. RESULTS: The solution structure of the AMP-DNA aptamer complex shows, surprisingly, that two AMP molecules are intercalated at adjacent sites within a rectangular widened minor groove. Complex formation involves adaptive binding where the asymmetric internal bubble of the free DNA aptamer zippers up through formation of a continuous six-base mismatch segment which includes a pair of adjacent three-base platforms. The AMP molecules pair through their Watson-Crick edges with the minor groove edges of guanine residues. These recognition G.A mismatches are flanked by sheared G.A and reversed Hoogsteen G.G mismatch pairs. CONCLUSIONS: The AMP-DNA aptamer and AMP-RNA aptamer complexes have distinct tertiary structures and binding stoichiometries. Nevertheless, both complexes have similar structural features and recognition alignments in their binding pockets. Specifically, AMP targets both DNA and RNA aptamers by intercalating between purine bases and through identical G.A mismatch formation. The recognition G.A mismatch stacks with a reversed Hoogsteen G.G mismatch in one direction and with an adenine base in the other direction in both complexes. It is striking that DNA and RNA aptamers selected independently from libraries of 10(14) molecules in each case utilize identical mismatch alignments for molecular recognition with micromolar affinity within binding-site pockets containing common structural elements.     

Solution conformation of the N-(deoxyguanosin-8-yl)-1-aminopyrene ([AP]dG) adduct opposite dC in a DNA duplex

Combined NMR-molecular mechanics computational studies were undertaken on the C8-deoxyguanosine adduct formed by the carcinogen 1-nitropyrene embedded in the d(C5-[AP]G6-C7).d(G16-C17-G18) sequence context in a 11-mer duplex, with dC opposite the modified deoxyguanosine. The exchangeable and nonexchangeable protons of the aminopyrene moiety and the nucleic acid were assigned following analysis of two-dimensional NMR data sets in H2O and D2O solution. There was a general broadening of several proton resonances for the three nucleotide d(G16-C17-G18) segment positioned opposite the [AP]dG6 lesion site resulting in weaker NOEs involving these protons in the adduct duplex. The solution conformation of the [AP]dG.dC 11-mer duplex has been determined by incorporating intramolecular and intermolecular proton-proton distances defined by upper and lower bounds deduced from NOESY spectra as restraints in molecular mechanics computations in torsion angle space. The aminopyrene ring of [AP]dG6 is intercalated into the DNA helix between intact Watson-Crick dC5.dG18 and dC7.dG16 base pairs. The modified deoxyguanosine ring of [AP]dG6 is displaced into the major groove and stacks with the major groove edge of dC5 in the adduct duplex. Both carbon and proton chemical shift data for the sugar resonances of the modified deoxyguanosine residue are consistent with a syn glycosidic torsion angle for the [AP]dG6 residue. The dC17 base on the partner strand is displaced from the center of the helix toward the major groove as a consequence of the aminopyrene ring intercalation into the helix. This base-displaced intercalative structure of the [AP]dG.dC 11-mer duplex exhibits several unusually shifted proton resonances which can be accounted for by the ring current contributions of the deoxyguanosinyl and pyrenyl rings of the [AP]dG6 adduct. In summary, intercalation of the aminopyrene moiety is accompanied by displacement of both [AP]dG6 and the partner dC17 into the major groove in the [AP]dG.dC 11-mer duplex.     

Solution conformation of the (-)-trans-anti-[BP]dG adduct opposite a deletion site in a DNA duplex: intercalation of the covalently attached benzo[a]pyrene into the helix with base displacement of the modified deoxyguanosine into the minor groove

A combined NMR-computational approach was employed to determine the solution structure of the (-)-trans-anti-[BP]dG adduct positioned opposite a -1 deletion site in the d(C1-C2-A3-T4-C5- [BP]G6-C7-T8-A9-C10-C11).d(G12-G13-T14-A15-G1 6-G17-A18-T19-G20-G21) sequence context. The (-)-trans-anti-[BP]dG moiety is derived from the binding of the (-)-anti-benzo[a]pyrene diol epoxide [(-)-anti-BPDE] to N2 of dG6 and has a 10R absolute configuration at the [BP]dG linkage site. The exchangeable and non-exchangeable protons of the benzo[a]pyrenyl moiety and the nucleic acid were assigned following analysis of two-dimensional NMR data sets in H2O and D2O solution. The solution conformation has been determined by incorporating intramolecular and intermolecular proton-proton distances defined by lower and upper bounds deduced from NOESY spectra as restraints in molecular mechanics computations in torsion angle space followed by restrained molecular dynamics calculations based on a NOE distance and intensity refinement protocol. Our structural studies establish that the aromatic BP ring system intercalates into the helix opposite the deletion site, while the modified deoxyguanosine residue is displaced into the minor groove with its face parallel to the helix axis. The intercalation site is wedge-shaped and the BP aromatic ring system stacks over intact flanking Watson-Crick dG.dC base pairs. The modified deoxyguanosine stacks over the minor groove face of the sugar ring of the 5'-flanking dC5 residue. The BP moiety is positioned with the benzylic ring oriented toward the minor groove and the distal pyrenyl aromatic ring directed toward the major groove. This conformation strikingly contrasts with the corresponding structure in the full duplex with the same 10R (-)-trans-anti-[BP]dG lesion positioned opposite a complementary dC residue [de los Santos et al. (1992) Biochemistry 31, 5245-5252); in this case the aromatic BP ring system is located in the minor groove, and there is no disruption of the [BP]dG.dC Watson-Crick base pairing alignment. The intercalation-base displacement features of the 10R (-)-trans-anti-[BP]dG adduct opposite a deletion site have features in common to those of the 10S (+)-trans-anti-[BP]dG adduct opposite a deletion site previously reported by Cosman et al. [(1994)(Biochemistry 33, 11507-11517], except that there is a nearly 180 degrees rotation of the BP residue about the axis of the helix at the base-displaced intercalation site and the modified deoxyguanosine is positioned in the opposite groove. In the 10S adduct, the benzylic ring is in the major groove and the aromatic ring systems point toward the minor groove. This work extends the theme of opposite orientations of adducts derived from chiral pairs of (+)- and (-)-anti-BPDE enantiomers; both 10S and 10R adducts can be positioned with opposite orientations either in the minor groove or at base displaced intercalation sites, depending on the presence or absence of the partner dC base in the complementary strand.     

Molecular modelling and footprinting studies of DNA minor groove binders: bisquaternary ammonium heterocyclic compounds

We report new quantitative footprinting data which reveal differences in binding constants of bisquaternary ammonium heterocyclic compounds (BQA) with AT-rich DNA sites depending on the ligand structure and on the size and sequence of the DNA binding site. In an attempt to understand the dependence of binding affinity on the ligand structure we have performed quantum-chemical AM1 calculations on the BQA compounds and on subunits to explore the conformational space and to calculate the electronic and structural features of individual ligand conformations. Due to the properties of the rotatable backbone bonds, there is a large number of possible conformations with almost equal energy. We present a new method for the calculation of the radius of curvature of molecular structures. Assuming that strong binders should have a shape complementary to the DNA minor groove, this measure is used to select the optimum conformations for DNA-drug binding. The approach yields the correct ligand conformation for SN6999, for which an X-ray DNA-drug structure is known. The curvature of the optimum conformations of all ligands is compared with the experimental binding constants. A correlation is found between curvature and binding constant provided other structural factors do not vary. Therefore, we conclude that within structurally similar BQA compounds the extent of curvature is the relevant quantity which modulates the binding affinity.     

The high resolution crystal structure of the DNA decamer d(AGGCATGCCT)

The crystal structure of the DNA decamer d(AGGCATGCCT) has been determined to a resolution of 1.3 A and R factor of 13.9%. The structure has a unique conformation with each of the decamer single strands forming base-pairing interactions with two symmetry-related strands. The central eight bases of the decamer form an A-DNA octamer duplex with one symmetry-related strand whilst the terminal 5'-A and T-3' bases are flipped out and away from the octamer helix axis to form base-pairing interactions with a second symmetry-related strand. These A.T base-pairs lie perpendicular to the crystallographic c axis and pack within the unit cell in conjunction with a symmetry-related A.T base-pair displaced by 3.4 A degrees along the c axis. A novel base triplet interaction of the type A*(G.C) is present in the structure with interaction from the major groove side of the terminal 5'-A base to the minor groove of the central A-DNA octamer. This structure reports the first example of cobalt hexammine binding to a right-handed DNA duplex. The crystallographic asymmetric unit contains two cobalt hexammine ligands with one site in the major groove coordinating via hydrogen bonds to the 5'-AGG bases, and the second site located between DNA molecules and interacting with the oxygen atoms of phosphate groups.     

Non-covalent complexes between DNA-binding drugs and double-stranded deoxyoligonucleotides: a study by ionspray mass spectrometry

The non-covalent complexes between some DNA-binding drugs and duplex oligodeoxynucleotides were studied by ionspray mass spectrometry, with the aim of evaluating the suitability of this technique to screen rapidly a series of drugs exerting their activity through non-covalent binding to specific base sequences of DNA. Two classes of drugs were considered, distamycins (which show affinity for the minor groove of DNA) and anthracyclines (which interact through intercalation between bases). For the former, d(CGCGAATTCGCG)2 was chosen as the model oligodeoxynucleotide. Following optimization of sample preparation and instrumental conditions, the complexes of different distamycins were observed; depending on the ligand considered, 1:1 or 2:1 complexes were formed preferentially. A semi-quantitative evaluation of the relative affinities was made by measuring the ratio of the complexes signals to those of the duplex, and also by competitive binding with equimolar amounts of distamycin. For anthracyclines, the daunorubicin-d(CGATCG)2 complex was chosen as the model for a preliminary mass spectrometric study; however, the signals of the duplex and the complex were very low compared with the monomer signal. Since the complex was known to be stable in solution, this was ascribed to gas-phase instability, probably caused by electrostatic repulsion between negatively charged phosphate groups.     

Interaction of minor groove ligands to an AAATT/AATTT site: correlation of thermodynamic characterization and solution structure

A combination of circular dichroism spectroscopy, titration calorimetry, and optical melting has been used to investigate the association of the minor groove ligands netropsin and distamycin to the central A3T2 binding site of the DNA duplex d(CGCAAATTGGC).d(GCCAATTTGCG). For the complex with netropsin at 20 degrees C, a ligand/duplex stoichiometry of 1:1 was obtained with Kb approximately 4.3 x 10(7) M-1, delta Hb approximately -7.5 kcal mol-1, delta Sb approximately 9.3 cal K-1 mol-1, and delta Cp approximately 0. Previous NMR studies characterized the distamycin complex with A3T2 at saturation as a dimeric side-by-side complex. Consistent with this result, we found a ligand/duplex stoichiometry of 2:1. In the current study, the relative thermodynamic contributions of the two distamycin ligands in the formation of this side-by-side complex (2:1 Dst.A3T2) were evaluated and compared with the thermodynamic characteristics of netropsin binding. The association of the first distamycin molecule of the 2:1 Dst.A3T2 complex yielded the following thermodynamic profile: Kb approximately 3.1 x 10(7) M-1, delta Hb = -12.3 kcal mol-1, delta Sb = -8 cal K-1 mol-1, and delta Cp = -42 cal K-1 mol-1. The binding of the second distamycin molecule occurs with a lower Kb of approximately 3.3 x 10(6) M-1, a more favorable delta Hb of -18.8 kcal mol-1, a more unfavorable delta Sb of -34 cal K-1 mol-1, and a higher delta Cp of -196 cal K-1 mol-1. The latter term indicates an ordering of electrostricted and structural water molecules by the complexes. These results correlate well with the NMR titrations and are discussed in context of the solution structure of the 2:1 Dst.A3T2 complex.     

Structural basis of DNA recognition and mechanism of quadruplex formation by the beta subunit of the Oxytricha telomere binding protein

Interactions of the beta subunit of the Oxytricha nova telomere binding protein with the telomeric DNA sequences, d(T4G4)2 and dT6(T4G4)2, have been investigated in vitro using Raman and fluorescence spectroscopies. Raman difference spectra show that the beta subunit binds to both d(T4G4)2 and dT6(T4G4)2 but promotes the formation of a parallel-stranded quadruplex only in dT6(T4G4)2, thus demonstrating the importance of the telomeric 5' tail for in vitro recognition and guanine quadruplex formation. While d(T4G4)2 is not a suitable substrate for quadruplex promotion by the beta subunit, the Raman spectra reveal other structural rearrangements of this DNA strand upon beta subunit binding, including changes in guanine glycosyl torsion angles from syn to anti and disruption of carbonyl hydrogen-bonding interactions. The conformation of d(T4G4)2 in the beta:d(T4G4)2 complex is suggested as a plausible intermediate along the pathway to formation of the parallel-stranded guanine quadruplex. Fluorescence band shifts indicate that at least one of the two tryptophans of the beta subunit is shielded from solvent as a consequence of DNA binding in both the beta:dT6(T4G4)2 and beta:d(T4G4)2 complexes. However, the Raman spectra of these complexes suggest no significant changes in the beta subunit secondary structure attendant with DNA binding. A model for beta subunit binding by Oxytricha telomeric DNA sequences and a mechanism for quadruplex formation are proposed. A key feature of this model is the use of a telomeric hairpin secondary structure as the recognition motif.     

Dicationic diarylfurans as anti-Pneumocystis carinii agents

Seven dicationic 2,5-diarylfurans have been synthesized, and their interactions with poly(dA-dT) and the duplex oligomer d(CGCCAATTCGCG)2 were evaluated by Tm measurements. The inhibition of topoisomerase II isolated from Giardia lamblia, the inhibition of growth of G. lamblia in cell culture by these furans, and the effectiveness of these compounds against Pneumocystis carinii in the immunosuppressed rat model have been assessed. Strong binding affinities to poly(dA-dT) and to the oligomer were observed for the dicationic furans, and the interaction strength is directly correlated to the biological activity of the compounds. An X-ray structure for the complex of the dicationic amidine derivative, 2,5-bis(4-guanylphenyl)furan (1), with the oligomer demonstrates the snug fit of these compounds with the AATT minor-groove binding site and hydrogen bonds to AT base pairs at the floor of the minor groove. The stronger DNA binding molecules are the most effective inhibitors of topoisomerase II and G. lamblia in cell culture, and there is a correlation for both DNA interaction and topoisomerase II inhibition with the biological activity of these compounds against G. lamblia. Compound 1 is the most effective against P. carinii, it is more active and less toxic than pentamidine on intravenous administration and it is also effective by oral dosage. The results presented here suggest a model for the biological action of these compounds in which the dication first binds in the minor groove of DNA and forms a complex that results in the inhibition of the microbial topoisomerase II enzyme.     

Designer DNA-binding drugs: the crystal structure of a meta-hydroxy analogue of Hoechst 33258 bound to d(CGCGAATTCGCG)2

An analogue of the DNA binding compound Hoechst 33258, which has the para hydroxyl group altered to be at the meta position, together with the replacement of one benzimidazole group by pyridylimidazole, has been cocrystallized with the dodecanucleotide sequence d(CGCGAATTCGCG)2. The X-ray structure has been determined at 2.2 A resolution and refined to an R factor of 20.1%. The ligand binds in the minor groove at the sequence 5'-AATTC with the bulky piperazine group extending over the CxG base pair. This binding is stabilised by hydrogen bonding and numerous close van der Waals contacts to the surface of the groove walls. The meta-hydroxyl group was found in two distinct orientations, neither of which participates in direct hydrogen bonds to the exocyclic amino group of a guanine base. The conformation of the drug differs from that found previously in other X-ray structures of Hoechst 33258-DNA complexes. There is significant variation between the minor groove widths in the complexes of Hoechst 33258 and the meta-hydroxyl derivative as a result of these conformational differences. Reasons are discussed for the inability of this derivative to actively recognise guanine.     

Crystal structure of a B-DNA dodecamer containing inosine, d(CGCIAATTCGCG), at 2.4 A resolution and its comparison with other B-DNA dodecamers

The crystal structure of the dodecamer, d(CGCIAATTCGCG), has been determined at 2.4 A resolution by molecular replacement, and refined to an R-factor of 0.174. The structure is isomorphous with that of the B-DNA dodecamer, d(CGCGAATTCGCG), in space group P2(1)2(1)2(1) with cell dimensions of a = 24.9, b = 40.4, and c = 66.4 A. The initial difference Fourier maps clearly indicated the presence of inosine instead of guanine. The structure was refined with 44 water molecules, and compared to the parent dodecamer. Overall the two structures are very similar, and the I:C forms Watson-Crick base pairs with similar hydrogen bond geometry to the G:C base pairs. The propeller twist angle is low for I4:C21 and relatively high for the I16:C9 base pair (-3.2 degrees compared to -23.0 degrees), and the buckle angles alter, probably due to differences in the contacts with symmetry related molecules in the crystal lattice. The central base pairs of d(CGCIAATTCGCG) show the large propeller twist angles, and the narrow minor groove that characterize A-tract DNA, although I:C base pairs cannot form the major groove bifurcated hydrogen bonds that are possible for A:T base pairs.     

Structural characterization of an N-acetyl-2-aminofluorene (AAF) modified DNA oligomer by NMR, energy minimization, and molecular dynamics

An N-acetyl-2-aminofluorene (AAF) modified deoxyoligonucleotide duplex, d(C1-C2-A3-C4-[AAF-G5]-C6-A7-C8-C9).d(G10-G11-T12-G13-C14-++ +G15-T16-G17-G18), was studied by one- and two-dimensional NMR spectroscopy. Eight of the nine complementary nucleotides form Watson-Crick base pairs, as shown by NOEs between the guanine imino proton and cytosine amino protons for G.C base pairs or by an NOE between the thymine imino proton and adenine H2 proton for A.T base pairs. The AAF-G5 and C14 bases show no evidence of complementary hydrogen bond formation to each other. The AAF-G5 base adopts a syn conformation, as indicated by NOEs between the G5 imino proton and the A3-H3' and A3-H2'/H2" protons and by NOEs between the fluorene-H1 proton of AAF and the G5-H1' or C6-H1' proton. The NOEs from the C4-H6 proton to C4 sugar protons are weak, and thus the glycosidic torsion angle in this nucleotide is not well defined by these NMR data. The remaining bases are in the anti conformation, as depicted by the relative magnitude of the H8/H6 to H2' NOEs when compared to the H8/H6 to H1' NOEs. The three base pairs on each end of the duplex exhibit NOEs characteristic of right-handed B-form DNA. Distance restraints obtained from NOESY data recorded at 32 degrees C using a 100-ms mixing time were used in conformational searches by molecular mechanics energy minimization studies. The final, unrestrained, minimum-energy conformation was then used as input for an unrestrained molecular dynamics simulation. Chemical exchange cross peaks are observed, and thus the AAF-9-mer exists in more than a single conformation on the NMR time scale. The NMR data, however, indicate the presence of a predominant conformation (> or = 70%). The structure of the predominant conformation of the AAF-9-mer shows stacking of the fluorene moiety on an adjacent base pair, exhibiting features of the base-displacement [Grunberger, D., Nelson, J. H., et al. (1970) Proc. Natl. Acad. Sci. U.S.A. 66, 488-494] and insertion-denaturation models [Fuchs, R.P.P., & Daune, M. (1971) FEBS Lett. 14, 206-208], while the distal ring of the fluorene moiety protrudes into the minor groove.     

On the importance of van der Waals interaction in the groove binding of DNA with ligands: restrained molecular dynamics study

Comparison of interaction energy between an oligonucleotide and a DNA-binding ligand in the minor and major groove modes was made by use of restrained molecular dynamics. Distortion in DNA was found for the major groove mode whereas less significant changes for both ligand and DNA were detected for the minor groove binding after molecular dynamics simulation. The conformation of the ligand obtained from the major groove modes resembles that computed with the ligand soaked in water. The van der Waals contact energy was found to be as significant as electrostatic energy and more important for difference in binding energy between these two binding modes. The importance of van der Waals force in groove binding was supported by computations on the complex formed by the repressor peptide fragment from the bacteriophage 434 and its operator oligonucleotide.     

A molecular mechanics and dynamics study of the minor adduct between DNA and the carcinogen 2-(acetylamino)fluorene (dG-N2-AAF)

Experimental studies involving the carcinogenic aromatic amine 2-(acetylamino)fluorene (AAF) have afforded two acetylated DNA adducts, the major one bound to C8 of guanine and a minor adduct bound to N2 of guanine. The minor adduct may be important in carcinogenesis because it persists, while the major adduct is rapidly repaired. Primer extension studies of the minor adduct have indicated that it blocks DNA synthesis, with some bypass and misincorporation of adenine opposite the lesion [Shibutani, S., and Grollman, A.P. (1993) Chem. Res. Toxicol. 6, 819-824]. No experimental structural information is available for this adduct. Extensive minimized potential energy searches involving thousands of trials and molecular dynamics simulations were used to study the conformation of this adduct in three sequences: I, d(C1-G2-C3-[AAF]G4-C5-G6-C7).d(G8-C9-G10-C11-G12-C13-G14+ ++); II, the sequence of Shibutani and Grollman, d(C1-T2-A3-[AAF]G4-T5-C6-A7).d(T8-G9-A10-C11-T12-A13-G14); and III, which is the same as II but with a mismatched adenine in position 11, opposite the lesion. AAF was located in the minor groove in the low-energy structures of all sequences. In the lowest energy form of the C3-[AAF]G4-C5 sequence I, the fluorenyl rings point in the 3' direction along the modified strand and the acetyl in the 5' direction. These orientations are reversed in the second lowest energy structure of this sequence, and the energy of this structure is 1.4 kcal/mol higher. Watson Crick hydrogen bonding is intact in both structures. In the two lowest energy structures of the A3-[AAF]G4-T5 sequence II, the AAF is also located in the minor groove with Watson-Crick hydrogen bonding intact. However, in the lowest energy form, the fluorenyl rings point in the 5' direction and the acetyl in the 3' direction. The energy of the structure with opposite orientation is 5.1 kcal/mol higher. In sequence III with adenine mismatched to the modified guanine, the lowest energy form also had the fluorenyl rings oriented 5' in the minor groove with intact Watson-Crick base pairing. However, the mispaired adenine adopts a syn orientation with Hoogsteen pairing to the modified guanine. These results suggest that the orientation of the AAF in the minor groove may be DNA sequence dependent. Mobile aspects of favored structures derived from molecular dynamics simulations with explicit solvent and salt support the essentially undistorting nature of this lesion, which is in harmony with its persistence in mammalian systems.     

Minor groove-directed and intercalative ligand-DNA interactions in the poisoning of human DNA topoisomerase I by protoberberine analogs

Spectroscopic, calorimetric, DNA cleavage, electrophoretic, and computer modeling techniques have been employed to characterize the DNA binding and topoisomerase poisoning properties of three protoberberine analogs, 8-desmethylcoralyne (DMC), 5,6-dihydro-8-desmethylcoralyne (DHDMC), and palmatine, which differ in the chemical structures of their B- and/or D-rings. DNA topoisomerase-mediated cleavage assays revealed that these compounds were unable to poison mammalian type II topoisomerase. By contrast, the three protoberberine analogs poisoned human topoisomerase I according to the following hierarchy: DHDMC > DMC > palmatine. DNA binding by all three protoberberine analogs induced negative flow linear dichroism signals as well as unwinding of the host duplex. These two observations are consistent with an intercalative mode of protoberberine binding to duplex DNA. However, a comparison of the DNA binding properties for DMC and DHDMC, which differ only by the state of saturation at the 5,6 positions of the B-ring, revealed that the protoberberine analogs do not "behave" like classic DNA intercalators. Specifically, saturation of the 5-6 double bond in the B-ring of DMC, thereby converting it to the DHDMC molecule, was associated with enhanced DNA unwinding as well as a reversal of DNA binding preference from a DNA duplex with an inaccessible or occluded minor groove {poly[d(G-C)]2} to DNA duplexes with accessible or unobstructed minor grooves {poly[d(A-T)]2 and poly[d(I-C)]2}. In addition, a comparison of the DNA binding properties for DHDMC and palmatine revealed that transferring the 11-methoxy moiety on the D-ring of DHDMC to the 9 position, thereby converting it to palmatine, was associated with a reduction in binding affinity for both duplexes with unobstructed minor grooves as well as for duplexes with occluded minor grooves. These DNA binding properties are consistent with a "mixed-mode" DNA binding model for protoberberines in which a portion of the ligand molecule intercalates into the double helix, while the nonintercalated portion of the ligand molecule protrudes into the minor groove of the host duplex, where it is thereby available for interactions with atoms lining the floor and/or walls of the minor groove. Furthermore, saturation at the 5,6 positions of the B-ring, which causes the A-ring to be tilted relative to the plane formed by rings C and D, appears to stabilize the interaction between the host duplex and the minor groove-directed portion of the protoberberine ligand. Computer modeling studies on the DHDMC-poly[d(A-T)]2 complex suggest that this interaction may involve van der Waals contacts between the ligand A-ring and backbone sugar atoms lining the minor groove of the host duplex. The hierarchy of topoisomerase I poisoning noted above suggests that this minor groove-directed interaction may play an important role in topoisomerase I poisoning by protoberberine analogs. In the aggregate, our results presented here, coupled with the recent demonstration of topoisomerase I poisoning by minor groove-binding terbenzimidazoles [Sun, Q., Gatto, B., Yu, C., Liu, A. , Liu, L. F., & LaVoie, E. J. (1995) J. Med. Chem. 38, 3638-3644], suggest that minor groove-directed ligand-DNA interactions may be of general importance in the poisoning of topoisomerase I.