Intermolecular interactions and water structure in a condensed phase B-DNA crystal

By controlled dehydration, the unit cells of dodecamer DNA–drug crystals have been shrunk from 68 000 (normal state) to 60 000 (partially dehydrated intermediate state) to 51 000 ų (fully dehydrated state), beyond which no further solvent loss occurs. The total solvent content in the normal crystals is ~40% by volume, reducing to ~20% in the fully dehydrated phase. The 25% reduction in cell volume induced a dramatic enhancement in the resolution of the X-ray diffraction data (from 2.6 to beyond 1.5 Å). We have determined the structures of the normal, partially dehydrated and fully dehydrated crystals. Details of the ligand binding have been presented in the preceding article. The present paper describes the unique features of the structure of the fully dehydrated phase. This structure was refined with 9015 unique observed reflections to R = 14.9%, making it one of the most reliable models of B-form DNA available. The crystals exist as infinite polymeric networks, in which neighbouring dodecamer duplexes are crosslinked through phosphate oxygens via direct bonding to magnesium cations. The DNA is packed so tightly that there is essentially only a single layer of solvent between adjacent molecules. The details of the crystal packing, magnesium bridging, DNA hydration and DNA conformation are described and compared with other experimental evidence related to DNA condensation.     

The 1.19 A X-ray structure of 2'-O-Me(CGCGCG)(2) duplex shows dehydrated RNA with 2-methyl-2,4-pentanediol in the minor groove

The crystal and molecular structure of 2′-O-Me(CGCGCG)₂ has been determined at 1.19 Å resolution, at 100 K, using synchrotron radiation. The structure in space group P3₂12 is a half-turn right-handed helix that includes two 2-methyl-2,4-pentanediol (MPD) molecules bound in the minor groove. The structure deviates from A-form RNA. The duplex is overwound with an average value of 9.7 bp per turn, characterised as having a C3′-endo sugar pucker, very low base pair rise and high helical twist and inclination angles. The structure includes 65 ordered water molecules. Only a single row of water molecules is observed in the minor groove due to the presence of hydrophobic 2′-O-methyl groups. As many as five magnesium ions are located in the structure. Two are in the major groove and interact with O⁶ and N⁷ of guanosine and N⁴ of cytidine residues through their hydration spheres. This work provides the first example of molecular interactions of nucleic acids with MPD, which was used as a precipitant, cryo-solvent and resolution enhancing agent. The two MPD molecules intrude into the hydration network in the minor groove, each forming hydrogen bonds between their secondary hydroxyl group and exo-amino functions of guanosine residues. Comparison of the 2′-O-Me(CGCGCG)₂ structure in the P3₂12 and P6₁22 crystals delineates stability of the water network within the minor groove to dehydration by MPD and is of interest for evaluating factors governing small molecule binding to RNA. Intrusion of MPD into the minor groove of 2′-O-Me(CGCGCG)₂ is discussed with respect to RNA dehydration, a prerequisite of Z-RNA formation.     

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.     

Potent inhibition of werner and bloom helicases by DNA minor groove binding drugs

Maintenance of genomic integrity is vital to all organisms. A number of human genetic disorders, including Werner Syndrome, Bloom Syndrome and Rothmund–Thomson Syndrome, exhibit genomic instability with some phenotypic characteristics of premature aging and cancer predisposition. Presumably the aberrant cellular and clinical phenotypes in these disorders arise from defects in important DNA metabolic pathways such as replication, recombination or repair. These syndromes are all characterized by defects in a member of the RecQ family of DNA helicases. To obtain a better understanding of how these enzymes function in DNA metabolic pathways that directly influence chromosomal integrity, we have examined the effects of non-covalent DNA modifications on the catalytic activities of purified Werner (WRN) and Bloom (BLM) DNA helicases. A panel of DNA-binding ligands displaying unique properties for interacting with double helical DNA was tested for their effects on the unwinding activity of WRN and BLM helicases on a partial duplex DNA substrate. The levels of inhibition by a number of these compounds were distinct from previously reported values for viral, prokaryotic and eukaryotic helicases. The results demonstrate that BLM and WRN proteins exhibit similar sensitivity profiles to these DNA-binding ligands and are most potently inhibited by the structurally related minor groove binders distamycin A and netropsin (K_i ≤1 µM). The distinct inhibition of WRN and BLM helicases by the minor groove binders suggest that these helicases unwind double-stranded DNA by a related mechanism.     

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.     

Solution structure of an arabinonucleic acid (ANA)/RNA duplex in a chimeric hairpin: comparison with 2′-fluoro-ANA/RNA and DNA/RNA hybrids

