C-C-A-G-G-C-m5C-T-G-G. Helical fine structure, hydration, and comparison with C-C-A-G-G-C-C-T-G-G

Methylation and hydroxymethylation of cytosine alter activity and fidelity of translesion DNA polymerases

Epigenetic cytosine methylation covers most of genomic CpG dinucleotides in human cells. In addition to common deamination-mediated mutagenesis at CpG sites, an alternative deamination-independent pathway associated with DNA polymerase activity was previously described. This mutagenesis is characterized by the TCG→TTG mutational signature and is believed to arise from dAMP misincorporation opposite 5-methylcytosine (mC) or its oxidized derivative 5-hydroxymethylcytosine (hmC) by B-family replicative DNA polymerases with disrupted proofreading 3→5'-exonuclease activity. In addition to being less stable and pro-mutagenic themselves, cytosine modifications also increase the risk of adjacent nucleotides damage, including the formation of 8-oxo-2'-deoxyguanosine (8-oxoG), a well-known mutagenic lesion. The effect of cytosine methylation on error-prone DNA polymerases lacking proofreading activity and involved in repair and DNA translesion synthesis remains unexplored. Here we analyze the efficiency and fidelity of translesion Y-family polymerases (Pol κ, Pol η, Pol ι and REV1) and primase-polymerase PrimPol opposite mC and hmC as well as opposite 8-oxoG adjacent to mC in the TCG context. We demonstrate that epigenetic cytosine modifications suppress Pol ι and REV1 activities and lead to increasing dAMP misincorporation by PrimPol, Pol κ and Pol ι in vitro. Cytosine methylation also increases misincorporation of dAMP opposite the adjacent 8-oxoG by PrimPol, decreases the TLS activity of Pol η opposite the lesion but increases dCMP incorporation opposite 8-oxoG by REV1. Altogether, these data suggest that methylation and hydroxymethylation of cytosine alter activity and fidelity of translesion DNA polymerases.

Diagnostic and Prognostic Characteristics of Circulating Free DNA Methylation Detected by the Electrochemical Method in Malignant Tumors

Simple Summary

Previous studies have established an electrochemical detection method for the rapid detection of cfDNA (circulating free DNA) methylation and found that this technology might be a potential method for cancer diagnosis. However, the underlying mechanism and its role in the diagnosis and prognosis of malignant tumors are not well-characterized. In present study, we utilized the electrochemical detection method to detect the DNA methylation status by using electron microscopies and infrared spectroscopy and found that DNA with different methylated levels adsorbed to the gold surface differently, which was likely mediated by hydrophobic bonds. In addition, after detection of the cfDNA methylation status from 505 normal individuals, 725 cancer patients before treatment, and 549 patients after treatment, we found that the cfDNA adsorption rate could be used as an indicator for the diagnosis and prognosis prediction of pan-cancer. Our research provides a novel method for the liquid biopsy of cancer for diagnosis and prognosis predictions.

Abstract

Prior research has established an electrochemical method based on the differential adsorption capacity of gold surfaces with different methylated DNA degrees and found that this method might be valuable for cancer diagnosis by detecting circulating free DNA methylation. However, further investigation on the underlying mechanism and validation of its diagnostic and prognostic values in a large cohort of malignant tumors was limited. We found that DNA with different methylation levels formed particles of diverse sizes on the gold surface. Hydrophobic bonds played a significant role in the binding process of methylated DNA to the gold surface. The detection condition of an adsorption time of 10 min and temperature of 20 °C was optimal. In a large cohort of plasma samples from the patients with different malignant tumors, as well as normal individuals, we found that the electrochemical detection method based on the differential adsorption capacity of methylated DNA degree on a gold surface could be used as a noninvasive tool for malignant tumor diagnosis and prognostic evaluation. The diagnostic efficiency of this method in malignant tumors was even slightly better than that of the current tumor biomarkers widely used in routine clinical practice (circulating free DNA (cfDNA) vs. carcinoembryonic antigen (CEA), 0.8131 vs. 0.7191 and cfDNA vs. CA19-9, 0.7687 vs. 0.6693).

A Conformational Switch in the Zinc Finger Protein Kaiso Mediates Differential Readout of Specific and Methylated DNA Sequences

Recognition of the epigenetic mark 5-methylcytosine (mC) at CpG sites in DNA has emerged as a novel function of many eukaryotic transcription factors (TFs). It remains unclear why the sequence specificity of these TFs differs for CpG-methylated motifs and consensus motifs. Here, we dissect the structural and dynamic basis for this differential DNA binding specificity in the human zinc finger TF Kaiso, which exhibits high affinity for two consecutive mCpG sites in variable contexts and also for a longer, sequence-specific Kaiso binding site (KBS). By integrating structural analysis and DNA binding studies with targeted protein mutagenesis and nucleotide substitutions, we identify distinct mechanisms for readout of methylated and KBS motifs by Kaiso. We show that a key glutamate residue (E535), critical for mCpG site recognition, adopts different conformations in complexes with specific and methylated DNA. These conformational differences, together with intrinsic variations in DNA flexibility and/or solvation at TpG versus mCpG sites, contribute to the different DNA affinity and sequence specificity. With methylated DNA, multiple direct contacts between E535 and the 5' mCpG site dominate the binding affinity, allowing for tolerance of different flanking DNA sequences. With KBS, Kaiso employs E535 as part of an indirect screen of the 5' flanking sequence, relying on key tyrosine-DNA interactions to stabilize an optimal DNA conformation and select against noncognate sites. These findings demonstrate how TFs use conformational adaptation and exploit variations in DNA flexibility to achieve distinct DNA readout outcomes and target a greater variety of regulatory and epigenetic sites than previously appreciated.

Symmetry in Nucleic-Acid Double Helices

