DNA Base Stack Research Articles

Kinetics of DNA force-induced melting

Once stretched beyond its B-form contour length, double-stranded DNA reveals a sudden increase in length over approximately constant force at about 65 pN. During this conformational transition, DNA base pairing and base stacking are disrupted, converting double stranded DNA (dsDNA) into single stranded DNA (ssDNA). While thermodynamic and recent chemical labeling and fluorescence imaging experiments indicate that this transition is analogous to thermal melting, the kinetics of DNA force-induced melting have not been characterized. We present a predictive model of force-induced melting in which thermal fluctuations induce local melting and re-annealing of DNA. These fluctuations are stabilized by the application of tension during the overstretching transition, favoring the conversion to ssDNA as the applied force is increased. This model quantitatively predicts small changes in the melting force as the pulling rate is varied. We verify that the DNA melting force varies with pulling rate, consistent with this model, and that DNA force-induced melting depends only weakly on pulling rate at slow pulling rates, as melting occurs cooperatively with a domain size of 100-200 base pairs. As the pulling rate is increased beyond the natural duplex opening rate, the melting force depends strongly on pulling rate and the melted domain size decreases to 5-10 base pairs, as the DNA is ripped sequentially from the free ends (or any boundary). The final strand separation occurs at much higher forces, representing the nonequilibrium ripping of the most stable regions that remain at the end of the low force transition. The results indicate that force only weakly enhances base pair opening, while strongly inhibiting base pair closing.

Characterization of The Open Complex of Yeast Mitochondrial RNA Polymerase

The mitochondrial RNA polymerase (mtRNAP) of Saccharomyces cerevisiae, consisting of a complex of Rpo41 and Mtf1, is homologous to the phage single polypeptide T7/T3/SP6 RNA polymerase and mammalian mitochondrial RNA polymerases. The yeast mtRNAP recognizes a conserved nonanucleotide sequence to initiate specific transcription. However, there is little clue about how mtRNAP recognizes the promoter, melts the double-stranded DNA and forms a transcription-competent open complex, which are steps essential to initiating RNA synthesis. In this work, we have defined the region of the nonanucleotide melted by the mtRNAP using 2-aminopurine fluorescence changes that are sensitive to changes in DNA base stacking interactions. We have also characterized DNA bending of the promoter induced by mtRNAP by fluorescence resonance energy transfer (FRET) measurements. We show that DNA melting requires the simultaneous presence of Rpo41 and Mtf1, whereas DNA bending can be induced by Rpo41 itself and regulated by MtF1. By monitoring DNA melting and bending under various conditions, we are able to dissect roles of Rpo41 and MtF1 in DNA recognition and forming the open complex. These results shed new insights into mechanisms of transcription initiation of mitochondrial RNA polymerases.

Induction of Cell Death by Ternary Copper(II) Complexes of l-Tyrosine and Diimines: Role of Coligands on DNA Binding and Cleavage and Anticancer Activity

