DNA minor groove electrostatic potential: influence of sequence-specific transitions of the torsion angle gamma and deoxyribose conformations

The structural adjustments of the sugar-phosphate DNA backbone (switching of the γ angle (O5′–C5′–C4′–C3′) from canonical to alternative conformations and/or C2′-endo → C3′-endo transition of deoxyribose) lead to the sequence-specific changes in accessible surface area of both polar and non-polar atoms of the grooves and the polar/hydrophobic profile of the latter ones. The distribution of the minor groove electrostatic potential is likely to be changing as a result of such conformational rearrangements in sugar-phosphate DNA backbone. Our analysis of the crystal structures of the short free DNA fragments and calculation of their electrostatic potentials allowed us to determine: (1) the number of classical and alternative γ angle conformations in the free B-DNA; (2) changes in the minor groove electrostatic potential, depending on the conformation of the sugar-phosphate DNA backbone; (3) the effect of the DNA sequence on the minor groove electrostatic potential. We have demonstrated that the structural adjustments of the DNA double helix (the conformations of the sugar-phosphate backbone and the minor groove dimensions) induce changes in the distribution of the minor groove electrostatic potential and are sequence-specific. Therefore, these features of the minor groove sizes and distribution of minor groove electrostatic potential can be used as a signal for recognition of the target DNA sequence by protein in the implementation of the indirect readout mechanism.     

DNA conformational analysis in solution by uranyl mediated photocleavage.

Uranyl mediated photocleavage of double stranded DNA is proposed as a general probing for DNA helix conformation in terms of minor groove width/electronegative potential. Specifically, it is found that A/T-tracts known to constitute strong distamycin binding sites are preferentially photocleaved by uranyl in a way indicating strongest uranyl binding at the center of the minor groove of the AT-region. The A-tracts of kinetoplast DNA show the highest reactivity at the 3'-end of the tract--as opposed to cleavage by EDTA/Fell--in accordance with the minor groove being more narrow at this end. Finally, uranyl photocleavage of the internal control region (ICR) of the 5S-RNA gene yields a cleavage modulation pattern fully compatible with that obtained by DNase I which also--in a more complex way--senses DNA minor groove width.     

DNA and its counterions: a molecular dynamics study

The behaviour of mobile counterions, Na⁺ and K⁺, was analysed around a B-DNA double helix with the sequence CCATGCGCTGAC in aqueous solution during two 50 ns long molecular dynamics trajectories. The movement of both monovalent ions remains diffusive in the presence of DNA. Ions sample the complete space available during the simulation time, although individual ions sample only about one-third of the simulation box. Ions preferentially sample electronegative sites around DNA, but direct binding to DNA bases remains a rather rare event, with highest site occupancy values of <13%. The location of direct binding sites depends greatly on the nature of the counterion. While Na⁺ binding in both grooves is strongly sequence-dependent with the preferred binding site in the minor groove, K⁺ mainly visits the major groove and binds close to the centre of the oligomer. The electrostatic potential of an average DNA structure therefore cannot account for the ability of a site to bind a given cation; other factors must also play a role. An extensive analysis of the influence of counterions on DNA conformation showed no evidence of minor groove narrowing upon ion binding. A significant difference between the conformations of the double helix in the different simulations can be attributed to extensive α/γ transitions in the phosphate backbone during the simulation with Na⁺. These transitions, with lifetimes over tens of nanoseconds, however, appear to be correlated with ion binding to phosphates. The ion-specific conformational properties of DNA, hitherto largely overlooked, may play an important role in DNA recognition and binding.     

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

Location of spermine and other polyamines on DNA as revealed by photoaffinity cleavage with polyaminobenzenediazonium salts