In nature and in the test tube, nucleic acids occur in many different forms. Apart from single-stranded, coiled molecules, DNA and RNA prefer to form helical arrangements, in which the bases are stacked to shield their hydrophobic surfaces and expose their polar edges. Focusing on double helices, we describe the crucial role played by symmetry in shaping DNA and RNA structure. The base pairs in nucleic-acid double helices display rotational pseudo-symmetry. In the Watson–Crick base pairs found in naturally occurring DNA and RNA duplexes, the symmetry axis lies in the base-pair plane, giving rise to two different helical grooves. In contrast, anti-Watson–Crick base pairs have a dyad axis perpendicular to the base-pair plane and identical grooves. In combination with the base-pair symmetry, the syn/anti conformation of paired nucleotides determines the parallel or antiparallel strand orientation of double helices. DNA and RNA duplexes in nature are exclusively antiparallel. Watson–Crick base-paired DNA or RNA helices display either right-handed or left-handed helical (pseudo-) symmetry. Genomic DNA is usually in the right-handed B-form, and RNA double helices adopt the right-handed A-conformation. Finally, there is a higher level of helical symmetry in superhelical DNA in which B-form double strands are intertwined in a right- or left-handed sense.

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.

Unmethylated and methylated CpG dinucleotides distinctively regulate the physical properties of DNA

In eukaryotic cells, DNA has to bend significantly to pack inside the nucleus. Physical properties of DNA such as bending flexibility and curvature are expected to affect DNA packaging and partially determine the nucleosome positioning patterns inside a cell. DNA CpG methylation, the most common epigenetic modification found in DNA, is known to affect the physical properties of DNA. However, its detailed role in nucleosome formation is less well-established. In this study, we evaluated the effect of defined CpG patterns (unmethylated and methylated) on DNA structure and their respective nucleosome-forming ability. Our results suggest that the addition of CpG dinucleotides, either as a (CG)n stretch or (CGX8 )n repeats at 10 bp intervals, lead to reduced hydrodynamic radius and decreased nucleosome-forming ability of DNA. This effect is more predominant for a DNA stretch ((CG)5) located in the middle of a DNA fragment. Methylation of CpG sites, surprisingly, seems to reduce the difference in DNA structure and nucleosome-forming ability among DNA constructs with different CpG patterns. Our results suggest that unmethylated and methylated CpG patterns can play very different roles in regulating the physical properties of DNA. CpG methylation seems to reduce the DNA conformational variations affiliated with defined CpG patterns. Our results can have significant bearings in understanding the nucleosome positioning pattern in living organisms modulated by DNA sequences and epigenetic features.

DNA methylation effects on tetra-nucleosome compaction and aggregation

DNA CpG methylation has been associated with chromatin compaction and gene silencing. Whether DNA methylation directly contributes to chromatin compaction remains an open question. In this study, we used fluorescence fluctuation spectroscopy (FFS) to evaluate the compaction and aggregation of tetra-nucleosomes containing specific CpG patterns and methylation levels. The compactness of both unmethylated and methylated tetra-nucleosomes is dependent on DNA sequences. Specifically, methylation of the CpG sites located in the central dyad and the major grooves of DNA seem to have opposite effects on modulating the compactness of tetra-nucleosomes. The interactions among tetra-nucleosomes, however, seem to be enhanced because of DNA methylation independent of sequence contexts. Our finding can shed light on understanding the role of DNA methylation in determining nucleosome positioning pattern and chromatin compactness.

Epi-fingerprinting and epi-interventions for improved crop production and food quality

Increasing crop production at a time of rapid climate change represents the greatest challenge facing contemporary agricultural research. Our understanding of the genetic control of yield derives from controlled field experiments designed to minimize environmental variance. In spite of these efforts there is substantial residual variability among plants attributable to Genotype × Environment interactions. Recent advances in the field of epigenetics have revealed a plethora of gene control mechanisms that could account for much of this unassigned variation. These systems act as a regulatory interface between the perception of the environment and associated alterations in gene expression. Direct intervention of epigenetic control systems hold the enticing promise of creating new sources of variability that could enhance crop performance. Equally, understanding the relationship between various epigenetic states and responses of the crop to specific aspects of the growing environment (epigenetic fingerprinting) could allow for a more tailored approach to plant agronomy. In this review, we explore the many ways in which epigenetic interventions and epigenetic fingerprinting can be deployed for the improvement of crop production and quality.

Crystal structure of interleukin-6 in complex with a modified nucleic acid ligand

IL-6 is a secreted cytokine that functions through binding two cell surface receptors, IL-6Rα and gp130. Because of its involvement in the progression of several chronic inflammatory diseases, IL-6 is a target of pharmacologic interest. We have recently identified a novel class of ligands called SOMAmers (S low Off-rate Modified Aptamers) that bind IL-6 and inhibit its biologic activity. SOMAmers exploit the chemical diversity of protein-like side chains assembled on flexible nucleic acid scaffolds, resulting in an expanded repertoire of intra- and intermolecular interactions not achievable with conventional aptamers. Here, we report the co-crystal structure of a high affinity SOMAmer (K_d = 0.20 nm) modified at the 5-position of deoxyuridine in a complex with IL-6. The SOMAmer, comprised of a G-quartet domain and a stem-loop domain, engages IL-6 in a clamp-like manner over an extended surface exhibiting close shape complementarity with the protein. The interface is characterized by substantial hydrophobic interactions overlapping the binding surfaces of the IL-6Rα and gp130 receptors. The G-quartet domain retains considerable binding activity as a disconnected autonomous fragment (K_d = 270 nm). A single substitution from our diversely modified nucleotide library leads to a 37-fold enhancement in binding affinity of the G-quartet fragment (K_d = 7.4 nm). The ability to probe ligand surfaces in this manner is a powerful tool in the development of new therapeutic reagents with improved pharmacologic properties. The SOMAmer·IL-6 structure also expands our understanding of the diverse structural motifs achievable with modified nucleic acid libraries and elucidates the nature with which these unique ligands interact with their protein targets.

Structural insights into how 5-hydroxymethylation influences transcription factor binding

