Sequence-selective metal ion binding to DNA oligonucleotides

Insights into DNA solvation found in protein-DNA structures

The proteins that bind double-helical DNA present various microenvironments that sense and/or induce signals in the genetic material. The high-resolution structures of protein-DNA complexes reveal the nature of both the microenvironments and the conformational responses in DNA and protein. Complex networks of interactions within the structures somehow tie the protein and DNA together and induce the observed spatial forms. Here we show how the cumulative buildup of amino acid atoms around the sugars, phosphates, and bases in different protein-DNA complexes produces a binding cloud around the double helix and how different types of atoms fill that cloud. Rather than focusing on the principles of molecular binding and recognition suggested by the arrangements of amino acids and nucleotides in the macromolecular complexes, we consider the proteins in contact with DNA as organized solvents. We describe differences in the mix of atoms that come in closest contact with DNA, subtle sequence-dependent features in the microenvironment of the sugar-phosphate backbone, a direct link between the localized buildup of ionic species and the electrostatic potential surfaces of the DNA bases, and sites of atomic buildup above and below the basepair planes that transmit the unique features of the base environments along the chain backbone. The inferences about solvation that can be drawn from the survey provide new stimuli for improvement of nucleic acid force fields and fresh ideas for exploration of the properties of DNA in solution.

Optimizing the Chemiluminescence of a Light-Producing Deoxyribozyme

Supernova is a chemiluminescent deoxyribozyme recently discovered in our group. It transfers the phosphate group from the 1,2-dioxetane substrate CDP-Star to its 5' hydroxyl group, which triggers a decomposition reaction and the production of light. Here we investigated the effects of reaction conditions on the ability of Supernova to generate a chemiluminescent signal (using a plate reader assay) and to phosphorylate itself (using a ligation assay). Our experiments indicate that multiple zinc ions are required for catalytic function, suggesting links between Supernova and protein enzymes that catalyze similar reactions. They also show how factors such as pH, potassium concentration, CDP-Star concentration, and DNA concentration affect the reaction. By combining information from different experiments, the rate enhancement of light production was increased by more than 1000-fold. These results should be useful for applications in which Supernova is used as a sensor.

Supernova: A Deoxyribozyme that Catalyzes a Chemiluminescent Reaction

Functional DNA molecules are useful components in nanotechnology and synthetic biology. To expand the toolkit of functional DNA parts, in this study we used artificial evolution to identify a glowing deoxyribozyme called Supernova. This deoxyribozyme transfers a phosphate from a 1,2-dioxetane substrate to its 5' hydroxyl group, which triggers a chemiluminescent reaction and a flash of blue light. An engineered version of Supernova is only catalytically active in the presence of an oligonucleotide complementary to its 3' end, demonstrating that light production can be coupled to ligand binding. We anticipate that Supernova will be useful in a wide variety of applications, including as a signaling component in allosterically regulated sensors and in logic gates of molecular computers.

The structural plasticity of nucleic acid duplexes revealed by WAXS and MD

Double-stranded DNA (dsDNA) and RNA (dsRNA) helices display an unusual structural diversity. Some structural variations are linked to sequence and may serve as signaling units for protein-binding partners. Therefore, elucidating the mechanisms and factors that modulate these variations is of fundamental importance. While the structural diversity of dsDNA has been extensively studied, similar studies have not been performed for dsRNA. Because of the increasing awareness of RNA's diverse biological roles, such studies are timely and increasingly important. We integrate solution x-ray scattering at wide angles (WAXS) with all-atom molecular dynamics simulations to explore the conformational ensemble of duplex topologies for different sequences and salt conditions. These tightly coordinated studies identify robust correlations between features in the WAXS profiles and duplex geometry and enable atomic-level insights into the structural diversity of DNA and RNA duplexes. Notably, dsRNA displays a marked sensitivity to the valence and identity of its associated cations.

Mechanisms of DNA damage induced by morin, an inhibitor of amyloid β-peptide aggregation

