Molecular structure of the netropsin-d(CGCGATATCGCG) complex: DNA conformation in an alternating AT segment

The molecular structure of the complex between a minor groove binding drug (netropsin) and the DNA dodecamer d(CGCGATATCGCG) has been solved and refined by single-crystal X-ray diffraction analysis to a final R factor of 20.0% to 2.4-A resolution. The crystal is similar to that of the other related dodecamers with unit cell dimensions of a = 25.48 A, b = 41.26 A, and c = 66.88 A in the space group P2(1)2(1)2(1). In the complex, netropsin binds to the central ATAT tetranucleotide segment in the narrow minor groove of the dodecamer B-DNA double helix as expected. However, in the structural refinement the drug is found to fit the electron density in two orientations equally well, suggesting the disordered model. This agrees with the results from solution studies (chemical footprinting and NMR) of the interactions between minor groove binding drugs (e.g., netropsin and distamycin A) and DNA. The stabilizing forces between drug and DNA are provided by a combination of ionic, van der Waals, and hydrogen-bonding interactions. No bifurcated hydrogen bond is found between netropsin and DNA in this complex due to the unique dispositions of the hydrogen-bond acceptors (N3 of adenine and O2 of thymine) on the floor of the DNA minor groove. Two of the four AT base pairs in the ATAT stretch have low propeller twist angles, even though the DNA has a narrow minor groove. Alternating helical twist angles are observed in the ATAT stretch with lower twist in the ApT steps than in the TpA step.     

Recognition of bulged DNA by a neocarzinostatin product via an induced fit mechanism

The binding of the wedge-shaped isostructural analogue of the biradical species of the chromphore of antitumor antibiotic neocarzinostatin to sequence-specific bulged DNAs results in alterations in ellipticity of the DNAs. Circular dichroism (CD) spectroscopic results suggest that the drug specifically recognizes bulges of DNA via a combination of conformational selection and induced fit, not by binding to a preorganized site. Analysis of circular dichroism spectra indicates that the degree of induced fit observed is primarily a consequence of optimising van der Waals contacts with the walls of the bulge cavity. The effective recognition of the bulge site on duplex DNA appears to depend to a significant extent on the bent groove space being flexible enough to be able to adopt the geometrically optimal conformation compatible with the wedge-shaped drug molecule, rather than involving 'lock and key' recognition. The spectroscopic results indicate a change of DNA conformation, consistent with an allosteric binding model. Spectroscopic studies with various bulged DNAs also reveal that the binding strength directly correlates with the stability of the bulge structures.     

Sequence-dependent variation in DNA minor groove width dictates orientational preference of Hoechst 33258 in A-tract recognition: solution NMR structure of the 2:1 complex with d(CTTTTGCAAAAG)(2)

The solution structure of the dodecamer duplex d(CTTTTGCAAAAG)₂ and its 2:1 complex with the bis-benzimidazole Hoechst 33258 has been investigated by NMR and NOE-restrained molecular dynamics (rMD) simulations. Drug molecules are bound in each of the two A-tracts with the bulky N-methylpiperazine ring of each drug located close to the central TG (CA) step, binding essentially to the narrow minor groove of each A-tract. MD simulations over 1 ns, using an explicit solvation model, reveal time-averaged sequence-dependent narrowing of the minor groove from the 3′-end towards the 5′-end of each TTTT sequence. Distinct junctions at the TpG (CpA) steps, characterised by large positive roll, low helical and propeller twists and rapid AT base pair opening rates, add to the widening of the groove at these sites and appear to account for the bound orientation of the two drug molecules with the N-methylpiperazine ring binding in the wider part of the groove close to the junctions. Comparisons between the free DNA structure and the 2:1 complex (heavy atom RMSD 1.55 Å) reveal that these sequence-dependent features persist in both structures. NMR studies of the sequence d(GAAAAGCTTTTC)₂, in which the A-tracts have been inverted with the elimination of the TpG junctions, results in loss of orientational specificity of Hoechst 33258 and formation of multiple bound species in solution, consistent with the drug binding in a number of different orientations.     

Binding of hairpin pyrrole and imidazole polyamides to DNA: relationship between torsion angle and association rate constants

