Sequence dependent effects of CpG cytosine methylation. A joint 1H-NMR and 31P-NMR study

Kinetic Referencing Allows Identification of Epigenetic Cytosine Modifications by Single-Molecule Hybridization Kinetics and Superresolution DNA-PAINT Microscopy

We develop a DNA origami-based internal kinetic referencing system with a colocalized reference and target molecule to provide increased sensitivity and robustness for transient binding kinetics. To showcase this, we investigate the subtle changes in binding strength of DNA oligonucleotide hybrids induced by cytosine modifications. These cytosine modifications, especially 5-methylcytosine but also its oxidized derivatives, have been increasingly studied in the context of epigenetics. Recently revealed correlations of epigenetic modifications and disease also render them interesting biomarkers for early diagnosis. Internal kinetic referencing allows us to probe and compare the influence of the different epigenetic cytosine modifications on the strengths of 7-nucleotide long DNA hybrids with one or two modified nucleotides by single-molecule imaging of their transient binding, revealing subtle differences in binding times. Interestingly, the influence of epigenetic modifications depends on their position in the DNA strand, and in the case of two modifications, effects are additive. The sensitivity of the assay indicates its potential for the direct detection of epigenetic disease markers.

Cytosine base modifications regulate DNA duplex stability and metabolism

DNA base modifications diversify the genome and are essential players in development. Yet, their influence on DNA physical properties and the ensuing effects on genome metabolism are poorly understood. Here, we focus on the interplay of cytosine modifications and DNA processes. We show by a combination of in vitro reactions with well-defined protein compositions and conditions, and in vivo experiments within the complex networks of the cell that cytosine methylation stabilizes the DNA helix, increasing its melting temperature and reducing DNA helicase and RNA/DNA polymerase speed. Oxidation of methylated cytosine, however, reverts the duplex stabilizing and genome metabolic effects to the level of unmodified cytosine. We detect this effect with DNA replication and transcription proteins originating from different species, ranging from prokaryotic and viral to the eukaryotic yeast and mammalian proteins. Accordingly, lack of cytosine methylation increases replication fork speed by enhancing DNA helicase unwinding speed in cells. We further validate that this cannot simply be explained by altered global DNA decondensation, changes in histone marks or chromatin structure and accessibility. We propose that the variegated deposition of cytosine modifications along the genome regulates DNA helix stability, thereby providing an elementary mechanism for local fine-tuning of DNA metabolism.

Combined Effects of Methylated Cytosine and Molecular Crowding on the Thermodynamic Stability of DNA Duplexes

Methylated cytosine within CpG dinucleotides is a key factor for epigenetic gene regulation. It has been revealed that methylated cytosine decreases DNA backbone flexibility and increases the thermal stability of DNA. Although the molecular environment is an important factor for the structure, thermodynamics, and function of biomolecules, there are few reports on the effects of methylated cytosine under a cell-mimicking molecular environment. Here, we systematically investigated the effects of methylated cytosine on the thermodynamics of DNA duplexes under molecular crowding conditions, which is a critical difference between the molecular environment in cells and test tubes. Thermodynamic parameters quantitatively demonstrated that the methylation effect and molecular crowding effect on DNA duplexes are independent and additive, in which the degree of the stabilization is the sum of the methylation effect and molecular crowding effect. Furthermore, the effects of methylation and molecular crowding correlate with the hydration states of DNA duplexes. The stabilization effect of methylation was due to the favorable enthalpic contribution, suggesting that direct interactions of the methyl group with adjacent bases and adjacent methyl groups play a role in determining the flexibility and thermodynamics of DNA duplexes. These results are useful to predict the properties of DNA duplexes with methylation in cell-mimicking conditions.

6-phenylpyrrolocytosine as a fluorescent probe to examine nucleotide flipping catalyzed by a DNA repair protein

Cellular exposure to tobacco-specific nitrosamines causes formation of promutagenic O⁶ -[4-oxo-4-(3-pyridyl)but-1-yl]guanine (O⁶ -POB-G) and O⁶ -methylguanine (O⁶ -Me-G) adducts in DNA. These adducts can be directly repaired by O⁶ -alkylguanine-DNA alkyltransferase (AGT). Repair begins by flipping the damaged base out of the DNA helix. AGT binding and base-flipping have been previously studied using pyrrolocytosine as a fluorescent probe paired to the O⁶ -alkylguanine lesion, but low fluorescence yield limited the resolution of steps in the repair process. Here, we utilize the highly fluorescent 6-phenylpyrrolo-2'-deoxycytidine (6-phenylpyrrolo-C) to investigate AGT-DNA interactions. Synthetic oligodeoxynucleotide duplexes containing O⁶ -POB-G and O⁶ -Me-G adducts were placed within the CpG sites of codons 158, 245, and 248 of the p53 tumor suppressor gene and base-paired to 6-phenylpyrrolo-C in the opposite strand. Neighboring cytosine was either unmethylated or methylated. Stopped-flow fluorescence measurements were performed by mixing the DNA duplexes with C145A or R128G AGT variants. We observe a rapid, two-step, nearly irreversible binding of AGT to DNA followed by two slower steps, one of which is base-flipping. Placing 5-methylcytosine immediately 5' to the alkylated guanosine causes a reduction in rate constant of nucleotide flipping. O⁶ -POB-G at codon 158 decreased the base flipping rate constant by 3.5-fold compared with O⁶ -Me-G at the same position. A similar effect was not observed at other codons.