Transcription factor binding and high resolution crystallographic studies (1.3 Å) of Dickerson-Drew duplexes with cytosine, methylcytosine and hydroxymethylcytosine bases provide evidence that C-5 cytosine modifications could regulate transcription by context dependent effects on DNA transcription factor interactions.

Differential effect of three base modifications on DNA thermostability revealed by high resolution melting

High resolution melting (HRM) can detect and quantify the presence of 5-methylcytosine (5mC) in DNA samples, but the ability of HRM to diagnose other DNA modifications remains unexplored. The DNA bases N6-methyladenine and 5-hydroxymethylcytosine occur across almost all phyla. While their function remains controversial, their presence perturbs DNA structure. Such modifications could affect gene regulation, chromatin condensation and DNA packaging. Here, we reveal that DNA containing N6-methyladenine or 5-hydroxymethylcytosine exhibits reduced thermal stability compared to cytosine-methylated DNA. These thermostability changes are sufficiently divergent to allow detection and quantification by HRM analysis. Thus, we report that HRM distinguishes between sequence-identical DNA differing only in the modification type of one base. This approach is also able to distinguish between two DNA fragments carrying both N6-methyladenine and 5-methylcytosine but differing only in the distance separating the modified bases. This finding provides scope for the development of new methods to characterize DNA chemically and to allow for low cost screening of mutant populations of genes involved in base modification. More fundamentally, contrast between the thermostabilizing effects of 5mC on dsDNA compared with the destabilizing effects of N6-methyladenine (m6A) and 5-hydroxymethylcytosine (5hmC) raises the intriguing possibility of an antagonistic relationship between modification types with functional significance.

B-DNA structure is intrinsically polymorphic: even at the level of base pair positions

Increasingly exact measurement of single crystal X-ray diffraction data offers detailed characterization of DNA conformation, hydration and electrostatics. However, instead of providing a more clear and unambiguous image of DNA, highly accurate diffraction data reveal polymorphism of the DNA atomic positions and conformation and hydration. Here we describe an accurate X-ray structure of B-DNA, painstakingly fit to a multistate model that contains multiple competing positions of most of the backbone and of entire base pairs. Two of ten base-pairs of CCAGGCCTGG are in multiple states distinguished primarily by differences in slide. Similarly, all the surrounding ions are seen to fractionally occupy discrete competing and overlapping sites. And finally, the vast majority of water molecules show strong evidence of multiple competing sites. Conventional resolution appears to give a false sense of homogeneity in conformation and interactions of DNA. In addition, conventional resolution yields an average structure that is not accurate, in that it is different from any of the multiple discrete structures observed at high resolution. Because base pair positional heterogeneity has not always been incorporated into model-building, even some high and ultrahigh-resolution structures of DNA do not indicate the full extent of conformational polymorphism.

Geometric properties of nucleic acids with potential for autobuilding

Algorithms and geometrical properties are described for the automated building of nucleic acids in experimental electron density.

Medium- to high-resolution X-ray structures of DNA and RNA molecules were investigated to find geometric properties useful for automated model building in crystallographic electron-density maps. We describe a simple method, starting from a list of electron-density ‘blobs’, for identifying backbone phosphates and nucleic acid bases based on properties of the local electron-density distribution. This knowledge should be useful for the automated building of nucleic acid models into electron-density maps. We show that the distances and angles involving C1′ and the P atoms, using the pseudo-torsion angles and that describe the …P—C1′—P—C1′… chain, provide a promising basis for building the nucleic acid polymer. These quantities show reasonably narrow distributions with asymmetry that should allow the direction of the phosphate backbone to be established.

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.

On the molecular basis of uracil recognition in DNA: comparative study of T-A versus U-A structure, dynamics and open base pair kinetics

Uracil (U) can be found in DNA as a mismatch paired either to adenine (A) or to guanine (G). Removal of U from DNA is performed by a class of enzymes known as uracil–DNA–glycosylases (UDG). Recent studies suggest that recognition of U–A and U–G mismatches by UDG takes place via an extra-helical mechanism. In this work, we use molecular dynamics simulations to analyze the structure, dynamics and open base pair kinetics of U–A base pairs relative to their natural T–A counterpart in 12 dodecamers. Our results show that the presence of U does not alter the local conformation of B-DNA. Breathing dynamics and base pair closing kinetics are only weakly dependent on the presence of U versus T, with open T–A and U–A pairs lifetimes in the nanosecond timescale. Additionally, we observed spontaneous base flipping in U–A pairs. We analyze the structure and dynamics for this event and compare the results to available crystallographic data of open base pair conformations. Our results are in agreement with both structural and kinetic data derived from NMR imino proton exchange measurements, providing the first detailed description at the molecular level of elusive events such as spontaneous base pair opening and flipping in mismatched U–A sequences in DNA. Based on these results, we propose that base pair flipping can occur spontaneously at room temperature via a 3-step mechanism with an open base pair intermediate. Implications for the molecular basis of U recognition by UDG are discussed.

Molecular diagnostics for personal medicine using a nanopore

Semiconductor nanotechnology has created the ultimate analytical tool: a nanopore with single molecule sensitivity. This tool offers the intriguing possibility of high-throughput, low cost sequencing of DNA with the absolute minimum of material and preprocessing. The exquisite single molecule sensitivity obviates the need for costly and error-prone procedures like polymerase chain reaction amplification. Instead, nanopore sequencing relies on the electric signal that develops when a DNA molecule translocates through a pore in a membrane. If each base pair has a characteristic electrical signature, then ostensibly a pore could be used to analyze the sequence by reporting all of the signatures in a single read without resorting to multiple DNA copies. The potential for a long read length combined with high translocation velocity should make resequencing inexpensive and allow for haplotyping and methylation profiling in a chromosome.

Conformations of the sugar-phosphate backbone in helical DNA crystal structures

The DNA backbone geometry was analyzed for 96 crystal structures of oligodeoxynucleotides. The ranges and mean values of the torsion angles for the observed subclasses of the A-, B-, and Z-DNA conformational types were determined by analyzing distributions of the torsion angles and scattergrams relating pairs of angles.