Hybrids of RNA and arabinonucleic acid (ANA) as well as the 2′-fluoro-ANA analog (2′F-ANA) were recently shown to be substrates of the enzyme RNase H. Although RNase H binds to double-stranded RNA, no cleavage occurs with such duplexes. Therefore, knowledge of the structure of ANA/RNA hybrids may prove helpful in the design of future antisense oligonucleotide analogs. In this study, we have determined the NMR solution structures of ANA/RNA and DNA/RNA hairpin duplexes and compared them to the recently published structure of a 2′F-ANA/RNA hairpin duplex. We demonstrate here that the sugars of RNA nucleotides of the ANA/RNA hairpin stem adopt the C3′-endo (north, A-form) conformation, whereas those of the ANA strand adopt a ‘rigid’ O4′-endo (east) sugar pucker. The DNA strand of the DNA/RNA hairpin stem is flexible, but the average DNA/RNA hairpin structural parameters are close to the ANA/RNA and 2′F-ANA/RNA hairpin parameters. The minor groove width of ANA/RNA, 2′F-ANA/RNA and DNA/RNA helices is 9.0 ± 0.5 Å, a value that is intermediate between that of A- and B-form duplexes. These results rationalize the ability of ANA/RNA and 2′F-ANA/RNA hybrids to elicit RNase H activity.     

Molecular dynamics simulations and free energy calculations of netropsin and distamycin binding to an AAAAA DNA binding site

Molecular dynamics simulations have been performed on netropsin in two different charge states and on distamycin binding to the minor groove of the DNA duplex d(CGCGAAAAACGCG)·d(CGCGTTTTTCGCG). The relative free energy of binding of the two non-covalently interacting ligands was calculated using the thermodynamic integration method and reflects the experimental result. From 2 ns simulations of the ligands free in solution and when bound to DNA, the mobility and the hydrogen-bonding patterns of the ligands were studied, as well as their hydration. It is shown that even though distamycin is less hydrated than netropsin, the loss of ligand–solvent interactions is very similar for both ligands. The relative mobilities of the ligands in their bound and free forms indicate a larger entropic penalty for distamycin when binding to the minor groove compared with netropsin, partially explaining the lower binding affinity of the distamycin molecule. The detailed structural and energetic insights obtained from the molecular dynamics simulations allow for a better understanding of the factors determining ligand–DNA binding.     

DNA bending induced by DNA (cytosine-5) methyltransferases

DNA bending induced by six DNA (cytosine-5) methyltransferases was studied using circular permutation gel mobility shift assay. The following bend angles were obtained: M.BspRI (GG^m5CC), 46–50°; M.HaeIII (GG^m5CC), 40–43°; M.SinI (GGW^m5CC), 34–37°; M.Sau96I (GGN^m5CC), 52–57°; M.HpaII (C^m5CGG), 30°; and M.HhaI (G^m5CGC), 13°. M.HaeIII was also tested with fragments carrying a methylated binding site, and it was found to induce a 32° bend. A phase-sensitive gel mobility shift assay, using a set of DNA fragments with a sequence-directed bend and a single methyltransferase binding site, indicated that M.HaeIII and M.BspRI bend DNA toward the minor groove. The DNA curvature induced by M.HaeIII contrasts with the lack of DNA bend observed for a covalent M.HaeIII–DNA complex in an earlier X-ray study. Our results and data from other laboratories show a correlation between the bending properties and the recognition specificities of (cytosine-5) methyltransferases: enzymes recognizing a cytosine 3′ to the target cytosine tend to induce greater bends than enzymes with guanine in this position. We suggest that the observed differences indicate different mechanisms employed by (cytosine-5) methyltransferases to stabilize the helix after the target base has flipped out.     

Similarities and differences in interaction of K+ and Na+ with condensed ordered DNA. A molecular dynamics computer simulation study

Four 20 ns molecular dynamics simulations have been performed with two counterions, K⁺ or Na⁺, at two water contents, 15 or 20 H₂O per nucleotide. A hexagonal simulation cell comprised of three identical DNA decamers [d(5′-ATGCAGTCAG) × d(5′-TGACTGCATC)] with periodic boundary condition along the DNA helix was used. The simulation setup mimics the DNA state in oriented DNA fibers or in crystals of DNA oligomers. Variation of counterion nature and water content do not alter averaged DNA structure. K⁺ and Na⁺ binding to DNA are different. K⁺ binds to the electronegative sites of DNA bases in the major and the minor grooves, while Na⁺ interacts preferentially with the phosphate groups. Increase of water causes a shift of both K⁺ and Na⁺ from the first hydration shell of O1P/O2P and of the DNA bases in the minor groove with lesser influence for the cation binding to the bases in the major groove. Mobility of both water and cations in the K–DNA systems is faster than in the Na–DNA systems: Na⁺ organizes and immobilizes water structure around itself and near DNA while for K⁺ water is less organized and more dynamic.     

