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.     

B-A and B-Z transitions in deoxyoligoduplexes containing 4- and 5- methylcytosine

Self-complementary oligodeoxynucleotides: GGACCCGGGTCC, GGA4mCCCGGGTCC, GGA5mCCCGGGTCC, CGCGCGCG, CG4mCGCGCG, CG5mCGCGCG were synthetized to study the contribution of methyl groups into the energetics of the three known cooperative transitions in DNA: helix-coil, B-A and B-Z With the use of circular dichroism and absorbtion methods the profiles of the above transitions were obtained by variation of temperature (helix-coil), trifluoroethanol fraction (B-A), NaCl and trifluorethanol contents (B-Z). On the basis of the transition widths and shifts of the transition points due to the methylations the energetics of the methyl groups was estimated. 5mC stabilizes the B form relatively the A form by 0.33 kcal/mol; while 4mC by 0.5 kcal/mol. In the B-Z transition 5 mC stabilizes the Z form by 0.28 kcal/mol relatively the B form; 4mC stabilizes also the Z form although by 0.14 kcal/mol only. Thus, these naturally occurring modifications could modulate substantially the ability of a DNA piece to shift into the A or Z form.     

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.     

Quantitative analysis of DNA secondary structure from solvent-accessible surfaces: the B- to Z-DNA transition as a model