DNA Modification Readers and Writers and Their Interplay

Genomic DNA is modified in a postreplicative manner and several modifications, the enzymes responsible for their deposition as well as proteins that read these modifications, have been described. Here, we focus on the impact of DNA modifications on the DNA helix and review the writers and readers of cytosine modifications and how they interplay to shape genome composition, stability, and function.

Oxidative damage to epigenetically methylated sites affects DNA stability, dynamics and enzymatic demethylation

DNA damage can affect various regulatory elements of the genome, with the consequences for DNA structure, dynamics, and interaction with proteins remaining largely unexplored. We used solution NMR spectroscopy, restrained and free molecular dynamics to obtain the structures and investigate dominant motions for a set of DNA duplexes containing CpG sites permuted with combinations of 5-methylcytosine (mC), the primary epigenetic base, and 8-oxoguanine (oxoG), an abundant DNA lesion. Guanine oxidation significantly changed the motion in both hemimethylated and fully methylated DNA, increased base pair breathing, induced BI→BII transition in the backbone 3′ to the oxoG and reduced the variability of shift and tilt helical parameters. UV melting experiments corroborated the NMR and molecular dynamics results, showing significant destabilization of all methylated contexts by oxoG. Notably, some dynamic and thermodynamic effects were not additive in the fully methylated oxidized CpG, indicating that the introduced modifications interact with each other. Finally, we show that the presence of oxoG biases the recognition of methylated CpG dinucleotides by ROS1, a plant enzyme involved in epigenetic DNA demethylation, in favor of the oxidized DNA strand. Thus, the conformational and dynamic effects of spurious DNA oxidation in the regulatory CpG dinucleotide can have far-reaching biological consequences.

A sequence environment modulates the impact of methylation on the torsional rigidity of DNA

We explain, with atomic-level detail, how DNA methylation may contribute to the regulation of biological processes. We show that DNA sequence environment modulates the impact of methylation on the physical properties of the double helix.

Computational Modeling Approach in Probing the Effects of Cytosine Methylation on the Transcription Factor Binding to DNA

INTRODUCTION

Cytosine methylation at CpG dinucleotides is a chief mechanism in epigenetic modification of gene expression patterns. Previous studies demonstrated that increased CpG methylation of Sp1 sites at -268 and -346 of protein kinase C ε promoter repressed the gene expression.

MATERIALS & METHODS

The present study investigated the impact of CpG methylation on the Sp1 binding via molecular modeling and electrophoretic mobility shift assay. Each of the Sp1 sites contain two CpGs. Methylation of either CpG lowered the binding affinity of Sp1, whereas methylation of both CpGs produced a greater decrease in the binding affinity. Computation of van der Waals (VDW) energy of Sp1 in complex with the Sp1 sites demonstrated increased VDW values from one to two sites of CpG methylation. Molecular modeling indicated that single CpG methylation caused underwinding of the DNA fragment, with the phosphate groups at C1, C4 and C5 reoriented from their original positions. Methylation of both CpGs pinched the minor groove and increased the helical twist concomitant with a shallow, hydrophobic major groove. Additionally, double methylation eliminated hydrogen bonds on recognition helix residues located at positions -1 and 1, which were essential for interaction with O6/N7 of G-bases. Bonding from linker residues Arg565, Lys595 and Lys596 were also reduced. Methylation of single or both CpGs significantly affected hydrogen bonding from all three Sp1 DNA binding domains, demonstrating that the consequences of cytosine modification extend beyond the neighboring nucleotides.

RESULTS

The results indicate that cytosine methylation causes subtle structural alterations in Sp1 binding sites consequently resulting in inhibition of side chain interactions critical for specific base recognition and reduction of the binding affinity of Sp1.

Quantitation of a low level coeluting impurity present in a modified oligonucleotide by both LC-MS and NMR

This paper describes the use of two complementary techniques, LC-MS and NMR, to quantify a low level mono phosphate substituted impurity in an oligonucleotide drug substance. This impurity is the result of a sulphurisation failure, leading to the production of a sequence where a phosphorothioate linkage is replaced by a phosphate. Few quantitative methods are possible to analyse these challenging molecules especially if reversed phase ion pair chromatography, one of the most commonly used techniques for the separation of oligonucleotides, is unable to resolve the impurity in question. With the use of a standard addition method it could be demonstrated that both analytical techniques show equivalency and furthermore, the LC-MS method alone with additional validation has the potential to perform this quantitative assay with a high degree of accuracy.