Determination of DNA minor groove width in distamycin-DNA complexes by solid-state NMR

We have performed solid-state ³¹P-¹⁹F REDOR nuclear magnetic resonance (NMR) experiments to monitor changes in minor groove width of the oligonucleotide d(CGCAAA^2′FUTGGC)·d(GCCAAT(pS)TT GCG) (A₃T₂) upon binding of the drug distamycin A at different stoichiometries. In the hydrated solid-state sample, the minor groove width for the unbound DNA, measured as the ^2′FU7–pS19 inter-label distance, was 9.4 ± 0.7 Å, comparable to that found for similar A:T-rich DNAs. Binding of a single drug molecule is observed to cause a 2.4 Å decrease in groove width. Subsequent addition of a second drug molecule results in a larger conformational change, expanding this minor groove width to 13.6 Å, consistent with the results of a previous solution NMR study of the 2:1 complex. These ³¹P-¹⁹F REDOR results demonstrate the ability of solid-state NMR to measure distances of 7–14 Å in DNA–drug complexes and provide the first example of a direct spectroscopic measurement of minor groove width in nucleic acids.     

Structural Ordering of Disordered Ligand-Binding Loops of Biotin Protein Ligase into Active Conformations as a Consequence of Dehydration

Mycobacterium tuberculosis (Mtb), a dreaded pathogen, has a unique cell envelope composed of high fatty acid content that plays a crucial role in its pathogenesis. Acetyl Coenzyme A Carboxylase (ACC), an important enzyme that catalyzes the first reaction of fatty acid biosynthesis, is biotinylated by biotin acetyl-CoA carboxylase ligase (BirA). The ligand-binding loops in all known apo BirAs to date are disordered and attain an ordered structure only after undergoing a conformational change upon ligand-binding. Here, we report that dehydration of Mtb-BirA crystals traps both the apo and active conformations in its asymmetric unit, and for the first time provides structural evidence of such transformation. Recombinant Mtb-BirA was crystallized at room temperature, and diffraction data was collected at 295 K as well as at 120 K. Transfer of crystals to paraffin and paratone-N oil (cryoprotectants) prior to flash-freezing induced lattice shrinkage and enhancement in the resolution of the X-ray diffraction data. Intriguingly, the crystal lattice rearrangement due to shrinkage in the dehydrated Mtb-BirA crystals ensued structural order of otherwise flexible ligand-binding loops L4 and L8 in apo BirA. In addition, crystal dehydration resulted in a shift of ∼3.5 Å in the flexible loop L6, a proline-rich loop unique to Mtb complex as well as around the L11 region. The shift in loop L11 in the C-terminal domain on dehydration emulates the action responsible for the complex formation with its protein ligand biotin carboxyl carrier protein (BCCP) domain of ACCA3. This is contrary to the involvement of loop L14 observed in Pyrococcus horikoshii BirA-BCCP complex. Another interesting feature that emerges from this dehydrated structure is that the two subunits A and B, though related by a noncrystallographic twofold symmetry, assemble into an asymmetric dimer representing the ligand-bound and ligand-free states of the protein, respectively. In-depth analyses of the sequence and the structure also provide answers to the reported lower affinities of Mtb-BirA toward ATP and biotin substrates. This dehydrated crystal structure not only provides key leads to the understanding of the structure/function relationships in the protein in the absence of any ligand-bound structure, but also demonstrates the merit of dehydration of crystals as an inimitable technique to have a glance at proteins in action.     

A novel form of intercalation involving four DNA duplexes in an acridine-4-carboxamide complex of d(CGTACG)(2)

