Crystal structure analysis of the B-DNA dodecamer CGTGAATTCACG.

The crystal structure of the DNA dodecamer C-G-T-G-A-A-T-T-C-A-C-G has been determined at a resolution of 2.5 A, with a final R factor of 15.8% for 1475 nonzero reflections measured at 0 degrees C. The structure is isomorphous with that of the Drew dodecamer, with the space group P2(1)2(1)2(1) and cell dimensions of a = 24.94 A, b = 40.78 A, and c = 66.13 A. The asymmetric unit contains all 12 base pairs of the B-DNA double helix and 36 water molecules. The structure of C-G-T-G-A-A-T-T-C-A-C-G is very similar to that of C-G-C-G-A-A-T-T-C-G-C-G, with no major alterations in helix parameters. Water peaks in the refined structure appear to represent a selection of peaks that were observed in the Drew dodecamer. The minor-groove spine of hydration at 2.5 A is fragmentary, but as Narendra et al. (1991) [Biochemistry (following paper in this issue)] have observed, lowering the temperature leads to a more complete representation of the spine.     

Placental 5-methylcytosine and 5-hydroxymethylcytosine patterns associate with size at birth

Altered placental function as a consequence of aberrant imprinted gene expression may be one mechanism mediating the association between low birth weight and increased cardiometabolic disease risk. Imprinted gene expression is regulated by epigenetic mechanisms, particularly DNA methylation (5mC) at differentially methylated regions (DMRs). While 5-hydroxymethylcytosine (5hmC) is also present at DMRs, many techniques do not distinguish between 5mC and 5hmC. Using human placental samples, we show that the expression of the imprinted gene CDKN1C associates with birth weight. Using specific techniques to map 5mC and 5hmC at DMRs controlling the expression of CDKN1C and the imprinted gene IGF2, we show that 5mC enrichment at KvDMR and DMR0, and 5hmC enrichment within the H19 gene body, associate positively with birth weight. Importantly, the presence of 5hmC at imprinted DMRs may complicate the interpretation of DNA methylation studies in placenta; future studies should consider using techniques that distinguish between, and permit quantification of, both modifications.     

Crystal structure of a B-DNA dodecamer containing inosine, d(CGCIAATTCGCG), at 2.4 A resolution and its comparison with other B-DNA dodecamers.

The crystal structure of the dodecamer, d(CGCIAATTCGCG), has been determined at 2.4 A resolution by molecular replacement, and refined to an R-factor of 0.174. The structure is isomorphous with that of the B-DNA dodecamer, d(CGCGAATTCGCG), in space group P2(1)2(1)2(1) with cell dimensions of a = 24.9, b = 40.4, and c = 66.4 A. The initial difference Fourier maps clearly indicated the presence of inosine instead of guanine. The structure was refined with 44 water molecules, and compared to the parent dodecamer. Overall the two structures are very similar, and the I:C forms Watson-Crick base pairs with similar hydrogen bond geometry to the G:C base pairs. The propeller twist angle is low for I4:C21 and relatively high for the I16:C9 base pair (-3.2 degrees compared to -23.0 degrees), and the buckle angles alter, probably due to differences in the contacts with symmetry related molecules in the crystal lattice. The central base pairs of d(CGCIAATTCGCG) show the large propeller twist angles, and the narrow minor groove that characterize A-tract DNA, although I:C base pairs cannot form the major groove bifurcated hydrogen bonds that are possible for A:T base pairs.     

Structural features of B-DNA dodecamer crystal structures: influence of crystal packing versus base sequence

