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.     

Temperature-dependent CD and NMR studies on a synthetic oligonucleotide containing a B-Z junction at high salt

It is now accepted that two or more conformations may exist within the same DNA molecule, thereby generating conformational junctions. The presence of B-Z junctions between right- and left-handed DNA conformations has been detected in plasmids by a number of techniques. Preliminary characterization of the first example of a B-Z junction is a short DNA oligonucleotide has recently been reported [Sheardy, R. D. (1988) Nucleic Acids Res. 16, 1153-1167]. We report additional CD and NMR data that support the existence of the junction in this model oligomer. These studies indicate that only three base pairs are involved in the junction and only one of these is dramatically distorted. Furthermore, the NMR saturation-transfer experiments suggest the junction's internal motion is temperature dependent.     

M.Smal is an N4-methylcytosine specific DNA-methylase.

An enzymatic activity rendering DNA immune to the action of the Smal restriction endonuclease in the presence of S-adenosyl-L-methionine has been detected in Serratia marcescens Sb. This methylase, M.Smal, modifies the second cytosine residue of the substrate sequence CCCGGG yielding N4-methylcytosine.     

DNA methylation in thermophilic bacteria: N4-methylcytosine, 5-methylcytosine, and N6-methyladenine.

While determining the minor and major base composition of the DNA from 17 types of thermophilic bacteria by high performance liquid chromatography (HPLC) of enzymatic digests, we have discovered a novel base, N4-methylcytosine (m4C). Its structure was proven by comparison of the DNA-derived nucleoside to the analogous authentic compound by HPLC, UV spectroscopy, and mass spectroscopy. Eight of the bacterial DNAs contained m4C. Only two contained the common minor base, 5-methylcytosine (m5C), and neither of these was from an extreme thermophile. The other prevalent modified base of bacterial DNA, N6-methyladenine (m6A), was found in nine of the DNAs. Restriction analysis revealed that four of the DNAs had dam-type (Gm6ATC) methylation patterns. Due to the propensity of m5C residues to be deaminated by heat to thymine residues and to inefficient repair of the resulting mismatched base pairs, thermophiles with optimal growth temperatures of greater than or equal to 60 degrees C generally may avoid having m5C in their genomes. Instead, some of them have deamination-resistant m4C residues.     

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.     

Effect of N4-methylcytosine and 5-methylcytosine on the stability of the DNA helix

The thermodynamic parameters (delta H, delta S) of the helix-coil transition of self-complementary oligonucleotides d(CGCGCGCG), d(CG5mCGCGCG), d(CG4mCGCGCG), d(GGACCCGGGTCC), d(GGA5mCCCGGGTCC), and d(GGA4mCCCGGGTCC) were determined. The substitution of 4mC for C was found to decrease the melting temperature of the oligonucleotides. The destabilization effect of the two substitutions is equivalent to the change of A.T for G.C pair. The free energy decrease of the helix-coil transition due to the introduction of two 4mC into an octanucleotide was estimated to be 1,24 kcal/mol.     

NMR studies on an oligodeoxynucleotide containing 2-aminopurine opposite adenine

A heteroduplex containing the mismatch 2-aminopurine (AP)-adenine has been synthesized and studied by proton NMR. The mismatch was incorporated into the sequence d[CGG(AP)GGC].d-(GCCACCG). One-dimensional nuclear Overhauser effect measurements in H2O and two-dimensional nuclear Overhauser effect spectra in D2O show AP.A base pairs in a wobble structure in which both bases are in the anti conformation. The adenine is stacked well in the helix, but the helix twist between the adenine and neighboring cytosine in the 3' direction is unusually small. As a result, the aminopurine on the opposite strand is somewhat pushed out of the helix. From the measurements of the imino proton line widths, the two adjacent G.C base pairs are not found to be significantly destabilized by the presence of the purine-purine wobble pair.     

N4-methylcytosine as a minor base in bacterial DNA.

The DNA base composition, including the minor base content, of 26 strains of bacteria was determined. The studied bacteria are sources of widely used restriction endonucleases. Approximately 35% of the bacterial DNAs contained N4-methylcytosine, about 60% contained 5-methylcytosine, and about 90% had N6-methyladenine.     

Structural model for an oligonucleotide containing a bulged guanosine by NMR and energy minimization