The mononuclear mixed ligand copper(II) complexes of the type [Cu(L-tyr)(diimine)](ClO(4)), where tyr is L-tyrosine and diimine is 2,2'-bipyridine (bpy) (1), 1,10-phenanthroline (phen) (2), 5,6-dimethyl-1,10-phenanthroline (5,6-dmp) (3), and dipyrido[3,2-d:2',3'-f]quinoxaline (dpq) (4), have been isolated and characterized by analytical and spectral methods. In the X-ray crystal structure 3 Cu(II) possesses a distorted square pyramidal coordination geometry with the two nitrogen atoms of 5,6-dmp ligand and the amine nitrogen and carboxylate oxygen atoms of L-tyrosine located at the equatorial sites and the coordinated water molecule present in the apical position. The electronic absorption and electron paramagnetic resonance (EPR) spectral parameters reveal that the complexes retain their square-based geometries even in solution. All of the complexes display a ligand field band in the visible region (600-700 nm) in Tris-HCl/NaCl buffer (5:50 mM) at pH 7.2 and also axial EPR spectra in acetonitrile at 77 K with g(parallel) > g(perpendicular) indicating a d(x(2)-y(2)) ground state. The g(parallel) and A(parallel) values of 2.230 and (170-180) x 10(-4) cm(-1), respectively, conform to a square-based CuN(3)O coordination chromophore, which is consistent with the X-ray crystal structure of 3. The interaction of the complexes with calf thymus DNA (CT DNA) has been explored by using physical methods to propose modes of DNA binding of the complexes. Absorption (K(b)) and emission spectral studies and viscosity measurements indicate that 4 interacts with DNA more strongly than all of the other complexes through partial intercalation of the extended planar ring of dpq with DNA base stack. Interestingly, complex 3 exhibits a DNA binding affinity that is higher than that of 2, which suggests the involvement of 5,6-dimethyl groups on the phen ring in hydrophobic interaction with DNA surface. In contrast with the increase in relative viscosities of DNA bound to 2-4, the viscosity of DNA bound to 1 decreases, indicating the shortening of the DNA chain length by means of the formation of kinks or bends. All complexes exhibit effective DNA (pUC19 DNA) cleavage at 100 microM complex concentrations, and the order of DNA cleavage ability varies as 3 > 2 > 4 > 1. Interestingly, 3 exhibits a DNA cleavage rate constant that is higher than that of the other complexes only at 100 microM concentration, whereas 4 exhibits the highest cleavage rate constant at 80 microM complex concentration. The oxidative DNA cleavage follows the order 4 > 3 > 2 > 1. Mechanistic studies reveal that the DNA cleavage pathway involves hydroxyl radicals. Interestingly, only 4 displays efficient photonuclease activity upon irradiation with 365 nm light, which occurs through double-strand DNA breaks involving hydroxyl radicals. Furthermore, cytotoxicity studies on the nonsmall lung cancer (H-460) cell line show that the IC(50) values of 2-4 are more or less equal to cisplatin for the same cell line, indicating that they have the potential to act as very effective anticancer drugs in a time-dependent manner. The study of cytological changes reveals the higher induction of apoptosis and mitotic catastrophe for 4 and 3, respectively. The alkaline single-cell gel electrophoresis (comet assay), DNA laddering, and AO/EB and Hoechst 33258 staining assays have also been employed in finding the extent of DNA damage. Flow cytometry analysis shows an increase in the percentage of cells with apoptotic morphological features in the sub-G(0)/G(1) phase for 4, whereas it shows mitotic catastrophe for 3.

Thermal melting studies of alkyne- and ferrocene-containing PNA bioconjugates

The preparation of new metal-containing Peptide Nucleic Acids (PNAs) is currently a field of research intensively studied for various purposes, e.g. DNA biosensors. The role played by the metal centre, notably on the stability of the PNA.DNA hybrid, is obviously crucial, but has not yet been fully investigated. In this work, UV-Vis spectroscopic measurements of solutions of DNA.PNA hybrids, whose 11/12-mer PNA oligomers contained either one or two alkyne- (1) or ferrocene-containing (2) PNA monomers, were carried out to determine the effect of these monomers on the thermal stability of the hybrids (PNA: H-Gly-X-gggtc-Y-agctt-X-Lys-NH2 with X = 1 or and Y = 1 or 2 or blank position). Supplementary CD spectroscopic measurements were performed to gain insight into the structures of the PNA.DNA duplexes formed. The effect of both modified monomers was found to depend on their actual positions within the PNA sequences. Insertions at the N- or C-termini of a PNA oligomer did not change the melting temperatures (T(m) values of about 72 degrees C) of the DNA.PNA hybrids significantly. Insertion of monomers 1 or 2 in the middle of a PNA sequence induced a substantial decrease in the T(m) of the hybrids (by about 23 degrees C) when bound to the same DNA oligomer. Interestingly, it was found that the type of modification, namely alkyne or ferrocene, did not significantly influence the T(m) values in these cases. However, the thermal stability of hybrids with the DNA oligomers containing one to four additional thymines and the PNA oligomers containing the ferrocene moiety in its middle, varied significantly with the number of thymines added compared to its alkyne analogues (DeltaT(m) up to -13 degrees C). The presence of the ferrocene moiety induced a significant decrease in thermal stability of the hybrids, probably due to its bulkiness. In order to assess the effect of PNA backbone rigidity on the stability of DNA.PNA hybrids, PNA oligomers with an internal amino acid, propargylglycine (Pgl) or the dipeptide glycine-propargylglycine (Gly-Pgl), were synthesised. It was assumed that the orientation of the alkyne moiety in the Pgl-containing PNA sequence is not identical to an alkyne-containing PNA sequence, as a significantly higher T(m) value (DeltaT(m) = +10 degrees C) was measured. It is anticipated that the alkyne moiety in Pgl is not facing the DNA base and therefore does not disturb as much the neighbouring nucleobases and base-stacking of the complementary DNA, in contrast to the alkyne moiety of 1. Interestingly, no significant differences in the thermal stability of the hybrids was observed between Pgl-containing and dipeptide-containing PNA oligomers, although the former contracts the PNA backbone by three atoms.