We have analyzed the set of inter and intra base pair parameters for each dinucleotide step in single crystal structures of dodecamers, solved at high and medium resolution and all crystallized in P2(1)2(1)2(1) space group. The objective was to identify whether all the structures which have either the Drew-Dickerson (DD) sequence d[CGCGAATTCGCG] with some base modification or related sequence (non-DD), would display the same sequence dependent structural variability about its palindromic sequence, despite the molecule being bent at one end because of similar crystal lattice packing effect. Most of the local doublet parameters for base pairs steps G2-C3 and G10-C11 positions, symmetrically situated about the lateral two-fold, were significantly correlated between themselves. In non-DD sequences, significant correlations between these positional parameters were absent. The different range of local step parameter values at each sequence position contributed to the gross feature of smooth helix axis bending in all structures. The base pair parameters in some of the positions, for medium resolution DD sequence, were quite unlike the high-resolution set and encompassed a higher range of values. Twist and slide are the two main parameters that show wider conformational range for the middle region of non-DD sequence structures in comparison to DD sequence structures. On the contrary, the minor and major groove features bear good resemblance between DD and non-DD sequence crystal structure datasets. The sugar-phosphate backbone torsion angles are similar in all structures, in sharp contrast to base pair parameter variation for high and low resolution DD and non-DD sequence structures, consisting of unusual (epsilon = g-, xi = t) BII conformation at the 10th position of the dodecamer sequence. Thus examining DD and non-DD sequence structures packed in the same crystal lattice arrangement, we infer that inter and intra base pair parameters are as symmetrically equivalent in its value as the symmetry related step for the palindromic DD sequence about lateral two-fold axis. This feature would lead us to agree with the conclusion that DNA conformation is not substantially affected by end-to-end or lateral inter-molecular interaction due to crystal lattice packing effect. Non-DD sequence structures acquire step parameter values which reflect the altered sequence at each of the dodecamer sequence position in the orthorhombic lattice while showing similar gross features of DD sequence structures.     

Identifying 5-methylcytosine and related modifications in DNA genomes.

Intense interest in the biological roles of DNA methylation, particularly in eukaryotes, has produced at least eight different methods for identifying 5-methylcytosine and related modifications in DNA genomes. However, the utility of each method depends not only on its simplicity but on its specificity, resolution, sensitivity and potential artifacts. Since these parameters affect the interpretation of data, they should be considered in any application. Therefore, we have outlined the principles and applications of each method, quantitatively evaluated their specificity,resolution and sensitivity, identified potential artifacts and suggested solutions, and discussed a paradox in the distribution of m5C in mammalian genomes that illustrates how methodological limitations can affect interpretation of data. Hopefully, the information and analysis provided here will guide new investigators entering this exciting field.     

Interaction of AluI, Cfr6I and PvuII restriction-modification enzymes with substrates containing either N4-methylcytosine or 5-methylcytosine

The cleavage specificity of R.Cfr6I, an isoschizomer of PvuII restriction endonuclease was determined to be 5'CAG decreases CTG and the methylation specificity of Cfr6I and PvuII methylases, 5'CAG4mCTG. Thus, M.Cfr6I and M.PvuII are new additions to the list of methylases with N4-methylcytosine specificity. Neither of the above RM enzymes acts on the substrates containing either N4-methylcytosine or 5-methylcytosine in a cognate methylation position.     

Transient reduction of 5-methylcytosine and 5-hydroxymethylcytosine is associated with active DNA demethylation during regeneration of zebrafish fin

Although dedifferentiation, transformation of differentiated cells into progenitor cells, is a critical step in the regeneration of amphibians and fish, the molecular mechanisms underlying this process, including epigenetic changes, remain unclear. Dot blot assays and immunohistochemical analyses revealed that, during regeneration of zebrafish fin, the levels of 5-methylcytosine (5mC) and 5-hydroxymethylcytosine (5hmC) are transiently reduced in blastema cells and cells adjacent to the amputation plane at 30 h post-amputation (hpa), and the level of 5mC, but not 5hmC, is almost restored by 72 hpa. We observed that the dedifferentiated cells showed reduced levels of 5mC and 5hmC independent of cell proliferation by 24 hpa. Furthermore, expressions of the proposed demethylation- and DNA repair-related genes were detected during fin regeneration. Taken together, our findings illustrate that the transient reduction of 5mC and 5hmC in dedifferentiated cells is associated with active demethylation during regeneration of zebrafish fin.     