Z-DNA crystallography

We review the effect of sequence on the structure of left-handed Z-DNA in single crystals. The various substituent groups that define a nucleotide base as guanine, cytosine,thymine, or adenine affect both the DNA conformation and the organization of solvent around the duplex. These are discussed in terms of their effect on the ability of sequences to adopt the unusual Z-DNA structure. In addition, the experimental and theoretical methods used to treat DNA hydration are discussed as they relate to the stability of Z-DNA . Finally, we argue that Z-DNA , as defined by the crystal conformation, is sufficient in itself to account for the physical properties of left-handed conformations observed in polymers and in genomic sequences

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.

ROS1 5-methylcytosine DNA glycosylase is a slow-turnover catalyst that initiates DNA demethylation in a distributive fashion

Arabidopsis ROS1 belongs to a family of plant 5-methycytosine DNA glycosylases that initiate DNA demethylation through base excision. ROS1 displays the remarkable capacity to excise 5-meC, and to a lesser extent T, while retaining the ability to discriminate effectively against C and U. We found that replacement of the C5-methyl group by halogen substituents greatly decreased excision of the target base. Furthermore, 5-meC was excised more efficiently from mismatches, whereas excision of T only occurred when mispaired with G. These results suggest that ROS1 specificity arises by a combination of selective recognition at the active site and thermodynamic stability of the target base. We also found that ROS1 is a low-turnover catalyst because it binds tightly to the abasic site left after 5-meC removal. This binding leads to a highly distributive behaviour of the enzyme on DNA substrates containing multiple 5-meC residues, and may help to avoid generation of double-strand breaks during processing of bimethylated CG dinucleotides. We conclude that the biochemical properties of ROS1 are consistent with its proposed role in protecting the plant genome from excess methylation.

Nanoelectromechanics of methylated DNA in a synthetic nanopore

Methylation of cytosine is a covalent modification of DNA that can be used to silence genes, orchestrating a myriad of biological processes including cancer. We have discovered that a synthetic nanopore in a membrane comparable in thickness to a protein binding site can be used to detect methylation. We observe a voltage threshold for permeation of methylated DNA through a <2 nm diameter pore, which we attribute to the stretching transition; this can differ by >1 V/20 nm depending on the methylation level, but not the DNA sequence.

Methylated DNA recognition during the reversal of epigenetic silencing is regulated by cysteine and serine residues in the Epstein-Barr virus lytic switch protein

Epstein-Barr virus (EBV) causes infectious mononucleosis and is associated with various malignancies, including Burkitt's lymphoma and nasopharyngeal carcinoma. Like all herpesviruses, the EBV life cycle alternates between latency and lytic replication. During latency, the viral genome is largely silenced by host-driven methylation of CpG motifs and, in the switch to the lytic cycle, this epigenetic silencing is overturned. A key event is the activation of the viral BRLF1 gene by the immediate-early protein Zta. Zta is a bZIP transcription factor that preferentially binds to specific response elements (ZREs) in the BRLF1 promoter (Rp) when these elements are methylated. Zta's ability to trigger lytic cycle activation is severely compromised when a cysteine residue in its bZIP domain is mutated to serine (C189S), but the molecular basis for this effect is unknown. Here we show that the C189S mutant is defective for activating Rp in a Burkitt's lymphoma cell line. The mutant is compromised both in vitro and in vivo for binding two methylated ZREs in Rp (ZRE2 and ZRE3), although the effect is striking only for ZRE3. Molecular modeling of Zta bound to methylated ZRE3, together with biochemical data, indicate that C189 directly contacts one of the two methyl cytosines within a specific CpG motif. The motif's second methyl cytosine (on the complementary DNA strand) is predicted to contact S186, a residue known to regulate methyl-ZRE recognition. Our results suggest that C189 regulates the enhanced interaction of Zta with methylated DNA in overturning the epigenetic control of viral latency. As C189 is conserved in many bZIP proteins, the selectivity of Zta for methylated DNA may be a paradigm for a more general phenomenon.

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.

Methylation of cytosine at C5 in a CpG sequence context causes a conformational switch of a benzo[a]pyrene diol epoxide-N2-guanine adduct in DNA from a minor groove alignment to intercalation with base displacement

It is well known that CpG dinucleotide steps in DNA, which are highly methylated at the 5-position of cytosine (meC) in human tissues, exhibit a disproportionate number of mutations within certain codons of the p53 gene. There is ample published evidence indicating that the reactivity of guanine with anti-B[a]PDE (a metabolite of the environmental carcinogen benzo[a]pyrene) at CpG mutation hot spots is enhanced by the methylation of the cytosine residue flanking the target guanine residue on the 5'-side. In this work we demonstrate that such a methylation can also dramatically affect the conformational characteristics of an adduct derived from the reaction of one of the two enantiomers of anti-B[a]PDE with the exocyclic amino group of guanine ([BP]G adduct). A detailed NMR study indicates that the 10R (-)-trans-anti-[BP]G adduct undergoes a transition from a minor groove-binding alignment of the aromatic BP ring system in the unmethylated C-[BP]G sequence context, to an intercalative BP alignment with a concomitant displacement of the modified guanine residue into the minor groove in the methylated meC-[BP]G sequence context. By contrast, a minor groove-binding alignment was observed for the stereoisomeric 10S (+)-trans-anti-[BP]G adduct in both the C-[BP]G and meC-[BP]G sequence contexts. This remarkable conformational switch resulting from the presence of a single methyl group at the 5-position of the cytosine residue flanking the lesion on the 5'-side, is attributed to the hydrophobic effect of the methyl group that can stabilize intercalated adduct conformations in an adduct stereochemistry-dependent manner. Such conformational differences in methylated and unmethylated CpG sequences may be significant because of potential alterations in the cellular processing of the [BP]G adducts by DNA transcription, replication, and repair enzymes.

How strong can the bend be on a DNA helix from cisplatin? DFT and MP2 quantum chemical calculations of cisplatin-bridged DNA purine bases

