Minor groove binding of a bis-quaternary ammonium compound: the crystal structure of SN 7167 bound to d(CGCGAATTCGCG)2
Minor groove binding of a bis-quaternary ammonium compound: the crystal structure of SN 7167 bound to d(CGCGAATTCGCG) 2 Christopher J. Squire*, George R. Clark and William A. Denny1
Chemistry Department and 1Cancer Research Laboratory, School of Medicine, University of Auckland, Auckland, New Zealand
Received June 20, 1997;Revised and Accepted August 26, 1997
NDB no. GDL053
ABSTRACT
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 Å 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.
INTRODUCTION
Compounds which bind reversibly in the minor groove of DNA possess a wide spectrum of biological behaviour, including antimicrobial, antiprotozoal and anticancer activities. It has been suggested that these effects are mediated by blocking transcription (1 -4 ), by inhibition of polyamine biosynthesis (5 ) and/or by inhibition of topoisomerase I (6 ) or II (7 ,8 ). Minor groove binders also show excellent molecular recognition properties, binding with high specificity to AT-rich sequences up to 6 bp in length. For this reason, they are also of interest as carriers of intrinsically less specific functionalities, such as nitrogen mustard alkylating agents, and the resultant `DNA-targeted' alkylators have been shown to possess significantly modified patterns of DNA adduct formation (9 -12 ). The minor groove binding motif is also evident in several classes of highly cytotoxic natural products which are DNA sequence-specific alkylators, such as the anthramycins (13 ) and duocarmycins (14 ).
For these reasons there is wide interest in the molecular level interactions of reversible minor groove binders with DNA, to determine if there are common structural features for this binding mode. Particular attention has been paid to the polypyrrole and bisbenzimidazole series, with many NMR (15 ,16 ) and X-ray (17 -19 ) studies of selected analogues with a range of oligodeoxynucleotides.
Another class of reversible minor groove binders which has received less attention in structural studies is that which can be loosely called the polybenzamides. An early example of this class, a terephthalanilide, was evaluated clinically as an anticancer drug (20 ). Later structure-activity studies on bis-quaternary salt analogues showed a positive correlation between experimental antileukaemic activity and selectivity of binding to AT-rich over GC-rich DNA sequences (21 ). The bis-quaternary cation 4-[4-[4-[4(N-methylquinolinium)amino]benzamido]anilino]-N-methylpyridinium (SN 6999, Fig. 1 ) and analogues have also been shown to be very potent inhibitors of multidrug-resistant strains of the malaria parasite Plasmodium falciparum (22 ). Recent work has shown that polybenzamides are suitable DNA targeting carriers for both DNA cleaving (23 ) and DNA alkylating (10 ,24 ) moieties.
Most structural work on polybenzamide DNA binding has been carried out on SN 6999 and an amino analogue SN 7167 (Fig. 1 ). NMR studies confirmed that SN 6999 binds preferentially to AT-rich oligonucleotides in the minor groove and showed the formation of specific hydrogen bonds (25 -27 ). An X-ray crystallographic study (28 ) of SN 6999 complexed to the modified oligodeoxynucleotide d(CGC[e6G]AATTCGCG)2 showed the ligand bound in the minor groove, but offset from the centre of the 5'-AATT tract. The structure was unusual in that the DNA was very bent and in a manner not seen before, even in other groove binder-[e6G] complexes. It appears that the semi-rigid nature of the SN 6999 molecule forced this extreme distortion of the DNA to occur, with the distortion being facilitated by the weaker hydrogen bonding in the [e6G]·C base pair. The [e6G] base contains an ethoxy substitution of O6. The hydrogen bond from O6 to cytosine N4 is eliminated by this substitution and the [e6G] base therefore makes only two hydrogen bonds to cytosine, instead of the normal three.
MATERIALS AND METHODS
Synthesis and crystallization
The DNA dodecamer d(CGCGAATTCGCG)2 was purchased from Oligos Etc. Inc. and annealed before use. The drug SN 7167 was prepared as previously reported (29 ). Crystals of the complex (yellow prisms) were grown from hanging drops at room temperature. The crystal used for data collection was grown at room temperature from a drop containing 4.8 µl 50% 2-methylpentane-2,4-diol, 4 µl 30 mM sodium cacodylate, pH 6.5, 4.8 µl 5 mM SN 7167, 3 µl 200 mM MgCl2, 2 µl 5 mM spermine hydrochloride and 4.8 µl 5 mM dodecamer. The droplet was equilibrated against a reservoir containing 1 ml 50% 2-methyl-2,4-pentanediol. An X-ray quality crystal was obtained after ~3 months from a gelatinous precipitate of the complex. The crystallization process was not highly reproducible because of formation of this precipitate.
Data collection
The single crystal used for data collection was of approximate dimensions 0.3 × 0.2 × 0.1 mm and was mounted in a 0.7 mm Lindemann quartz capillary with a small amount of mother liquor. Diffraction data were collected at 293 K using a Rigaku R-axis 2C diffractometer equipped with two image plate area detectors, a rotating anode X-ray generator and a graphite monochromator. Oscillation data was collected with a crystal to detector distance of 10 cm and an oscillation angle of 2o. Forty four frames were collected through a total angle of 88o of rotation to obtain data to a maximum resolution of 2.6 Å. The crystal showed no significant decay during data collection. Data processing was performed using DENZO (30 ). After merging, the data set comprised 1955 unique reflections (91.4% complete) with a merging R value of 9.5%.
Structure refinement
The unit cell dimensions of the crystal are a = 24.56, b = 40.15 and c = 65.68 Å, in the orthorhombic space group P 212121. This cell is close to that reported for related native dodecamers and groove binder complexes and the coordinates for the native structure (PDB entry 9BNA; 31 ) were used as the initial model for refinement. The refinement process was carried out using the program X-PLOR v.3.1 (32 ). Rigid body refinement of the whole duplex model was carried out with the resolution range of the data gradually increased from 8.0-4.0 Å (496 reflections) to 8.0-3.0 Å (1198 reflections). The DNA duplex was then divided into 24 rigid groups, comprising the 22 nucleotide and 2 nucleoside units, and the upper resolution limit gradually increased to the maximum of 2.6 Å (1888 reflections) during rigid body refinement. The R factor converged at 36.0%. Positional and temperature factor refinement further reduced the R factor to 29.3%. Electron density maps were calculated and displayed using the graphics package Alberta/Caltech TOM v.3.1. The DNA molecule showed a good fit to the electron density in all regions of the model and a continuous lobe of density located within the minor groove (Fig. 2 ) could be assigned to the SN 7167 molecule. Coordinates of the ligand SN 6999 were obtained from the published structure of the SN 6999- d(CGC[e6G]AATTCGCG)2 complex and were altered using the computer modelling package InsightII (33 ) to create a molecule of SN 7167. The electron density in the groove showed sufficient detail to allow the ligand to be placed in a unique direction and position in the groove. There was some uncertainty in the rotational conformation of the BQ-amido moiety and two possible ligand forms were refined further. The first rotamer has N37 pointing directly out of the groove (this can be called the `amino out' conformation). The second rotamer has the same amino group facing the floor of the groove, the BQ-amido moiety having rotated 180o into the groove (`amino in'). Partial charges were calculated for the SN 7167 molecule by empirical calculations using the cvff (default) forcefield of InsightII. Planar restraints were applied to each individual aromatic ring system. Positional and temperature factor refinements were continued for both models with the SN 7167 molecules included in the calculations. Both `amino in' and `amino out' structures refined to the same R factor value. At this stage the ligand conformations and DNA binding were examined in more detail to try to elucidate the correct structure. The `amino out' structure was chemically and stereochemically sensible and showed a good fit to the experimental electron density. The `amino in' rotamer appeared to fit the electron density less well and the molecule looked strained into a distorted conformation. An unfavourable ligand conformation in itself is not enough evidence to say that this structure is not correct, as the nature of the DNA and the groove binding process invariably result in conformations that are different from the minimum energy conformations favoured in solution. However, the `amino in' conformation was excluded from further refinement because of the presence of a repulsive contact of 2.6 Å between the ligand atom O19 and O2-T19. Modelling of this section of the structure using the Amber forcefield of InsightII indicated that this contact is strongly repulsive in nature and that these atoms would be forced apart to at least 3 Å. The resulting ligand conformation would not fit the experimental electron density well. Because of these considerations, the `amino in' rotamer was discarded as a likely conformation and only the `amino out' rotamer complex was refined further.
RESULTS AND DISCUSSION
Structure of the complex
The DNA duplex adopts a B-DNA conformation, with the SN 7167 ligand lying in the narrow central region of the minor groove but offset towards one end (Fig. 3 ). The overall crystal packing is analogous to that observed for the native sequence and the various groove binder-dodecamer complexes. In all these structures the two terminal GC base pairs of each helix interact with a neighbouring duplex with the minor grooves slotting into each other and being involved in inter-duplex hydrogen bonding. The ligand binds in the minor groove over the 5'-AATTC site with the quinolinium ring system (Q) lying towards the 5'-end of the sequence. The quinolinium extends beyond the 5'-A site towards the adjacent guanine and the pyridinium (P) extends to the C9 base at the 3'-end, and so the SN 7167 molecule binds over a site ~5[1/2] bp in length. There is no evidence of disorder or mobility along the sequence and the ligand adopts a unique position and orientation in the groove.
Ligand-DNA interactions
The nitrogen of the SN 7167 amido linker faces the floor of the minor groove and hydrogen bonds to O2-T19. The hydrogen bonding distance at 3.4 Å is probably too long to be considered a true hydrogen bond but it does indicate a weak interaction. The two amine linkers of the ligand have their hydrogens pointing directly out of the groove and do not hydrogen bond to the DNA. There are 15 contacts between the ligand and the minor groove which are shorter than 3.5 Å (Table 1 ) and all are long range at 3.1-3.4 Å. The ligand does not show a strong preference to interact with one of the DNA strands over the other. The shortest contacts involve the two rings on the ends of the molecule, although the pyridinium ring (P) has only one close contact with the DNA. There is also a possible charge-induced dipolar interaction between the quinolinium atom N1 and O4'-A6, but with a distance of 3.3 Å this interaction can also be considered weak.
The SN 7167 molecule is approximately planar overall with relatively small dihedral angles between ring planes (Table 2 ). The largest dihedral angle found, between the two central rings (BQ and BP) at 17o, is interesting in that co-planarity of the two ring systems linked by the trans amide would normally be favoured. The other two ring systems at the ends of the molecule have more free rotation about their amine linkers but require relatively small twists in order to follow the groove. Average temperature factors for the individual ring systems of the SN 7167 molecule show that the BQ ring is held somewhat more rigidly in place than the remainder of the molecule. The BQ ring has an average atomic temperature factor of 24 Å2, while the remaining systems have values of 43, 41 and 51 Å2 for the Q, BP and P rings respectively. This low value of 24 Å2 could be indicative of a lack of rotational mobility for the BQ ring because of the exocyclic amino group substitution.
29 Atwell,G.J. and Cain,B.F. (1974) J. Med. Chem., 17, 930-934.MEDLINE Abstract
30 Zbyszek,O. (1993) In Sawyer,L., Isaacs,N. and Bailey,S. (eds), Proceedings of the CCP4 Study Weekend: `Data Collection and Processing'. SERC Laboratory, Daresbury, UK, pp. 56-62.