Synthesis and physical characterization of DNA fragments containing N4-methylcytosine and 5-methylcytosine.

The synthesis of N4-methyl-2'-deoxycytidine and its fully protected mononucleotide, suitable for the oligonucleotide synthesis by phosphotriester method is described. A set of octanucleotides - d(CGCGCGCG), d(CG5mCGCGCG), d(CG4mCGCGCG) and dodecanucleotides - d(GGACCCGGGTCC), d(GGA5mCCCGGGTCC), d(GGA4mCCCGGGTCC) has been synthesized in a solution. Physical characterization of the oligonucleotide duplexes by means of UV and CD spectrometry provides the evidence that 4mC similarly to 5mC favours the B--greater than Z transition, although both of these methylated cytosines inhibit the B--greater than A conformational change. N4-Methylcytosine in contrast to 5-methylcytosine reduces the DNA double helix thermal stability.     

Crystal and solution structures of the B-DNA dodecamer d(CGCAAATTTGCG) probed by Raman spectroscopy: heterogeneity in the crystal structure does not persist in the solution structure

The self-complementary dodecamer d(CGCAAATTTGCG) crystallizes as a double helix of the B form and manifests a Raman spectrum with features not observed in Raman spectra of either DNA solutions or wet DNA fibers. A number of Raman bands are assigned to specific nucleoside sugar and phosphodiester conformations associated with this model B-DNA crystal structure. The Raman bands proposed as markers of the crystalline B-DNA structure are compared and contrasted with previously proposed markers of Z-DNA and A-DNA crystals. The results indicate that the three canonical forms of DNA can be readily distinguished by Raman spectroscopy. However, unlike Z-DNA and A-DNA, which retain their characteristic Raman fingerprints in aqueous solution, the B-DNA Raman spectrum is not completely conserved between crystal and solution states. The Raman spectra reveal greater heterogeneity of nucleoside conformations (sugar puckers) in the DNA molecules of the crystal structure than in those of the solution structure. The results are consistent with conversion of one-third of the dG residues from the C2'-endo/anti conformation in the solution structure to another conformation, deduced to be C1'-exo/anti, in the crystal. The dodecamer crystal also exhibits unusually broad Raman bands at 790 and 820 cm-1, associated with the geometry of the phosphodiester backbone and indicating a wider range of (alpha, zeta) backbone torsion angles in the crystal than in the solution structure. The results suggest that backbone torsion angles in the CGC and GCG sequences, which flank the central AAATTT sequence, are significantly different for crystal and solution structures, the former containing the greater diversity.(ABSTRACT TRUNCATED AT 250 WORDS)     

TET2 mutations are associated with specific 5-methylcytosine and 5-hydroxymethylcytosine profiles in patients with chronic myelomonocytic leukemia

Chronic myelomonocytic leukemia (CMML) has recently been associated with a high incidence of diverse mutations in genes such as TET2 or EZH2 that are implicated in epigenetic mechanisms. We have performed genome-wide DNA methylation arrays and mutational analysis of TET2, IDH1, IDH2, EZH2 and JAK2 in a group of 24 patients with CMML. 249 genes were differentially methylated between CMML patients and controls. Using Ingenuity pathway analysis, we identified enrichment in a gene network centered around PLC, JNK and ERK suggesting that these pathways, whose deregulation has been recently described in CMML, are affected by epigenetic mechanisms. Mutations of TET2, JAK2 and EZH2 were found in 15 patients (65%), 4 patients (17%) and 1 patient (4%) respectively while no mutations in the IDH1 and IDH2 genes were identified. Interestingly, patients with wild type TET2 clustered separately from patients with TET2 mutations, showed a higher degree of hypermethylation and were associated with higher risk karyotypes. Our results demonstrate the presence of aberrant DNA methylation in CMML and identifies TET2 mutant CMML as a biologically distinct disease subtype with a different epigenetic profile.     

