Rational Design of Pyrrolo[1,2-a]benzimidazole-Based Antitumor Agents Targeting the DNA Major Groove

Models have been developed for the interaction of the pyrrolo[1,2-a]benzimidazole (PBI) antitumor agents with the two-electron activating enzyme DT-diaphorase and the DNA major groove. The DT-diaphorase model and experimental results indicate that the S-enantiomer of 3-carbamido PBI can enantioselect ovarian cancers. The reduced PBI interacts with the DNA major groove at AT base pairs by forming Hoogsteen-like hydrogen bonds. The reduced 3-amino PBI forms three hydrogen bonds in the major groove with the amino group acting as an H-bond donor to the thymine carbonyl. The DNA-binding model will permit the design of major groove recognition agents.     

Ultrasonic cleavage of DNA in complexes with Ag(I), Cu(II), Hg(II)

We investigated a phenomenon of ultrasonic cleavage of DNA complexed with transition metal cations Ag(I), Cu(II) and Hg(II). We found the statistically significant dependence of relative intensity of cleavage on cation type and concentration. Each cation may cause two different types of distortion in the DNA double-helix depending on whether it binds to major or minor DNA groove. The intensity of ultrasonic cleavage decreases where the cation binds to the major DNA groove; the intensity of cleavage increases where the cation binds to the minor DNA groove and disturbs the hydrogen bonds of complementary base pairs or where it intercalates between bases. Both types of DNA distortion can affect the intensity of N?S intercon-version of deoxyribose.     

Theoretical approaches for evaluation of DNA-protein interactions

The notion of the DNA-recognizing segment of a protein for three-dimensional structures of DNA-protein complexes was formalized. Algorithms for calculating two parameters, hydrogen bonds and hydrophobic contacts, that characterize the interaction of the DNA-recognizing segment of the protein with the groove DNA major were proposed. DNA-recognizing protein segments in three-dimensional structures (of complexes) were classified according to these two parameters. These data were compared with the classification of the corresponding DNA-recognizing domains according to structure. The contribution of each pair amino acid residue-base to the interaction of protein with the DNA major groove was calculated.     

A study of the interaction of Cr(III) complexes and their selective binding with B-DNA: a molecular modeling approach

Molecular modeling and energy minimisation calculations have been used to investigate the interaction of chromium(III) complexes in different ligand environments with various sequences of B-DNA. The complexes are [Cr(salen)(H(2)O)(2)](+); salen denotes 1, 2 bis-salicylideneaminoethane, [Cr(salprn)(H(2)O)(2)](+); salprn denotes 1, 3 bis- salicylideneaminopropane, [Cr(phen)(3)](3+); phen denotes 1, 10 phenanthroline and [Cr(en)(3)](3+); en denotes ethylenediamine. All the chromium(III) complexes are interacted with the minor groove and major groove of d(AT)(12), d(CGCGAATTCGCG)(2) and d(GC)(12) sequences of DNA. The binding energy and hydrogen bond parameters of DNA-Cr complex adduct in both the groove have been determined using molecular mechanics approach. The binding energy and formation of hydrogen bonds between chromium(III) complex and DNA has shown that all complexes of chromium(III) prefer minor groove interaction as the favourable binding mode.     

Systematic Analysis of Possible Hydrogen Bonds between Amino Acid Side Chains and B-form DNA

Abstract

The ways in which amino acid side chains could make a pair of hydrogen bonds within the major groove of B DNA are systematically analyzed. Hydrogen bond donors within the major groove are characterized by determining the idealized position of the hydrogen bond acceptors that they might bond with. It appears that an amino acid side chain could, at most, contact a pair of base pairs. The ten possible pairs of base pairs are analyzed to determine how they could be recognized by the amino acid side chains.     

Double helix conformation, groove dimensions and ligand binding potential of a G/C stretch in B-DNA.

The self-complementary DNA fragment CCGGCGCCGG crystallizes in the rhombohedral space group R3 with unit cell parameters a = 54.07 A and c = 44.59 A. The structure has been determined by X-ray diffraction methods at 2.2 A resolution and refined to an R value of 16.7%. In the crystal, the decamer forms B-DNA double helices with characteristic groove dimensions: compared with B-DNA of random sequence, the minor groove is wide and deep and the major groove is rather shallow. Local base pair geometries and stacking patterns are within the range commonly observed in B-DNA crystal structures. The duplex bears no resemblance to A-form DNA as might have been expected for a sequence with only GC base pairs. The shallow major groove permits an unusual crystal packing pattern with several direct intermolecular hydrogen bonds between phosphate oxygens and cytosine amino groups. In addition, decameric duplexes form quasi-infinite double helices in the crystal by end-to-end stacking. The groove geometries and accessibilities of this molecule as observed in the crystal may be important for the mode of binding of both proteins and drug molecules to G/C stretches in DNA.     