Solvent structure and its interactions have been suggested to play a critical role in defining the conformation of polynucleotides and other macromolecules. In this work, we attempt to quantitate solvent effects on the well-studied conformational transition between right-handed B- and left-handed Z-DNA. The solvent-accessible surfaces of the hexamer sequences d(m5CG)3, d(CG)3, d(CA)3, and d(TA)3 were calculated in their B- and Z-DNA conformations. The difference in hydration free energies between the Z and the B conformations (delta delta GH(Z-B] was determined from these surfaces to be -0.494 kcal/mol for C-5 methylated d(CG), 0.228 kcal/mol for unmethylated d(CG), 0.756 kcal/mol for d(CA)-d(TG), and 0.896 kcal/mol for d(TA) dinucleotides. These delta delta GH(Z-B) values were compared to the experimental B- to Z-DNA transition energies of -0.56 kcal/mol that we measured for C-5 methylated d(CG), 0.69-1.30 kcal/mol reported for unmethylated d(CG), 1.32-1.48 kcal/mol reported for d(CA)-d(TG), and 2.3-2.4 kcal/mol for d(TA) dinucleotides. From this comparison, we found that the calculated delta delta GH(Z-B) of these dinucleotides could account for the previous observation that the dinucleotides were ordered as d(m5CG) greater than d(CG) greater than d(CA)-d(TG) greater than d(TA) in stability as Z-DNA. Furthermore, we predicted that one of the primary reasons for the inability of d(TA) sequences to form Z-DNA results from a decrease in exposed hydrophilic surfaces of adjacent base pairs due to the C-5 methyl group of thymine; thus, d(UA) dinucleotides should be more stable as Z-DNA than the analogous d(TA) dinucleotides.(ABSTRACT TRUNCATED AT 250 WORDS)     

The distribution of 5-methylcytosine in the nuclear genome of plants.

We have determined the 5-methylcytosine (5mC) content in high molecular weight DNA, from two dicot (tobacco and pea) and two monocot (wheat and maize) plant species, fractionated according to base composition. The results show that the proportion of 5mC in the genomic fractions increases linearly with their guanine + cytosine (G + C) content while the proportion of non-methylated cytosine remains almost constant. This can be interpreted as a consequence of a difference in mutation pressure related to spontaneous deamination of 5mC to thymine between the different compartments of plant genomes.     

Thermal stability of the Z conformation of the hexanucleoside pentaphosphate d(br5CGbr5CGbr5CG): evidence for a conformational transition before melting

The thermal stability of the hexanucleoside pentaphosphate d(br5CGbr5CGbr5CG) has been studied at two nucleotide concentrations, in the presence of 1 M NaClO4. At low nucleotide concentration (7 X 10(-5) M), circular dichroism experiments show a conformational transition from the Z conformation to another conformation, named X, which is not the B conformation, as the temperature is increased from 0 to 35 degrees C. Between 40 and 65 degrees C, another transition is observed which corresponds to the melting of the X conformation. At higher nucleotide concentration (2 X 10(-3) M), circular dichroism and 31P nuclear magnetic resonance experiments show that at low temperature (br5dC-dG)3 adopts the Z conformation. There are associations between the oligonucleotides which progressively disappear as the temperature increases. In the range 35-60 degrees C a transition from the Z conformation to another conformation is observed. This new conformation is the X conformation detected at low nucleotide concentration.     

A parallel stranded linear DNA duplex incorporating dG.dC base pairs

DNA oligonucleotides with appropriately designed complementary sequences can form a duplex in which the two strands are paired in a parallel orientation and not in the conventional antiparallel double helix of B-DNA. All parallel stranded (ps) molecules reported to date have consisted exclusively of dA.dT base pairs. We have substituted four dA.dT base pairs of a 25-nt parallel stranded linear duplex (ps-D1.D2) with dG.dC base pairs. The two strands still adopt a duplex structure with the characteristic spectroscopic properties of the ps conformation but with a reduced thermodynamic stability. Thus, the melting temperature of the ps duplex with four dG.dC base pairs (ps-D5.D6) is 10-16 degrees C lower and the van't Hoff enthalpy difference delta HvH for the helix-coil transition is reduced by 20% (in NaCl) and 10% (in MgCl2) compared to that of ps-D1.D2. Based on energy minimizations of a ps-[d(T5GA5).d(A5CT5)] duplex using force field calculations we propose a model for the conformation of a trans dG.dC base pair in a ps helix.     

Competing B-Z and helix-coil conformational transitions in supercoiled plasmid DNA.

The formation of melted regions from A + T-rich sequences and left-handed Z-DNA by alternating purine-pyrimidine sequences will both be facilitated by negative supercoiling, and thus if the sequences are present within the same plasmid molecule they will compete for the free energy of supercoiling. We have studied a series of plasmids that contain either (CG)8 or (TG)12 sequences in either G + C or A + T-rich contexts, by means of two-dimensional gel electrophoresis and chemical modification. We observe both B-Z and helix-coil transitions in all plasmids at elevated temperatures and low ionic strength. The plasmids fall into a number of different classes, in terms of the conformational behavior. As the superhelix density is increased, pCG8/vec ((CG)8 in G + C-rich context) undergoes an initial B-Z transition, followed by melting transitions in sequences remote from the (CG)8 sequence. The two transitions are coupled through the topology of the molecule but are otherwise independent. When the (CG)8 sequence was placed in an A + T-rich context (pCG8/col), the helix-coil transition was perturbed by the presence of the Z-DNA segment. Replacement of the (CG)8 tracts by (TG)12 sequences resulted in a further level of interaction between the transitions. Statistical mechanical modeling of the transitions suggested that at intermediate levels of negative supercoiling the Z-DNA formed by the (TG)12 sequence has a lowered probability due to the helix-coil transition in the A + T-rich sequences. These studies illustrate the complexities of competing conformational equilibria in supercoiled DNA molecules.     

Hydroxyl-radical-induced oxidation of 5-methylcytosine in isolated and cellular DNA

The methylation and oxidative demethylation of cytosine in CpG dinucleotides plays a critical role in the regulation of genes during cell differentiation, embryogenesis and carcinogenesis. Despite its low abundance, 5-methylcytosine (5mC) is a hotspot for mutations in mammalian cells. Here, we measured five oxidation products of 5mC together with the analogous products of cytosine and thymine in DNA exposed to ionizing radiation in oxygenated aqueous solution. The products can be divided into those that arise from hydroxyl radical (•OH) addition at the 5,6-double bond of 5mC (glycol, hydantoin and imidazolidine products) and those that arise from H-atom abstraction from the methyl group of 5mC including 5-hydroxymethylcytosine (5hmC) and 5-formylcytosine (5fC). Based on the analysis of these products, we show that the total damage at 5mC is about 2-fold greater than that at C in identical sequences. The formation of hydantoin products of 5mC is favored, compared to analogous reactions of thymine and cytosine, which favor the formation of glycol products. The distribution of oxidation products is sequence dependent in specific ODN duplexes. In the case of 5mC, the formation of 5hmC and 5fC represents about half of the total of •OH-induced oxidation products of 5mC. Several products of thymine, cytosine, 5mC, as well as 8-oxo-7,8-dihydroguanine (8oxoG), were also estimated in irradiated cells.     

Volumetric properties of the formation of double stranded DNA: a nearest-neighbor analysis

The kinetics of the helix-coil transition have been studied by performing UV-monitored melting and reannealing curves of DNA and analyzing the resultant hysteresis between these curves. The analysis assumes a single-step bimolecular transition with duplex formation defined as the forward reaction. Volume parameters of the helix-coil transition were obtained by measuring the pressure dependence of the rate constants from 5-200 MPa. The data were interpreted in terms of several possible nearest-neighbor models, ranging from one to eleven parameters. Twenty-four oligonucleotide duplexes 22 base pairs in length were used to solve for individual nearest-neighbor activation volumes and transition volumes. Statistically, the most valid fit of the volumetric data was obtained with a six-parameter model in which the directionality of the dinucleotide steps is not considered, for example, 5'AG/CT is the same as 5'GA/TC. The resultant transition volumes at 48 degrees C ranged from -7.1 +/- 0.8 mL/mol (GC/CG) to +2.9 +/- 0.3 mL/mol (AA/TT). The success of the six-parameter model suggests that the relative size of the nearest-neighbor dinucleotides is the most important factor determining the magnitude of the volumetric parameters. The finding that the magnitude of the volumetric parameters correlates with the change in the solvent accessible surface area of the bases during the helix-coil transition corroborates this hypothesis.     

The recognition domain of the methyl-specific endonuclease McrBC flips out 5-methylcytosine

DNA cytosine methylation is a widespread epigenetic mark. Biological effects of DNA methylation are mediated by the proteins that preferentially bind to 5-methylcytosine (5mC) in different sequence contexts. Until now two different structural mechanisms have been established for 5mC recognition in eukaryotes; however, it is still unknown how discrimination of the 5mC modification is achieved in prokaryotes. Here we report the crystal structure of the N-terminal DNA-binding domain (McrB-N) of the methyl-specific endonuclease McrBC from Escherichia coli. The McrB-N protein shows a novel DNA-binding fold adapted for 5mC-recognition. In the McrB-N structure in complex with methylated DNA, the 5mC base is flipped out from the DNA duplex and positioned within a binding pocket. Base flipping elegantly explains why McrBC system restricts only T4-even phages impaired in glycosylation [Luria, S.E. and Human, M.L. (1952) A nonhereditary, host-induced variation of bacterial viruses. J. Bacteriol., 64, 557-569]: flipped out 5-hydroxymethylcytosine is accommodated in the binding pocket but there is no room for the glycosylated base. The mechanism for 5mC recognition employed by McrB-N is highly reminiscent of that for eukaryotic SRA domains, despite the differences in their protein folds.     

Examination of the specificity of DNA methylation profiling techniques towards 5-methylcytosine and 5-hydroxymethylcytosine

DNA cytosine-5 methylation is a well-studied epigenetic pathway implicated in gene expression control and disease pathogenesis. Different technologies have been developed to examine the distribution of 5-methylcytosine (5mC) in specific sequences of the genome. Recently, substantial amounts of 5-hydroxymethylcytosine (5hmC), most likely derived from enzymatic oxidation of 5mC by TET1, have been detected in certain mammalian tissues. Here, we have examined the ability of several commonly used DNA methylation profiling methods to distinguish between 5mC and 5hmC. We show that techniques based on sodium bisulfite treatment of DNA are incapable of distinguishing between the two modified bases. In contrast, techniques based on immunoprecipitation with anti-5mC antibody (methylated DNA immunoprecipitation, MeDIP) or those based on proteins that bind to methylated CpG sequences (e.g. methylated-CpG island recovery assay, MIRA) do not detect 5hmC and are specific for 5mC unless both modified bases occur in the same DNA fragment. We also report that several methyl-CpG binding proteins including MBD1, MBD2 and MBD4 do not bind to sequences containing 5hmC. Selective mapping of 5hmC will require the development of unique tools for the detection of this modified base.     

Sequence-dependent effects on DNA stability resulting from guanosine replacements by inosine.

The effect of replacing a G.C base-pair with an I.C base-pair on DNA stability was investigated for a related set of 14-mers. DNA melting analysis of the 14-mers was used to determine delta Hzero, delta Szero and delta G(zero)37 of the double to single stranded transition. All 14mers were shown to have B-DNA character by circular dichroism analysis. 14mers substituted with a single inosine in place of guanosine at different positions showed that consequences on DNA stability are sequence-dependent. Large changes in delta Hzero and delta Szero result when inosine is substituted within the trinucleotide sequence d(TCG).d(CGA) while substitution within d(TCC).d(GGA) causes minor changes in enthalpy and entropy. Moreover, some 14-mers with two inosine substitutions five base-pairs apart showed non-additive free energy changes for the double to single stranded transition. These results clearly indicate that the structural consequences of replacing a single guanosine with an inosine are transmitted over a significant distance.     

Infrared spectral studies of the non regularly alternating purine-pyrimidine hexamers d(m5CGGCM5CG), d(CBr8GGCCBr8G) and d(CGCGGC)

The oligonucleotides d(m5CGGCm5CG), d(CBr8GGCCBr8G) and d(CGCGGC) have been prepared and studied by infrared spectroscopy. The three sequences contain two GC pairs which are out of purine-pyrimidine alternation with the rest of the sequence. From the IR data of the d(m5CGGCm5CG) hexamer, it is shown that all of the dG residues adopt a syn conformation. The marker IR bands for the C3' endo syn conformation are at 1410, 1354, 1320 and 925 cm-1 whereas those for the C2' endo anti conformation at 1420, 1374 and 890 cm-1 are clearly absent. This result implies that the two adjacent guanines of the d(m5CGGCm5CG) sequence are in syn conformation. It is suggested that duplex formation occurs in d(CGCGGC) films and that all of the guanines are in syn conformation. In contrast, the central non-brominated guanine of the d(CBr8GGCCBr8G) hexamer is found in anti conformation, as expected in a Z type structure of the non-alternating region.     

B to Z transitions of the short DNA hairpins formed from the oligomer sequences: d[(CG)3X4(CG)3] (X = A, T, G, C).

Circular Dichroism (CD) spectra were collected as a function of sodium perchlorate concentration [NaClO4] for the set of DNA hairpins formed from the oligomer sequences d[(CG)3X4(CG)3] where X = A, T, G or C. Over the range in salt concentration from 0 to 4.0 M NaClO4, the CD spectra invert in a manner characteristic of the B to Z transition. A factor analysis routine is described and employed to determine the least number of basis spectra required to fit the measured spectra of each hairpin over the entire salt range examined. In every case, linear combinations of only two sub-spectra fit the experimental spectra of the hairpins with greater than 98% accuracy, indicating the spectrally monitored structural transitions are two-state. From the relative weights of the individual sub-spectra, B-Z transition curves are constructed. The transitions are analyzed in terms of a simple two-state equilibrium model which yields an evaluation of the transition free-energy, delta GB-Z, as a function of NaClO4 concentration. At 1.0 M NaClO4 and 21 degrees C, delta GB-Z = 5.4, 4.9, 3.6 and 2.3 kcal/mole for the G4, T4, A4 and C4 loop hairpins, respectively.     

Effect of P-chirality of internucleotide bonds on B-Z conversion of stereodefined self-complementary phosphorothioate oligonucleotides of the [PS]-d(CG)4 and [PS]-d(GC)4 series

Diastereomerically pure, partially modified (in selected positions) or fully modified phosphorothioate oligomers of the [PS]-d(CG)(4) and [PS]-d(GC)(4) series were investigated with respect to their ability to adopt the left-handed conformation at high sodium chloride concentration. NaCl induces the B-Z transition of [All-S(P)R(P)-PS]-d(CG)(4) with a midpoint of transition at ca. 2 M, which is approximately 1 M less than for unmodified d(CG)(4). Also, [All-R(P)S(P)-PS]-d(GC)(4) at 5 M NaCl converts to the Z form to the extent of ca. 55%, while the unmodified d(GC)(4) counterpart does not convert at all. This enhanced ability of stereodefined phosphorothioate oligomers to adopt the Z conformation is discussed in terms of already known structural factors (hydrogen bonding and water bridges) facilitating the B-Z transition, identified for unmodified d(CG)(n) oligonucleotides. By CD spectroscopy, the [All-S(P)-PS]-d(CG)(4) oligomer at a NaCl concentration higher than 0.01 M adopts a unique conformation as assessed from the presence of an additional negative band centered at 282 nm.     

Calculation of the stability of mini-duplexes of DNA in an electrolyte solution using a molecular mechanics method in approximating a statistical model of the environment

Approximation of the statistical model of environment (SME) to estimate the energy of the macromolecule in electrolyte solution has been developed and used for calculating the conformational energy of nucleic acids by means of the molecular mechanics method. Calculation of base pairs opening delta Hcalop enthalpies and enthalpies of DNA miniduplexes dissociation delta Hcaldis were performed for 10 types of diduplexes. The approximation SME enable to perform calculations of the absolute base-dependent values of delta Hcalop which coincide in the range of 1 kcal/mol with the experimental base-dependent values of the helix-coil transition enthalpies. Values of dissociation enthalpies delta Hcalop greater than delta Hcaldis for all miniduplexes, the difference of delta Hcalop--delta Hcaldis determine the base-dependent energy of the helix-coil boundary. The values of activation barriers for the strands dissociation delta H d not equal to congruent to 8 kcal/mol and association delta H not equal to as congruent to 4 kcal/mol were obtained for GG/CC and AA/TT duplexes. It is concluded that the approximation SME enables to increase substantially the accuracy of the calculations of the macromolecule conformational rearrangement enthalpies in the electrolytes solution.     

Base-resolution detection of N4-methylcytosine in genomic DNA using 4mC-Tet-assisted-bisulfite- sequencing

Restriction-modification (R-M) systems pose a major barrier to DNA transformation and genetic engineering of bacterial species. Systematic identification of DNA methylation in R-M systems, including N⁶-methyladenine (6mA), 5-methylcytosine (5mC) and N⁴-methylcytosine (4mC), will enable strategies to make these species genetically tractable. Although single-molecule, real time (SMRT) sequencing technology is capable of detecting 4mC directly for any bacterial species regardless of whether an assembled genome exists or not, it is not as scalable to profiling hundreds to thousands of samples compared with the commonly used next-generation sequencing technologies. Here, we present 4mC-Tet-assisted bisulfite-sequencing (4mC-TAB-seq), a next-generation sequencing method that rapidly and cost efficiently reveals the genome-wide locations of 4mC for bacterial species with an available assembled reference genome. In 4mC-TAB-seq, both cytosines and 5mCs are read out as thymines, whereas only 4mCs are read out as cytosines, revealing their specific positions throughout the genome. We applied 4mC-TAB-seq to study the methylation of a member of the hyperthermophilc genus, Caldicellulosiruptor, in which 4mC-related restriction is a major barrier to DNA transformation from other species. In combination with MethylC-seq, both 4mC- and 5mC-containing motifs are identified which can assist in rapid and efficient genetic engineering of these bacteria in the future.     

Effect of 5-methylcytosine on the stability of triple-stranded DNA--a thermodynamic study.

We have previously shown that the pyrimidine oligonucleotide 5'CTTCCTCCTCT (Y11) recognizes the double-helical stem of hairpin 5'GAAGGAGGAGA-T4-TCTCCTCCTTC (h26) by triple-helix formation (1). In this paper, we report the effect on triplex formation of substituting the cytosine residues of Y11 with 5-methylcytosines (5meY11). In addition, we have studied the thermodynamics of the interaction between h26 and 5meY11. The results can be summarised as follows: (i) gel electrophoresis shows that at T = 5 degrees C and pH 5, both Y11 and 5meY11 form DNA triple helices with h26, whereas at pH 6.8 only the methylated strand binds to h26; (ii) pH-stability curves of the DNA triplexes formed from h26 + Y11 and h26 + 5meY11 show that Y11 and 5meY11 are semi-protonated at pH 5.7 and 6.7, respectively. Thus, it is concluded that cytosine methylation expands the pH range compatible with triplex formation by one pH unit; (iii) as the unmethylated triplex (h26:Y11), the methylated one (h26:5meY11) denatures in a biphasic manner, in which the low temperature transition results from the dissociation of 5meY11 from h26. The Tm of the triplex to h26 plus 5meY11 transition is strongly enhanced (about 10 degrees C) by cytosine methylation. A van 't Hoff analysis of denaturation curves is presented; (iv) DSC experiments show that triplex formation between 5meY11 and h26 is characterized by delta H = -237 +/- 25 kJ/mol and delta S = -758 +/- 75 J/Kmol, corresponding to an average delta H of -21 kJ/mol and delta S of -69 J/Kmol per Hoogsteen base pair; (v) the thermodynamic analysis indicates that the extra stability imparted to the triplex by methylcytosine is entropic in origin.