Morin is a potential inhibitor of amyloid β-peptide aggregation. This aggregation is involved in the pathogenesis of Alzheimer's disease. Meanwhile, morin has been found to be mutagenic and exhibits peroxidation of membrane lipids concurrent with DNA strand breaks in the presence of metal ions. To clarify a molecular mechanism of morin-induced DNA damage, we examined the DNA damage and its site specificity on ³²P-5'-end-labeled human DNA fragments treated with morin plus Cu(II). The formation of 8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxodG), an indicator of oxidative DNA damage, was also determined in calf thymus DNA treated with morin plus Cu(II). Morin-induced DNA strand breaks and base modification in the presence of Cu(II) were dose dependent. Morin plus Cu(II) caused piperidine-labile lesions preferentially at thymine and guanine residues. The DNA damage was inhibited by methional, catalase and Cu(I)-chelator bathocuproine. The typical •OH scavengers ethanol, mannitol and sodium formate showed no inhibitory effect on DNA damage induced by morin plus Cu(II). When superoxide dismutase was added to the solution, DNA damage was not inhibited. In addition, morin plus Cu(II) increased 8-oxodG formation in calf thymus DNA fragments. We conclude that morin undergoes autoxidation in the presence of Cu(II) via a Cu(I)/Cu(II) redox cycle and H₂O₂ generation to produce Cu(I)-hydroperoxide, which causes oxidative DNA damage.

Experimental and first-principles investigation of the adsorption and entrapping of guanine with SiO2 clusters of sol–gel silicate material for understanding DNA photodamage

A sol–gel silicate matrix containing entrapped guanine was prepared. The SiO₂ matrix provides UVA protection by reducing the light penetration to the entrapped guanine molecules.

Ion distributions around left- and right-handed DNA and RNA duplexes: a comparative study

The ion atmosphere around nucleic acids is an integral part of their solvated structure. However, detailed aspects of the ionic distribution are difficult to probe experimentally, and comparative studies for different structures of the same sequence are almost non-existent. Here, we have used large-scale molecular dynamics simulations to perform a comparative study of the ion distribution around (5'-CGCGCGCGCGCG-3')2 dodecamers in solution in B-DNA, A-RNA, Z-DNA and Z-RNA forms. The CG sequence is very sensitive to ionic strength and it allows the comparison with the rare but important left-handed forms. The ions investigated include Na(+), K(+) and Mg(2 +), with various concentrations of their chloride salts. Our results quantitatively describe the characteristics of the ionic distributions for different structures at varying ionic strengths, tracing these differences to nucleic acid structure and ion type. Several binding pockets with rather long ion residence times are described, both for the monovalent ions and for the hexahydrated Mg[(H2O)6](2+) ion. The conformations of these binding pockets include direct binding through desolvated ion bridges in the GpC steps in B-DNA and A-RNA; direct binding to backbone oxygens; binding of Mg[(H2O)6](2+) to distant phosphates, resulting in acute bending of A-RNA; tight 'ion traps' in Z-RNA between C-O2 and the C-O2' atoms in GpC steps; and others.

Mg2+ 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 the tetradecanucleotide d(CCCCGGTACCGGGG)2 as an A-DNA duplex

The crystal structure of the tetradecanucleotide sequence d(CCCCGGTACCGGGG)(2) has been determined at 2.5 Å resolution in the tetragonal space group P4(1). This sequence was designed with the expectation of a four-way junction. However, the sequence crystallized as an A-DNA duplex and represents more than one full turn of the A-helix. The crystallographic asymmetric unit consists of one tetradecanucleotide duplex. The structural parameters of the A-type DNA duplex structure and the crystal-packing arrangement are described. One Mn(2+) ion was identified with direct coordination to the N7 position of G(13) and a water molecule at the major-groove side of the C(2)·G(13) base pair.

A DNA-based assay for toxic chemicals in wastewater