Synthesis, DNA binding and cleavage studies of Ni(II) complexes with fused aromatic N-containing ligands

The three Ni(II) complexes of fused aromatic N-containing ligands such as [Ni(bnp) 3](PF 6) 2 ( 1), [Ni(phen) 2(bnp)](PF 6) 2 ( 2) and [Ni(bpy) 2(bnp)](PF 6) 2 ( 3) (where bnp = dibenzo(b)1,8-naphthpyridine, phen = 1,10-phenanthroline and bpy = bipyridine) were synthesized and structurally characterized. Elemental analysis, magnetic and spectroscopic data suggested octahedral geometry for all the complexes. Binding of these complexes with (ds)DNA were analyzed by absorption spectra, viscosity and thermal denaturation studies. Detailed analysis revealed that the metal complexes intercalates into the DNA base stack as intercalator. The oxidative cleavage activities of the complexes were studied with supercoiled (SC)pUC19 DNA by using gel electrophoresis, and the results show that complexes have potent nuclease activity.

Ping-pong electron transfer through DNA.

Either way: An IrIII complex strongly coupled to a DNA base stack can inject a hole or an electron in DNA upon irradiation (see picture). With strands containing the modified bases CPA and CPC/BrU, which serve as oxidative and reductive probes, respectively, photolysis of Ir–DNA conjugates leads to a “ping-pong” electron transfer, with both hole and electron migration.

DNA binding property, nuclease activity and cytotoxicity of Zn(II) complexes of terpyridine derivatives