Can A Denaturant Stabilize DNA? Pyridine Reverses DNA Denaturation in Acidic pH

The stability of DNA is highly dependent on the properties of the surrounding solvent, such as ionic strength, pH, and the presence of denaturants and osmolytes. Addition of pyridine is known to unfold DNA by replacing π-π stacking interactions between bases, stabilizing conformations in which the nucleotides are solvent exposed. We show here experimental and theoretical evidences that pyridine can change its role and in fact stabilize the DNA under acidic conditions. NMR spectroscopy and MD simulations demonstrate that the reversal in the denaturing role of pyridine is specific, and is related to its character as pseudo groove binder. The present study sheds light on the nature of DNA stability and on the relationship between DNA and solvent, with clear biotechnological implications.

C-5 Propynyl Modifications Enhance the Mechanical Stability of DNA

Increased thermal or mechanical stability of DNA duplexes is desired for many applications in nanotechnology or -medicine where DNA is used as a programmable building block. Modifications of pyrimidine bases are known to enhance thermal stability and have the advantage of standard base-pairing and easy integration during chemical DNA synthesis. Through single-molecule force spectroscopy experiments with atomic force microscopy and the molecular force assay we investigated the effect of pyrimidines harboring C-5 propynyl modifications on the mechanical stability of double-stranded DNA. Utilizing these complementary techniques, we show that propynyl bases significantly increase the mechanical stability if the DNA is annealed at high temperature. In contrast, modified DNA complexes formed at room temperature and short incubation times display the same stability as non-modified DNA duplexes.

Parallel force assay for protein-protein interactions

Quantitative proteome research is greatly promoted by high-resolution parallel format assays. A characterization of protein complexes based on binding forces offers an unparalleled dynamic range and allows for the effective discrimination of non-specific interactions. Here we present a DNA-based Molecular Force Assay to quantify protein-protein interactions, namely the bond between different variants of GFP and GFP-binding nanobodies. We present different strategies to adjust the maximum sensitivity window of the assay by influencing the binding strength of the DNA reference duplexes. The binding of the nanobody Enhancer to the different GFP constructs is compared at high sensitivity of the assay. Whereas the binding strength to wild type and enhanced GFP are equal within experimental error, stronger binding to superfolder GFP is observed. This difference in binding strength is attributed to alterations in the amino acids that form contacts according to the crystal structure of the initial wild type GFP-Enhancer complex. Moreover, we outline the potential for large-scale parallelization of the assay.

Comparison of the structural and dynamic effects of 5-methylcytosine and 5-chlorocytosine in a CpG dinucleotide sequence

Inflammation-mediated reactive molecules can result in an array of oxidized and halogenated DNA-damage products, including 5-chlorocytosine ((Cl)C). Previous studies have shown that (Cl)C can mimic 5-methylcytosine ((m)C) and act as a fraudulent epigenetic signal, promoting the methylation of previously unmethylated DNA sequences. Although the 5-halouracils are good substrates for base-excision repair, no repair activity has yet been identified for (Cl)C. Because of the apparent biochemical similarities of (m)C and (Cl)C, we have investigated the effects of (m)C and (Cl)C substitution on oligonucleotide structure and dynamics. In this study, we have constructed oligonucleotide duplexes containing C, (Cl)C, and (m)C within a CpG dinucleotide. The thermal and thermodynamic stability of these duplexes were found to be experimentally indistinguishable. Crystallographic structures of duplex oligonucleotides containing (m)C or (Cl)C were determined to 1.2 and 1.9 Å resolution, respectively. Both duplexes are B-form and are superimposable on a previously determined structure of a cytosine-containing duplex with a rmsd of approximately 0.25 Å. NMR solution studies indicate that all duplexes containing cytosine or the cytosine analogues are normal B-form and that no structural perturbations are observed surrounding the site of each substitution. The magnitude of the base-stacking-induced upfield shifts for nonexchangeable base proton resonances are similar for each of the duplexes examined, indicating that neither (m)C nor (Cl)C significantly alter base-stacking interactions. The (Cl)C analogue is paired with G in an apparently normal geometry; however, the G-imino proton of the (Cl)C-G base pair resonates to higher field relative to (m)C-G or C-G, indicating a weaker imino hydrogen bond. Using selective ¹⁵N-enrichment and isotope-edited NMR, we observe that the amino group of (Cl)C rotates at roughly half of the rate of the corresponding amino groups of the C-G and (m)C-G base pairs. The altered chemical shifts of hydrogen-bonding proton resonances for the (Cl)C-G base pair as well as the slower rotation of the (Cl)C amino group can be attributed to the electron-withdrawing inductive property of the 5-chloro substituent. The apparent similarity of duplexes containing (m)C and (Cl)C demonstrated here is in accord with results of previous biochemical studies and further suggests that (Cl)C is likely to be an unusually persistent form of DNA damage.

Effects of DNA methylation on nucleosome stability