Chemical toxicants, particularly metal ions, are a major contaminant in global waterways. Live-organism bioassays used to monitor chemical toxicants commonly involve measurements of activity or survival of a freshwater cladoceran (Ceriodaphnia dubia) or light emitted by the marine bacterium Vibrio fischeri, used in the commercial Microtox® bioassay. Here we describe a novel molecule-based assay system employing DNA as the chemical biosensor. Metals bind to DNA, causing structural changes that expel a bound (intercalated) fluorescent reporter dye. Analyses of test data using 48 wastewater samples potentially contaminated by metal ions show that the DNA-dye assay results correlate with those from C. dubia and Microtox bioassays. All three assays exhibit additive, antagonistic, and synergistic responses that cannot be predicted by knowing individual metal concentrations. Analyses of metals in these samples imply the presence of chemical toxicants other than metal ions. The DNA-dye assay is robust, has a 12-month shelf life, and is only slightly affected by sample pH in the range 4 to 9. The assay is completed in a matter of minutes, and its portability makes it well suited as a screening assay for use in the field. We conclude that the DNA-dye test is a surrogate bioassay suitable for screening chemical toxicity.

Computational exploration of mobile ion distributions around RNA duplex

Atomically detailed distributions of ions around an A-form RNA are computed. Different mixtures of monovalent and divalent ions are considered explicitly. Studies of tightly bound and of diffusive (but bound) ions around 25 base pairs RNA are conducted in explicit solvent. Replica exchange simulations provide detailed equilibrium distributions with moderate computing resources (20 ns of simulation using 64 replicas). The simulations show distinct behavior of single and double charged cations. Binding of Mg(2+) ion includes tight binding to specific sites while Na(+) binds only diffusively. The tight binding of Mg(2+) is with a solvation shell while Na(+) can bind directly to RNA. Negative mobile ions can be found near the RNA but must be assisted by proximate and mobile cations. At distances larger than 16 A from the RNA center, a model of RNA as charged rod in a continuum of ionic solution provides quantitative description of the ion density (the same as in atomically detailed simulation). At shorter distances, the structure of RNA (and ions) has a significant impact on the pair correlation functions. Predicted binding sites of Mg(2+) at the RNA surface are in accord with structures from crystallography. Electric field relaxation is investigated. The relaxation due to solution rearrangements is completed in tens of picoseconds, while the contribution of RNA tumbling continues to a few nanoseconds.

Quantification of single-stranded nucleic acid and oligonucleotide interactions with metal ions by affinity capillary electrophoresis--Part II

The interactions between tetranucleotides or heptanucleotides and inorganic cations have been measured by affinity CZE. The variation in migration behavior with increasing concentrations of divalent cations (Ca(2+), Mg(2+) and Ni(2+)) in the running buffer was investigated and quantified by the calculation of binding constants for mononuclear and multinuclear interactions. In addition to these fundamental studies of binding equilibria, the effect of sequence and the position of the guanine transition metal-binding site in the oligonucleotide have been investigated.

Internal versus terminal metalation of double-helical oligodeoxyribonucleotides

The formation of adducts between cis-[Pt(NH(3))(2)Cl(2)], Zn(II), and Mn(II) and double-stranded oligodeoxynucleotides was studied by 1D and 2D (1)H, (31)P, and (15)N NMR spectroscopy. For labile adducts involving Zn(II) and Mn(II), both (1)H chemical shifts (Zn(II)) and (1)H line-broadening effects (Mn(II)) showed that in the hexamer [d(GGCGCC)](2) I, the terminal G(1)-N7 is the exclusive binding site, while for the dodecamer [d(GGTACCGGTACC)](2) II, which contains both a terminal and internal GG pair, the preference for metal binding is the internal guanine G(7). Zn(II) binding to II was confirmed by natural-abundance 2D [(1)H,(15)N] HMBC NMR spectroscopy, which unambiguously showed that G(7)-N7 is the preferred binding site. The long duplex [d(GGTATATATACCGGTATATATACC)](2) III was expected to have a more pronounced accumulation of electrostatic potential towards the central part of the sequence (vs the terminal part) than does II. However, the Zn(II) titration of III showed no increase in coordination with the internal Gs (vs the terminal Gs), compared with what was observed for II. The reaction between the nonlabile metal complex cis-[PtCl(2)((15)NH(3))(2)] (cisplatin) and II showed a slight preference for the internal GG pair over the terminal GG pair. However, when the diaqua form of cisplatin cis-[Pt((15)NH(3))(2)(H(2)O)(2)] was reacted with II a more pronounced binding preference for the internal GG pair was observed.