Water mediated Dickerson-Drew-type crystal of DNA dodecamer containing 2'-deoxy-5-formyluridine

To investigate the role of divalent cations in crystal packing, a Dickerson-Drew-type dodecamer with the sequence d(CGCGAATXCGCG), containing 2'-deoxy-5-formyluridine at X, was crystallized under several conditions with Ba(2+) ion instead of Mg(2+) ion. The crystal structure is isomorphous with the original Dickerson-type crystal containing Mg(2+) ion. In the Mg(2+)-free crystals, however, a five-membered ring of water molecules occupies the same position as the magnesium site found in the Mg(2+)-containing crystals, and connects the two duplexes similarly to the hydrated Mg(2+) ion. It has been concluded that the five-membered water molecules can take the place of the hydrated magnesium cation in crystallization. The 5-formyluracil residues form the canonical Watson-Crick pair with the opposite adenine residues.     

The primary mode of binding of cisplatin to a B-DNA dodecamer: C-G-C-G-A-A-T-T-C-G-C-G

When cisplatin [cis- diamminodichloroplatinum (II)] is diffused into pre-grown crystals of the B-DNA double-helical dodecamer C-G-C-G-A-A-T-T-C-G-C-G, it binds preferentially to the N7 positions of guanines, with what probably is an aquo bridge between Pt and the adjacent O6 atom of the same guanine. The entire guanine ring moves slightly toward the platinum site, into the major groove. Only three of the eight potential cisplatin binding sites on guanines actually are occupied, and this differential reactivity can be explained in terms of the relative freedom of motion of guanines toward the major groove. This shift of guanines upon ligation may weaken the glycosyl bond and assist in the depurination that leads to mismatch SOS repair and G.C. to T.A. transversion.     

Sequence-based 5-mers highly correlated to epigenetic modifications in genes interactions

One of the main concerns in biology is extracting sophisticated features from DNA sequence for gene interaction determination, receiving a great deal of researchers’ attention. The epigenetic modifications along with their patterns have been intensely recognized as dominant features affecting on gene expression. However, studying sequenced-based features highly correlated to this key element has remained limited. The main objective in this research was to propose a new feature highly correlated to epigenetic modifications capable of classification of genes. In this paper, classification of 34 genes in PPAR signaling pathway associated with muscle fat tissue in human was performed. Using different statistical outlier detection methods, we proposed that 5-mers highly correlated to epigenetic modifications can correctly categorize the genes involved in the same biological pathway or process. Thirty-four genes in PPAR signaling pathway were classified via applying a proposed feature, 5-mers strongly associated to 17 different epigenetic modifications. For this, diverse statistical outlier detection methods were applied to specify the group of thoroughly correlated genes. The results indicated that these 5-mers can appropriately identify correlated genes. In addition, our results corresponded to GeneMania interaction information, leading to support the suggested method. The appealing findings imply that not only epigenetic modifications but also their highly correlated 5-mers can be applied for reconstructing gene regulatory networks as supplementary data as well as other applications like physical interaction, genes prioritization, indicating some sort of data fusion in this analysis.     

Cleavage of methylated CCCGGG sequences containing either N4-methylcytosine or 5-methylcytosine with MspI, HpaII, SmaI, XmaI and Cfr9I restriction endonucleases.

The cleavage specificity of R.Cfr9I was determined to be C decreases CCGGG whereas the methylation specificity of M.Cfr9I was C4mCCGGG. The action of MspI, HpaII, SmaI, XmaI and Cfr9I restriction endonucleases on an unmethylated parent d(GGACCCGGGTCC) dodecanucleotide duplex and a set of oligonucleotide duplexes, containing all possible substitutions of either 4mC or 5mC for C in the CCCGGG sequence, was investigated. It was found that 4mC methylation, in contrast to 5mC, renders the CCCGGG site resistant to practically all the investigated endonucleases. The cleavage of methylated substrates with restriction endonucleases is discussed.