CG base pair recognition within DNA triple helices using N-methyl-3H-pyrrolo[2,3-d]pyrimidin-2(7H)-one nucleoside analogues

Triplex-mediated recognition of Py.Pu base pairs in DNA is a greater challenge than for Pu.Py base pairs as fewer hydrogen bonds are presented for binding in the major groove. Initial studies on m-aminophenyl-modified analogues of the bicyclic nucleoside N-methyl-3H-pyrrolo[2,3-d]pyrimidin-2(7 H)-one suggest that selective recognition of the CG base pair is possible.     

Extra-helical guanine interactions in DNA

We review the extra-helical guanine interactions present in many oligonucleotide crystals. Very often terminal guanines interact with other guanines in the minor groove of neighboring oligonucleotides through N2 x N3 hydrogen bonds. In other cases the interaction occurs with the help of Ni2+ ions. Guanine/netropsin stacking in the minor groove has also been found. From these studies we conclude that guanine may have multiple extra-helical interactions. In particular it may be considered a very effective minor groove binder, which could be used in the design of sequence selective binding drugs. Interactions through the major groove are seldom encountered, but might be present when DNA is stretched. Such interactions are also analyzed, since they might be important for homologous chromosome pairing during meiosis.     

The propeller DNA conformation of poly(dA).poly(dT).

Physical properties of the DNA duplex, poly(dA).poly(dT) differ considerably from the alternating copolymer poly(dAT). A number of molecular models have been used to describe these structures obtained from fiber X-ray diffraction data. The recent solutions of single crystal DNA dodecamer structures with segments of oligo-A.oligo-T have revealed the presence of a high propeller twist in the AT regions which is stabilized by the formation of bifurcated (three-center) hydrogen bonds on the floor of the major groove, involving the N6 amino group of adenine hydrogen bonding to two O4 atoms of adjacent thymine residues on the opposite strand. Here we show that it is possible to incorporate the features of the single crystal analysis, specifically high propeller twist, bifurcated hydrogen bonds, and a narrow minor groove, as well as the close interstrand NMR signal between adenine HC2 and ribose HC1' of the opposite strand, into a model that is fully compatible with the diffraction data obtained from poly(dA).poly(dT).     

Mg(2+) in the Major Groove Modulates B-DNA Structure and Dynamics

This study investigates the effect of Mg(2+) bound to the DNA major groove on DNA structure and dynamics. The analysis of a comprehensive dataset of B-DNA crystallographic structures shows that divalent cations are preferentially located in the DNA major groove where they interact with successive bases of (A/G)pG and the phosphate group of 5'-CpA or TpG. Based on this knowledge, molecular dynamics simulations were carried out on a DNA oligomer without or with Mg(2+) close to an ApG step. These simulations showed that the hydrated Mg(2+) forms a stable intra-strand cross-link between the two purines in solution. ApG generates an electrostatic potential in the major groove that is particularly attractive for cations; its intrinsic conformation is well-adapted to the formation of water-mediated hydrogen bonds with Mg(2+). The binding of Mg(2+) modulates the behavior of the 5'-neighboring step by increasing the BII (ε-ζ>0°) population of its phosphate group. Additional electrostatic interactions between the 5'-phosphate group and Mg(2+) strengthen both the DNA-cation binding and the BII character of the 5'-step. Cation binding in the major groove may therefore locally influence the DNA conformational landscape, suggesting a possible avenue for better understanding how strong DNA distortions can be stabilized in protein-DNA complexes.     

Structure of a bis-amidinium derivative of hoechst 33258 complexed to dodecanucleotide d(CGCGAATTCGCG)2: the role of hydrogen bonding in minor groove drug-DNA recognition.

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.     

Wing 1 of protein HOP2 is as important as helix 3 in DNA binding by MD simulation

