Preferred positions of AA and TT dinucleotides in aligned nucleosomal DNA sequences.

Multiple alignment of 118 nucleosomal DNA sequences by maximizing simultaneously match of AA dinucleotides and match of TT dinucleotides results in a pattern of the dinucleotide distributions which is characteristic of the nucleosomal DNA sequences. The AA dinucleotides are found to be distributed symmetrically relative to the TT dinucleotide distribution, around the middle point of the nucleosomal DNA sequence. The distances between major peaks of the distributions are multiples of about 10.4 bases. The peaks of the TT distribution are shifted by 6 bases downstream from the peaks of the AA distribution.     

The relationship between periodic dinucleotides and the nucleosomal DNA deformation revealed by normal mode analysis

Nucleosomes, which contain DNA and proteins, are the basic unit of eukaryotic chromatins. Polymers such as DNA and proteins are dynamic, and their conformational changes can lead to functional changes. Periodic dinucleotide patterns exist in nucleosomal DNA chains and play an important role in the nucleosome structure. In this paper, we use normal mode analysis to detect significant structural deformations of nucleosomal DNA and investigate the relationship between periodic dinucleotides and DNA motions. We have found that periodic dinucleotides are usually located at the peaks or valleys of DNA and protein motions, revealing that they dominate the nucleosome dynamics. Also, a specific dinucleotide pattern CA/TG appears most frequently.     

How many numbers are required to specify sequence-dependent properties of polynucleotides?

There are 10 unique dinucleotides of double-stranded DNA, but only 8 independent nearest-neighbor energies that occur in circular DNA, as shown by D. M. Gray and I. Tinoco [(1970) Biopolymers 9, 223-244]. We extend that analysis to include end effects, and show that the number of unique dinucleotide pairs (including ends) is 14, but there are only 12 independent energies. We discuss how these 12 energies (or spectra or any other pairwise additive property) can be measured and displayed, and how they should and should not be compared between experimenters. As an example, we analyzed the recently reported melting curves [M.J. Doktycz et al. (1992) Biopolymers, 32, 849-864.] of 16 DNA dumbbells in two different Na+ environments. This analysis reveals a new means for evaluating end effects and the emergence of longer than nearest-neighbor interactions at low salt concentration.     

Preparation and properties of 3-halopyridine--adenine dinucleotides, NAD+ analogues and model compounds.

The preparation of model compounds 1-(2',6'-dichlorobenzyl)-3-halogenopyridinium and the study of their properties were achieved. Their chemical reduction to the corresponding 1,4-dihydropyridines is proved by spectroscopic analysis. 3-Iodopyridine--adenine dinucleotide was prepared by enzymic transglycosidation while the 3-chloro, 3-bromo and 3-iodo pyridine--adenine dinucleotides were synthesized from 3-amino-pyridine--adenine dinucleotide. The 3-halogenopyridine--adenine dinucleotides were proved to be active as hydrogen acceptors with alcohol as a substrate. The absorption band at 290 nm of cinnamaldehyde appeared to be a very sensitive tool for studying the enzymic reaction. With the alcohol dehydrogenase from yeast, only slight activity was detected. 3-Halogenopyridine--adenine dinucleotides are competitive inhibitors with respect to nicotinamide--adenine dinucleotide with alcohol dehydrogenase from yeast, lactate dehydrogenase and malate dehydrogenase. The use of 3-iodopyridine--adenine dinucleotide as a heavy-atom derivative for X-ray structure determination is proposed.     

Theoretical studies of cis-Pt(II)-diammine binding to duplex DNA