We present three-dimensional structural models for a DNA oligomer containing a bulged guanosine based on proton NMR data and energy minimization computations. The nonexchangeable proton resonances of the duplex 5'd(GATGGGCAG).d(CTGCGCCATC) are assigned by nuclear Overhauser effect spectroscopy (NOESY) and correlated spectroscopy connectivities, and the NMR spectrum is compared with that of a regular 8-mer of similar sequence, 5'd(GATGGCAG).d(CTGCCATC). Experimental proton-proton distances are obtained from NOESY spectra acquired with mixing times of 100, 150, and 200 ms. A refined three-dimensional structure for the bulge-containing duplex is calculated from regular B DNA starting coordinates by using the AMBER molecular mechanics program [Weiner, S. J., Kollman, P. A., Case, D. A., Singh, U. C., Ghio, C., Alagona, G., Profeta, S., & Weiner, P. (1984) J. Am. Chem. Soc. 106, 765-784]. We compare structures obtained by building the helix in three and four base pair increments with structures obtained by direct minimization of the entire nine base sequence, with and without experimental distance constraints. The general features of all the calculated structures are very similar. The helix is of the B family, with the extra guanine stacked into the helix, and the helix axis is bent by 18-23 degrees, in agreement with gel mobility data for bulge-containing sequences [Rice, J. A. (1987) Ph.D. Thesis, Yale University].     

Synthesis of decadeoxyribonucleotides containing N6-methyladenine, N4-methylcytosine, and 5-methylcytosine: recognition and cleavage by restriction endonucleases (nucleosides and nucleotides part 74).

The naturally-occurring modified bases, N6-methyladenine, N4-methylcytosine, and 5-methylcytosine were chemically introduced in place of the adenine or cytosine in the decadeoxyribonucleotides containing recognition sequences of Bgl II, Sau 3AI, Mbo I and Mfl I. The modified oligomers bind to the enzymes but the rates of cleavage by the enzymes are variable.     

NMR and CD studies of triple-helical peptides.

Triple-helix formation of the peptide (Pro-Hyp-Gly)10 was monitored by nmr and CD spectroscopy. The two-dimensional nmr spectra indicated that the Gly C alpha H and Pro C delta H proton resonances shift upfield in going from the nonhelical to helical form, while hydroxy-proline resonances are unchanged. The integrated areas of the helical and nonhelical resonances could be monitored in the one-dimensional nmr spectrum, and indicate that in the (Pro-Hyp-Gly)10 about 90% of the residues are in a defined triple-helical conformation. The introduction of a glycine to alanine substitution or the deletion of a single hydroxyproline residue in the stable triple-helical peptide (Pro-Hyp-Gly)10 still allows trimers to be formed, but the trimers show a substantial loss of triple helix and decreased thermal stability compared with (Pro-Hyp-Gly)10. Two computer models were generated for the Gly----Ala peptide, one with the Ala side chains packed inside the helix and the other with the region containing the alanines forming a beta-bend that loops out from the helix. The nmr data is more consistent with the latter model.     

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.     

Covalent carcinogenic lesions in DNA: NMR studies of O6-methylguanosine containing oligonucleotide duplexes

We report on proton and phosphorus high resolution NMR investigations of the self-complementary dodecanucleotide d(C1-G2-N3-G4-A5-A6-T7-T8-C9-O6meG10-C11-G12) duplexes (henceforth called O6 meG.N 12-mers), N = C, T, A and G, which contain N3.O6meG10 interactions in the interior of the helix. These sequences containing a single modified O6meG per strand were prepared by phosphoamidite synthesis and provide an excellent model for probing the structural basis for covalent carcinogenic lesions in DNA. Distance dependent nuclear Overhauser effect (NOE) measurements and line widths of imino protons demonstrate that the N3 and O6meG.10 bases stack into the duplex and are flanked by stable Watson-Crick base pairs at low temperature for all four O6meG.N 12-mer duplexes. The imino proton of T3 in the O6meG.T 12-mer and G3 in the O6meG.N 12-mer helix, which are associated with the modification site, resonate at unusually high field (8.5 to 9.0 ppm) compared to imino protons in Watson-Crick base pairs (12.5 to 14.5 ppm). The nonexchangeable base and sugar protons have been assigned from two dimensional correlated (COSY) and nuclear Overhauser effect (NOESY) measurements on the O6meG.N 12-mer helices. The directionality of the distance dependent NOEs establish all O6meG.N duplexes to be right-handed helices in solution. The glycosidic torsion angles are in the anti range at the N3.O6meG10 modification site except for O6meG10 in the O6meG.G 12-mer duplex which adopts a syn configuration. This results in altered NOEs between the G3 (anti).O6meG10 (syn) pair and flanking G2.C11 and G4.C9 base pairs in the O6meG.G 12-mer duplex. We observe pattern reversal for cross peaks in the COSY spectrum linking the sugar H1' protons with the H2',2" protons at the G2 and O6meG10 residues in the O6meG.N 12-mer duplexes with the effect least pronounced for the O6meG.T 12-mer helix. The proton chemical shift and NOE data have been analyzed to identify regions of conformational perturbations associated with N3.O6meG10 modification sites in the O6meG.N 12-mer duplexes. The proton decoupled phosphorus spectrum of O6meG.T 12-mer duplex exhibits an unperturbed phosphodiester backbone in contrast to the phosphorus spectra of the O6meG.C 12-mer, O6meG.G 12-mer and O6meG.A 12-mer duplexes which exhibit phosphorus resonances dispersed over 2 ppm characteristic of altered phosphodiester backbones at the modification site. Tentative proposals are put forward for N3.O6meG10 pairing models based on the available NMR data and serve as a guide for the design of future experiments.     