The repair of programmed DNA double-strand breaks through recombination is required for proper association and disjunction of the meiotic homologous chromosomes. Meiosis-specific protein HOP2 plays essential roles in recombination by promoting recombinase activities. The N-terminal domain of HOP2 interacts with DNA through helix 3 (H3) and wing 1 (W1). Mutations in wing 1 (Y65A/K67A/Q68A) slightly weakened the binding but mutations in helices 2 and 3 (Q30A/K44A/K49A) nearly abolished the binding. To better understand such differential effects at atomic level, molecular dynamics simulations were employed. Despite losing some hydrogen bonds, the W1-mutant DNA complex was rescued by stronger hydrophobic interactions. For the wild type and W1-mutant, the protein was found to slide along the DNA grooves as the DNA rolls along its double-helix axis. This motion could be functionally important to facilitate the precise positioning of the single-stranded DNA with the homologous double-stranded DNA. The sliding motion was reduced in the W1-mutant. The H-mutant nearly lost all intermolecular interactions. Moreover, an additional mutation in wing 1 (Y65A/K67A/Q68A/K69A) also caused complete complex dissociation. Therefore, both wing 1 and helix 3 make important contribution to the DNA binding, which could be important to the strand invasion function of HOP2 homodimer and HOP2-MND1 heterodimer. Similar to cocking a medieval crossbow with the archer’s foot placed in the stirrup, wing 1 may push the minor groove to cause distortion while helix 3 grabs the major groove.     

Structural basis of replication origin recognition by the DnaA protein

Escherichia coli DnaA binds to 9 bp sequences (DnaA boxes) in the replication origin, oriC, to form a complex initiating chromosomal DNA replication. In the present study, we determined the crystal structure of its DNA-binding domain (domain IV) complexed with a DnaA box at 2.1 Å resolution. DnaA domain IV contains a helix–turn–helix motif for DNA binding. One helix and a loop of the helix– turn–helix motif are inserted into the major groove and 5 bp (3′ two-thirds of the DnaA box sequence) are recognized through base-specific hydrogen bonds and van der Waals contacts with the C5-methyl groups of thymines. In the minor groove, Arg399, located in the loop adjacent to the motif, recognizes three more base pairs (5′ one-third of the DnaA box sequence) by base-specific hydrogen bonds. DNA bending by ~28° was also observed in the complex. These base-specific interactions explain how DnaA exhibits higher affinity for the strong DnaA boxes (R1, R2 and R4) than the weak DnaA boxes (R3 and M) in the replication origin.     

C-H.O hydrogen bonds in minor groove of A-tracts in DNA double helices

Analysis of available B-DNA type oligomeric crystal structures as well as protein-bound DNA fragments (solved using data with resolution <2.6 A) indicates that in both data sets, a majority of the (3'-Ade) H2..O2(3'-Thy/Cyt) distances in AA.TT and GA.TC dinucleotide steps, are considerably shorter than their values in a uniform fibre model, and are smaller than their optimum separation distance. Since the electropositive C2-H2 group of adenine is in close proximity of the electronegative keto oxygen atoms of both pyrimidine bases in the antiparallel strand of the double-helical DNA structures, it suggests the possibility of intra-base-pair as well as cross-strand C-H..O hydrogen bonds in the minor groove. The C2-H2..O2 hydrogen bonds within the A.T base-pairs could be a natural consequence of Watson-Crick pairing. However, the close cross-strand interactions between the bases at the 3'-ends of the AA.TT and GA.TC steps arise due to the local sequence-dependent geometry of these steps. While the base-pair propeller twist in these steps is comparable to the fibre model, some of the other local parameters such as base-pair opening angle and inter-base-pair slide show coordinated changes, leading to these shorter C2-H2..O2 distances. Hence, in addition to the well-known minor groove hydration, it appears that favourable C2-H2..O2 cross-strand interactions may play a role in imparting a characteristic geometry to AA.TT and GA.TC steps, as well as An.Tn and GAn.TnC tracts, which leads to a narrow minor groove in these regions.     

Functional hydrogen-bonding map of the minor groove binding tracks of six DNA polymerases

Recent studies have identified amino acid side chains forming several hydrogen bonds in the DNA minor groove as potentially important in polymerase replication of DNA. Few studies have probed these interactions on the DNA itself. Using non-hydrogen-bonding nucleoside isosteres, we have now studied effects in both primer and template strands with several polymerases to investigate the general importance of these interactions. All six polymerases show differences in the H-bonding effects in the minor groove. Two broad classes of activity are seen, with a first group of DNA polymerases (KF(-), Taq, and HIV-RT) that efficiently extends nonpolar base pairs containing nucleoside Q (9-methyl-1H-imidazo[4,5-b]pyridine) but not the analogue Z (4-methylbenzimidazole), implicating a specific minor groove interaction at the first extension site. A second group of polymerases (Pol alpha, Pol beta, and T7(-)) fails to extend all non-H-bonding base pairs, indicating that these enzymes may need minor groove hydrogen bonds at both minor groove sites or that they are especially sensitive to noncanonical DNA structure or stability. All DNA polymerases examined use energetically important minor groove interactions to probe newly synthesized base pairs before extending them. The positions of these interactions vary among the enzymes, and only a subset of the interactions identified structurally appears to be functionally important. In addition, polymerases appear to be differently sensitive to small changes in base pair geometry.     