Methylation of DNA at CpG dinucleotides represents one of the most important epigenetic mechanisms involved in the control of gene expression in vertebrate cells. In this report, we conducted nucleosome reconstitution experiments in conjunction with high-throughput sequencing on 572 KB of human DNA and 668 KB of mouse DNA that was unmethylated or methylated in order to investigate the effects of this epigenetic modification on the positioning and stability of nucleosomes. The results demonstrated that a subset of nucleosomes positioned by nucleotide sequence was sensitive to methylation where the modification increased the affinity of these sequences for the histone octamer. The features that distinguished these nucleosomes from the bulk of the methylation-insensitive nucleosomes were an increase in the frequency of CpG dinucleotides and a unique rotational orientation of CpGs such that their minor grooves tended to face toward the histones in the nucleosome rather than away. These methylation-sensitive nucleosomes were preferentially associated with exons as compared to introns while unmethylated CpG islands near transcription start sites became enriched in nucleosomes upon methylation. The results of this study suggest that the effects of DNA methylation on nucleosome stability in vitro can recapitulate what has been observed in the cell and provide a direct link between DNA methylation and the structure and function of chromatin.

Effects of cytosine hydroxymethylation on DNA strand separation

Cytosine hydroxymethylation is an epigenetic control factor in higher organisms. New discoveries of the biological roles of hydroxymethylation serve to raise questions about how this epigenetic modification exerts its functions and how organisms discriminate cytosine hydroxymethylation from methylation. Here, we report investigations that reveal an effect of cytosine hydroxymethylation on mechanical properties of DNA under load. The findings are based on molecular force assay measurements and steered molecular dynamics simulations. Molecular force assay experiments identified significant effects of hydroxymethylation on stretching-induced strand separation; the underlying physical mechanism has been revealed by steered molecular dynamics simulations. We find that hydroxymethylation can either upregulate or downregulate DNA's strand separation propensity, suggesting that hydroxymethylation can control gene expression by facilitating or obstructing the action of transcription machinery or the access to chromosomal DNA.

Cytosine methylation alters DNA mechanical properties

DNA methylation plays an essential role in transcriptional control of organismal development in epigenetics, from turning off a specific gene to inactivation of entire chromosomes. While the biological function of DNA methylation is becoming increasingly clear, the mechanism of methylation-induced gene regulation is still poorly understood. Through single-molecule force experiments and simulation we investigated the effects of methylation on strand separation of DNA, a crucial step in gene expression. Molecular force assay and single-molecule force spectroscopy revealed a strong methylation dependence of strand separation. Methylation is observed to either inhibit or facilitate strand separation, depending on methylation level and sequence context. Molecular dynamics simulations provided a detailed view of methylation effects on strand separation, suggesting the underlying physical mechanism. According to our study, methylation in epigenetics may regulate gene expression not only through mechanisms already known but also through changing mechanical properties of DNA.

Direct detection and quantification of methylation in nucleic acid sequences using high-resolution melting analysis

High-resolution melting (HRM) analysis exploits the reduced thermal stability of DNA fragments that contain base mismatches to detect single nucleotide polymorphisms (SNPs). However, the capacity of HRM to reveal other features of DNA chemistry remains unexplored. DNA methylation plays a key role in regulating gene expression and is essential for normal development in many higher organisms. The presence of methylated bases perturbs the double-stranded DNA structure, although its effect on DNA thermal stability is largely unknown. Here, we reveal that methylated DNA has enhanced thermal stability and is sufficiently divergent from nonmethylated DNA to allow detection and quantification by HRM analysis. This approach reliably distinguishes between sequence-identical DNA differing only in the methylation of one base. The method also provides accurate discrimination between mixes of methylated and nonmethylated DNAs, allowing discrimination between DNA that is 1% and 0% methylated and also between 97.5% and 100% methylated. Thus, the method provides a new means of adjusting thermal optima for DNA hybridization and PCR-based techniques and to empirically measure the impact of DNA methylation marks on the thermostability of regulatory regions. In the longer term, it could enable the development of new techniques to quantify methylated DNA.

Cytosine methylation effects on the repair of O6-methylguanines within CG dinucleotides