N-methylpyrrole (Py)-N-methylimidazole (Im) polyamides are small organic molecules that bind to DNA with sequence specificity and can be used as synthetic DNA-binding ligands. In this study, five hairpin eight-ring Py–Im polyamides 1–5 with different number of Im rings were synthesized, and their binding behaviour was investigated with surface plasmon resonance assay. It was found that association rate (k_a) of the Py–Im polyamides with their target DNA decreased with the number of Im in the Py–Im polyamides. The structures of four-ring Py–Im polyamides derived from density functional theory revealed that the dihedral angle of the Py amide carbonyl is 14∼18°, whereas that of the Im is significantly smaller. As the minor groove of DNA has a helical structure, planar Py–Im polyamides need to change their conformation to fit it upon binding to the minor groove. The data explain that an increase in planarity of Py–Im polyamide induced by the incorporation of Im reduces the association rate of Py–Im polyamides. This fundamental knowledge of the binding of Py–Im polyamides to DNA will facilitate the design of hairpin Py–Im polyamides as synthetic DNA-binding modules.     

Out-of-shape DNA minor groove binders: induced fit interactions of heterocyclic dications with the DNA minor groove

DB921 and DB911 are benzimidazole-biphenyl isomers with terminal charged amidines. DB911 has a central meta-substituted phenyl that gives it a shape similar to those of known minor groove binding compounds. DB921 has a central para-substituted phenyl with a linear conformation that lacks the appropriate radius of curvature to match the groove shape. It is thus expected that DB911, but not DB921, should be an effective minor groove binder, but we find that DB921 not only binds in the groove but also has an unusually high binding constant in SPR experiments (2.9 x 10(8) M(-)(1), vs 2.1 x 10(7) M(-)(1) for DB911). ITC thermodynamic analysis with an AATT sequence shows that the stronger binding of DB921 is due to a more favorable binding enthalpy relative to that of DB911. CD results support minor groove binding for both compounds but do not provide an explanation for the binding of DB921. X-ray crystallographic analysis of DB921 bound to AATT shows that an induced fit structural change in DB921 reduces the twist of the biphenyl to complement the groove, and places the functional groups in position to interact with bases at the floor of the groove. The phenylamidine of DB921 forms indirect contacts with the bases through a bound water. The DB921-water pair forms a curved binding module that matches the shape of the minor groove and provides a number of strong interactions that are not possible with DB911. This result suggests that traditional views of compound curvature required for minor groove complex formation should be reevaluated.     

Groove Width and Depth of B-DNA Structures Depend on Local Variation in Slide

Abstract

The groove widths of DNA helix, especially minor groove width, are generally believed to be important for recognition of DNA by various types of ligands. It has been postulated earlier that large negative propeller twist, in the AT rich regions compresses the minor groove of duplex DNA A systematic study has now been carried out by generating models with different values of local doublet and intra-basepair parameters and calculating their minor groove widths. It is found that several local doublet parameters affect the minor groove width but it depends most strongly on the local step parameters roll and slide when each parameter is considered individually. However, a detailed analysis of the various local parameters within the B-DNA family of crystal structures indicates that propeller twist and slide are most strongly correlated with the observed values of minor groove width. The groove depth is also strongly correlated with slide. Thus the local base sequence dependent variations in slide can modify both the groove width and depth and consequently determine the ligand binding properties of DNA.     

Groove width and depth of B-DNA structures depend on local variation in slide.

The groove widths of DNA helix, especially minor groove width, are generally believed to be important for recognition of DNA by various types of ligands. It has been postulated earlier that large negative propeller twist, in the AT rich regions compresses the minor groove of duplex DNA. A systematic study has now been carried out by generating models with different values of local doublet and intra-basepair parameters and calculating their minor groove widths. It is found that several local doublet parameters affect the minor groove width but it depends most strongly on the local step parameters roll and slide when each parameter is considered individually. However, a detailed analysis of the various local parameters within the B-DNA family of crystal structures indicates that propeller twist and slide are most strongly correlated with the observed values of minor groove width. The groove depth is also strongly correlated with slide. Thus the local base sequence dependent variations in slide can modify both the groove width and depth and consequently determine the ligand binding properties of DNA.     

Variations in duplex DNA conformation detected by the binding of monoclonal autoimmune antibodies.

Four monoclonal antibodies (Jel 229, 239, 241, 242) which bound to duplex DNA were prepared from two autoimmune female NZB/NZW mice. Their binding to various nucleic acids was investigated by a competitive solid phase radioimmune assay which allows the estimation of relative binding constants. None of the antibodies showed any consistent variation of binding constant with base composition and thus they must recognize features of the DNA backbone. Jel 241 binds across the major groove but the interaction with poly(pyrimidine) X poly(purine) DNAs was barely detectable. This antibody appears to recognize the "alternating-B" conformation which is promoted by methylation of pyrimidines in alternating sequences. The other three antibodies bind in the minor groove. In particular, for Jel 229 the preferred antigen was poly(dG) X poly(dC) with only weak binding to poly(dA) X poly(dT). This suggests a requirement for a wide minor groove. Thus autoimmune antibodies provide examples of "analogue" recognition and can be used to detect structural variations in the grooves of duplex DNA.     

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.     