The binding of cis-Pt(II) diammine (cis-DP) to double-stranded DNA was studied with several kinked conformations that can accommodate the formation of a square planar complex. Molecular mechanics (MM) calculations were performed to optimize the molecular fit. These results were combined with quantum mechanical (QM) calculations to ascertain the relative energetics of ligand binding through water vs direct binding of the phosphate to the ammine and platinum, and to guide the selection of DNA conformations to model complex formation. Based on QM and MM calculations, models are proposed that may be characterized by several general features. A structure involving hydrogen bonding between each ammine and distinct adjacent phosphate groups, referred to as closed conformation (CC), has already been reported. This is also found in the crystal structure of small dimers. We report alternative conformations that may be important in platination of duplex DNA. They are characterized by an intermediate conformation (IC), involving hydrogen bonding between one ammine and phosphate group, and an open conformation (OC), without ammine phosphate hydrogen bonding. The IC and OC can be stabilized by water bridges in the space between the ammine and the phosphate groups. Sugar puckers alternate from the type C(2')-endo or C(1')-exo (S), to the type C(3')-endo or C(2')-exo (N), with intermediate types near O(1')-endo (O). In general, the sugar puckers alternate from S to N to S through the platinated region (3'-TpG*pG*p-5'), with the complexed strand exhibiting, (3')-S*-N*-S-(5') alternation, while the complementary strand shows either (3')-S*-N*-S-(5') or (3')-S*-N*-O-(5') alternation. In both the OC and IC, a hydrogen bond is found between the ammine and O4(T) on thymine (T) at the (3') end, adjacent to the complex site. There is a continuous range of backbone conformations through the platinated region which relate the OC to the IC. The models presented suggest that the dynamics of the binding of the cis-Pt(II)-diammines to adjacent N7(G) in double-stranded DNA may encompass several conformational possibilities, and that water bridges may play a roll in supporting open and intermediate conformations. Proton-proton distances are reported to assist in the experimental determination of conformations.     

Interaction between nucleic acids and glycans: In silico and in vitro

Quantum-chemical (PM3) calculations prove the possibility of hydrogen bonding between all nucleic bases and carboxyls of acidic or hydroxymethyls of neutral sugars, with energy in the former case comparable to that for canonical GC or AT pairing. However, there is appreciable energy preference for carboxyl H-bonding with purines, and for hydroxymethyl H-bonding with pyrimidines. Simulation reveals that the H-bonds formed in purine-uronide and pyrimidine-hexose pairs can give rise to double-stranded nucleic acid-polysaccharide complexes. Indeed, dot hybridization with radioactive probes and hypochromic effects in solution testify that a hexose homoglycan (amylose) selectively forms a complex with polypyrimidine, whereas uronic acid homoglycans selectively complex with polypurine. Complexing is also observed between heteropolymers such as thymus DNA and hyaluronic acid. These results are consistent with the idea of template-directed synthesis of the polysaccharide moiety of glycosaminoglycans involving nucleic acids.     

Purification,characterization and docking studies of the HIN domain of human myeloid nuclear differentiation antigen (MNDA)

The HIN domain of myeloid nuclear differentiation antigen (MNDA) was expressed and purified as a monomer using E. coli JM109 as host. The protein interacted with double-stranded DNA at a K_d of 3.15 μM and did not recognize the termini of double-stranded DNA. Isothermal titration calorimetry indicated that the interaction between the protein and double-stranded DNA is mainly mediated by electrostatic attractions and hydrogen bonding. We developed a model to analyze the potential DNA binding site of the MNDA HIN domain. Based on the model, molecular docking and mutation studies suggest that the double-stranded DNA binding site of the protein is different from other HIN–DNA structures. This work facilitates the design of specific drugs against pathogens detected by human MNDA.     

The mechanism of the supercoiled DNA recognition by eukaryotic type I topoisomerases. I. The enzyme interaction with nonspecific oligonucleotides

Interaction of the DNA type I topoisomerases from the murine and human placenta cells with nonspecific oligonucleotides was analyzed. The contributions of strong and week nonspecific electrostatic, van der Waals's, and hydrophobic interactions, and hydrogen bonding of the enzymes to the complex formation with the single- and double-stranded DNAs were determined. The factors that determine the top-priority recognition of the topologically stressed DNA were revealed. The results were interpreted in comparison with the X-ray analysis data for human DNA topoisomerase I.     

The Mechanism of the Supercoiled DNA Recognition by the Eukaryotic Type I Topoisomerases: I. The Enzyme Interaction with Nonspecific Oligonucleotides

Interaction of the DNA type I topoisomerases from the murine and human placenta cells with nonspecific oligonucleotides was analyzed. The contributions of strong and week nonspecific electrostatic, van der Waals's, and hydrophobic interactions, and hydrogen bonding of the enzymes to the complex formation with the single- and double-stranded DNAs were determined. The factors that determine the top-priority recognition of the topologically stressed DNA were revealed. The results were interpreted in comparison with the X-ray analysis data for human DNA topoisomerase I.     

Sequence-dependent contribution of distal binding domains to CAP protein-DNA binding affinity.