O(6)-alkyldeoxyguanine adducts induced by tobacco-specific nitrosamines are repaired by O(6)-alkylguanine DNA alkyltransferase (AGT), which transfers the O(6)-alkyl group from the damaged base to a cysteine residue within the protein. In the present study, a mass spectrometry-based approach was used to analyze the effects of cytosine methylation on the kinetics of AGT repair of O(6)-methyldeoxyguanosine (O(6)-Me-dG) adducts placed within frequently mutated 5'-CG-3' dinucleotides of the p53 tumor suppressor gene. O(6)-Me-dG-containing DNA duplexes were incubated with human recombinant AGT protein, followed by rapid quenching, acid hydrolysis, and isotope dilution high pressure liquid chromatography-electrospray ionization tandem mass spectrometry analysis of unrepaired O(6)-methylguanine. Second-order rate constants were calculated in the absence or presence of the C-5 methyl group at neighboring cytosine residues. We found that the kinetics of AGT-mediated repair of O(6)-Me-dG were affected by neighboring 5-methylcytosine ((Me)C) in a sequence-dependent manner. AGT repair of O(6)-Me-dG adducts placed within 5'-CG-3' dinucleotides of p53 codons 245 and 248 was hindered when (Me)C was present in both DNA strands. In contrast, cytosine methylation within p53 codon 158 slightly increased the rate of O(6)-Me-dG repair by AGT. The effects of (Me)C located immediately 5' and in the base paired position to O(6)-Me-dG were not additive as revealed by experiments with hypomethylated sequences. Furthermore, differences in dealkylation rates did not correlate with AGT protein affinity for cytosine-methylated and unmethylated DNA duplexes or with the rates of AGT-mediated nucleotide flipping, suggesting that (Me)C influences other kinetic steps involved in repair, e.g. the rate of alkyl transfer from DNA to AGT.

An FTIR investigation of flanking sequence effects on the structure and flexibility of DNA binding sites

Fourier transform infrared (FTIR) spectroscopy and a library of FTIR marker bands have been used to examine the structure and relative flexibilities conferred by different flanking sequences on the EcoRI binding site. This approach allowed us to examine unique peaks and subtle changes in the spectra of d(AAAGAATTCTTT)(2), d(TTCGAATTCGAA)(2), and d(CGCGAATTCGCG)(2) and thereby identify local changes in base pairing, base stacking, backbone conformation, glycosidic bond rotation, and sugar puckering in the studied sequences. The changes in flanking sequences induce differences in the sugar puckers, glycosidic bond rotation, and backbone conformations. Varying levels of local flexibility are observed within the sequences in agreement with previous biological activity assays. The results also provide supporting evidence for the presence of a splay in the G(4)-C(9) base pair of the EcoRI binding site and a potential pocket of flexibility at the G(4) cleavage site that have been proposed in the literature. In sum, we have demonstrated that FTIR is a powerful methodology for studying the effect of flanking sequences on DNA structure and flexibility, for it can provide information about the local structure of the nucleic acid and the overall relative flexibilities conferred by different flanking sequences.

DNA methylation impacts the cleavage activity of Chlorella virus topoisomerase II

Topoisomerase II from Paramecium bursaria chlorella virus-1 (PBCV-1) and chlorella virus Marburg-1 (CVM-1) displays an extraordinarily high in vitro DNA cleavage activity that is 30-50 times higher than that of human topoisomerase IIalpha. This remarkable scission activity may reflect a unique role played by the type II enzyme during the viral life cycle that extends beyond the normal control of DNA topology. Alternatively, but not mutually exclusively, it may reflect an adaptation to some aspect of the viral environment that differs from the in vitro conditions. To this point, the genomes of many chlorella viruses contain high levels of N6-methyladenine (6mA) and 5-methylcytosine (5mC), but the DNA employed in vitro is unmodified. Therefore, to determine whether methylation impacts the ability of chlorella virus topoisomerase II to cleave DNA, the effects of 6mA and 5mC on the PBCV-1 and CVM-1 enzymes were examined. Results indicate that 6mA strongly inhibits DNA scission mediated by both enzymes, while 5mC has relatively little effect. At levels of 6mA and 5mC methylation comparable to those found in the CVM-1 genome (10% 6mA and 42% 5mC), the level of DNA cleavage decreased approximately 4-fold. As determined using a novel rapid quench pre-equilibrium DNA cleavage system in conjunction with oligonucleotide binding and ligation assays, this decrease appears to be caused primarily by a slower forward rate of DNA scission. These findings suggest that the high DNA cleavage activity of chlorella virus topoisomerase II on unmodified nucleic acid substrates may reflect, at least in part, an adaptation to act on methylated genomic DNA.

Sequence preference for BI/BII conformations in DNA: MD and crystal structure data analysis

Deciphering sequence information from sugar-phosphate backbone is finely tuned through the conformational substates of DNA. BII conformation, one of the conformational substates of B-DNA, is known to play a key role in DNA-protein recognition. BI and BII are identified by the epsilon-zeta difference, which is negative in BI and positive in BII. Our analysis of MD and crystal structures shows that BII conformation is sequence specific and dinucleotides GC, CG, CA, TG, TA show high preference to take up BII conformation, while TT, TC, CT, CC dinucleotides rarely take up this conformation. Significant changes were observed in the dinucleotide parameters viz. twist, roll, and slide for the steps having BII conformation. Interestingly, the magnitude of variation in the dinucleotide parameters is seen to depend mainly on two factors, the magnitude of epsilon-zeta difference and the presence or absence of BII conformation in the second strand, across the WC base-paired dinucleotide step. Based on these two factors, the conformational substate of a dinucleotide step can be further classified as BI.BI (BI conformation in both strands), BI.BII (BI conformation in one strand and BII conformation in the other), and BII.BII (BII conformation in both strands). The occurrence of BII in both strands was found to be quite rare and thus, it can be concluded that BI.BI and BI.BII hybrid steps are more favorable than a BII.BII step. In conformity with the sequence preference seen for dinucleotides in each strand, BII.BII combination of backbone conformation was observed only for GC, CG, CA, and TG containing dinucleotide steps. We further classified BII.BII step as strong BII and weak BII depending on the magnitude of the average epsilon-zeta difference. The dinucleotide steps which belong to the category of strong BII, have large twist, high positive slide and negative roll values, while those in the weak BII group have roll, twist, and slide values similar to that of hybrid BI.BII steps. This conformational property could be contributing to the groove opening/closing and thus can modulate protein-DNA interaction.

