A new D-DNA form of poly(dA-dT).poly(dA-dT): an A-DNA type structure with reversed Hoogsteen pairing

The D-DNA double helix model of poly(dA-dT).poly(dA-dT) proposed in the literature is not in accordance with some notable experimental facts and physicochemical conditions to which it is related. Thus, the fibre X-ray diffraction pattern of D-DNA obtained at a relative humidity lower than that giving the A-DNA form is singularly not taken into account when one assumes that there is only one D structure of B-DNA type. We rather suggest that there are actually two different forms of D-DNA, namely D(A) which partakes in the D-A-B transitions and D(B) associated with the D-B change of conformation. Although these two DNA structures have the same helical parameters (pitch and number of residues per turn), in agreement with X-ray data, their detailed conformations are considerably different. Whereas D(B) is indeed the structure generally defined as D-DNA, a critical analysis based on a comparison between different possible DNA double helices leads us to propose dihedral angles, a set of atomic coordinates and a stereo view of another new form of D-DNA, the D(A) structural model. It is a right-handed double helix with a dinucleotide as the repeat unit. The furanose rings are of the A-DNA type (C3' endo) and the bases are hydrogen bonded according to the reversed Hoogsteen pairing. Such a disposition renders the D(A) model unsuitable for poly(dI-dC).poly(dI-dC), the other alternating polynucleotide observed in the D(B) structure. The consistency of these two different D-DNA structures of poly(dA-dT).poly(dA-dT) with the general aspects of hydration and helix-helix transitions of DNA, as well as with the conformational variability of AT base sequences, is discussed.     

Determination of the DNA-binding characteristics of ethidium bromide, proflavine, and cisplatin by flow injection analysis: usefulness in studies on antitumor drugs

Flow injection analysis was used to study the reactions occurring between DNA and certain compounds that bind to its double helix, deforming this and even breaking it, such that some of them (e.g., cisplatin) are endowed with antitumoral activity. Use of this technique in the merging zones and stopped-flow modes afforded data on the binding parameters and the kinetic characteristics of the process. The first compound studied was ethidium bromide (EtdBr), used as a fluorescent marker because its fluorescence is enhanced when it binds to DNA. The DNA-EtdBr binding parameters, the apparent intrinsic binding constant (0.31+/-0.02 microM(-1)), and the maximum number of binding sites per nucleotide (0.327+/-0.009) were determined. The modification introduced in these parameters by the presence of proflavine (Prf), a classic competitive inhibitor of the binding of EtdBr to the DNA double helix, was also studied, determining the value of the intrinsic binding constant of Prf (K(Prf) = 0.119+/-9x10(-3) microM(-1)). Finally, we determined the binding parameters between DNA and EtdBr in the presence of the antitumor agent cisplatin, a noncompetitive inhibitor of such binding. This provided information about the binding mechanism as well as the duration and activity of the binding of the compound in its pharmacological use.     

Phased psoralen cross-links do not bend the DNA double helix

Although the chemical reaction of psoralens with nucleic acids is well understood, the structure of psoralen-DNA cross-linked products is still not clear. Model building studies base on the crystal structure of the psoralen-thymine monoadduct suggest that each cross-link bends the DNA double helix by 46.5 degrees [Pearlman, D. A., Holbrook, S. R., Pirkle, D. H., & Kim, S.-H. (1985) Science (Washington, D.C.) 227, 1304-1308]. On the other hand, Sinden and Hagerman [Sinden, R. R., & Hagerman, P. J. (1984) Biochemistry 23, 6299-6303] find that, in solution, psoralen cross-linked DNA is not bent. Here we use gel electrophoresis to test the validity of the current models. We have synthesized a series of DNA fragments (21-24 base pairs in length), each containing one unique T-A site for 4'-(hydroxymethyl)-4,5',8-trimethylpsoralen (HMT) cross-linking. Because of an estimated 28 degrees unwinding of the helix by HMT [Wiesehahn, G., & Hearst, J. E. (1978) Proc. Natl. Acad. Sci. U.S.A. 75, 2703-2707], one expects that the 22-bp cross-linked fragment will be repeated nearly in phase with the average helical screw when multimerized. In that sequence ligation will maximally amplify any deformation to the double helix. We find that the ligated multimers of cross-linked DNA migrate close to the multimers of non-cross-linked DNA on polyacrylamide gels. Our observations place an upper limit of 10 degrees on DNA bending induced by psoralen cross-linking and indicate unwinding by about 1 bp, as well as stiffening of the double helix. These properties are not unexpected for classical intercalators.     

The role of molecular structure of sugar‐phosphate backbone and nucleic acid bases in the formation of single‐stranded and double‐stranded DNA structures

Our previous DFT computations of deoxydinucleoside monophosphate complexes with Na⁺‐ions (dDMPs) have demonstrated that the main characteristics of Watson‐Crick (WC) right‐handed duplex families are predefined in the local energy minima of dDMPs. In this work, we study the mechanisms of contribution of chemically monotonous sugar‐phosphate backbone and the bases into the double helix irregularity. Geometry optimization of sugar‐phosphate backbone produces energy minima matching the WC DNA conformations. Studying the conformational variability of dDMPs in response to sequence permutation, we found that simple replacement of bases in the previously fully optimized dDMPs, e.g. by constructing Pyr‐Pur from Pur‐Pyr, and Pur‐Pyr from Pyr‐Pur sequences, while retaining the backbone geometry, automatically produces the mutual base position characteristic of the target sequence. Based on that, we infer that the directionality and the preferable regions of the sugar‐phosphate torsions, combined with the difference of purines from pyrimidines in ring shape, determines the sequence dependence of the structure of WC DNA. No such sequence dependence exists in dDMPs corresponding to other DNA conformations (e.g., Z‐family and Hoogsteen duplexes). Unlike other duplexes, WC helix is unique by its ability to match the local energy minima of the free single strand to the preferable conformations of the duplex. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 640–650, 2014.     

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

Lambda phage cro repressor interaction with its operator DNA: 2'-deoxy-5-fluorouracil OR3 analogues

The experiments here show that chemically synthesized DNA containing fluorine at selected sites can be used to test specific predictions of a model for cro repressor--operator interaction. This is done by observation of the perturbation to the fluorine-19 NMR spectra of analogues of OR3 synthesized with 2'-deoxy-5-fluorouracil at specific positions in the DNA helix. Although the three-dimensional structure of the cro repressor from phage lambda has been determined by Matthews and co-workers [Anderson, W., Ohlendorf, D., Takeda, Y., & Matthews, B. (1981) Nature (London) 290, 754-758], direct structural observations on the complex of the protein with its specific DNA recognition sequence, OR3, are limited. From that structure of the protein, alone, a model of its complex to DNA was built by fitting B-form DNA, with some distortion [Ohlendorf, D., Anderson, W., Fisher, R., Takeda, Y., & Matthews, B. (1982) Nature (London) 298, 718-723]. That model proposes that the cro repressor contacts only one side of this DNA double helix and a number of specific protein--DNA contacts. To test the model, 2'-deoxy-5-fluorouracil was used to place the fluorine-19 nuclear spin-label on the side of the DNA contacting the cro repressor and on the opposite side facing away from the cro repressor. The results presented here are consistent with the prediction that lambda phage cro repressor contacts only one side of the DNA double helix.     

Structure of the Rad50 DNA double‐strand break repair protein in complex with DNA

The Mre11–Rad50 nuclease–ATPase is an evolutionarily conserved multifunctional DNA double‐strand break (DSB) repair factor. Mre11–Rad50's mechanism in the processing, tethering, and signaling of DSBs is unclear, in part because we lack a structural framework for its interaction with DNA in different functional states. We determined the crystal structure of Thermotoga maritima Rad50^NBD (nucleotide‐binding domain) in complex with Mre11^HLH (helix‐loop‐helix domain), AMPPNP, and double‐stranded DNA. DNA binds between both coiled‐coil domains of the Rad50 dimer with main interactions to a strand‐loop‐helix motif on the NBD. Our analysis suggests that this motif on Rad50 does not directly recognize DNA ends and binds internal sites on DNA. Functional studies reveal that DNA binding to Rad50 is not critical for DNA double‐strand break repair but is important for telomere maintenance. In summary, we provide a structural framework for DNA binding to Rad50 in the ATP‐bound state.     

Potency and property of DNA conformational damage determine ATM activation: two new DNA topological probes

We have shown previously the intercalation geometry of a series of acenaphtho [1,2-b] pyrrole derivatives with DNA double helix in vitro. In this report we chose a couple of intercalating analogues and a Chinese traditional medicine Tanshinone IIA as probes to investigate the response of DNA damage sensor ataxia-telangiectasia mutated (ATM) protein toward the DNA topological change in vivo. The two analogues (1)a (3-(4-Methyl-piperazin)-8-oxo-8H-acenaphtho [1,2-b]pyrrole -9-carbonitrile) and (3)a (3-(3-Dimethylamino-propylamino)-8-oxo-8H-acenaphtho[1,2-b]pyrrole-9- carbonitrile) could unwind double helix to different extents, whereas Tanshinone IIA could wind the double helix. Using a combination of circular dichroism (CD) studies and immunoflurescence assays, we found for the first time that the ATM protein kinase can respond to the unwinding chromatin conformational damage caused by (1)a and (3)a, while it could not be activated by the winding effects caused by Tanshinone IIA. Moreover, the amount of ATM protein phosphorylation is consistent with the degree of unwinding conformational damage. The average number of ATM foci in an MCF-7 cell is 32 +/- 1.5 at 6 microM (1)a, which is significantly higher than the 8 microM (3)a exposure (15 +/- 0.5, p < 0.5). A new couple of DNA topological probes, (1)a and (3)a have been found for the future semi-quantitative investigation of factors involved in the DNA damage pathway.     

Mitonuclear discordance in wolf spiders: Genomic evidence for species integrity and introgression

Systematists and taxonomists have benefited greatly from the emergence of molecular methods. Species identification has become straightforward through DNA barcoding and the rapid build‐up of massive DNA barcode reference libraries. In animals, mitonuclear discordance can significantly complicate the process of species identification and delimitation. The causes of mitonuclear discordance are either biological (e.g., introgression, incomplete lineage sorting, horizontal gene transfer androgenesis) or induced by operational factors (e.g., human error with specimen misidentification or incorrect species delimitation). Moreover, endosymbionts may play an important role in promoting fixation of mitochondrial genomes. Here, we study the mitonuclear discordance of wolf spiders species (Lycosidae) (independent cases from Alopecosa aculeata and Pardosa pullata groups) that share identical COI DNA barcodes. We approached the case utilizing double‐digest restriction site‐associated DNA sequencing (ddRADseq) to obtain and analyse genomic‐scale data. Our results suggest that the observed cases of mitonuclear discordance are not due to operational reasons but result from biological processes. Further analysis indicated introgression and that incomplete lineage sorting is unlikely to have been responsible for the observed discrepancy. Additional survey of endosymbionts provided ideas on further research and their role in shaping mitochondrial DNA distribution patterns. Thus, ddRADseq grants an efficient way to study the taxonomy of problematic groups with insight into underlying evolutionary processes.     

Prokaryotic DNA Methyltransferases: The Structure and the Mechanism of Interaction with DNA

The review considers current views on the function of DNA methyltransferases (MTases) that belong to prokaryotic type II restriction–modification systems. A commonly accepted classification of MTases is described along with their primary and tertiary structures and molecular mechanisms of their specific interaction with DNA (including methylation). MTase inhibitors are also considered. Special emphasis is placed on the flipping of the target heterocyclic base out of the double helix and on the methods employed in its analysis. Base flipping is a fundamentally new type of DNA conformational changes and is also of importance in the case of other DNA-operating enzymes. MTases show unique sequence homology, and are similar in structure of functional centers and in the mechanism of methylation. These data contribute to the understanding of the general biological significance of methylation, since prokaryotic and eukaryotic MTases are structurally and functionally similar.     

The structure of duplex DNA in the nucleosome core particle: An alternative model

Many studies indirectly indicate that the conformation ofin vivo duplex DNA is the double helix. The most direct view, from the X-ray analysis of the nucleosome core particle, has also been interpreted in terms of the double helix structure. However, an alternative possibility exists; that the duplex adopts a metastable side-by-side conformation which readily converts to the double helix on removal of protein. Evidence for the existence of this conformation has been obtained from a reanalysis of the electron density map for the nucleosome particle.     

Nonhistone proteins HMG1 and HMG2 unwind DNA double helix.

In a previous communication we have shown that both HMG1 and HMG2 nonhistone proteins change the DNA helical structure and the binding of HMG1 and HMG2 to DNA induces a net unwinding equivalent of DNA double helix (Javaherian, K., Liu, L. F. and Wang, J. C. (1978) Science, 199, 1345-1346). Employing melting absorption technique, we now show that in the presence of salt HMG1 and HMG2 destabilize DNA whereas in the absence of salt, they both stabilize DNA molecules. Consequently the folded structure of HMG must play an important role in melting DNA. Furthermore, by measuring topological winding number using competition unwinding experiments, we conclude that HMG1 has a higher affinity for a single-stranded DNA relative to double-stranded DNA. These results together suggest that HMG1 and HMG2 unwind DNA double helix by local denaturation of the DNA base pairs. The net unwinding angles have been measured to be 22 degrees and 26 degrees per molecule of HMG1 and HMG2 respectively.     

Strong sequence-dependent polymorphism in adduct-induced DNA structure: analysis of single N-2-acetylaminofluorene residues bound within the NarI mutation hot spot

We have used a set of chemical probes to characterize and to compare the structural deformation of double-stranded oligomers bearing a single N-2-acetylaminofluorene (AAF) adduct covalently bound to each of the three guanine residues located within the frameshift mutation hot spot sequence -G1G2CG3CC-(NarI site). Two classes of chemical probes have been used, probes that sense the geometry of the helix, giving rise to cuts at every nucleotide (for example, 1,10-phenanthroline-copper), and probes that react with specific bases depending on their conformation (e.g., diethyl pyrocarbonate). For all probes that were tested, a distinct pattern of reactivity was observed according to the position of the adduct within the DNA sequence, revealing an important polymorphism in the adduct-induced DNA structure. With 1,10-phenanthroline-copper at least three base pairs 3' of the AAF-modified guanine were reactive on each strand, showing that the deformation of the DNA helix extends over a region of 4-6 bases pairs centered around the adduct and sensed by the probe in both strands. With the base-specific probes, reactivities were limited to the base complementary to the modified guanine and to adjacent bases. Within this sequence context, the three possible AAF adducts have previously been shown to exhibit strong differences in biological responses such as excision repair [Seeberg, E., & Fuchs, R. P. P. (1990) Proc. Natl. Acad. Sci. U.S.A. 87, 191-194] and mutagenesis [Burnouf, D., Koehl, P., & Fuchs, R. P. P. (1989) Proc. Natl. Acad. Sci. U.S.A. 86, 4147-4151].(ABSTRACT TRUNCATED AT 250 WORDS)     

DNA binding and alkylation by the "left half" of azinomycin B

Azinomycin B (also known as carzinophilin A) contains two electrophilic functional groups-an epoxide and an aziridine residue-that react with nucleophilic sites in duplex DNA to form cross-links at 5'-dGNT and 5'-dGNC sequences. Although the aziridine residue of azinomycin is undoubtedly required for cross-link formation, analogues containing an intact epoxide group but no aziridine residue retain significant biological activity. Azinomycin epoxide analogues (e.g., 5 and 6) are of interest due to their potent biological activity and because there is evidence that azinomycin may decompose in vivo to yield such compounds. To investigate the chemical events underlying the toxicity of azinomycin epoxides, DNA binding and alkylation by synthetic analogues of azinomycin B (6, 8, and 9) that comprise the naphthalene-containing "left half" of the antibiotic have been investigated. The epoxide-containing analogue of azinomycin (6) efficiently alkylates guanosine residues in duplex DNA. DNA alkylation by 6 is facilitated by noncovalent binding of the compound to the double helix. The results of UV-vis absorbance, fluorescence spectroscopy, DNA winding, viscometry, and equilibrium dialysis experiments indicate that the naphthalene group of azinomycin binds to DNA via intercalation. Equilibrium dialysis experiments provide an estimated binding constant of (1.3 +/- 0.3) x 10(3) M(-)(1) for the association of a nonalkylating azinomycin analogue (9) with duplex DNA. The DNA-binding and alkylating properties of the azinomycin epoxide 6 provide a basis for understanding the cytotoxicity of azinomycin analogues which contain an epoxide residue but no aziridine group and may provide insight into the mechanisms by which azinomycin forms interstrand DNA cross-links.     

Melting of Cross-Linked DNA V. Cross-Linking Effect Caused by Local Stabilization of the Double Helix

Abstract

DNA interstrand cross-links are usually formed due to bidentate covalent or coordination binding of a cross-linking agent to nucleotides of different strands. However interstrand linkages can be also caused by any type of chemical modification that gives rise to a strong local stabilization of the double helix. These stabilized sites conserve their helical structure and prevent local and total strand separation at temperatures above the melting of ordinary AT and GC base pairs. This local stabilization makes DNA melting fully reversible and independent of strand concentration like ordinary covalent interstrand cross-links. The stabilization can be caused by all the types of chemical modifications (interstrand cross-links, intrastrand cross-links or monofunctional adducts) if they give rise to a strong enough local stabilization of the double helix. Our calculation demonstrates that an increase in stability by 25 to 30 kcal in the free energy of a single base pair of the double helix is sufficient for this “cross-linking effect” (i.e. conserving the helicity of this base pair and preventing strand separation after melting of ordinary base pairs). For the situation where there is more then one stabilized site in a DNA duplex (e.g., 1 stabilized site per 1000 bp), a lower stabilization per site is sufficient for the “cross-linking effect” (18–20 kcal). A substantial increase in DNA stability was found in various experimental studies for some metal-based anti-tumor compounds. These compounds may give rise to the effect described above. If ligand induced stabilization is distributed among several neighboring base pairs, a much lower minimum increase per stabilized base pair is sufficient to produce the cross-linking effect (1 bp- 24.4 kcal; 5 bp- 5.3 kcal; 10 bp- 2.9 kcal, 25 bp- 1.4 kcal; 50 bp- 1.0 kcal). The relatively weak non-covalent binding of histones or protamines that cover long regions of DNA (20–40 bp) can also cause this effect if the salt concentration of the solution is sufficiently low to cause strong local stabilization of the double helix. Stretches of GC pairs more than 25 bp in length inserted into poly(AT) DNA also exhibit properties of stabilizing interstrand cross-links.     

Conformational and chiral selection of oligonucleotides

In view of a better understanding of chiral selection of oligonucleotides, we have studied the hybridization of D- and L-CNA (cyclohexane nucleic acids) and D- and L-DNA, with chiral D-beta-homo-DNA and achiral PNA (peptide nucleic acids). PNA hybridizes as well with D-DNA, L-DNA as with D-beta-homo-DNA. The structure of the PNA x D-beta-homo-DNA complex is different from the PNA x DNA duplexes. D-CNA prefers D-DNA as hybridization partner, while L-CNA prefers D-beta-homo-DNA as hybridization partner. The conformation of the enantiomeric oligonucleotides D-CNA and L-CNA in the supramolecular complex with D-DNA and D-beta-homo-DNA, respectively, is different. These data may contribute to the confirmation of a hypothesis of the existence of achiral informative polymers as RNA predecessor, and to the understanding of homochirality of nucleic acids.     

Sequence-dependent anisotropic flexibility of B-DNA. A conformational study

Bending flexibility of the six tetrameric duplexes was investigated d(AAAA):d(TTTT), d(AATT)2, d(TTAA)2, d(GGGG):d(CCCC), d(GGCC)2 and d(CCGG)2,. The tetramers were extended in the both directions by regular double helices. The stiffness of the B-DNA double helix when bent into the both grooves proved to be less than that in the perpendicular direction by an order of magnitude. Such an anisotropy is a property of the sugar-phosphate backbone structure. The calculated fluctuations of the DNA bending along the dyad axis, 5-7 degree, are in agreement with experimental value of the DNA persistence length. Anisotropy of the double helix is sequence-dependent: most easily bent into the minor groove are the tetramers with purine-pyrimidine dimer (RY) in the middle. In contrast, YR dinucleotides prefer bending into the major groove. Moreover, they have an equilibrium bend of 6-12 degree into this groove. The above inequality is caused by stacking interaction of the bases. The bend in the central dimer is distributed to some extent between the adjacent links, though the main fraction of the bend remains within the central link. Variation of the sugar-phosphate geometry in the bent helix is inessential, so that DNA remains within the B-family of forms: namely, when the helical axis is bent by 20 degree. the backbone dihedral angles vary by no more than 15 degree. The obtained results are in accord with x-ray structure of the B-DNA dodecamer; they further substantiate our early model of DNA wrapping in the nucleosome by means of "mini-kinks" separated by a half-pitch of the double helix, i.e. by 5-6 b.p. Sequence-dependent anisotropy of DNA presumably dictates the three-dimensional structure of DNA in solution as well. We have found that nonrandom allocation of YR dimers leads to the systematic bends in equilibrium structure of certain DNA fragments.     

基于三链核酸的DNA计算

一种研究DNA计算的新模型——三链DNA计算模型在本文中提出。此模型是在近年三链核酸的研究成果的基础上建立的。并应用于求解可满足性问题（SAT）,这是一个困难的NP-完全问题。不同于以住的DNA计算方法,基于三链核酸的分子算法通过寡聚脱氧核苷酸（ODN）在RecA蛋白的介导下与同源的双链DNA匹配成三链DNA进行基本的运算,这样可以有效减少计算中的错误。依据分子生物学的实验方法,该算法切实可行并且有效。     

Quantitative prediction of 3D solution shape and flexibility of nucleic acid nanostructures

DNA nanotechnology enables the programmed synthesis of intricate nanometer-scale structures for diverse applications in materials and biological science. Precise control over the 3D solution shape and mechanical flexibility of target designs is important to achieve desired functionality. Because experimental validation of designed nanostructures is time-consuming and cost-intensive, predictive physical models of nanostructure shape and flexibility have the capacity to enhance dramatically the design process. Here, we significantly extend and experimentally validate a computational modeling framework for DNA origami previously presented as CanDo [Castro,C.E., Kilchherr,F., Kim,D.-N., Shiao,E.L., Wauer,T., Wortmann,P., Bathe,M., Dietz,H. (2011) A primer to scaffolded DNA origami. Nat. Meth., 8, 221-229.]. 3D solution shape and flexibility are predicted from basepair connectivity maps now accounting for nicks in the DNA double helix, entropic elasticity of single-stranded DNA, and distant crossovers required to model wireframe structures, in addition to previous modeling (Castro,C.E., et al.) that accounted only for the canonical twist, bend and stretch stiffness of double-helical DNA domains. Systematic experimental validation of nanostructure flexibility mediated by internal crossover density probed using a 32-helix DNA bundle demonstrates for the first time that our model not only predicts the 3D solution shape of complex DNA nanostructures but also their mechanical flexibility. Thus, our model represents an important advance in the quantitative understanding of DNA-based nanostructure shape and flexibility, and we anticipate that this model will increase significantly the number and variety of synthetic nanostructures designed using nucleic acids.