The structures of the complexes formed between 9-amino-[N-(2-dimethyl-amino)butyl]acridine-4-carboxamide and d(CG^5BrUACG)₂ and d(CGTACG)₂ have been solved by X-ray crystallography using MAD phasing methodology and refined to a resolution of 1.6 Å. The complexes crystallised in space group C222. An asymmetric unit in the brominated complex comprises two strands of DNA, one disordered drug molecule, two cobalt (II) ions and 19 water molecules (31 in the native complex). Asymmetric units in the native complex also contain a sodium ion. The structures exhibit novel features not previously observed in crystals of DNA/drug complexes. The DNA helices stack in continuous columns with their central 4 bp adopting a B-like motif. However, despite being a palindromic sequence, the terminal GC base pairs engage in quite different interactions. At one end of the duplex there is a CpG dinucleotide overlap modified by ligand intercalation and terminal cytosine exchange between symmetry-related duplexes. A novel intercalation complex is formed involving four DNA duplexes, four ligand molecules and two pairs of base tetrads. The other end of the DNA is frayed with the terminal guanine lying in the minor groove of the next duplex in the column. The structure is stabilised by guanine N7/cobalt (II) coordination. We discuss our findings with respect to the effects of packing forces on DNA crystal structure, and the potential effects of intercalating agents on biochemical processes involving DNA quadruplexes and strand exchanges. NDB accession numbers: DD0032 (brominated) and DD0033 (native).     

Unusual DNA duplex and hairpin motifs

Single-stranded DNA or double-stranded DNA has the potential to adopt a wide variety of unusual duplex and hairpin motifs in the presence (trans) or absence (cis) of ligands. Several principles for the formation of those unusual structures have been established through the observation of a number of recurring structural motifs associated with different sequences. These include: (i) internal loops of consecutive mismatches can occur in a B-DNA duplex when sheared base pairs are adjacent to each other to confer extensive cross- and intra-strand base stacking; (ii) interdigitated (zipper-like) duplex structures form instead when sheared G·A base pairs are separated by one or two pairs of purine·purine mismatches; (iii) stacking is not restricted to base, deoxyribose also exhibits the potential to do so; (iv) canonical G·C or A·T base pairs are flexible enough to exhibit considerable changes from the regular H-bonded conformation. The paired bases become stacked when bracketed by sheared G·A base pairs, or become extruded out and perpendicular to their neighboring bases in the presence of interacting drugs; (v) the purine-rich and pyrimidine-rich loop structures are notably different in nature. The purine-rich loops form compact triloop structures closed by a sheared G·A, A·A, A·C or sheared-like G_anti·C_syn base pair that is stacked by a single residue. On the other hand, the pyrimidine-rich loops with a thymidine in the first position exhibit no base pairing but are characterized by the folding of the thymidine residue into the minor groove to form a compact loop structure. Identification of such diverse duplex or hairpin motifs greatly enlarges the repertoire for unusual DNA structural formation.     

Kinetic Complexity, Homogeneity, and Copy Number of Chloroplast DNA from the Marine Alga Olisthodiscus luteus

The kinetic complexity of chloroplast DNA isolated from the chromophytic alga Olisthodiscus luteus has been determined. Using optical reassociation techniques, it was shown that the plastid DNA of this alga reacted as a single component with a second order rate constant of 4.1 molar⁻¹ and second⁻¹ (Cot_½ 0.24 molar second) under conditions equivalent to 180 millimolar Na⁺ and 60°C. Given the 92 × 10⁵ dalton complexity calculated for this chloroplast genome, an Olisthodiscus cell contains 650 plastome copies. Although this complement remains constant throughout the growth cycle of the organism, the ploidy level of an individual chloroplast shows significant plasticity and is dependent upon the number of chloroplasts present per cell. Experiments with the DNA fluorochrome Hoechst dye 33258 (bisbenzimide) demonstrate that plastids isolated from all phases of cell growth each possess a ring-shaped nucleoid containing detectable DNA. Olisthodiscus chloroplast DNA showed no sequence mismatch when thermal denaturation profiles of reassociated chloroplast DNA were examined, thus all plastome copies are essentially identical. Finally, reassociation studies demonstrated that no foldback (short inverted repeat) sequences were present in the plastid genome although significant hairpin loop structures were observed in control nuclear DNA samples.     

Intracellular Freezing in the Infective Juveniles of Steinernema feltiae: An Entomopathogenic Nematode

Taking advantage of their optical transparency, we clearly observed the third stage infective juveniles (IJs) of Steinernema feltiae freezing under a cryo-stage microscope. The IJs froze when the water surrounding them froze at −2°C and below. However, they avoid inoculative freezing at −1°C, suggesting cryoprotective dehydration. Freezing was evident as a sudden darkening and cessation of IJs'' movement. Freeze substitution and transmission electron microscopy confirmed that the IJs of S. feltiae freeze intracellularly. Ice crystals were found in every compartment of the body. IJs frozen at high sub-zero temperatures (−1 and −3°C) survived and had small ice crystals. Those frozen at −10°C had large ice crystals and did not survive. However, the pattern of ice formation was not well-controlled and individual nematodes frozen at −3°C had both small and large ice crystals. IJs frozen by plunging directly into liquid nitrogen had small ice crystals, but did not survive. This study thus presents the evidence that S. feltiae is only the second freeze tolerant animal, after the Antarctic nematode Panagrolaimus davidi, shown to withstand extensive intracellular freezing.     