New insights into the structure of An tracts and B'-B' bends in DNA

Energy calculations suggest that the currently available NOE distance constraints for An tracts in DNA are incapable of distinguishing between structures with a narrowed minor groove arising from a large propeller twist with a small inclination or from a small propeller twist with a large negative inclination. Furthermore, analysis of published data, together with energy estimations, strongly argue against bifurcated hydrogen bonding between A and T residues being the cause of the anomalous structural properties of An tracts. A conformational analysis of the B'-B' junction has been performed in which a single variable base pair has been inserted between two regions of B' structure. We have calculated low-energy structures for AnGAn,AnCAn,AnTAn,AnCTn, and TnCAn duplexes, where the An and Tn tracts were fixed in the anomalous B' conformation. Upon optimization, all these structures were found to contain a pronounced roll-like bending into the major groove at the site of the insertion. The important factors in the formation of these B'-B' bends are the destruction of the B' conformation and the concomitant widening of the minor groove at the junction region in order to reduce minor groove interstrand base clashes and improve interstrand stacking energy. If the B' conformation has strong negative inclination, the improved intrastrand stacking energy also contributes to the bending. In calculations of duplexes with An and Tn tracts in the B conformation instead of B', the bending disappears.     

A molecular dynamics study of the effect of G.T mispairs on the conformation of DNA in solution

The effect of G.T mispair incorporation into a double-helical environment was examined by molecular dynamics simulation. The 60-ps simulations performed on the two hexanucleotide duplexes d (G3C3)2 and d(G3TC2)2 included 10 Na+ counterions and first hydration shell waters. The resulting backbone torsional angle trajectories were analyzed to select time spans representative of conformational domains. The average backbone angles and helical parameters of the last time span for both duplexes are reported. During the simulation the hexamers retained B-type DNA structures that differed from typical A- or B-DNA forms. The overall helical structures for the two duplexes are vary similar. The presence of G.T mispairs did not alter the overall helical structure of the oligonucleotide duplex. Large propeller twist and buckle angles were obtained for both duplexes. The purine/pyrimidine crossover step showed a large decrease in propeller twist in the normal duplex but not in the mismatch duplex. Upon the formation of wobble mispairs in the mismatched duplex, the guanines moved into the minor groove and the thymines moved into the major groove. This helped prevent purine/purine clash and created a deformation in the relative orientation of the glycosidic bonds. It also exposed the free O4 of the thymines in the major groove and N2 of the guanines in the minor groove to interactions with solvent and counterions. These factors seemed to contribute to the apparently higher rigidity of the mismatched duplex during the simulation.     

Inhibition of DNA binding by human estrogen-related receptor 2 and estrogen receptor alpha with minor groove binding polyamides

Human estrogen-related receptor 2 (hERR2, ESRRB, ERRbeta, NR3B2) belongs to a class of nuclear receptors that bind DNA through sequence-specific interactions with a 5'-AGGTCA-3' estrogen response element (ERE) half-site in the major groove and an upstream 5'-TNA-3' site in the minor groove. This minor groove interaction is mediated by a C-terminal extension (CTE) of the DNA binding domain and is unique to the estrogen-related receptors. We have used synthetic pyrrole-imidazole polyamides, which bind specific sequences in the minor groove, to demonstrate that DNA binding by hERR2 is sensitive to the presence of polyamides in both the upstream minor groove CTE site and the minor groove of the ERE half-site. Thus, polyamides can inhibit hERR2 by two mechanisms, by direct steric blockage of minor groove DNA contacts mediated by the CTE and by changing the helical geometry of DNA such that major groove interactions are weakened. To confirm the generality of the latter approach, we show that the dimeric human estrogen receptor alpha (hERalpha, ESR1, NR3A1), which binds in the major groove of the ERE, can be inhibited by a polyamide bound in the opposing minor groove of the ERE. These results highlight two mechanisms for inhibition of protein-DNA interactions and extend the repertoire of DNA recognition motifs that can be inhibited by polyamides. These molecules may thus be useful for controlling expression of hERR2- or hERalpha-responsive genes.     

Conformational effects of DNA stretching