Structural effects of cytosine methylation on DNA sugar pucker studied by FTIR

This FTIR investigation concerns structural consequences of 5-methylation of cytosine in a DNA decamer in solution. Methylation of DNA is an important functional signal in transcription, but its effect on DNA structure is variable and not fully understood. Here, single and multiple 5-methylcytosine substitutions are introduced into the self-complementary sequence d(CCGGCGCCGG)(2). No major structural effect of methylation on the DNA duplex in solution is seen in the IR spectra: The overall B-form character of the backbone and S-type of sugar puckering are maintained in all the studied sequences, in agreement with previous literature. However, certain significant effects are detected in the IR regions sensitive to sugar pucker and glycosidic torsional angle. A single or multiple 5-methylcytosine substitution in d(CCGGCGCCGG)(2) leads to a doublet splitting of the S-type 840-820 cm(-1) sugar conformational band. The results suggest the coexistence of two different major sugar puckers within the S-conformational family, with an increased relative contribution of the C2'-endo type of sugar in the methylated sequences. In addition, a partial or full downshift of the guanosine/anti marker band at 1,375 cm(-1) in the methylated sequences reflects a change in the value of the dihedral angle chi of guanosine upon methylation. The IR spectra are interpreted in terms of localized transitions between the BI and BII subconformational states of the B-DNA backbone caused by the methylation. An increased amount of the BII subconformer in the methylated sequences should give rise to a structurally more rigid conformation, in agreement with earlier observations on DNA backbone dynamics and bending flexibility in methylated DNA.

The influence of the thymine C5 methyl group on spontaneous base pair breathing in DNA

Sequences of four or more AT base pairs without a 5'-TA-3' step, so-called A-tracts, influence the global properties of DNA by causing curvature of the helix axis if phased with the helical repeat and also influence nucleosome packaging. Hence it is interesting to understand this phenomenon on the molecular level, and numerous studies have been devoted to investigations of dynamical and structural features of A-tract DNA. It was early observed that anomalously slow base pair-opening kinetics were a striking physical property unique to DNA A-tracts (Leroy, J. L., Charretier, E., Kochoyan, M., and Gueron, M. (1988) Biochemistry 27, 8894-8898). Furthermore, a strong correlation between DNA curvature and anomalously slow base pair-opening dynamics was found. In the present work it is shown, using imino proton exchange measurements by NMR spectroscopy that the main contribution to the dampening of the base pair-opening fluctuations in A-tracts comes from the C5 methylation of the thymine base. Because the methyl group has been shown to have a very limited effect on the DNA curvature as well as the structure of the DNA helix, the thymine C5 methyl group stabilizes the helix directly. Empirical potential energy calculations show that methylation of the tract improves the stacking energy of a base pair with its neighbors in the tract by 3-4 kcal/mol.

Quantitative assessment of the use of modified nucleoside triphosphates in expression profiling: differential effects on signal intensities and impacts on expression ratios

BACKGROUND

The power of DNA microarrays derives from their ability to monitor the expression levels of many genes in parallel. One of the limitations of such powerful analytical tools is the inability to detect certain transcripts in the target sample because of artifacts caused by background noise or poor hybridization kinetics. The use of base-modified analogs of nucleoside triphosphates has been shown to increase complementary duplex stability in other applications, and here we attempted to enhance microarray hybridization signal across a wide range of sequences and expression levels by incorporating these nucleotides into labeled cRNA targets.

RESULTS

RNA samples containing 2-aminoadenosine showed increases in signal intensity for a majority of the sequences. These results were similar, and additive, to those seen with an increase in the hybridization time. In contrast, 5-methyluridine and 5-methylcytidine decreased signal intensities. Hybridization specificity, as assessed by mismatch controls, was dependent on both target sequence and extent of substitution with the modified nucleotide. Concurrent incorporation of modified and unmodified ATP in a 1:1 ratio resulted in significantly greater numbers of above-threshold ratio calls across tissues, while preserving ratio integrity and reproducibility.

CONCLUSIONS