The B3LYP/6-31G(d) level of theory was used for the optimization of [Pt(NH(3))(4)](2+), [Pt(NH(3))(3)(H(2)O)](2+), cis-[Pt(NH(3))(2)(H(2)O)(2)](2+), and related platinum complexes. In addition, water or ammonium ligands were replaced by DNA purine bases so that finally cis-diammineplatinum with two bases (Pt-bridged complexes) is obtained. Single point calculations using the MP2/6-31+G(d) method were performed on the obtained reference geometries and were utilized for estimating bond dissociation energies (BDEs) and stabilization energies, and for electron density analyses. After reoptimization, IR spectra were determined from HF second derivatives. It was found that replacement of both water and ammonium by the DNA base is an exothermic process (20-50 kcal/mol depending on the ligands present in the complex). Asymmetric structures with one interbase H-bond were obtained for cis-diammine[bond](N(7),N(7)'-diadenine)[bond]platinum and mixed cis-diammine[bond](N(7)-adenine)[bond](N(7)-guanine)[bond]platinum complexes. In the case of the diguanine Pt-bridge, a symmetrical complex with two ammonium...O(6) H-bonds was found. The higher stabilization energy of the di-guanine complex is linked to a larger component of the Coulombic interaction. However, the BDE of Pt[bond]N(7)(G) is smaller in this complex than the BDE of Pt[bond]N(7)(G) from the mixed Pt[bond]AG complex. Also, steric repulsion of the ligands is about 10 kcal/mol smaller for the asymmetrical Pt[bond]AA and Pt[bond]AG bridges. The influence of the trans effect on DBE can be clearly seen. Adenine exhibits the largest trans effect, followed by guanine, ammonium, and water. The strength of the H-bond can be determined from the IR spectra. The strongest H-bond is the interbase H-bridge between adenine and guanine in the mixed Pt[bond]AG complex; otherwise, the H-bonds of adenine complexes are weaker than in guanine complexes. BDE can be traced in the guanine-containing complexes. The nature of the covalent bonding is analyzed in terms of partial charges and MO. A general explanation of the lower affinity of transition metals to oxygen than nitrogen can be partially seen in the less favorable geometrical orientation of lone electron pairs of oxygen.

Endogenous 5-methylcytosine protects neighboring guanines from N7 and O6-methylation and O6-pyridyloxobutylation by the tobacco carcinogen 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone

All CG dinucleotides along exons 5-8 of the p53 tumor suppressor gene contain endogenous 5-methylcytosine (MeC). These same sites (e.g., codons 157, 158, 245, 248, and 273) are mutational hot spots in smoking-induced lung cancer. Several groups used the UvrABC endonuclease incision assay to demonstrate that methylated CG dinucleotides of the p53 gene are the preferred binding sites for the diol epoxides of bay region polycyclic aromatic hydrocarbons (PAH). In contrast, effects of endogenous cytosine methylation on the distribution of DNA lesions induced by tobacco-specific nitrosamines, e.g., 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK), have not been elucidated. In the work presented here, a stable isotope labeling HPLC-ESI-MS/MS approach was employed to analyze the reactivity of the N7 and O6 positions of guanines within hemimethylated and fully methylated CG dinucleotides toward NNK-derived methylating and pyridyloxobutylating species. 15N3-labeled guanine bases were placed within synthetic DNA sequences representing endogenously methylated p53 codons 154, 157, and 248, followed by treatment with acetylated precursors to NNK diazohydroxides. HPLC-ESI-MS/MS analysis was used to determine the relative yields of N7- and O6-guanine adducts at the 15N3-labeled position. In all cases, the presence of MeC inhibited the formation of N7-methylguanine, O6-methylguanine, and O6-pyridyloxobutylguanine at a neighboring G, with the greatest decrease observed in fully methylated dinucleotides and at guanines preceded by MeC. Furthermore, the O6-Me-dG/N7-Me-G molar ratios were decreased in the presence of the 5'-neighboring MeC, suggesting that the observed decline in O6-alkylguanine adduct yields is, at least partially, a result of an altered reactivity pattern in methylated CG dinucleotides. These results indicate that, unlike N2-guanine adducts of PAH diol epoxides, NNK-induced N7- and O6-alkylguanine adducts are not preferentially formed at the endogenously methylated CG sites within the p53 tumor suppressor gene.

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.

A-tract DNA disfavours triplex formation

Optimisation of DNA triplex stability is of fundamental importance in the anti-gene strategy. In the present work, thermal denaturation studies by UV-spectrophotometry and structural and dynamical characterizations by NMR spectroscopy have been used systematically to investigate the effects on triplex stability of isolated insertions of different base triplets into an otherwise homogeneous 15-mer dT x dA-dT oligo-triplex. It is found that insertion of a single central C(+) x G-C or T x D-T triplet (D=2,6-diaminopurine) leads to a pronounced stabilization (up to 20 deg. C if the cytosine base is C5 methylated) at acidic as well as neutral pH. To a smaller degree, this is the case also for a C(+) x I-C triplet insertion. Using imino proton exchange measurements, it is shown that insertion of a DT base-pair in the underlying duplex perturbs the intrinsic A-tract structure in the same way as has been shown for a GC insert. We propose that the intrinsic properties of A-tract duplex DNA (e. g. high propeller twist and rigidity) are unfavourable for triplex formation and that GC- or DT-inserts stabilize the triplex by interfering with the A-tract features of the underlying duplex. The C(+) x I-C triplet without the N2 amino group in the minor groove is readily accommodated within the typical, highly propeller-twisted A-tract structure. This might be related to its smaller effect on the stability of the corresponding triplex. These results may be valuable for understanding DNA triplex formation in vivo as well as for the design of efficient triplex-forming oligonucleotides and in choosing suitable target sequences in the anti-gene strategy.

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.

Influence of the dynamic positions of cations on the structure of the DNA minor groove: sequence-dependent effects

