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.     

GC base sequence recognition by oligo(imidazolecarboxamide) and C-terminus-modified analogues of distamycin deduced from circular dichroism, proton nuclear magnetic resonance, and methidiumpropylethylenediaminetetraacetate-iron(II) footprinting studies

The DNA binding properties of a series of imidazole-containing and C-terminus-modified analogues 4-7 of distamycin are described. These analogues contain one to four imidazole units, respectively. Data from the ethidium displacement assay showed that these compounds bind in the minor groove of DNA, with the relative order of binding constants of 6 (Im3) > 7 (Im4) > 5 (Im2) > 4 (Im1). The reduced binding constants of these compounds for poly(dA-dT) relative to distamycin, while they still interact strongly with poly(dG-dC), provided evidence of GC sequence acceptance. The preferences for GC-rich sequences by these compounds were established from a combination of circular dichroism (CD) titration, proton nuclear magnetic resonance (1H-NMR), and methidiumpropylethylenediaminetetraacetate-iron(II) [MPE.Fe-(II)] footprinting studies. In the CD studies, these compounds produced significantly larger DNA-induced ligand bands with poly(dG-dC) than poly(dA-dT) at comparable ligand concentrations. 1H-NMR studies of the binding of 5 to d-[CATGGCCATG]2 provided further evidence of the recognition of GC sequences by these compounds, and suggested that the ligand was located on the underlined sequence in the minor groove with the C-terminus oriented over the T residue. MPE footprinting studies on a GC-rich BamHI/SalI fragment of pBR322 provided unambiguous evidence for the GC sequence selectivity for some of these compounds. Compounds 4 and 7 produced poor footprints on the gels; however, analogues 5 and 6 gave strong footprints.(ABSTRACT TRUNCATED AT 250 WORDS)     

Mono- and dysfunctional nitrogen mustard analogues of the DNA minor groove binder pibenzimol. Synthesis, cytotoxicity and interaction with DNA

Two series of mono- and dysfunctional aniline mustards linked to a bisbenzimidazole minor groove binder have been prepared using a new method (polyphosphate ester-mediated direct coupling of appropriate mustard acids with a preformed advanced phenylenediamine intermediate). As the linker chain attaching the mustard was lengthened the binding site size of the compounds to calf thymus DNA remained essentially constant at 2.6 nucleotides, but reversible binding strength declined by a factor of 2. Analogues with longer linker chains alkylated DNA much more rapidly than those with shorter chains, consistent with the electronic factors. The short chain analogues also failed to alkylate a 120 bp HindIII to Bg/II fragment of the gpt gene, as measured by gel electrophoresis cleavage assays. The longer chain analogues (both mono- and dysfunctional mustards) showed patterns of DNA alkylation that varied with chain length. In particular, while most compounds showed substantial N7 alkylation at many guanine residues, the analogue with a (CH2)3 linker chain showed strong alkylation at adenine sites in poly-AT regions. For the longer chain analogues, the bifunctional mustards were substantially (10- to 20-fold) more cytotoxic than the corresponding monofunctional analogues.     

(+)-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.     

A peptide interaction in the major groove of RNA resembles protein interactions in the minor groove of DNA

A 17-amino acid arginine-rich peptide from the bovine immunodeficiency virus Tat protein has been shown to bind with high affinity and specificity to bovine immunodeficiency virus transactivation response element (TAR) RNA, making contacts in the RNA major groove near a bulge. We show that, as in other peptide-RNA complexes, arginine and threonine side chains make important contributions to binding but, unexpectedly, that one isoleucine and three glycine residues also are critical. The isoleucine side chain may intercalate into a hydrophobic pocket in the RNA. Glycine residues may allow the peptide to bind deeply within the RNA major groove and may help determine the conformation of the peptide. Similar features have been observed in protein-DNA and drug-DNA complexes in the DNA minor groove, including hydrophobic interactions and binding deep within the groove, suggesting that the major groove of RNA and minor groove of DNA may share some common recognition features.     

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.     

DNA-Directed alkylating agents. 7. Synthesis, DNA interaction, and antitumor activity of bis(hydroxymethyl)- and bis(carbamate)-substituted pyrrolizines and imidazoles