Cheminformatics methods for novel nanopore analysis of HIV DNA termini

BACKGROUND

Channel current feature extraction methods, using Hidden Markov Models (HMMs) have been designed for tracking individual-molecule conformational changes. This information is derived from observation of changes in ionic channel current blockade "signal" upon that molecule's interaction with (and occlusion of) a single nanometer-scale channel in a "nanopore detector". In effect, a nanopore detector transduces single molecule events into channel current blockades. HMM analysis tools described are used to help systematically explore DNA dinucleotide flexibility, with particular focus on HIV's highly conserved (and highly flexible/reactive) viral DNA termini. One of the most critical stages in HIV's attack is the binding between viral DNA and the retroviral integrase, which is influenced by the dynamic-coupling induced high flexibility of a CA/TG dinucleotide positioned precisely two base-pairs from the blunt terminus of the duplex viral DNA. This suggests the study of a family of such CA/TG dinucleotide molecules via nanopore measurement and cheminformatics analysis.

RESULTS

HMMs are used for level identification on the current blockades, HMM/EM with boosted variance emissions are used for level projection pre-processing, and time-domain FSAs are used to parse the level-projected waveform for kinetic information. The observed state kinetics of the DNA hairpins containing the CA/TG dinucleotide provides clear evidence for HIV's selection of a peculiarly flexible/interactive DNA terminus.

Preferential binding and structural distortion by Fe2+ at RGGG-containing DNA sequences correlates with enhanced oxidative cleavage at such sequences

Certain DNA sequences are known to be unusually sensitive to nicking via the Fe2+-mediated Fenton reaction. Most notable are a purine nucleotide followed by three or more G residues, RGGG, and purine nucleotides flanking a TG combination, RTGR. Our laboratory previously demonstrated that nicking in the RGGG sequences occurs preferentially 5' to a G residue with the nicking probability decreasing from the 5' to 3'end of these sequences. Using 1H NMR to characterize Fe2+ binding within the duplex CGAGTTAGGGTAGC/GCTACCCTAACTCG and 7-deazaguanine-containing (Z) variants of it, we show that Fe2+ binds preferentially at the GGG sequence, most strongly towards its 5' end. Substitutions of individual guanines with Z indicate that the high affinity Fe2+ binding at AGGG involves two adjacent guanine N7 moieties. Binding is accompanied by large changes in specific imino, aromatic and methyl proton chemical shifts, indicating that a locally distorted structure forms at the binding site that affects the conformation of the two base pairs 3' to the GGG sequence. The binding of Fe2+ to RGGG contrasts with that previously observed for the RTGR sequence, which binds Fe2+ with negligible structural rearrangements.

EDTA-derivatized deoxythymidine as a tool for rapid determination of protein binding polarity to DNA by intermolecular paramagnetic relaxation enhancement

EDTA-derivatized deoxythymidine (dT-EDTA), incorporated into DNA and complexed to Fe2+ in the presence of dithiothreitol, is a widely used reagent for sequence-specific cleavage of duplex DNA. Using HPLC/electrospray mass spectrometry, we show that cleavage is specific to Fe2+, and no cleavage occurs when DNA-EDTA is complexed to other metal ions such as Ca2+, Mn2+, and Fe3+ even after many days. Because dT-EDTA can be incorporated at any desired position of a synthetic oligonucleotide, DNA-EDTA is ideally suited for the measurement of intermolecular paramagnetic relaxation enhancement effects between a paramagnetic ion chelated to DNA-EDTA and a bound protein. Measurements on the SRY/DNA-EDTA complex using two double-stranded oligonucleotides bearing dT-EDTA at opposite ends of the sequence indicate that intermolecular 1HN-T2 enhancement by chelated Mn2+ can be used to readily ascertain the polarity of protein binding to DNA and to derive quantitative long-range distance information for structure refinement. In the case of the SRY-DNA complex, excellent agreement between observed and calculated 1HN-T2 paramagnetic relaxation enhancement data can be achieved with insignificant shifts in the atomic coordinates ( approximately 0.25 A for all heavy atoms) while simultaneously satisfying all other experimental restraints. A unique feature of DNA-EDTA is that the relaxation enhancement effect can be tuned by judicious choice of the paramagnetic metal ion, thereby permitting a wide range of long-range intermolecular electron-proton distances, ranging from approximately 9 to 35 A, to be probed.

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