Incorporation of 2-aminoadenosine triphosphate into cRNA targets is a promising method for increasing signal detection in microarrays. Furthermore, this approach can be optimized to minimize impact on yield of amplified material and to increase the number of expression changes that can be detected.

Dynamic impact of methylation at the M. Hhai target site: a solid-state deuterium NMR study

Base methylation plays an important role in numerous biological functions of DNA, from inhibition of cleavage by endonucleases to inhibition of transcription factor binding. Studies of nucleic acid structure have shown little differences in unmethylated DNAs and the identical sequence containing methylated analogues. We have investigated changes in the local dynamics of DNA upon substitution of a methylated cytosine analogue for cytosine using solid-state deuterium NMR. In particular, we have observed changes in the local dynamics at the target site of the M. HhaI restriction system. These studies observe changes in the amplitudes of the local backbone dynamics at the actual target site of the HhaI methyltransferase. This conclusion is another indication that the significant result of base methylation is to perturb the local dynamics, and therefore the local conformational flexibility, of the DNA helix, inhibiting or restricting the protein's ability to manipulate the DNA helix in order to perform its chemical alterations.

Impact of CpG methylation on structure, dynamics and solvation of cAMP DNA responsive element

Methylation of CpG motifs in DNA is involved in the control of gene expression and in several other epigenic effects. It suppresses also the immuno-stimulation properties of bacterial or viral DNAs that contain CPGS: However, effects of methylation on the DNA structure and dynamics are not clear. Here we carried out a 10 ns MD simulation, confronted to an NMR analysis, of a hexadecanucleotide with the cAMP responsive element (CRE) DNA methylated at its center: d(GAGATGAmCGTCATCTC)(2) (CREmet). Methylation does not introduce significant structure modification but reduces the dynamics. Molecular mechanics and generalized Born solvation energy calculations showed that the stiffness of CREmet arises from both a restriction of the conformational space by the bulky methyl groups and a folding of DNA around the hydrophobic methyls. The latter effect is favored when the GpA steps belonging to the TGA binding half-sites adopt the BII conformation. The inability of the methylated DNAs to interact with their protein partners-either transcription factors for gene regulation or a Toll-like receptor for immunostimulation-could result from both the obstacle created by methyls, preventing crucial interactions, and the loss of DNA flexibility, reducing its adaptability. Results are discussed in the light of NMR and crystallographic data.

The dynamic impact of CpG methylation in DNA

Solid-state deuterium NMR is used to investigate perturbations of the local, internal dynamics in the EcoRI restriction binding site, -GAATTC- induced by cytidine methylation. Methylation of the cytidine base in this sequence is known to suppress hydrolysis by the EcoRI restriction enzyme. Previous solid-state deuterium NMR studies have detected large amplitude motions of the phosphate-sugar backbone at the AT-CG junction of the unmethylated DNA sequence. This study shows that methylation of the cytidine base in a CpG dinucleotide reduces the amplitudes of motions of the phosphate-sugar backbone. These observations suggest a direct link between suppression of the amplitudes of localized, internal motions of the sugar-phosphate backbone of the DNA and inhibition of restriction enzyme cleavage.

An A-type double helix of DNA having B-type puckering of the deoxyribose rings

DNA usually adopts structure B in aqueous solution, while structure A is preferred in mixtures of trifluoroethanol (TFE) with water. However, the octamer d(CCCCGGGG) and other d(C(n)G(n)) fragments of DNA provide CD spectra that suggest that the base-pairs are stacked in an A-like fashion even in aqueous solution. Yet, d(CCCCGGGG) undergoes a cooperative TFE-induced transition into structure A, indicating that an important part of the aqueous duplex retains structure B. NMR spectroscopy shows that puckering of the deoxyribose rings is of the B-type. Hence, combination of the information provided by CD spectroscopy and NMR spectroscopy suggests an unprecedented double helix of DNA in which A-like base stacking is combined with B-type puckering of the deoxyribose rings. In order to determine whether this combination is possible, we used molecular dynamics to simulate the duplex of d(CCCCGGGG). Remarkably, the simulations, completely unrestrained by the experimental data, provided a very stable double helix of DNA, exhibiting just the intermediate B/A features described above. The double helix contained well-stacked guanine bases but almost unstacked cytosine bases. This generated a hole in the double helix center, which is a property characteristic for A-DNA, but absent from B-DNA. The minor groove was narrow at the double helix ends but wide at the central CG step where the Watson-Crick base-pairs were buckled in opposite directions. The base-pairs stacked tightly at the ends but stacking was loose in the duplex center. The present double helix, in which A-like base stacking is combined with B-type sugar puckering, is relevant to replication and transcription because both of these phenomena involve a local B-to-A transition.

Impact of C5-cytosine methylation on the solution structure of d(GAAAACGTTTTC)2. An NMR and molecular modelling investigation