Membrane Organization of the Desiccation-Tolerant Moss Tortula ruralis in Dehydrated States

Membrane organization of the desiccation tolerant moss Tortula ruralis was studied in several intensely dehydrated states (75% relative humidity [RH], 90% RH, plasmolysis in molar salt, freezing to −20°C) by ³¹P nuclear magnetic resonance and ultrastructural analyses. Both methods revealed that even at 75% RH (−400 bars), the moss cellular membranes retained extended phospholipid bilayers. Ultrastructural analyses of the fully hydrated moss showed an extensive proliferation of membrane vesicles in the endoplasmic reticulum. During dehydration, these vesicles form layers of membrane under the plasmalemma and in some cases appear to fuse with the surface membrane. This suggests that these vesicles may serve as a reservoir of membranes to accommodate for membrane surface area changes during desiccation and subsequent rehydration.     

Recognition of RNA duplex by a neomycin–Hoechst 33258 conjugate

DNA minor groove binding drugs such as Hoechst 33258 have been shown to bind to a number of RNA structures. Similarly, RNA binding ligands such as neomycin have been shown by us to bind to a number of A-form DNA structures. A neomycin–Hoechst 33258 conjugate was recently shown to bind B-DNA, where Hoechst exhibits high affinity for the minor groove of A/T tract DNA and neomycin docks into the major groove. Further studies now indicate that the Hoechst moiety of the conjugate can be driven to bind RNA duplex as a consequence of neomycin binding in the RNA major groove. This is the first example of Hoechst 33258 binding to RNA duplex not containing bulges or loop motifs.     

Synthesis,antiproliferative and antitrypanosomal activities,and DNA binding of novel 6-amidino-2-arylbenzothiazoles

A series of 6-amidinobenzothiazoles, linked via phenoxymethylene or directly to the 1,2,3-triazole ring with a p-substituted phenyl or benzyl moiety, were synthesised and evaluated in vitro against four human tumour cell lines and the protozoan parasite Trypanosoma brucei. The influence of the type of amidino substituent and phenoxymethylene linker on antiproliferative and antitrypanosomal activities was observed, showing that the imidazoline moiety had a major impact on both activities. Benzothiazole imidazoline 14a, which was directly connected to N-1-phenyl-1,2,3-triazole, had the most potent growth-inhibitory effect (IC₅₀ = 0.25 µM) on colorectal adenocarcinoma (SW620), while benzothiazole imidazoline 11b, containing a phenoxymethylene linker, exhibited the best antitrypanosomal potency (IC₉₀ = 0.12 µM). DNA binding assays showed a non-covalent interaction of 6-amidinobenzothiazole ligands, indicating both minor groove binding and intercalation modes of DNA interaction. Our findings encourage further development of novel structurally related 6-amidino-2-arylbenzothiazoles to obtain more selective anticancer and anti-HAT agents.     

Crystal structures of Escherichia coli exonuclease I in complex with single-stranded DNA provide insights into the mechanism of processive digestion

Escherichia coli Exonuclease I (ExoI) digests single-stranded DNA (ssDNA) in the 3′-5′ direction in a highly processive manner. The crystal structure of ExoI, determined previously in the absence of DNA, revealed a C-shaped molecule with three domains that form a central positively charged groove. The active site is at the bottom of the groove, while an extended loop, proposed to encircle the DNA, crosses over the groove. Here, we present crystal structures of ExoI in complex with four different ssDNA substrates. The structures all have the ssDNA bound in essentially the predicted manner, with the 3′-end in the active site and the downstream end under the crossover loop. The central nucleotides of the DNA form a prominent bulge that contacts the SH3-like domain, while the nucleotides at the downstream end of the DNA form extensive interactions with an ‘anchor’ site. Seven of the complexes are similar to one another, but one has the ssDNA bound in a distinct conformation. The highest-resolution structure, determined at 1.95 Å, reveals an Mg²⁺ ion bound to the scissile phosphate in a position corresponding to Mg^B in related two-metal nucleases. The structures provide new insights into the mechanism of processive digestion that will be discussed.