We report measurements of the relative binding affinity of CAP for DNA sequences which have been systematically mutated in the region flanking the consensus binding site. Our experiments focus on the locus one helical turn from the dyad axis where DNA bending toward the minor groove is induced upon C-AP binding. The binding free energy and extent of bending are moderately well correlated for the set of 56 sequences. Changes in binding affinity spanning a factor of about 50 could be accounted for by additive contributions of dinucleotides; with a few exceptions, the relative ranking of dinucleotide contributions to binding and bending are similar. We conclude that dinucleotides are the smallest independent unit required for quantitative interpretation of CAP-induced DNA bending and binding in the distal domains of the CAP consensus binding site. The imperfect correlation between binding strength and extent of bending implies that sequence changes affect protein binding strength not only by altering the DNA deformation energy required to form the complex, but also by affecting directly the free energy of interaction between protein and DNA.     

Characterization of inherent curvature in DNA lacking polyadenine runs.

Sequence-directed DNA curvature is most commonly associated with AA dinucleotides in the form of polyadenine runs. We demonstrate inherent curvature in DNA which lacks AA/TT dinucleotides using the criteria of polyacrylamide gel mobility and efficiency of DNA cyclization. These studies are based upon two 21-base pair synthetic DNA fragments designed to exhibit fixed curvature according to deflections made to the helical axis by non-AA dinucleotide stacks. Repeats of these sequences display anomalously slow migration in polyacrylamide gels. Moreover, both sequences describe helical conformations that are closed into circles by DNA ligase at much smaller sizes than is typical of nondeformed DNA. Chemical cleavage of these DNA molecules with hydroxyl radical is also consistent with local variation in helical conformation at specific dinucleotide steps.     

A Streptomyces glaucescens endodeoxyribonuclease which shows a strong preference for CC dinucleotide.

We have characterized a deoxyribonuclease from Streptomyces glaucescens that cleaves double-stranded DNA preferably between the dinucleotide 5'-CC-3'. The cleavage specificity was demonstrated by both analysis of the terminal nucleotides of the generated fragments and DNA sequencing of partially digested DNA. Digestion of lambda DNA with this enzyme resulted in the production of double-stranded fragments with 5' and/or 3'-protruding single-stranded tails. DNase I footprinting experiments indicated that the nuclease specifically binds to its cleavage sites on the DNA under non-catalytic conditions. The enzyme is not affected by cytosine methylation in hemimethylated DNA.     

Recognition of methylated DNA through methyl-CpG binding domain proteins

DNA methylation is a key regulatory control route in epigenetics, involving gene silencing and chromosome inactivation. It has been recognized that methyl-CpG binding domain (MBD) proteins play an important role in interpreting the genetic information encoded by methylated DNA (mDNA). Although the function of MBD proteins has attracted considerable attention and is well characterized, the mechanism underlying mDNA recognition by MBD proteins is still poorly understood. In this article, we demonstrate that the methyl-CpG dinucleotides are recognized at the MBD-mDNA interface by two MBD arginines through an interplay of hydrogen bonding and cation-π interaction. Through molecular dynamics and quantum-chemistry calculations we investigate the methyl-cytosine recognition process and demonstrate that methylation enhances MBD-mDNA binding by increasing the hydrophobic interfacial area and by strengthening the interaction between mDNA and MBD proteins. Free-energy perturbation calculations also show that methylation yields favorable contribution to the binding free energy for MBD-mDNA complex.     

Repertoires of the Nucleosome-Positioning Dinucleotides

It is generally accepted that the organization of eukaryotic DNA into chromatin is strongly governed by a code inherent in the genomic DNA sequence. This code, as well as other codes, is superposed on the triplets coding for amino acids. The history of the chromatin code started three decades ago with the discovery of the periodic appearance of certain dinucleotides, with AA/TT and RR/YY giving the strongest signals, all with a period of 10.4 bases. Every base-pair stack in the DNA duplex has specific deformation properties, thus favoring DNA bending in a specific direction. The appearance of the corresponding dinucleotide at the distance 10.4 xn bases will facilitate DNA bending in that direction, which corresponds to the minimum energy of DNA folding in the nucleosome. We have analyzed the periodic appearances of all 16 dinucleotides in the genomes of thirteen different eukaryotic organisms. Our data show that a large variety of dinucleotides (if not all) are, apparently, contributing to the nucleosome positioning code. The choice of the periodical dinucleotides differs considerably from one organism to another. Among other 10.4 base periodicities, a strong and very regular 10.4 base signal was observed for CG dinucleotides in the genome of the honey bee A. mellifera. Also, the dinucleotide CG appears as the only periodical component in the human genome. This observation seems especially relevant since CpG methylation is well known to modulate chromatin packing and regularity. Thus, the selection of the dinucleotides contributing to the chromatin code is species specific, and may differ from region to region, depending on the sequence context.     