Different models for minor groove structures predict that the conformation is essentially fixed by sequence and has an influence on local ion distribution or alternatively that temporal positions of ions around the minor groove can affect the structure if they neutralize cross-strand phosphate charges. Our previous studies show that the minor groove in an AATT dodecamer responds to local sodium ion positions and is narrow when ions neutralize cross-strand phosphate-phosphate charges [J. Am. Chem. Soc. 2000, 122, 10513-10520]. Previous results from a number of laboratories have shown that G-tracts often have a wider minor groove than A-tracts, but they do not indicate whether this is due to reduced flexibility or differences in ion interactions. We have undertaken a molecular dynamics study of a d(TATAGGCCTATA) duplex to answer this question. The results show that the G-tract has the same amplitude of minor groove fluctuations as the A-tract sequence but that it has fewer ion interactions that neutralize cross-strand phosphate charges. These results demonstrate that differences in time-average groove width between A- and G-tracts are due to differences in ion interactions at the minor groove. When ions neutralize the cross-strand phosphates, the minor groove is narrow. When there are no neutralizing ion interactions, the minor groove is wide. The population of structures with no ion interactions is larger with the GGCC than with the AATT duplex, and GGCC has a wider time-average minor groove in agreement with experiment.

DNA-interactions and nuclear localisation of the chromosomal HMG domain protein SSRP1 from maize

The structure-specific recognition protein 1 (SSRP1) is a member of the protein family containing a high mobility group (HMG) domain DNA-binding motif. We have functionally characterised the 71.4 kDa Zm-SSRP1 protein from maize. The chromatin-associated Zm-SSRP1 is detected by immunoblot analysis in maize leaves, kernels and suspension culture cells, but not in roots. Mediated by its HMG domain, recombinant Zm-SSRP1 interacts structure-specifically with supercoiled DNA and DNA minicircles when compared with linear DNA. In linear duplex DNA, the protein does not recognise a specific sequence, but it binds preferentially to sequences containing the deformable dinucleotide TG, as demonstrated by a random oligonucleotide selection experiment. Zm-SSRP1 modulates DNA structure by bending the target sequence, since it promotes the circularisation of short DNA fragments in the presence of DNA ligase. Moreover, Zm-SSRP1 facilitates the formation of nucleoprotein structures, as measured using the bacterial site-specific beta-mediated recombination reaction. Analysis of the subcellular localisation of various SSRP1-GFP fusions revealed that, in contrast to HMG domain transcription factors, the nuclear localisation sequence of Zm-SSRP1 is situated within a 20-amino acid residue region adjacent to the HMG domain rather than within the DNA-binding domain. The results are discussed in the context of the likely function of SSRP1 proteins in transcription and replication.

Effect of neighboring bases on base-pair stacking orientation: a molecular dynamics study

It is generally believed that base-pair stacking interaction in DNA double helix is one of the strongest interactions that governs sequence directed structural variability. However, X-ray crystal structures of some base-paired doublet sequences have been seen to adopt different structures when flanked by different base-pairs. DNA crystal database, however, is still too small to make good statistical inference about effect of such flanking residues. Influence of neighboring residue on the local helical geometry of a base-paired doublet in B-DNA has been investigated here using molecular dynamics simulation. We have generated ensembles of structures for d(CA).d(TG) and d(AA).d(TT) base-paired doublets located at the centers of d(CGCGCAAAGCG).d(CGCTTTGCGCG) and d(CGCGAAAACGCG).d(CGCGTTTTCGCG) sequences along with their analogs by varying the bases either at 5'- or 3'- position to the central doublet. Comparison of base paired doublet parameters for the ensembles of structures show that stacking geometry of d(CA).d(TG) doublet depends on some of the flanking base-pairs. On the other hand d(AA).d(TT) doublet remains nearly unperturbed when the flanking residues are altered.

Crystal studies of B-DNA: the answers and the questions

The history of DNA crystallography is reviewed and is followed by discussion of the methods used for structure determination. The features of B-DNA molecular and crystal structures are described.

C5 cytosine methylation at CpG sites enhances sequence selectivity of mitomycin C-DNA bonding

We have established that UvrABC nuclease is equally efficient in cutting mitomycin C (MC)-DNA monoadducts formed at different sequences and that the degree of UvrABC cutting represents the extent of drug-DNA bonding. Using this method we determined the effect of C5 cytosine methylation on the DNA monoalkylation by MC and the related analogues N-methyl-7-methoxyaziridinomitosene (MS-NMA) and 10-decarbamoylmitomycin C (DC-MC). We have found that C5 cytosine methylation at CpG sites greatly enhances MC and MS-NMA DNA adduct formation at those sites while reducing adduct formation at non-CpG sequences. In contrast, although DC-MC DNA bonding at CpG sites is greatly enhanced by CpG methylation, its bonding at non-CpG sequences is not appreciably affected. These cumulative results suggest that C5 cytosine methylation at CpG sites enhances sequence selectivity of drug-DNA bonding. We propose that the methylation pattern and status (hypo- or hypermethylation) of genomic DNA may determine the cells' susceptibility to MC and its analogues, and these effects may, in turn, play a crucial role in the antitumor activities of the drugs.

Effects of 5-[S-(2,4-dinitrophenyl)-thio]-2'-deoxyuridine analog incorporation on the structure and stability of DNA hybrids: implications for the design of nucleic acid probes

Labeled nucleic acid probes are used as diagnostic tools by detecting changes in gene expression upon hybridization to target RNAs or DNAs that are related to specific disease genes. 5-[S-(2, 4-Dinitrophenyl)-thio]-2'-deoxyuridine analog represents an excellent nucleic acid label, containing the DNP group which functions both as a probe and as a precursor for the introduction of a variety of fluorescent groups. This study describes thermal denaturation hybridization experiments with oligonucleotides containing the 5-[S-(2,4-dinitrophenyl)-thio]-2'-deoxyuridine analog. Using molecular modeling techniques, the effects of this analog on the hybrid structure and stability were examined, including (i) analog conformation, (ii) hydrogen bonding, (iii) stacking interactions and (iv) hybrid helical geometry. This analog does not prohibitively affect the hybrid thermal stability and incorporation of the analog does not compromise the structural integrity of the double helix. In particular, the sequence-dependence of the analog effects and the dependence on the modification site relative to the end(s) of the helix were investigated. Findings described here should provide guidelines in the rational design of nucleic acid probes.

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.