Two zinc(II) terpyridine complexes Zn(atpy)(2)(PF(6))(2) (1) (atpy = 4'-p-N9'-adeninylmethylphenyl-2,2':6,2''-terpyridine) and Zn(ttpy)(2)(PF(6))(2) (2) (ttpy = 4'-p-tolyl-2,2':6,2''-terpyridine) have been synthesized and characterized by elemental analysis, (1)H NMR and electrospray mass spectroscopy. The structure of complex 2 was also determined by X-ray crystallography, which revealed a ZnN(6) coordination in an octahedral geometry with two terpyridine acting as equatorial ligands. The circular dichroism data showed that complex 1 exhibited an ICD signal at around 300 nm and induced more evident disturbances on DNA base stacking than complex 2, reflecting the impact of the adenine moiety on DNA binding modes. Complex 1 exhibited higher cleavage activity to supercoiled pUC 19 DNA than complex 2 under aerobic conditions, suggesting a promotional effect of adenine moiety in DNA nuclease ability. Interestingly, both complexes demonstrated potent in vitro cytotoxicity against a series human tumor cell lines such as human cervix carcinoma cell line (HeLa), human liver carcinoma cell line (HepG2), human galactophore carcinoma cell line (MCF-7) and human prostate carcinoma cell line (pc-3). The cytotoxicity is averagely 10 times more active than the anticancer drug cisplatin.

Abundance and length of simple repeats in vertebrate genomes are determined by their structural properties

Microsatellites are abundant in vertebrate genomes, but their sequence representation and length distributions vary greatly within each family of repeats (e.g., tetranucleotides). Biophysical studies of 82 synthetic single-stranded oligonucleotides comprising all tetra- and trinucleotide repeats revealed an inverse correlation between the stability of folded-back hairpin and quadruplex structures and the sequence representation for repeats > or =30 bp in length in nine vertebrate genomes. Alternatively, the predicted energies of base-stacking interactions correlated directly with the longest length distributions in vertebrate genomes. Genome-wide analyses indicated that unstable sequences, such as CAG:CTG and CCG:CGG, were over-represented in coding regions and that micro/minisatellites were recruited in genes involved in transcription and signaling pathways, particularly in the nervous system. Microsatellite instability (MSI) is a hallmark of cancer, and length polymorphism within genes can confer susceptibility to inherited disease. Sequences that manifest the highest MSI values also displayed the strongest base-stacking interactions; analyses of 62 tri- and tetranucleotide repeat-containing genes associated with human genetic disease revealed enrichments similar to those noted for micro/minisatellite-containing genes. We conclude that DNA structure and base-stacking determined the number and length distributions of microsatellite repeats in vertebrate genomes over evolutionary time and that micro/minisatellites have been recruited to participate in both gene and protein function.

Uracil DNA glycosylase uses DNA hopping and short-range sliding to trap extrahelical uracils

The astonishingly efficient location and excision of damaged DNA bases by DNA repair glycosylases is an especially intriguing problem in biology. One example is the enzyme uracil DNA glycosylase (UNG), which captures and excises rare extrahelical uracil bases that have emerged from the DNA base stack by spontaneous base pair breathing motions. Here, we explore the efficiency and mechanism by which UNG executes intramolecular transfer and excision of two uracil sites embedded on the same or opposite DNA strands at increasing site spacings. The efficiency of intramolecular site transfer decreased from 41 to 0% as the base pair spacing between uracil sites on the same DNA strand increased from 20 to 800 bp. The mechanism of transfer is dominated by DNA hopping between landing sites of approximately 10 bp size, over which rapid 1D scanning likely occurs. Consistent with DNA hopping, site transfer at 20- and 56-bp spacings was unaffected by whether the uracils were placed on the same or opposite strands. Thus, UNG uses hopping and 3D diffusion through bulk solution as the principal pathways for efficient patrolling of long genomic DNA sequences for damage. Short-range sliding over the range of a helical turn allows for redundant inspection of very local DNA sequences and trapping of spontaneously emerging extrahelical uracils.

Intercalation of Daunomycin into Stacked DNA Base Pairs. DFT Study of an Anticancer Drug

We have computationally studied the intercalation of the antitumor drug daunomycin into six stacks of Watson-Crick DNA base pairs i.e., AT-AT, AT-TA, GC-AT, CG-TA, GC-GC, GC-CG) using density functional theory (DFT). The proton affinity of the DNA intercalater daunomycin in water was computed to be 159.2 kcal/mol at BP86/TZ2P, which is in line with the experimental observation that daunomycin is protonated under physiological conditions. The intercalation interaction of protonated daunomycin with two stacked DNA base pairs was studied through a hybrid approach in which intercalation is treated at LDA/TZP while the molecular structure of daunomycin and hydrogen-bonded Watson-Crick pairs is computed at BP86/TZ2P We find that the affinity of the drug for the six considered base pair dimers decreases in the order AT-AT > AT-TA > GC-AT > GC-TA > GC-CG > GC-GC, in excellent agreement with experimental data on the thermodynamics of the interaction between daunomycin and synthetic polynucleotides in aqueous solution. Our analyses show that the overall stability of the intercalation complexes comes mainly from π-π stacking but an important contribution to the computed and experimentally observed sequence specificity comes from hydrogen bonding between daunomycin and hetero atoms in the minor groove of AT base pairs.

Thermodynamics of single strand DNA base stacking