DNA is an extensible molecule, and an extended conformation of DNA is involved in some biological processes. We have examined the effect of elongation stress on the conformational properties of DNA base pairs by conformational analysis. The calculations show that stretching does significantly affect the conformational properties and flexibilities of base pairs. In particular, we have found that the propeller twist in base pairs reverses its sign upon stretching. The energy profile analysis indicates that electrostatic interactions make a major contribution to the stabilization of the positive-propeller-twist configuration in stretched DNA. This stretching also results in a monotonic decrease in the helical twist angle, tending to unwind the double helix. Fluctuations in most variables initially increase upon stretching, because of unstacking of base pairs, but then the fluctuations decrease as DNA is stretched further, owing to the formation of specific interactions between base pairs induced by the positive propeller twist. Thus, the stretching of DNA has particularly significant effects upon DNA flexibility. These changes in both the conformation and flexibility of base pairs probably have a role in functional interactions with proteins.     

Conformational Effects of DNA Stretching

DNA is an extensible molecule, and an extended conformation of DNA is involved in some biological processes. We have examined the effect of elongation stress on the conformational properties of DNA base pairs by conformational analysis. The calculations show that stretching does significantly affect the conformational properties and flexibilities of base pairs. In particular, we have found that the propeller twist in base pairs reverses its sign upon stretching. The energy profile analysis indicates that electrostatic interactions make a major contribution to the stabilization of the positive-propeller-twist configuration in stretched DNA. This stretching also results in a monotonic decrease in the helical twist angle, tending to unwind the double helix. Fluctuations in most variables initially increase upon stretching, because of unstacking of base pairs, but then the fluctuations decrease as DNA is stretched further, owing to the formation of specific interactions between base pairs induced by the positive propeller twist. Thus, the stretching of DNA has particularly significant effects upon DNA flexibility. These changes in both the conformation and flexibility of base pairs probably have a role in functional interactions with proteins.     

Unusually strong binding to the DNA minor groove by a highly twisted benzimidazole diphenylether: induced fit and bound water

RT29 is a dicationic diamidine derivative that does not obey the classical "rules" for shape and functional group placement that are expected to result in strong binding and specific recognition of the DNA minor groove. The compound contains a benzimidazole diphenyl ether core that is flanked by the amidine cations. The diphenyl ether is highly twisted and gives the entire compound too much curvature to fit well to the shape of the minor groove. DNase I footprinting, fluorescence intercalator displacement studies, and circular dichroism spectra, however, indicate that the compound is an AT specific minor groove binding agent. Even more surprisingly, quantitative biosensor-surface plasmon resonance and isothermal titration calorimetric results indicate that the compound binds with exceptional strength to certain AT sequences in DNA with a large negative enthalpy of binding. Crystallographic results for the DNA complex of RT29 compared to calculated results for the free compound show that the compound undergoes significant conformational changes to enhance its minor groove interactions. In addition, a water molecule is incorporated directly into the complex to complete the compound-DNA interface, and it forms an essential link between the compound and base pair edges at the floor of the minor groove. The calculated DeltaCp value for complex formation is substantially less than the experimentally observed value, which supports the idea of water being an intrinsic part of the complex with a major contribution to the DeltaCp value. Both the induced fit conformational changes of the compound and the bound water are essential for strong binding to DNA by RT29.     

DNA nicking by HinP1I endonuclease: bending, base flipping and minor groove expansion

HinP1I recognizes and cleaves the palindromic tetranucleotide sequence G↓CGC in DNA. We report three structures of HinP1I–DNA complexes: in the presence of Ca²⁺ (pre-reactive complex), in the absence of metal ion (binary complex) and in the presence of Mg²⁺ (post-reactive complex). HinP1I forms a back-to-back dimer with two active sites and two DNA duplexes bound on the outer surfaces of the dimer facing away from each other. The 10 bp DNA duplexes undergo protein-induced distortions exhibiting features of A-, B- and Z-conformations: bending on one side (by intercalation of a phenylalanine side chain into the major groove), base flipping on the other side of the recognition site (by expanding the step rise distance of the local base pair to Z-form) and a local A-form conformation between the two central C:G base pairs of the recognition site (by binding of the N-terminal helix in the minor groove). In the pre- and post-reactive complexes, two metals (Ca²⁺ or Mg²⁺) are found in the active site. The enzyme appears to cleave DNA sequentially, hydrolyzing first one DNA strand, as seen in the post-reactive complex in the crystalline state, and then the other, as supported by the observation that, in solution, a nicked DNA intermediate accumulates before linearization.     

Critical effect of the N2 amino group on structure, dynamics, and elasticity of DNA polypurine tracts