Although polyamines interact strongly with nucleic acids, X-ray and NMR studies have not revealed much structural information about spermine-DNA complexes. Therefore, it was of interest to look at the binding of polyamines to 32P-labeled DNA restriction fragments by sequencing gel electrophoresis of the photoaffinity cleavage products induced by polyaminobenzendiazonium salts. The shift of cleavage patterns observed on opposite strands as well as competition experiments with distamycin shows polyamines to be located in the minor groove of B-DNA and to depend on the nucleic acid polymorphism, jumping to the major groove in the A-form. The sequence selectivities of various polycations (spermine, putrescine, and cobalt (III) hexaammine) are similar and slightly favor A,T-rich regions. Taken together, these results show that polycations which are not point charges are guided by the electronegative potential along the nucleic acid and suggest fast crawling of the polyamine within the minor groove, due to individual NH2+ jumping between multiple equidistant and isoenergetic bidentate hydrogen-bonding sites. Such a picture could be the clue to the unexpected NMR and to the frequently silent X-ray behavior of polyamines when bound to DNA.     

Sequence-specific interaction of the Salmonella Hin recombinase in both major and minor grooves of DNA.

The Hin recombinase of Salmonella catalyzes a site-specific recombination event which leads to flagellar phase variation. Starting with a fully symmetrical recombination site, hixC, a set of 40 recombination sites which vary by pairs of single base substitutions was constructed. This set was incorporated into the Salmonella-specific bacteriophage P22 based challenge phage selection and used to define the DNA sequence determinants for the binding of Hin to DNA in vivo. The critical sequence-specific contacts between a Hin monomer and a 13 bp hix half-site are at two T:A base pairs in the major groove of the DNA which are separated by one base pair, and two consecutive A:T contacts in the minor groove. The base substitutions in the major groove recognition portion which were defective in binding Hin still retained residual binding capability in vivo, while the base pair substitutions affecting the minor groove recognition region lost all in vivo binding. Using in vitro binding assays, Hin was found to bind to hix symmetrical sites with A:T base pairs or I:C base pairs in the minor groove recognition sequences, but not to G:C base pairs. In separate in vitro binding assays, Hin was equally defective in binding to either a G:C or a I:C contact in a major groove recognition sequence. Results from in vitro binding assays to hix sites in which 3-deazaadenine was substituted for adenine are consistent with Hin making a specific contact to either the N3 of adenine or O2 of thymine in the minor groove within the hix recombination site on each symmetric half-site. These results taken with the results of previous studies on the DNA binding domain of Hin suggest a sequence-specific minor groove DNA binding motif.     

Putative binding modes of Ku70-SAP domain with double strand DNA: a molecular modeling study

The channel structure of the Ku protein elegantly reveals the mechanistic basis of sequence-independent DNA-end binding, which is essential to genome integrity after exposure to ionizing radiation or in V(D)J recombination. However, contradicting evidence indicates that this protein is also involved in the regulation of gene expression and in other regulatory processes with intact chromosomes. This computational study predicts that a putative DNA binding domain of this protein, the SAP domain, can form DNA-bound complexes with relatively high affinities (ΔG ≈ -20 kcal mol^-1). The binding modes are searched by low frequency vibration modes driven by the fully flexible docking method while binding affinities are calculated by the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method. We find this well defined 5 kDa domain with a helix-extended loop-helix structure is suitable to form favorable electrostatic and hydrophobic interactions with either the major groove or the minor groove of DNA. The calculation also reveals the sequence specified binding preference which may relate to the observed pause sites when Ku translocates along DNA and the perplex binding of Ku with circular DNA.     

Specific minor groove solvation is a crucial determinant of DNA binding site recognition

The DNA sequence preferences of nearly all sequence specific DNA binding proteins are influenced by the identities of bases that are not directly contacted by protein. Discrimination between non-contacted base sequences is commonly based on the differential abilities of DNA sequences to allow narrowing of the DNA minor groove. However, the factors that govern the propensity of minor groove narrowing are not completely understood. Here we show that the differential abilities of various DNA sequences to support formation of a highly ordered and stable minor groove solvation network are a key determinant of non-contacted base recognition by a sequence-specific binding protein. In addition, disrupting the solvent network in the non-contacted region of the binding site alters the protein''s ability to recognize contacted base sequences at positions 5–6 bases away. This observation suggests that DNA solvent interactions link contacted and non-contacted base recognition by the protein.     