The thermodynamics of the stacking to unstacking transitions of 24 single-stranded DNA sequences (ssDNA), 10-12 bases in length, in sodium phosphate buffer were determined from 10 to 95 degrees C, using differential scanning calorimetry (DSC). An additional 22 ssDNA sequences did not exhibit an S<=>U transition in this temperature range. The transition properties of the ssDNA sequences with <or=60% self-complementarily in the reverse direction were independent of concentration with transition temperatures ranging from 15 to 70 degrees C, van't Hoff transition enthalpies from 92 to 201 kJ mol(-1) and transition enthalpies from 5 to 75% of the corresponding van't Hoff transition enthalpies. Since all the 16 doublets and 60 of the 64 triplets are present in both the transition and the non-transition ssDNA sequences, it is unlikely that the nucleation subset initiating stacking of the sequence is a specific doublet or triplet subset. Of the 141 quadruplet subsets of the 46 sequences, each transition ssDNA sequence contained at least one or more quadruplets not found in the non-transition ssDNA sequences. It could be concluded that the thermal stability of the stacked conformation was dependent on the presence of a possible nucleation quadruplet and the length of the ssDNA sequence and not on the G or C content of the ssDNA sequence, nor on the number of purine bases in the sequence.

Investigations of Oligonucleotide Usage Variance Within and Between Prokaryotes

Oligonucleotide usage in archaeal and bacterial genomes can be linked to a number of properties, including codon usage (trinucleotides), DNA base-stacking energy (dinucleotides), and DNA structural conformation (di- to tetranucleotides). We wanted to assess the statistical information potential of different DNA ‘word-sizes’ and explore how oligonucleotide frequencies differ in coding and non-coding regions. In addition, we used oligonucleotide frequencies to investigate DNA composition and how DNA sequence patterns change within and between prokaryotic organisms. Among the results found was that prokaryotic chromosomes can be described by hexanucleotide frequencies, suggesting that prokaryotic DNA is predominantly short range correlated, i.e., information in prokaryotic genomes is encoded in short oligonucleotides. Oligonucleotide usage varied more within AT-rich and host-associated genomes than in GC-rich and free-living genomes, and this variation was mainly located in non-coding regions. Bias (selectional pressure) in tetranucleotide usage correlated with GC content, and coding regions were more biased than non-coding regions. Non-coding regions were also found to be approximately 5.5% more AT-rich than coding regions, on average, in the 402 chromosomes examined. Pronounced DNA compositional differences were found both within and between AT-rich and GC-rich genomes. GC-rich genomes were more similar and biased in terms of tetranucleotide usage in non-coding regions than AT-rich genomes. The differences found between AT-rich and GC-rich genomes may possibly be attributed to lifestyle, since tetranucleotide usage within host-associated bacteria was, on average, more dissimilar and less biased than free-living archaea and bacteria.

DNA binding shifts the redox potential of the transcription factor SoxR

Electrochemistry measurements on DNA-modified electrodes are used to probe the effects of binding to DNA on the redox potential of SoxR, a transcription factor that contains a [2Fe-2S] cluster and is activated through oxidation. A DNA-bound potential of +200 mV versus NHE (normal hydrogen electrode) is found for SoxR isolated from Escherichia coli and Pseudomonas aeruginosa. This potential value corresponds to a dramatic shift of +490 mV versus values found in the absence of DNA. Using Redmond red as a covalently bound redox reporter affixed above the SoxR binding site, we also see, associated with SoxR binding, an attenuation in the Redmond red signal compared with that for Redmond red attached below the SoxR binding site. This observation is consistent with a SoxR-binding-induced structural distortion in the DNA base stack that inhibits DNA-mediated charge transport to the Redmond red probe. The dramatic shift in potential for DNA-bound SoxR compared with the free form is thus reconciled based on a high-energy conformational change in the SoxR-DNA complex. The substantial positive shift in potential for DNA-bound SoxR furthermore indicates that, in the reducing intracellular environment, DNA-bound SoxR is primarily in the reduced form; the activation of DNA-bound SoxR would then be limited to strong oxidants, making SoxR an effective sensor for oxidative stress. These results more generally underscore the importance of using DNA electrochemistry to determine DNA-bound potentials for redox-sensitive transcription factors because such binding can dramatically affect this key protein property.