Differences between sites of binding to DNA and strand cleavage for complexes of bleomycin with iron or cobalt

The sequence specificity of bleomycin A5 and of its light-activated cobalt complex were compared by examining the relative cleavage of each strand of two DNA fragments by either species. Significant differences between the two metallobleomycins were observed. The iron-bleomycin (Fe-BLM) complex cleaved the DNA molecules preferentially at dinucleotides GpT and GpC, whereas the light-activated cobalt-bleomycin complex (Co-BLM) showed a preference for cutting at the dinucleotide GpA in addition to cleavage at every GpT dinucleotide. Further, new sites of preferential cleavage were noted for Co-BLM in regions of the DNA where enhanced reaction with DNAaseI can be observed in the presence of the antibiotic. No differences in the cutting behaviour of the Fe-BLM were evident upon irradiation of the reaction mixture. A reduction in the relative efficiency of cutting at GpC sequences by Co-BLM is responsible for the previously observed diminution of double-strand breaks under conditions of photoactivated cleavage. The results are discussed in terms of the likely production of highly reactive, diffusible cutting elements in the light activated reaction which cause cleavage of the DNA in regions where the antibiotic is not bound.     

CpG Dinucleotide Frequencies Reveal the Role of Host Methylation Capabilities in Parvovirus Evolution

Parvoviruses are rapidly evolving viruses that infect a wide range of hosts, including vertebrates and invertebrates. Extensive methylation of the parvovirus genome has been recently demonstrated. A global pattern of methylation of CpG dinucleotides is seen in vertebrate genomes, compared to “fractional” methylation patterns in invertebrate genomes. It remains unknown if the loss of CpG dinucleotides occurs in all viruses of a given DNA virus family that infect host species spanning across vertebrates and invertebrates. We investigated the link between the extent of CpG dinucleotide depletion among autonomous parvoviruses and the evolutionary lineage of the infected host. We demonstrate major differences in the relative abundance of CpG dinucleotides among autonomous parvoviruses which share similar genome organization and common ancestry, depending on the infected host species. Parvoviruses infecting vertebrate hosts had significantly lower relative abundance of CpG dinucleotides than parvoviruses infecting invertebrate hosts. The strong correlation of CpG dinucleotide depletion with the gain in TpG/CpA dinucleotides and the loss of TpA dinucleotides among parvoviruses suggests a major role for CpG methylation in the evolution of parvoviruses. Our data present evidence that links the relative abundance of CpG dinucleotides in parvoviruses to the methylation capabilities of the infected host. In sum, our findings support a novel perspective of host-driven evolution among autonomous parvoviruses.     

Characteristics of nucleosome core DNA and their applications in predicting nucleosome positions

By analyzing dinucleotide position-frequency data of yeast nucleosome-bound DNA sequences, dinucleotide periodicities of core DNA sequences were investigated. Within frequency domains, weakly bound dinucleotides (AA, AT, and the combinations AA-TT-TA and AA-TT-TA-AT) present doublet peaks in a periodicity range of 10-11 bp, and strongly bound dinucleotides present a single peak. A time-frequency analysis, based on wavelet transformation, indicated that weakly bound dinucleotides of core DNA sequences were spaced smaller (∼10.3 bp) at the two ends, with larger (∼11.1 bp) spacing in the middle section. The finding was supported by DNA curvature and was prevalent in all core DNA sequences. Therefore, three approaches were developed to predict nucleosome positions. After analyzing a 2200-bp DNA sequence, results indicated that the predictions were feasible; areas near protein-DNA binding sites resulted in periodicity profiles with irregular signals. The effects of five dinucleotide patterns were evaluated, indicating that the AA-TT pattern exhibited better performance. A chromosome-scale prediction demonstrated that periodicity profiles perform better than previously described, with up to 59% accuracy. Based on predictions, nucleosome distributions near the beginning and end of open reading frames were analyzed. Results indicated that the majority of open reading frames’ start and end sites were occupied by nucleosomes.     