Cytosine methylation in a CpG sequence leads to enhanced reactivity with Benzo[a]pyrene diol epoxide that correlates with a conformational change

Benzo[a]pyrene (B[a]P) is a widespread environmental carcinogen that must be activated by cellular metabolism to a diol epoxide form (BPDE) before it reacts with DNA. It has recently been shown that BPDE preferentially modifies the guanine in methylated 5'-CpG-3' sequences in the human p53 gene, providing one explanation for why these sites are mutational hot spots. Using purified duplex oligonucleotides containing identical methylated and unmethylated CpG sequences, we show here that BPDE preferentially modified the guanine in hemimethylated or fully methylated CpG sequences, producing between 3- and 8-fold more modification at this site. Analysis of this reaction using shorter duplex oligonucleotides indicated that it was the level of the (+)-trans isomer that was specifically increased. To determine if there were conformational differences between the methylated and unmethylated B[a]P-modified DNA sequences that may be responsible for this enhanced reactivity, a native polyacrylamide gel electrophoresis analysis was carried out using DNA containing isomerically pure B[a]P-DNA adducts. These experiments showed that each adduct resulted in an altered gel mobility in duplex DNA but that only the presence of a (+)-trans isomer and a methylated C 5' to the adduct resulted in a significant gel mobility shift compared with the unmethylated case.

Atomic-resolution crystal structures of B-DNA reveal specific influences of divalent metal ions on conformation and packing

Crystal structures of B-form DNA have provided insights into the global and local conformational properties of the double helix, the solvent environment, drug binding and DNA packing. For example, structures of the duplex with sequence CGCGAATTCGCG, the Dickerson-Drew dodecamer (DDD), established a unique geometry of the central A-tract and a hydration spine in the minor groove. However, our knowledge of the various interaction modes between metal ions and DNA is very limited and almost no information exists concerning the origins of the different effects on DNA conformation and packing exerted by individual metal ions. Crystallization of the DDD duplex in the presence of Mg(2+)and Ca(2+)yields different crystal forms. The structures of the new Ca(2+)-form and isomorphous structures of oligonucleotides with sequences GGCGAATTCGCG and GCGAATTCGCG were determined at a maximum resolution of 1.3 A. These and the 1.1 A structure of the DDD Mg(2+)-form have revealed the most detailed picture yet of the ionic environment of B-DNA. In the Mg(2+)and Ca(2+)-forms, duplexes in the crystal lattice are surrounded by 13 magnesium and 11 calcium ions, respectively.Mg(2+)and Ca(2+)generate different DNA crystal lattices and stabilize different end-to-end overlaps and lateral contacts between duplexes, thus using different strategies for reducing the effective repeat length of the helix to ten base-pairs. Mg(2+)crystals allow the two outermost base-pairs at either end to interact laterally via minor groove H-bonds, turning the 12-mer into an effective 10-mer. Ca(2+)crystals, in contrast, unpair the outermost base-pair at each end, converting the helix into a 10-mer that can stack along its axis. This reduction of a 12-mer into a functional 10-mer is followed no matter what the detailed nature of the 5'-end of the chain: C-G-C-G-A-ellipsis, G-G-C-G-A-ellipsis, or a truncated G-C-G-A-ellipsis Rather than merely mediating close contacts between phosphate groups, ions are at the origin of many well-known features of the DDD duplex structure. A Mg(2+)coordinates in the major groove, contributing to kinking of the duplex at one end. While Ca(2+)resides in the minor groove, coordinating to bases via its hydration shell, two magnesium ions are located at the periphery of the minor groove, bridging phosphate groups from opposite strands and contracting the groove at one border of the A-tract.

Reactivity of guanine at m5CpG steps in DNA: evidence for electronic effects transmitted through the base pairs

BACKGROUND

Mitomycin C (MC), a DNA cross-linking and alkylating agent, targets guanines in the m5CpG sequence with 2-3-fold preference over guanines in unmethylated CpG. Benzo[a]pyrenediolepoxide (BPDE) and several other aromatic carcinogens form guanine adducts with an identical selectivity for m5CpG, and in certain cancers G to T transversion mutation 'hotspots' in the p53 tumor suppressor gene are more frequent at this sequence than at guanines in other sequences. MC appears suitable to probe the general mechanism of this selectivity.

RESULTS

A 162-bp DNA fragment containing C, m5C or f5C (5-fluoro cytosine) at all cytosine positions was cross-linked by MC at guanines in CpG steps. The extent of cross-linking increased in the order f5C < C < m5C. Monoalkylation or cross-linking of duplex 12-mer oligonucleotides containing a single CpG, f5CpG or m5CpG step gave yields of adducts that increased in the same order. The rates showed a correlation with the Hammett sigma constant of the methyl and fluoro substituents of the cytosine. Only the base-pair cytosine substituent influenced reactivity of guanine.

CONCLUSIONS

The 2-amino group of guanine in the m5CpG sequence of DNA has a greater nucleophilic reactivity with mitomycin than CpG. Evidence is presented for a novel mechanism: transmission of the electron-donating effect of the 5-methyl substituent of the cytosine to guanine through H-bonding of the m5C.G base pair. The results explain the enhanced reaction of BPDE at m5CpG in DNA and the origin of G-T mutational hotspots in the p53 gene in cancer.

High base pair opening rates in tracts of GC base pairs