A series of bis(hydroxymethyl)-substituted imidazoles, thioimidazoles, and pyrrolizines and related bis(carbamates), linked to either 9-anilinoacridine (intercalating) or 4-(4-quinolinylamino)benzamide (minor groove binding) carriers, were synthesized and evaluated for sequence-specific DNA alkylation and cytotoxicity. The imidazole and thioimidazole analogues were prepared by initial synthesis of [(4-aminophenyl)alkyl]imidazole-, thioimidazole-, or pyrrolizine dicarboxylates, coupling of these with the desired carrier, and reduction to give the required bis(hydroxymethyl) alkylating moiety. The pyrrolizines were the most reactive alkylators, followed by the thioimidazoles, while the imidazoles were unreactive. The pyrrolizines and some of the thioimidazoles cross-linked DNA, as measured by agarose gel electrophoresis. Strand cleavage assays showed that none of the compounds reacted at purine N7 or N3 sites in the gpt region of the plasmid gpt2Eco, but the polymerase stop assay showed patterns of G-alkylation in C-rich regions. The corresponding thioimidazole bis(carbamates) were more selective than the bis(hydroxymethyl) pyrrolizines, with high-intensity bands at 5'-NCCN, 5'-NGCN and 5'-NCGN sequences in the PCR stopping assay ( indicates block sites). The data suggest that these targeted compounds, like the known thioimidazole bis(carbamate) carmethizole, alkylate exclusively at guanine residues via the 2-amino group, with little or no alkylation at N3 and N7 guanine or adenine sites. The cytotoxicities of the compounds correlated broadly with their reactivities, with the bis(hydroxymethyl)imidazoles being the least cytotoxic (IC50s >1 microM; P388 leukemia) and with the intercalator-linked analogues being more cytotoxic than the corresponding minor-groove-targeted ones. This was true also for the more reactive thioimidazole bis(carbamates) (IC50s 0.8 and 11 microM, respectively), but both were more active than the analogous "untargeted" carmethizole (IC50 20 microM). The bis(hydroxymethyl)pyrrolizine analogues were the most cytotoxic, with IC50s as low as 0.03 microM.     

Targeting of human retrotransposon integration is directed by the specificity of the L1 endonuclease for regions of unusual DNA structure

L1 elements are polyA retrotransposons which inhabit the human genome. Recent work has defined an endonuclease (L1 EN) encoded by the L1 element required for retrotransposition. We report the sequence specificity of this nicking endonuclease and the physical basis of its DNA recognition. L1 endonuclease is specific for the unusual DNA structural features found at the TpA junction of 5'(dTn-dAn) x 5'(dTn-dAn) tracts. Within the context of this sequence, substitutions which generate a pyrimidine-purine junction are tolerated, whereas purine-pyrimidine junctions greatly reduce or eliminate nicking activity. The A-tract conformation of the DNA substrate 5' of the nicked site is required for L1 EN nicking. Chemical or physical unwinding of the DNA helix enhances L1 endonuclease activity, while disruption of the adenine mobility associated with TpA junctions reduces it. Akin to the protein-DNA interactions of DNase I, L1 endonuclease DNA recognition is likely mediated by minor groove interactions. Unlike several of its homologues, however, L1 EN exhibits no AP endonuclease activity. Finally, we speculate on the implications of the specificity of the L1 endonuclease for the parasitic relationship between retroelements and the human genome.     

Asymmetric interactions of hexameric bacteriophage T7 DNA helicase with the 5'- and 3'-tails of the forked DNA substrate

Bacteriophage T7 DNA helicase requires two noncomplementary single-stranded DNA (ssDNA) tails next to a double-stranded DNA (dsDNA) region to initiate DNA unwinding. The interactions of the helicase with the DNA were investigated using a series of forked DNAs. Our results show that the helicase interacts asymmetrically with the two tails of the forked DNA. When the helicase was preassembled on the forked DNA before the start of unwinding, a DNA with 15-nucleotide (nt) 3'-tail and 35-nt 5'-tail was unwound with optimal rates close to 60 base pairs/s at 18 degrees C. When the helicase was not preassembled on the DNA, a >65-nt long 5'-tail was required for maximal unwinding rates of 12 base pairs/s. We show that the helicase interacts specifically with the ssDNA region and maintains contact with both ssDNA strands during DNA unwinding, since conversion of the two ssDNA tails to dsDNA structures greatly inhibited unwinding, and the helicase was unable to unwind past a nick in the dsDNA region. These studies have provided new insights into the mechanism of DNA unwinding. We propose an exclusion model of DNA unwinding in which T7 helicase hexamer interacts mainly with the ssDNA strands during DNA unwinding, encircling the 5'-strand and excluding the 3'-strand from the hole.     

Inhibition of restriction endonuclease activity by DNA binding fluorochromes