The electronic structure of the four nucleotide bases in DNA, of their stacks, and of their homopolynucleotides in the absence and presence of water

Using the ab initio Hartree-Fock crystal orbital method in its linear combination of atomic orbital form, the energy band structure of the four homo-DNA-base stacks and those of poly(adenilic acid), polythymidine, and polycytidine were calculated both in the absence and presence of their surrounding water molecules. For these computations Clementi's double zeta basis set was applied. To facilitate the interpretation of the results, the calculations were supplemented by the calculations of the six narrow bands above the conduction band of poly(guanilic acid) with water. Further, the sugar-phosphate chain as well as the water structures around poly(adenilic acid) and polythymidine, respectively, were computed. Three important features have emerged from these calculations. (1) The nonbase-type or water-type bands in the fundamental gap are all close to the corresponding conduction bands. (2) The very broad conduction band (1.70 eV) of the guanine stack is split off to seven narrow bands in the case of poly(guanilic acid) (both without and with water) showing that in the energy range of the originally guanine-stack-type conduction band, states belonging to the sugar, to PO(4)(-), to Na(+), and to water mix with the guanine-type states. (3) It is apparent that at the homopolynucleotides with water in three cases the valence bands are very similar (polycytidine, because it has a very narrow valence band, does not fall into this category). We have supplemented these calculations by the computation of correlation effects on the band structures of the base stacks by solving the inverse Dyson equation in its diagonal approximation taken for the self-energy the MP2 many body perturbation theory expression. In all cases the too large fundamental gap decreased by 2-3 eV. In most cases the widths of the valence and conduction bands, respectively, decreased (but not in all cases). This unusual behavior is most probably due to the rather large complexity of the systems. From all this emerges the following picture for the charge transport in DNA: There is a possibility in short segments of the DNA helix of a Bloch-type conduction of holes through the nucleotide base stacks of DNA combined with hopping (and in a lesser degree with tunneling). The motivation of this large scale computation was that recently in Zurich (ETH) they have performed high resolution x-ray diffraction experiments on the structure of the nucleosomes. The 8 nucleohistones in them are wrapped around by a DNA superhelix of 147 base pairs in the DNA B form. The most recent investigations have shown that between the DNA superhelix (mostly from its PO(4) (-) groups) there is a charge transfer to the positively charged side chains (first of all arginines and lysines) of the histones at 120 sites of the superhelix. This would cause a hole conduction in DNA and an electronic one in the proteins.

Mixed ligand ruthenium(ii) complexes of bis(pyrid-2-yl)-/bis(benzimidazol-2-yl)-dithioether and diimines: Study of non-covalent DNA binding and cytotoxicity