Interchange of DNA-binding modes in the deformed and ultrabithorax homeodomains: a structural role for the N-terminal arm

The deformed (Dfd) and ultrabithorax (Ubx) homeoproteins regulate developmental gene expression in Drosophila melanogaster by binding to specific DNA sequences within its genome. DNA binding is largely accomplished via a highly conserved helix-turn-helix DNA-binding domain that is known as a homeodomain (HD). Despite nearly identical DNA recognition helices and similar target DNA sequence preferences, the in vivo functions of the two proteins are quite different. We have previously revealed differences between the two HDs in their interactions with DNA. In an effort to define the individual roles of the HD N-terminal arm and recognition helix in sequence-specific binding, we have characterized the structural details of two Dfd/Ubx chimeric HDs in complex with both the Dfd and Ubx-optimal-binding site sequences. We utilized hydroxyl radical cleavage of DNA to assess the positioning of the proteins on the binding sites. The effects of missing nucleosides and purine methylation on HD binding were also analyzed. Our results show that both the Dfd and Ubx HDs have similar DNA-binding modes when in complex with the Ubx-optimal site. There are subtle but reproducible differences in these modes that are completely interchanged when the Dfd N-terminal arm is replaced with the corresponding region of the Ubx HD. In contrast, we showed previously that the Dfd-optimal site sequence elicits a very different binding mode for the Ubx HD, while the Dfd HD maintains a mode similar to that elicited by the Ubx-optimal site. Our current methylation interference studies suggest that this alternate binding mode involves interaction of the Ubx N-terminal arm with the minor groove on the opposite face of DNA relative to the major groove that is occupied by the recognition helix. As judged by hydroxyl radical footprinting and the missing nucleoside experiment, it appears that interaction of the Ubx recognition helix with the DNA major groove is reduced. Replacing the Dfd N-terminal arm with that of Ubx does not elicit a complete interchange of the DNA-binding mode. Although the position of the chimera relative to DNA, as judged by hydroxyl radical footprinting, is similar to that of the Dfd HD, the missing nucleoside and methylation interference patterns resemble those of the Ubx HD. Repositioning of amino acid side-chains without wholesale structural alteration in the polypeptide appears to occur as a function of N-terminal arm identity and DNA-binding site sequence. Complete interchange of binding modes was achieved only by replacement of the Dfd N-terminal arm and the recognition helix plus 13 carboxyl-terminal residues with the corresponding residues of Ubx. The position of the N-terminal arm in the DNA minor groove appears to differ in a manner that depends on the two base-pair differences between the Dfd and Ubx-optimal-binding sites. Thus, N-terminal arm position dictates the binding mode and the interaction of the recognition helix with nucleosides in the major groove.     

A negative electrostatic determinant mediates the association between the Escherichia coli trp repressor and its operator DNA.

The electrostatic potential surfaces were characterized for trp repressor models that bind to DNA with sequence specificity, without specificity, and not at all. Comparisons among the surfaces were used to isolate protein surface features likely to be important in DNA binding. Models that differ in protein conformation and tryptophan-analogue binding consistently showed positive potential associated with the protein surfaces that interact with the DNA major groove. However, negative potential is associated with the trp repressor surface that contacts the DNA minor groove. This negative potential is significantly neutralized in the protein conformation that is bound to DNA. Positive potential is also associated with the tryptophan binding-site surface, a consequence of the tryptophan- or tryptophan analogue-induced allosteric change. This protein region is complementary to the strongest negative potential associated with the DNA phosphate backbone and is also present in the isolated protein structure from the protein-DNA complex. The effects of charge-change mutation, pH dependence, and salt dependence on the electrostatic potential surfaces were also examined with regard to their effects on protein-DNA binding constants. A consistent model is formed that defines a role for long-range electrostatics early in the protein-DNA association process and complements previous structural, molecular association, and mutagenesis studies.     