Unrestrained 5-20-ns explicit-solvent molecular dynamics simulations using the Cornell et al. force field have been carried out for d[GCG(N)11GCG]2 (N, purine base) considering guanine*cytosine (G*C), adenine*thymine (A*T), inosine*5-methyl-cytosine (I*mC), and 2-amino-adenine*thymine (D*T) basepairs. The simulations unambiguously show that the structure and elasticity of N-tracts is primarily determined by the presence of the amino group in the minor groove. Simulated A-, I-, and AI-tracts show almost identical structures, with high propeller twist and minor groove narrowing. G- and D-tracts have small propeller twisting and are partly shifted toward the A-form. The elastic properties also differ between the two groups. The sequence-dependent electrostatic component of base stacking seems to play a minor role. Our conclusions are entirely consistent with available experimental data. Nevertheless, the propeller twist and helical twist in the simulated A-tract appear to be underestimated compared to crystallographic studies. To obtain further insight into the possible force field deficiencies, additional multiple simulations have been made for d(A)10, systematically comparing four major force fields currently used in DNA simulations and utilizing B and A-DNA forms as the starting structure. This comparison shows that the conclusions of the present work are not influenced by the force field choice.     

A unified model for the origin of DNA sequence-directed curvature

The fine structure of the DNA double helix and a number of its physical properties depend upon nucleotide sequence. This includes minor groove width, the propensity to undergo the B-form to A-form transition, sequence-directed curvature, and cation localization. Despite the multitude of studies conducted on DNA, it is still difficult to appreciate how these fundamental properties are linked to each other at the level of nucleotide sequence. We demonstrate that several sequence-dependent properties of DNA can be attributed, at least in part, to the sequence-specific localization of cations in the major and minor grooves. We also show that effects of cation localization on DNA structure are easier to understand if we divide all DNA sequences into three principal groups: A-tracts, G-tracts, and generic DNA. The A-tract group of sequences has a peculiar helical structure (i.e., B*-form) with an unusually narrow minor groove and high base-pair propeller twist. Both experimental and theoretical studies have provided evidence that the B*-form helical structure of A-tracts requires cations to be localized in the minor groove. G-tracts, on the other hand, have a propensity to undergo the B-form to A-form transition with increasing ionic strength. This property of G-tracts is directly connected to the observation that cations are preferentially localized in the major groove of G-tract sequences. Generic DNA, which represents the vast majority of DNA sequences, has a more balanced occupation of the major and minor grooves by cations than A-tracts or G-tracts and is thereby stabilized in the canonical B-form helix. Thus, DNA secondary structure can be viewed as a tug of war between the major and minor grooves for cations, with A-tracts and G-tracts each having one groove that dominates the other for cation localization. Finally, the sequence-directed curvature caused by A-tracts and G-tracts can, in both cases, be explained by the cation-dependent mismatch of A-tract and G-tract helical structures with the canonical B-form helix of generic DNA (i.e., a cation-dependent junction model).     

Sequence-specific recognition of DNA minor groove by an NIR-fluorescence switch-on probe and its potential applications

In molecular biology, understanding the functional and structural aspects of DNA requires sequence-specific DNA binding probes. Especially, sequence-specific fluorescence probes offer the advantage of real-time monitoring of the conformational and structural reorganization of DNA in living cells. Herein, we designed a new class of D2A (one-donor-two-acceptor) near-infrared (NIR) fluorescence switch-on probe named quinone cyanine–dithiazole (QCy–DT) based on the distinctive internal charge transfer (ICT) process for minor groove recognition of AT-rich DNA. Interestingly, QCy–DT exhibited strong NIR-fluorescence enhancement in the presence of AT-rich DNA compared to GC-rich and single-stranded DNAs. We show sequence-specific minor groove recognition of QCy–DT for DNA containing 5′-AATT-3′ sequence over other variable (A/T)4 sequences and local nucleobase variation study around the 5′-X(AATT)Y-3′ recognition sequence revealed that X = A and Y = T are the most preferable nucleobases. The live cell imaging studies confirmed mammalian cell permeability, low-toxicity and selective staining capacity of nuclear DNA without requiring RNase treatment. Further, Plasmodium falciparum with an AT-rich genome showed specific uptake with a reasonably low IC₅₀ value (<4 µM). The ease of synthesis, large Stokes shift, sequence-specific DNA minor groove recognition with switch-on NIR-fluorescence, photostability and parasite staining with low IC₅₀ make QCy–DT a potential and commercially viable DNA probe.