The majority of methylated deoxycytidines in human DNA are not in the CpG dinucleotide

The only natural postsynthetic modification known to occur in mammalian DNA is the methylation in the 5 position of deoxycytidines. Of the four 5'-CpN-3' dinucleotides (ie. CpG, CpC, CpA, and CpT), the dinucleotide which contains the highest proportion of deoxycytidines methylated is CpG, with 40 to 80% methylation in different mammalian genomes. It has also been shown that CpA, CpT, and CpC are methylated as well but to a much lower extent. Here we report the result of a full nearest neighbour analysis (together with quantitation of methylation levels in the 4 CpN dinucleotides) for DNA from human spleen. Using the values we have calculated the overall frequencies for all the methylated dinucleotides in the human genome. Because of the relative underrepresentation (by 7 to 10 fold) of the CpG dinucleotide, only 45.5% of total mC was present in mCpG, with 54.5% in mCpA, mCpT plus mCpC. These calculations have implications for studies into the function and significance of DNA methylation in mammalian cells.     

In vitro selection of DNAs with an increased propensity to form small circles

A protocol was devised to select for DNA molecules that efficiently form circles from a library of 126 base pair DNAs containing 90 randomized base pairs. After six rounds of selection, individual molecules from the library showed 20‐ to 100‐fold greater j‐factors compared with the starting library, validating the selection protocol. High‐throughput sequencing revealed a sinusoidal pattern of enrichment and de‐enrichment of A/T dinucleotides in the random region with a 10.4 base pair period associated with the helicity of DNA. A similar, but more moderate pattern of C/G dinucleotides was offset by precisely half a helical turn. While C/G dinucleotide enrichments were evenly distributed, A/T dinucleotide enrichments displayed a preference to cluster in individual DNA molecules. The most highly enriched 10 base pair sequences in the random region contained adjacent blocks of A/T and C/G trinucleotides present in some, but not all, rapidly cyclizing molecules. The phased dinucleotide enrichments closely match those present in accurately mapped yeast nucleosomes, confirming the importance of DNA bending in nucleosome formation. However, at certain sites the nucleosomal DNAs show dinucleotide enrichments that differ substantially from the cyclization data. These discrepancies can often be correlated with sequence specific contacts that form between histones and DNA. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 303–320, 2015.     

Deoxyribonucleic acid polymerase III of Escherichia coli. Characterization of associated exonuclease activities.

Purified DNA polymerase III has two distinct exonuclease activities: one initiates hydrolsis at the 3 termini, and the other at the 5 termini of single-stranded DNA. Both exonucleases have the same relative mobility on polyacrylamide gels as the polymerase activity. Molecular identity of the three activities is further indicated by their comparative rates of thermal inactivation and their sensitivity to ionic strength. The 3-5 exonuclease activity hydrolyzes only single-standed DNA. The rate of hydrolysis is twice the optimal rate of polymerization. The products are 5-mononucleotides, but the 3-5 activity is unable to cleave free dinucleotides or the 5-terminal dinucleotide of a polydeoxynucleotide chain. The 3-5 activity will not degrade 3-phosphoryl-terminated oligonucleotides such as d(pTpTpTp). The 5-3 activity catalyzes the hydrolysis of single-stranded DNA at 1/15 the rate of the 3-5 exonuclease. The 5-3 exonuclease requires the presence of a 5 single-stranded terminus in order to initiate hydrolysis, but will thereafter proceed into a double-stranded region. Although the limit products found during hydrolysis of substrates designed to assay specifically the 5-3 activity are predominantly mono- and dinucleotides, these products probably arise from the subsequent hydrolysis of oligonucleotides by the 3-5 hydrolytic activity. This interpretation is supported by (a) the relatively greater activity of the 3-5 exonuclease, (b) the inability of the enzyme to degrade d(pTpTpTp), and (c) the release of the 5 terminus of a single-stranded DNA molecule as an oligonucleotide. The 5-3 exonuclease attacks ultraviolet-irradiated duplex DNA which has first been incised by the Micrococcus luteus endonuclease specific for thymine dimers in DNA.