Quantitative structure of a complex between a minor-groove-specific drug and a bent DNA decamer duplex: use of 2D NMR data and NOESY constrained energy minimization

Two-dimensional nuclear magnetic resonance (2D NMR) studies on d(GA4T4C)2 and d(GT4A4C)2 [Sarma, M.H., et al. (1988) Biochemistry 27, 3423-3432; Gupta G., et al. (1988) Biochemistry 27, 7909-7919] showed that A.T pairs are propeller twisted. As a result, A/T tracts form a straight rigid structural block with an array of bifurcated inter base pair H bonds in the major groove. It was demonstrated (previous paper) that replacement of methyl group by hydrogen (changing from T to U) in the major groove does not disrupt the array of bifurcated H bonds in the major groove. In this article, we summarize results of 2D NMR and molecular mechanic studies on the effect of a minor-groove-binding A.T-specific drug on the structure d(GA4T4C)2. A distamycin analogue (Dst2) was used for this study. It is shown that Dst2 binds to the minor groove of d(GA4T4C)2 mainly driven by van der Waals interaction between A.T pairs and the drug; as a consequence, an array of bifurcated H bonds can be formed in the minor groove between amide/amino protons of Dst2 and A.T pairs of DNA. NOESY data suggest that Dst2 predominantly binds at the central 5 A.T pairs. NOESY data also reveal that, upon drug binding, d(GA4T4C)2 does not undergo any significant change in conformation from the free state; i.e., propeller-twisted A.T pairs are still present in DNA and hence the array of bifurcated H bonds must be preserved in the major groove. NOESY data for the A5-T6 sequence also indicate that there is little change in junction stereochemistry upon drug binding.     

Molecular Dynamics Simulation in Solvent of the Estrogen Receptor Protein DNA Binding Domain in Complex with a Non-Consensus Estrogen Response Element DNA Sequence

Abstract

We investigated protein/DNA interactions, using molecular dynamics simulations computed between a 10 Angstom water layer model of the estrogen receptor (ER) protein DNA binding domain (DBD) amino acids and DNA of a non-consensus estrogen response element (ERE) consisting of 29 nucleotide base pairs. This ERE nucleotide sequence occurs naturally upstream of the Xenopus laevis Vitelligenin AI gene. The ER DBD is encoded by three exons. Namely, exons 2 and 3 which encode the two zinc binding motifs and a sequence of exon 4 which encodes a predicted alpha helix. We generated a computer model of the ER DBD using atomic coordinates derived from the average of 30 nuclear magnetic resonance (NMR) spectroscopy coordinate sets. Amino acids on the carboxyl end of the ER DBD were disordered in both X-ray crystallography and NMR determinations and no coordinates were reported. This disordered region includes 10 amino acids of a predicted alpha helix encoded in exon 4 at the exon 3/4 splice junction. These amino acids are known to be important in DNA binding and are also believed to function as a nuclear translocation signal sequence for the ER protein. We generated a computer model of the predicted alpha helix consisting of the 10 amino acids encoded in exon 4 and attached this helix to the carboxyl end of the ER DBD at the exon 3/4 splice junction site. We docked the ER DBD model within the DNA major groove halfsites of the 29 base pair non-consensus ERE and flanking nucleotides. We constructed a solvated model with the ER DBD/ERE complex surrounded by a ten Angstrom water layer and conducted molecular dynamics simulations. Hydrogen bonding interactions were monitored. In addition, van der Waals and electrostatic interaction energies were calculated. Amino acids of the ER DBD DNA recognition helix formed both direct and water mediated hydrogen bonds at cognate codon-anticodon nucleotide base and backbone sites within the ERE DNA right major groove halfsite. Amino acids of the ER DBD exon 4 encoded predicted alpha helix formed direct and water mediated H-bonds with base and backbone sites of their cognate codon-anticodon nucleotides within the minor grooves flanking the ERE DNA major groove halfsites. These interactions together induced bending of the DNA into the protein.