Control of DNA minor groove width and Fis protein binding by the purine 2-amino group

The width of the DNA minor groove varies with sequence and can be a major determinant of DNA shape recognition by proteins. For example, the minor groove within the center of the Fis–DNA complex narrows to about half the mean minor groove width of canonical B-form DNA to fit onto the protein surface. G/C base pairs within this segment, which is not contacted by the Fis protein, reduce binding affinities up to 2000-fold over A/T-rich sequences. We show here through multiple X-ray structures and binding properties of Fis–DNA complexes containing base analogs that the 2-amino group on guanine is the primary molecular determinant controlling minor groove widths. Molecular dynamics simulations of free-DNA targets with canonical and modified bases further demonstrate that sequence-dependent narrowing of minor groove widths is modulated almost entirely by the presence of purine 2-amino groups. We also provide evidence that protein-mediated phosphate neutralization facilitates minor groove compression and is particularly important for binding to non-optimally shaped DNA duplexes.     

DNA shape, genetic codes, and evolution

Although the three-letter genetic code that maps nucleotide sequence to protein sequence is well known, there must exist other codes that are embedded in the human genome. Recent work points to sequence-dependent variation in DNA shape as one mechanism by which regulatory and other information could be encoded in DNA. Recent advances include the discovery of shape-dependent recognition of DNA that depends on minor groove width and electrostatics, the existence of overlapping codes in protein-coding regions of the genome, and evolutionary selection for compensatory changes in nucleotide composition that facilitate nucleosome occupancy. It is becoming clear that DNA shape is important to biological function, and therefore will be subject to evolutionary constraint.     

Comparative thermodynamics for monomer and dimer sequence-dependent binding of a heterocyclic dication in the DNA minor groove

Phenylamidine cationic groups linked by a furan ring (furamidine) and related symmetric diamidine compounds bind as monomers in the minor groove of AT sequences of DNA. DB293, an unsymmetric derivative with one of the phenyl rings of furamidine replaced with a benzimidazole, can bind to AT sequences as a monomer but binds more strongly to GC-containing minor-groove DNA sites as a stacked dimer. The dimer-binding mode has high affinity, is highly cooperative and sequence selective. In order to develop a better understanding of the correlation between structural and thermodynamic aspects of DNA molecular recognition, DB293 was used as a model to compare the binding of minor-groove agents with AT and mixed sequence DNA sites. Isothermal titration calorimetry and surface plasmon resonance results clearly show that the binding of DB293 and other related compounds into the minor groove of AT sequences is largely entropy-driven while the binding of DB293 as a dimer into the minor groove of GC-containing sequences is largely enthalpy-driven. At 25 degrees C, for example, the AT binding has DeltaG degrees, DeltaH degrees and TDeltaS degrees values of -9.6, -3.6 and 6.0 kcal/mol while the values for dimer binding to a GC-containing site are -9.0, -10.9 and -1.9 kcal/mol (per mol of bound compound), respectively. These results show that the thermodynamic components for binding of compounds of this type to DNA are very dependent on the structure, solvation and sequence of the DNA binding site.     

Selective Binding of Synthetic Polypeptides to DNA of Varying Composition and Sequence: Effect of Minor Groove Binding Drugs

Abstract

Repetitive basic polypeptides containing lysine or arginine as every third amino acid were shown to cause DNA condensation at physiological salt concentration connected with selective DNA binding with respect to DNA composition and sequence. This selectivity is very similar to that existing in the case of histone H1 and other basic proteins and does not depend on polypeptide chain conformation. The effect of the minor groove binding drugs netropsin and distamycin was tested to elucidate the origin of the binding selectivity. The results suggest that the binding preferences are due to the variations in the conformation in various types of B-DNA that depend on DNA composition and sequence. The most important factor affecting the selectivity is probably the value of the negative electrostatic potential in the minor groove.