Activity of type II restriction endonuclease is affected by many common factors including buffer composition and sequences flanking the recognition site (Brabec et al., Eur.J. Biochem. 216, 183, 1993). The successful development of Optical Mapping (Schwartz et al., Science, 262, 110, 1993; Meng et al., Nature Genet. 9, 432, 1995; Wang and Schwartz, PNAS, 1995 Cai et al., PNAS, 92, 5164, 1995) relied on optimization of light microscope-based imaging of fluorescently labeled DNA molecules during restriction endonuclease digestion. Little was known about the effects of commonly used DNA-fluorochromes on restriction endonuclease activity. Thus, we developed an enzyme activity assay using lambda bacteriophage DNA or adenovirus-2 DNA to evaluate the effects of five DNA binding fluorochromes (4'-6-daimidine-2-phenylindole (DAPI), ethidium bromide (EtdBr), ethidium bromide homodimer (EthD-1), bis-benzimide (H33258) and benzothiazolium-4-quinolinium dimer (TOTO-1)) on the enzymatic activities of eleven type II restriction endonucleases (Asc I, Csp I, Dra I, EcoR I, Hha I, Hind III, Not I, Rsr II, Sfi I, SgrA I and Sma I). We found that the minor groove binding fluorochrome, DAPI, did not measurably inhibit activity of this group, with the exception of Dra I. Similarly, another minor groove binding fluorochrome H33258 inhibited Dra I and Not I (slightly). The three intercalating fluorochromes EtdBr, EthD-1 and TOTO-1, however, variably inhibited the other enzymes. Since Beta-mercaptoethanol (Beta-ME) is used to discourage photodamage of stained DNA molecules, we also assessed its effect on restriction endonuclease activity. Interestingly, Dra I, Hind III, Sfi I and Sma I retained full activities at high concentration of Beta-ME (5%), but Asc I, Csp I, Not I, Rsr II and SgrA I showed varying sensitivities to the Beta-ME. Isoschizomers Csp I and Rsr II behaved differently to both fluorochromes and Beta-ME. The results presented here should provide a basis for further development of new Optical Mapping-based techniques requiring fluorescence labeling of other actively imaged enzymatic reactions.     

GC rich DNA oligonucleotides with narrow minor groove width

Investigation of the width of the minor groove using 500 MHz NMR spectroscopy in three closely related 11-mer B-DNA duplexes shows that the minor groove is narrow in a GC rich oligonucleotide, and that a narrow minor groove is not something endemic to DNAs with persistent repetitions of adenine nucleotides (A-tract DNA). The width of the groove is dictated by local sequence contexts and independent of neighboring A-tract DNA.     

A two-site mechanism for ATP hydrolysis by the asymmetric Rep dimer P2S as revealed by site-specific inhibition with ADP-A1F4

The Escherichia coli Rep helicase is a dimeric motor protein that catalyzes the transient unwinding of duplex DNA to form single-stranded (ss) DNA using energy derived from the binding and hydrolysis of ATP. In an effort to understand this mechanism of energy transduction, we have used pre-steady-state methods to study the kinetics of ATP binding and hydrolysis by an important intermediate in the DNA unwinding reaction--the asymmetric Rep dimer state, P2S, where ss DNA [dT(pT)15] is bound to only one subunit of the Rep dimer. To differentiate between the two potential ATPase active sites inherent in the dimer, we constructed dimers with one subunit covalently cross-linked to ss DNA and where one or the other of the ATPase sites was selectively complexed to the tightly bound transition state analog ADP-A1F4. We found that when ADP-A1F4 is bound to the Rep subunit in trans from the subunit bound to ss DNA, steady-state ATPase activity of 18 s(-1) per dimer (equivalent to wild-type P2S) was recovered. However, when the ADP-A1F4 and ss DNA are both bound to the same subunit (cis), then a titratable burst of ATP hydrolysis is observed corresponding to a single turnover of ATP. Rapid chemical quenched-flow techniques were used to resolve the following minimal mechanism for ATP hydrolysis by the unligated Rep subunit of the cis dimer: E + ATP <==> E-ATP <==> E'-ATP <==> E'-ADP-Pi <==> E-ADP-Pi <==> E-ADP + Pi <==> E + ADP + Pi, with K1 = (2.0 +/- 0.85) x 10(5) M(-1), k2 = 22 +/- 3.5 s(-1), k(-2) < 0.12 s(-1), K3 = 4.0 +/- 0.4 (k3 > 200 s(-1)), k4 = 1.2 +/- 0.14 s(-1), k(-4) < 1.2 s(-1), K5 = 1.0 +/- 0.2 mM, and K6 = 80 +/- 8 microM. A salient feature of this mechanism is the presence of a kinetically trapped long-lived tight nucleotide binding state, E'-ADP-Pi. In the context of our "subunit switching" model for Rep dimer translocation during processive DNA unwinding [Bjornson, K. B., Wong, I., & Lohman, T. M. (1996) J. Mol. Biol. 263, 411-422], this state may serve an energy storage function, allowing the energy from the binding and hydrolysis of ATP to be harnessed and held in reserve for DNA unwinding.     

Simultaneous determination of helical unwinding angles and intrinsic association constants in ligand-DNA complexes: the interaction between DNA and calichearubicin B

We present a helical unwinding assay for reversibly binding DNA ligands that uses closed circular DNA, topoisomerase I (Topo I), and two-dimensional agarose gel electrophoresis. Serially diluted Topo I relaxation reactions at constant DNA/ligand ratio are performed, and the resulting apparent unwinding of the closed circular DNA is used to calculate both ligand unwinding angle (phi) and intrinsic association constant (Ka). Mathematical treatment of apparent unwinding is formally analogous to that of apparent extinction coefficient data for optical binding titrations. Extrapolation to infinite DNA concentration yields the true unwinding angle of a given ligand and its association constant under Topo I relaxation conditions. Thus this assay delivers simultaneous structural and thermodynamic information describing the ligand-DNA complex. The utility of this assay has been demonstrated by using calichearubicin B (CRB), a synthetic hybrid molecule containing the anthraquinone chromophore of (DA) and the carbohydrate domain of calicheamicin gamma1I. The unwinding angle for CRB calculated by this method is -5. 3 +/- 0.5 degrees. Its Ka value is 0.20 x 10(6) M-1. For comparison, the unwinding angles of ethidium bromide and DA have been independently calculated, and the results are in agreement with canonical values for these compounds. Although a stronger binder to selected sites, CRB is a less potent unwinder than its parent compound DA. The assay requires only small amounts of ligand and offers an attractive option for analysis of DNA binding by synthetic and natural compounds.     

The herpes simplex virus type 1 origin binding protein. Specific recognition of phosphates and methyl groups defines the interacting surface for a monomeric DNA binding domain in the major groove of DNA

The UL9 gene of herpes simplex virus type 1 (HSV-1) encodes an origin binding protein (OBP). It is an ATP-dependent DNA helicase and a sequence-specific DNA-binding protein. The latter function is carried out by the C-terminal domain of OBP (DeltaOBP). We have now performed a quantitative analysis of the interaction between DeltaOBP and its recognition sequence, GTTCGCAC, in oriS. Initially optimal conditions for binding were carefully determined. We observed that complexes with different electrophoretic mobilities were formed. A cross-linking experiment demonstrated that nonspecific complexes containing 2 or more protein monomers per DNA molecule were formed at high protein concentrations. The specific complex formed at low concentrations of DeltaOBP had an electrophoretic mobility corresponding to a 1:1 complex. We then demonstrated that the methyl groups of thymine in the major groove were essential for high affinity binding. Changes in the minor groove had considerably smaller effects. Ethylation interference experiments indicated that specific contacts were made between OBP and three phosphates in the recognition sequence. Finally, these observations were used to present a model of the surface of DNA that interacts with DeltaOBP in a sequence-specific manner.     

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.     

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 new type of DNA minor-groove complex: carbazole dication-DNA interactions

The effect of opportunistic infections (OI) on immune-compromised populations has been known for decades, but the recent AIDS epidemic has sparked renewed interest in the development of new anti-OI agents. The mechanism of action of a series of cationic unfused-aromatic anti-OI drugs is believed to involve binding of the drug to AT sequences in the minor groove of DNA. Some new anti-OI drug candidates have been synthesized with fused aromatic ring systems (e.g. carbazoles) that do not resemble the classical paradigm for minor-groove interactions at AT sequences in DNA. To characterize the DNA interactions of these compounds, we have used UV-vis absorbance, fluorescence, kinetic measurements, and circular dichroism in conjunction with NMR spectroscopy to evaluate the structure of the complexes formed between the carbazoles and DNA. Application of these methods to carbazoles substituted at either the 3,6 or 2,7 positions with cationic imidazoline groups gave conclusive, but very surprising, evidence that both compounds bind strongly in the minor groove at AT DNA sequences. NMR and molecular modeling of the complexes formed between the 3,6- and 2,7-carbazoles and the self-complementary oligomer d(GCGAATTCGC) have been used to establish structural details for the minor-groove complex. These results have been used as constraints for molecular modeling calculations to construct models of the minor-groove-carbazole complexes and to draw conclusions regarding the molecular basis for the effects of substituent position on carbazole-DNA affinities. The surprising result is that the 2,7 carbazole binds in AT sequences with hydrogen bonds involving one imidazoline group and the carbazole NH. The 3,6-carbazole compound binds in a more "classical" model that uses both imidazoline groups for H-bonding while the carbazole NH points out of the minor groove. The carbazoles thus form a new type of DNA minor groove complex and their excellent biological activities indicate that a variety of fused-ring minor-groove binding agents should be investigated.