Probing structural factors stabilizing antisense oligonucleotide duplexes: NMR studies of a DNA.DNA duplex containing a formacetal linkage.

The duplex formed by annealing the formacetal backbone modified dodecamer d-(CGCGTTOCH2OTTGCGC) to its complementary strand, d(GCGCAAAACGCG) (duplex I), has been studied by NMR techniques and analyzed with reference to its unmodified counterpart (duplex II). Comparison of parameters such as 2D cross-peak intensities, coupling constants, and spectral patterns indicates that structural perturbations caused by the incorporation of the formacetal linkage are minimal and localized to the central T4.A4 block. Duplex I adopts a B-type helical conformation with regular Watson-Crick base pairing and normal minor groove width. The methylene group is accommodated along the phosphate backbone in a conformation similar to that of the PO2 group found in the B-form DNA family. The central T6-T7 base pairs of duplex I melt simultaneously with the duplex, indicating a cooperative transition to single strands. Although the formacetal linkage affects global melting, as evidenced by a 3 degree C reduction in Tm for duplex I with respect to duplex II, the present study indicates that this is not the result of localized premelting at the formacetal site of duplex I but rather reflects the subtle interplay of several structural and energy factors which need to be further explored.     

NMR studies of a DNA containing 8-hydroxydeoxyguanosine.

The effects of hydroxylation at the C8 of a deoxyguanosine residue in DNA were studied by NMR analysis of a self-complementary dodecanucleotide, d(C1-G2-C3-oh8G4-A5-A6-T7-T8-C9-G10-C11-G12), which has an 8-hydroxy-2'-deoxyguanosine (oh8dG) residue at the 4th position. NMR data indicate that the 8-hydroxyguanine (oh8G) base takes a 6,8-diketo tautomeric form and is base-paired to C with Watson-Crick type hydrogen bonds in a B-form structure. The thermal stability of the duplex is reduced, but the overall structure is much the same as that of the unmodified d(CGCGAATTCGCG) duplex. The structural changes caused by 8-hydroxylation of the deoxyguanosine, if any, are localized near the modification site.     

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.     

Comparison of the conformation of an oligonucleotide containing a central G-T base pair with the non-mismatch sequence by proton NMR.

We have recorded NOESY spectra of two non-selfcomplementary undecanucleotide duplexes. From the observed NOEs we do not detect any significant distortion of the helix when a G-C pair is replaced by a G-T pair and the normal interresidue connectivities can be followed through the mismatch site. We conclude that the 2D spectra of the non-exchangeable protons do not allow differentiation between a wobble or rare tautomer form for the mismatch. NOE measurements in H2O, however, clearly show that the mismatch adopts a wobble structure and give information on the hydration in the minor groove for the G-T base pair which is embedded between two A-T base pairs in the sequence.     

Conformational solution studies of neuropeptide gamma using CD and NMR spectroscopy.

Neuropeptide gamma is one of the largest members of the tachykinin family of peptides, exhibiting strong agonistic activity towards the NK-2 tachykinin receptor. This peptide was synthesized by the solid-phase method using the Fmoc chemistry. Circular-dichroism spectroscopy (CD) investigations of this peptide were performed in phosphate buffer, in the presence of sodium dodecylsulphate (SDS) micelles and trifluoroethanol (TFE) solutions and in DMSO-d6 using the 2D NMR technique in conjunction with two different theoretical approaches. The first assumes multiconformational equilibrium of the peptide studied characterized by the values of statistical weights of low-energy conformations. These calculations were performed using three different force fields ECEPP/3, AMBER4.1 and CHARMM (implemented in the X-PLOR program). The second method incorporates interproton distance and dihedral angle constraints into the starting conformation using the Simulated Annealing algorithm (X-PLOR program). The CD experiments revealed that although the peptide studied is flexible in polar solvents, a tendency to adopt a helical structure was observed in the hydrophobic environment. The NMR data (NOE effects) indicate a helical or reverse structure in the Ile7-His12 fragment of the peptide studied in DMSO-d6 solution. The results obtained cannot be interpreted in terms of a single conformation. Most of the conformations obtained with the ECEPP/3 force field possess a high content of a helical structure. None of the conformers, obtained with the AMBER4.1 and CHARMM force fields, can be considered as the dominant one. In all conformations several beta-turns were detected and in some cases gamma-turns were also found. But in fact, it is rather difficult to select the position of the secondary element(s) present in the structure of NPgamma in solution. All conformers calculated with the X-PLOR program (with using NMR derived distance and torsion angle constraints) are stabilized by several beta-turns. Common structural motives are a type IV beta-turn in the Gln6-His12 fragment. All conformations obtained using two approaches adopt very similar turn shapes in the middle region of molecule and a random structure on the N- and C-terminal fragments.