Sequence-dependent structural features of the DNA double helix have a strong influence on the base pair opening dynamics. Here we report a detailed study of the kinetics of base pair breathing in tracts of GC base pairs in DNA duplexes derived from 1H NMR measurements of the imino proton exchange rates upon titration with the exchange catalyst ammonia. In the limit of infinite exchange catalyst concentration, the exchange times of the guanine imino protons of the GC tracts extrapolate to much shorter base pair lifetimes than commonly observed for isolated GC base pairs. The base pair lifetimes in the GC tracts are below 5 ms for almost all of the base pairs. The unusually rapid base pair opening dynamics of GC tracts are in striking contrast to the behavior of AT tracts, where very long base pair lifetimes are observed. The implication of these findings for the structural principles governing spontaneous helix opening as well as the DNA-binding specificity of the cytosine-5-methyltransferases, where flipping of the cytosine base has been observed, are discussed.

Conversions of the left-handed form and the protonated form of DNA back to the bound right-handed form by sanguinarine and ethidium: A comparative study

The interaction of sanguinarine and ethidium with right-handed (B-form), left-handed (Z-form) and left-handed protonated (designated as H(L)-form) structures of poly(dG-dC).poly(dG-dC) and poly(dG-me5dC).poly(dG-me5dC) was investigated by measuring the circular dichroism and UV absorption spectral analysis. Both sanguinarine and ethidium bind strongly to the B-form DNA and convert the Z-form and the H(L)-form back to the bound right-handed form. Circular dichroic data also show that the conformation at the binding site is right-handed, even though adjacent regions of the polymer have a left-handed conformation either in Z-form or in H(L)-form. Both the rate and extent of B-form to Z-form transition were decreased by sanguinarine and ethidium under ionic conditions that otherwise favour the left-handed conformation of the polynucleotides. The rate of decrease is faster in the case of ethidium as compared to that of sanguinarine. Scatchard analysis of the spectrophotometric data shows that sanguinarine binds strongly to both the polynucleotides in a non-cooperative manner under B-form conditions, in sharp contrast to the highly-cooperative binding under Z-form and H(L)-form conditions. Correlation of binding isotherms with circular dichroism data indicates that the cooperative binding of sanguinarine under the Z-form and the H(L)-form conditions is associated with a sequential conversion of the polymer from a left-handed to a bound right-handed conformation. Determination of bound alkaloid concentration by spectroscopic titration technique and the measurement of circular dichroic spectra have enabled us to calculate the number of base pairs of Z-form and H(L)-form that adopt a right-handed conformation for each bound alkaloid. Analysis reveals that 2-3 base pairs (bp) of Z-form of poly(dG-dC).poly(dG-dC) and poly(dG-me5dC).poly(dG-me5dC) switch to the right-handed form for each bound sanguinarine, while approximately same number of base pairs switch to the bound right-handed form in complexes with H(L)-form of these polynucleotides. Comparative binding analysis shows that ethidium also converts approximately 2 bp of Z-form or H(L)-form to bound right-handed form under same experimental conditions. Since sanguinarine binds preferentially to alternating GC sequences, which are capable of undergoing the B to Z or B to H(L) transition, these effects may be an important part in understanding its extensive biological activities.

Hydration of the phosphate group in double-helical DNA

Water distributions around phosphate groups in 59 B-, A-, and Z-DNA crystal structures were analyzed. It is shown that the waters are concentrated in six hydration sites per phosphate and that the positions and occupancies of these sites are dependent on the conformation and type of nucleotide. The patterns of hydration that are characteristic of the backbone of the three DNA helical types can be attributed in part to the interactions of these hydration sites.

Sensitivity of NMR internucleotide distances to B-DNA conformation: underlying mechanics

Nuclear magnetic resonance (NMR) spectroscopy, combining correlated spectroscopy (COSY) coupling constant measurements with nuclear Overhauser effect spectroscopy (NOESY) interatomic distances, should make it possible to determine an averaged solution structure for DNA oligomers. However, even if such data could be obtained with high accuracy, it is not clear which structural parameters of DNA would be determined. Here, the relationships between measurable internucleotide distances and helical parameters are systematically studied through molecular modelling. Investigations are carried out using four representative sequences, (ACGT)n, (TCGA)n, (AGCT)n and (TGCA)n, composed of repeated tetranucleotides belonging to oligomers previously studied by NMR. Correlations between interatomic distances become evident and strong connections between distances and inter-base helical parameters are observed. Results imply that twist, roll, shift and slide values can be accurately determined from NMR data. Sequence independent mechanical coupling which link backbone and sugar conformations to helical twist are also described.

A DNA decamer with a sticky end: the crystal structure of d-CGACGATCGT

The crystal structure of d-CGACGATCGT has been determined to a resolution of 2.6 A. The molecule was synthesized by standard phosphoramidite procedures, and purified by anion-exchange HPLC. Crystals are monolclinic, space group P2(1), with unit cell dimensions, a = 26.45 A, b = 34.66 A, c = 32.17 A, beta = 113.45 degrees and Z = 4, containing a B-DNA double helix in each crystallographic asymmetric unit. The structure was solved using molecular replacement, aided by an isomorphous derivative, in which a bromine atom was attached to the 5 position of cytosine 8. Problems of fit between the search model and the structure ultimately obtained necessitated the use of Patterson correlation procedures between the determination of the orientation and the translation of the molecule. In all, 69 solvent molecules have been identified, and the structure has been refined to an R-factor of 0.214, using the 1421 reflections with F > 2sigma(F), collected at -120 degrees C. The sequence produces a molecule containing eight Watson-Crick base-pairs and a two-nucleotide 5'-sticky end at each end of the duplex. The sticky ends cohere with one another, so the molecules form continuous 10-fold double helices throughout the crystal, with each strand being interrupted by inherent staggered nicks. The relative angular relationships between helices in the structure differ from each other; most of the arrangements differ from Holliday junctions, whose rotational orientations are phased by a crossover and which are modeled to contain double helices that are exactly parallel or antiparallel. However, one helical juxtaposition in this crystal is similar to the alignment of double helices in parallel Holliday junctions. A survey of DNA decamers that also form infinite helices in crystals reveals relationships that approximate both parallel and antiparallel Holliday junction alignments.