Interaction between macrocyclic nickel complexes and the nucleotides GMP, AMP and ApG

Reactions between the nucleotides GMP, AMP and ApG and the complexes Ni(tren), Ni(cyclam) and NiCR in aqueous solution have been monitored by (1)H, (15)N NMR and UV spectroscopy. The three nickel complexes display different properties in reactions with nucleotides. Ni(tren) which has a pseudo-octahedral coordination geometry was shown to bind to all three nucleotides. Ni(cyclam) and NiCR, both with four nitrogen atoms in a square planar arrangement are not able to bind to nucleotides efficiently because of steric hindrance. Oxidation of Ni(cyclam) by KHSO(5) to produce trivalent Ni(III)(cyclam) improves the coordination capacity, while oxidation of NiCR does not produce a similar effect. The nucleotides interact with trivalent nickel complexes to different extent. Ni(III)CR is seen to oxidize GMP gradually but does not affect AMP significantly. Ni(III)(cyclam), on the other hand, does not oxidize either GMP or AMP at the 1:1 concentration of oxidant used. This result is probably due to the lower redox potential of Ni(cyclam). ApG binds less efficiently to the Ni complexes but is easier oxidized than the mononucleotides.

Simultaneous N7,O6-binding of guanine to two zinc centers and its possible biological significance