A series of mixed ligand ruthenium(II) complexes [Ru(pdto)(diimine)](ClO4)2/(PF6)2 1-3 and [Ru(bbdo)(diimine)](ClO4), 4-6, where pdto is 1,8-bis(pyrid-2-yl)-3,6-dithiooctane, bbdo is 1,8-bis(benzimidazol-2-yl)-3,6-dithiooctane and diimine is 1,10-phenanthroline (phen), dipyrido-[3,2-d:2',3'-f]-quinoxaline (dpq) and dipyrido[3,2-a:2',3'-c]phenazine (dppz), have been isolated and characterised by analytical and spectral methods. The complexes [Ru(pdto)(phen)](PF6)2 la, [Ru(pdto)(dpq)(Cl](PF6) 2a, [Ru(bbdo)(phen)](PF6)2 4a and [Ru(bbdo)(dpq)](ClO4)2 5 have been structurally characterized and their coordination geometries around ruthenium(II) are described as distorted octahedral. In la, 4a and 5 the two thioether sulfur and two py/bzim nitrogen atoms of the tetradentate pdto/bbdo ligand are folded around Ru(II) to give predominantly a "cis-alpha" configuration. (I)H NMR spectral data of the complexes support this configuration in solution. In [Ru(pdto)(dpq)Cl](PF6) 2a with a distorted octahedral coordination geometry, one of the two py nitrogens of pdto is not coordinated. The DNA binding constants (Kb: 2, 2.00 +/- 0.02 x 10(4) M(-1), s = 1.0; 3, 3.00 +/- 0.01 x 10(6) M(-1), s = 1.3) determined by absorption spectral titrations of the complexes with CT DNA reveal that 3 interacts with DNA more tightly than 2 through partial intercalation of the extended planar ring of coordinated dppz with the DNA base stack. The DNA binding affinities of the complexes increase with increase in the number of planar aromatic rings in the co-ligand, and on replacing both the py moieties in pdto complexes (1-3) by bzim moieties to give bbdo complexes (4-6). Upon interaction with CT DNA the complexes 1, 2, 5 and 6 show a decrease in anodic current in the cyclic voltammograms. On the other hand, interestingly, 3 and 4 show an increase in anodic current suggesting their involvement in electrocatalytic guanine oxidation. Interestingly, of all the complexes, only 6 alters the superhelicity of DNA upon binding with supercoiled pBR322 DNA. The cytotoxicities of the dppz complexes 3 and 6, which avidly bind to DNA, have been examined by screening them against cell lines of different cancer origins. It is noteworthy that 6 exhibits selectivity with higher cytotoxicity against the melanoma cancer cell line (A375) than other cell lines, potency approximately twice that of cisplatin and toxicity to normal cells 3 and 90 times less than cisplatin and adriamycin respectively.

Pi-pi stacking tackled with density functional theory.

Through comparison with ab initio reference data, we have evaluated the performance of various density functionals for describing π-π interactions as a function of the geometry between two stacked benzenes or benzene analogs, between two stacked DNA bases, and between two stacked Watson–Crick pairs. Our main purpose is to find a robust and computationally efficient density functional to be used specifically and only for describing π-π stacking interactions in DNA and other biological molecules in the framework of our recently developed QM/QM approach "QUILD". In line with previous studies, most standard density functionals recover, at best, only part of the favorable stacking interactions. An exception is the new KT1 functional, which correctly yields bound π-stacked structures. Surprisingly, a similarly good performance is achieved with the computationally very robust and efficient local density approximation (LDA). Furthermore, we show that classical electrostatic interactions determine the shape and depth of the π-π stacking potential energy surface.FigureAdditivity approximation for the π-π interaction between two stacked Watson–Crick base pairs in terms of pairwise interactions between individual basesElectronic supplementary materialThe online version of this article (doi:10.1007/s00894-007-0239-y) contains supplementary material, which is available to authorized users.

Mixed-Ligand Copper(II)-phenolate Complexes: Effect of Coligand on Enhanced DNA and Protein Binding, DNA Cleavage, and Anticancer Activity

The copper(II) complex [Cu(tdp)(ClO4)].0.5H2O (1), where H(tdp) is the tetradentate ligand 2-[(2-(2-hydroxyethylamino)ethylimino)methyl]phenol, and the mixed ligand complexes [Cu(tdp)(diimine)]+ (2-5), where diimine is 2,2'-bipyridine (bpy) (2), 1,10-phenanthroline (phen) (3), 3,4,7,8-tetramethyl-1,10-phenanthroline (tmp) (4), and dipyrido-[3,2-d:2',3'-f]-quinoxaline (dpq) (5), have been isolated and characterized by analytical and spectral methods. Complexes 1 and [Cu(tdp)(phen)]ClO4 (3) have been structurally characterized, and their coordination geometries around copper(II) are described as distorted octahedral. The equatorially coordinated ethanolic oxygen in 1 is displaced to an axial position upon incorporating the strongly chelating phen, as in 3. The solution structures of all the complexes have been assessed to be square-based using electronic absorption and electron paramagnetic resonance (EPR) spectroscopy. The interaction of the complexes with calf thymus DNA (CT DNA) has been explored by using absorption, emission, and circular dichroic spectral and viscometric studies, and modes of DNA binding for the complexes have been proposed. Absorption spectral (Kb = 0.071 +/- 0.005 (2), 0.90 +/- 0.03 (3), 7.0 +/- 0.2 (4), 9.0 +/- 0.1 x 10(5) M(-1) (5)), emission spectral (Kapp = 4.6 (1), 7.8 (2), 10.0 (3), 12.5 (4), 25.0 x 10(5) M(-1) (5)), and viscosity measurements reveal that 5 interacts with DNA more strongly than the other complexes through partial intercalation of the extended planar ring of the coordinated dpq with the DNA base stack. Interestingly, only complex 4 causes a B to A conformational change upon binding DNA. All the complexes hydrolytically cleave pBR322 supercoiled DNA in 10% DMF/5 mM Tris-HCl/50 mM NaCl buffer at pH 7.1 in the absence of an activating agent, and the cleavage efficiency varies in the order 5 > 3 > 2 > 4 > 1 with 5 displaying the highest Kcat value (5.47 +/- 0.10 h(-1)). The same order of cleavage is observed for the oxidative cleavage of DNA in the presence of ascorbic acid as a reducing agent. Interestingly, of all the complexes, only 5 displays efficient photonuclease activity through double-strand DNA breaks upon irradiation with 365 nm light through a mechanistic pathway involving hydroxyl radicals. The protein binding ability of 1-5 has been also monitored by using the plasma protein bovine serum albumin (BSA), and 4 exhibits a protein binding higher than that of the other complexes. Further, the anticancer activity of the complexes on human cervical epidermoid carcinoma cell line (ME180) has been examined. Interestingly, the observed IC50 values reveal that complex 4, which effects conformational change on DNA and binds to BSA more strongly, exhibits a cytotoxicity higher than the other complexes. It also exhibits approximately 100 and 6 times more potency than cisplatin and mitomycin C for 24 and 48 h incubation times, respectively, suggesting that 4 can be explored further as a potential anticancer drug. Complexes 4 and 5 mediate the arrest of S and G2/M phases in the cell cycle progression at 24 h harvesting time, which progress into apoptosis.

Ultrafast Energy Delocalization and Electron Transfer Dynamics in 2‐Aminopurine‐Containing Trinucleotides†

The fate of electronically excited states in DNA base stacks is of tremendous importance for subsequent photochemical damage reactions in the genome. In this study we present a femtosecond broadband pump-probe study on the adenine isomer 2-aminopurine (Ap) incorporated into trinucleotides. After selective excitation of Ap we can monitor energy delocalization between neighboring Ap moieties as well as excited state electron transfer, depending on the sequence of the trinucleotide. Our results establish the time scale for intrastand excimer formation and reveal the lifetime of the excimer state.

Electronic coupling in the excited electronic state of stacked DNA base homodimers

The nature of the electronic coupling of stacked nucleic acid bases adenine (A), thymine (T), and cytosine (C), in A-A, T-T, and C-C complexes in their excited states was investigated; a different character of the electronic coupling for the T-T complex was shown.

The catalytic power of uracil DNA glycosylase in the opening of thymine base pairs.

Uracil DNA glycosylase (UNG) locates uracil and its structural congener thymine in the context of duplex DNA using a base flipping mechanism. NMR imino proton exchange measurements were performed on free and UNG-bound DNA duplexes in which a single thymine (T) was paired with a series of adenine analogues (X) capable of forming one, two, or three hydrogen bonds. The base pair opening equilibrium for the free DNA increased 55-fold as the number of hydrogen bonds decreased, but the opening rate constants were nearly the same in the absence and presence of UNG. In contrast, UNG was found to slow the base pair closing rate constants (kcl) compared to each free duplex by a factor of 3- to 23-fold. These findings indicate that regardless of the inherent thermodynamic stability of the TX pair, UNG does not alter the spontaneous opening rate. Instead, the enzyme holds the spontaneously expelled thymine (or uracil) in a transient extrahelical sieving site where it may partition forward into the enzyme active site (uracil) or back into the DNA base stack (thymine).