The solution structures of d(GAAAACGTTTTC)2 and of its methylated derivative d(GAAAAMe5CGTTTTC)2 have been determined by NMR and molecular modelling in order to examine the impact of cytosine methylation on the central CpG conformation. Detailed 1H NMR and 31P NMR investigation of the two oligomers includes quantitative NOESY, 2D homonuclear Hartmann-Hahn spectroscopy, double-quantum-filtered COSY and heteronuclear 1H-31P correlation. Back-calculations of NOESY spectra and simulations of double-quantum-filtered COSY patterns were performed to gain accurate information on interproton distances and sugar phase angles. Molecular models under experimental constraints were generated by energy minimization by means of the molecular mechanics program JUMNA. The MORASS software was used to iteratively refine the structures obtained. After methylation, the oligomer still has a B-DNA conformation. However, there are differences in the structural parameters and the thermal stability as compared to the unmethylated molecule. Careful structural analysis shows that after methylation CpG departs from the usual conformation observed in other ACGT tetramers with different surroundings. Subtle displacements of bases, sugars and backbone imposed by the steric interaction of the two methyl groups inside the major groove are accompanied by severe pinching of the minor groove at the C-G residues.

Influence of pre-existing methylation on the de novo activity of eukaryotic DNA methyltransferase

Aberrant de novo methylation of CpG island DNA sequences has been observed in cultured cell lines or upon malignant transformation, but the mechanisms underlying this phenomenon are poorly understood. Using eukaryotic DNA (cytosine-5)-methyltransferase (of both human and murine origin), we have studied the in vitro methylation pattern of three CpG islands. Such sequences are intrinsically poor substrates of the enzyme, yet are efficiently methylated when a small amount of 5-methylcytosine is randomly introduced by the M.SssI prokaryotic DNA (cytosine-5)-methyltransferase prior to in vitro methylation by the eukaryotic enzyme. A stimulation was also found with several other double-stranded DNA substrates, either natural or of synthetic origin, such as poly(dG-dC).poly(dG-dC). An A + T-rich plasmid, pHb beta 1S, showed an initial stimulation, followed by a severe inhibition of the activity of DNA (cytosine-5)-methyltransferase. Methylation of poly(dI-dC).poly(dI-dC) was instead inhibited by pre-existing 5-methylcytosines. The extent of stimulation observed with poly(dG-dC).poly(dG-dC) depends on both the number and the distribution of the 5-methylcytosine residues, which probably must not be too closely spaced for the stimulatory effect to be exerted. The activity of the M.SssI prokaryotic DNA methyltransferase was not stimulated, but was inhibited by pre-methylation on either poly(dG-dC).poly(dG-dC) or poly(dI-dC).poly(dI-dC). The prokaryotic and eukaryotic DNA methyltransferases also differed in sensitivity to poly(dG-m5dC).poly(dG-m5dC), which is highly inhibitory for eukaryotic enzymes and almost ineffective on prokaryotic enzymes.

Solution structure of the CpG containing d(CTTCGAAG)2 oligonucleotide: NMR data and energy calculations are compatible with a BI/BII equilibrium at CpG

We report the analysis of the solution structure of the DNA duplex d(CTTCGAAG)2 compared to that of d(CATCGATG)2, the two oligonucleotides being related by the permutation of residues 2 and 7. An earlier study has demonstrated the malleability of CpG in the tetrad TCGA of d(CATCGATG)2 [Lefebvre et al. (1995) Biochemistry 34, 12019-12028]. Conformations of d(CTTCGAAG)2 were evaluated by (a) two-dimensional NMR, including proton and phosphorus experiments, (b) adiabatic mapping of the conformational space, (c) restrained molecular mechanics undertaken with sugar phase angle, epsilon-zeta difference angle, and NOE distances as input, and (d) back-calculation-refinement against NOE spectra at various mixing times. d(CTTCGAAG)2 like d(CATCGATG)2 exhibits a B-DNA conformation. However, significant differences are noted between the two oligonucleotides, extending up to the central CpG step, although this step resides in the same TCGA tetrad in both sequences. In structures obtained with refined NMR data, CpG adopts, for instance, a greater twist and a higher guanine phase within d(CTTCGAAG)2 compared to d(CATCGATG)2. In the former oligonucleotide, the structure of CpG resembles strikingly that found in the ACGT tetrad of the cAMP responsive element [Mauffret et al. (1992) J. Mol. Biol. 227, 852-875]. Moreover, two conformers with CpG either in the BII state (epsilon, zeta = g-, t) or in the BI state (epsilon, zeta = t, g-) are found equally stable for d(CTTCGAAG)2. The energy barrier from BI to BII comes to only 5.7 kcal/mol, and the path of the transition is very short. When calculations on d(CTTCGAAG)2 are performed taking the BI/BII equilibrium into account, the agreement with both the 1H and 31P data is found better than in the case with a single conformation taken alone. The BI/BII equilibrium may also occur in d(CATCGATG)2, but the amount of BII conformer is now found weaker compared to its analogue. The ability of the CpG phosphate groups to adopt the BII conformation could provide a satisfying explanation for the high mutation rates observed at these sites.