The reaction of ZnCl(2) with 9-ethylguanine (9-EtGH) produced a novel dinuclear Zn(II) complex, [Zn(2)Cl(4)(H(2)O)(mu-9-EtGH-N7,O6)(9-EtGH-N7), 1. The X-ray structure analysis (monoclinic, P2(1) (No. 4), a = 11.0636(6) A, b = 6.6546(4) A, c = 15.9630(9) A, beta = 101.069(1) degrees, V = 1153.4(1) A(3), Z = 2) revealed that one of the tetrahedrally coordinated Zn(II) atoms binds to the N7 site of 9-EtGH and to the exocyclic O6 atom of another 9-EtGH molecule. The remaining Zn(II) atom binds to the N7 site of the second 9-EtGH moiety.

DNA-dependent divalent cation binding in the nucleosome core particle

Sequence-specific binding of divalent cations to nucleosomal DNA can potentially influence nucleosome position and mobility, as well as modulate interactions with nuclear factors. We define the bonding and specificity of divalent cation interaction with nucleosomal DNA by characterizing Mn2+ binding in the x-ray structure of the nucleosome core particle at 1.9-A resolution. Manganese ions are found ligated with high occupancy in the major groove to 12 of 22 GG and GC base pair steps. The specific location and mode of metal binding is the consequence of unambiguous conformational differences between dinucleotide sites, owing to their sequence context and orientation in the nucleosome core.

Characterization of divalent cation localization in the minor groove of the A(n)T(n) and T(n)A(n) DNA sequence elements by (1)H NMR spectroscopy and manganese(II)

The localization of Mn(2+) in A-tract DNA has been studied by (1)H NMR spectroscopy using a series of self-complementary dodecamer oligonucleotides that contain the sequence motifs A(n)(n) and T(n)A(n), where n = 2, 3, or 4. Mn(2+) localization in the minor groove is observed for all the sequences that have been studied, with the position and degree of localization being highly sequence-dependent. The site most favored for Mn(2+) localization in the minor groove is near the 5'-most ApA step for both the T(n)A(n) and the A(n)T(n) series. For the T(n)A(n) series, this results in two closely spaced symmetry-related Mn(2+) localization sites near the center of each duplex, while for the A(n)T(n) series, the two symmetry-related sites are separated by as much as one half-helical turn. The degree of Mn(2+) localization in the minor groove of the T(n)A(n) series decreases substantially as the AT sequence element is shortened from T(4)A(4) to T(2)A(2). The A(n)T(n) series also exhibits length-dependent Mn(2+) localization; however, the degree of minor groove occupancy by Mn(2+) is significantly less than that observed for the T(n)A(n) series. For both A(n)T(n) and T(n)A(n) sequences, the 3'-most AH2 resonance is the least broadened of the AH2 resonances. This is consistent with the observation that the minor groove of A-tract DNA narrows in the 5' to 3' direction, apparently becoming too narrow after two base pairs for the entry of a fully hydrated divalent cation. The results that are reported illustrate the delicate interplay that exists between DNA nucleotide sequence, minor groove width, and divalent cation localization. The proposed role of cation localization in helical axis bending by A-tracts is also discussed.

Localization of Fe(2+) at an RTGR sequence within a DNA duplex explains preferential cleavage by Fe(2+) and H2O2

Nicking of duplex DNA by the iron-mediated Fenton reaction occurs preferentially at a limited number of sequences. Of these, purine-T-G-purine (RTGR) is of particular interest because it is a required element in the upstream regulatory regions of many genes involved in iron and oxidative-stress responses. In order to study the basis of this preferential nicking, NMR studies were undertaken on the RTGR-containing duplex oligonucleotide, d(CGCGATATGACACTAG)/d(CTAGTGTCATATCGCG). One-dimensional and two-dimensional 1H NMR measurements show that Fe(2+) interacts preferentially and reversibly at the ATGA site within the duplex at a rate that is rapid relative to the chemical-shift timescale, while selective paramagnetic NMR line-broadening of the ATGA guanine H8 suggests that Fe(2+) interacts with the guanine N7 moiety. Localization at this site is supported by Fe(2+) titrations of a duplex containing a 7-deazaguanine substitution in place of the guanine in the ATGA sequence. The addition of a 100-fold excess of Mg(2+) over Fe(2+) does not affect the Fe(2+)-dependent broadening. When the ATGA site in the duplex is replaced by ATGT, an RTGR site (GTGA) is created on the opposite strand. Preferential iron localization then takes place at the 3' guanine in GTGA but no longer at the guanine in ATGT, consistent with the lack of preferential cleavage of ATGT sites relative to ATGA sites.

Selective intercalation of charge neutral intercalators into GG and CG steps: implication of HOMO-LUMO interaction for sequence-selective drug intercalation into DNA

We have synthesized naphthopyranone epoxide 4 from D-isoascorbic acid together with its three diastereoisomers. DNA alkylation of ODNs containing 5'XGT3' and 5'TGY3' by 4 (11R, 13R), where X and Y are any nucleotide bases, occurred at all G residues except at G of the 5'TGC3' sequence. In contrast, the three other diastereoisomers of 4 showed only weak G alkylation activity. Differential (1)H NMR NOE of the 4-G adduct confirmed the G-N7 alkylation at the epoxide carbon of 4 with concomitant S(N)2 ring opening of the epoxide. Quantitative HPLC analysis of G alkylation efficiency for 4 showed the order of G alkylation susceptibility as TGGT approximately CGT >> TGA > AGT > TGT >> TGC. The order was fully consistent with those reported for aflatoxin B(1) oxide and kapurimycin A(3), suggesting that the sequence selectivity observed for these DNA alkylating agents is not structure dependent but most likely due to the intrinsic property of DNA sequences. We found that the order of G alkylation susceptibility obtained for 4 completely matched the calculated HOMO energy level of G-containing sequences. These results underscore that 4 is a unique molecular probe for ranking the HOMO level of G-containing sequences by well-known G alkylation chemistry and suggests that the intercalation of charge neutral intercalators is a HOMO-controlled process.

Superoxide dismutases enhance H2O2-induced DNA damage and alter its site specificity

Superoxide dismutases (SODs) are involved in the protection of cells from oxygen toxicity. However, several papers have reported that the overexpression of CuZn-SOD causes oxidative damage to cells. We investigated a mechanism by which an excess of SODs accelerates oxidative stress. The presence of CuZn-SOD, Mn-SOD or Mn(II) enhanced the frequency of DNA damage induced by hydrogen peroxide (H2O2) and Cu(II), and altered the site specificity of the latter: H2O2 induced Cu(II)-dependent DNA damage with high frequency at the 5'-guanine of poly G sequences; when SODs were added, the frequency of cleavages at thymine and cytosine residues increased. SODs also enhanced the formation of 8-oxo-7,8-dihydro-2'-deoxyguanosine by H2O2 and Cu(II). We conclude that SODs may increase carcinogenic risks, e.g. of tumors in Down syndrome.

Mechanism of guanine-specific DNA damage by oxidative stress and its role in carcinogenesis and aging

Reactive species generated by chemicals and UV radiation can cause sequence-specific DNA damage and play important roles in mutagenesis, carcinogenesis and aging. We have investigated sequence specificity of oxidative stress-mediated DNA damage by using 32P-labeled DNA fragments obtained from the human c-Ha-ras-1 and p53 genes. Free hydroxyl radical causes DNA damage with no marked site specificity. Reactive nitrogen species, sulfate radicals, nitrogen-centered radicals, benzoyloxyl radical and alkoxyl radical show different sequence specificity. Benzoyloxyl radical specifically causes damage to the 5'-G in GG sequence. UVA radiation also causes DNA damage at this site through electron transfer in the presence of certain photosensitizers. The 5'-G in GG sequence is easily oxidized because a large part of the highest occupied molecular orbital is distributed on this site. On the basis of these findings, the sequence specificity of DNA damage is presumably determined by (a) redox potential of reactive species; (b) ionization potential of DNA bases; and (c) site-specific binding of metal ion to DNA. Here we discuss the mechanisms of sequence-specific DNA damage in relation to carcinogenesis and aging.

Parallel intramolecular DNA triple helix with G and T bases in the third strand stabilized by Zn(2+) ions

We present evidence of formation of an intramolecular parallel triple helix with T*A.T and G*G.C base triplets (where * represents the hydrogen bonding interaction between the third strand and the duplex while. represents the Watson-Crick interactions which stabilize the duplex). The third GT strand, containing seven GpT/TpG steps, targets the polypurine sequence 5'-AGG-AGG-GAG-GAG-3'. The triple helix is obtained by the folding back twice of a 36mer, formed by three dodecamers tethered by hydroxyalkyl linkers (-L-). Due to the design of the oligonucleotide, the third strand orientation is parallel with respect to the polypurine strand. Triple helical formation has been studied in concentration conditions in which native gel electrophoresis experiments showed the absence of intermolecular structures. Circular dichroism (CD) and UV spectroscopy have been used to evidence the triplex structure. A CD spectrum characteristic of triple helical formation as well as biphasic UV and CD melting curves have been obtained in high ionic strength NaCl solutions in the presence of Zn(2+) ions. Specific interactions with Zn(2+) ions in low water activity conditions are necessary to stabilize the parallel triplex.

Stacking interaction of guanine with netropsin in the minor groove of d(CGTATATACG)2

We present the structure of the decanucleotide d(CGTATATACG) determined by single crystal X-ray diffraction at 1.58 A resolution. A netropsin drug is found in the minor groove with guanine stacked on a pyrrole ring of the drug, a feature described here for the first time. The stacked guanine is an extra-helical base coming from the end of a neighbour oligonucleotide. This observation may open the way to the development of minor groove binding drugs with a higher sequence selectivity. The oligonucleotide is in the B-conformation, but the terminal base-pairs are disrupted: the cytosine residues are disordered while the guanine residues penetrate into the minor groove of neighbouring duplexes. Four hydrated Ni ions with octahedral co-ordination are found associated with the N7 atoms of each guanine. The high affinity of these ions with guanine suggests that they may be used as probes for specific guanine residues.

Selective association between a macrocyclic nickel complex and extrahelical guanine residues

Nickel-dependent recognition and oxidation of guanine have been linked in part through the paramagnetic effects of nickel on the NMR of model oligonucleotide duplexes. Direct interaction between nickel and guanine N7 had originally been postulated from correlations between the efficiency of guanine oxidation and the environment surrounding its N7 position. (1)H and (31)P NMR spectra of DNA containing a single, isolated extrahelical guanine are consistent with selective binding of nickel to the N7 of this unique base over a background of nonspecific association to the phosphate backbone. The presence of a macrocyclic complex or simple salt of nickel did not detectably alter the structure of the duplex or extrahelical residue. Accordingly, nickel appeared to bind the extrahelical guanine N7 within the major groove as indicated by paramagnetic effects on the proton signals of nucleotides on the 5' but not 3' side of the nickel binding site. Similar (1)H NMR analysis of DNA containing a dynamic equilibrium of extrahelical guanine residues also suggested that the nickel complex did not affect the native distribution of structures. Oxidation of these sites by a nickel-mediated pathway consequently reflected their solvent accessibility in a general and metal-independent manner. The close proximity of the extrahelical guanines produced a composite of paramagnetic effects on each adjacent nucleotide resulting from both direct and proximal coordination of nickel.

Localization of ammonium ions in the minor groove of DNA duplexes in solution and the origin of DNA A-tract bending

Monovalent cation binding sites on nucleic acids in solution can be localized using the isotopically labeled ammonium ion (15NH4+) as a probe in high resolution NMR spectroscopy experiments. The application of this technique to a series of DNA duplexes reveals a preference for the binding of ammonium cations in the minor groove of A-tract sequences. These results are consistent with a recent report which indicates that some solvent electron densities previously identified as water molecules in DNA X-ray crystal structures are partially occupied by sodium ions. The sequence-specific nature of monovalent cation binding sites demonstrated here for A-tract DNA provides an explanation for the origin of sequence-directed bending.

Sequence-dependent binding of metal ion to DNA oligomeres. A comparison of molecular electrostatic potentials with NMR data

Experimentally observed sequence-selective binding of metal ion to DNA oligonucleotides have been compared with variations of electrostatic potential (EP) along the helix. Calculations of EP have been performed for three atomic models of the oligonucleotide duplex [d(CGCGAATTCGCG)2] using several variants of EP calculations, including a solution of non-linear Poisson-Boltzmann equation (NPBE). N7 atom of guanine adjacent to adenine base was identified as a region with the most negative electrostatic potential in the major groove. The EP value for the Me ion binding site surpasses the value for N7 of other guanines by 10-26% depending on particular duplex conformation. Qualitatively, the sequence dependent variations of EP near guanine N7 atoms are in agreement with the sequence-selective behavior of Mn(II) and Zn(II) ions as revealed by NMR experiments. But the difference in EP between the two most negative regions near guanine N7 atoms does not exceed 1.25 kT/e. Simple model suggests that metal ions are capable to form ion-hydrate complexes with G-Pu steps of DNA duplex. These complexes are formed via one Me...G and five Me...water coordination bonds with water molecules hydrogen bonded to two adjacent purine bases in the same chain. We suppose that such a stereospecific structural possibility is the main factor which control the sequence-selectivity in the metal ion binding. A combination of both mechanisms allows to explain sequence specific Mn(II) and Zn(II) binding to a set of oligonucleotides.

Interaction of Co, Mn, Mg and Al with d(GCGTACGC): a spectroscopic study

Spectroscopic study on the interactions of trace elements Co., Mn, Mg and Al with d(GCGTACGC) indicated the following: Al and Mg did not alter Tm values. Mn enhanced Tm at lower concentration and decreased it at higher concentrations. Interestingly Co at higher concentration elevated the Tm. These studies also showed lower concentrations of Mn displaced EtBr, whereas Al could displace it at higher ionic strength. Mg and Co displaced EtBr fluorescence at moderate concentrations. The binding constant values and CD spectra clearly indicated strong binding of these elements to DNA.