M.HhaI binds tightly to substrates containing mismatches at the target base

Examining DNA breathing with pyDNA-EPBD

MOTIVATION

The two strands of the DNA double helix locally and spontaneously separate and recombine in living cells due to the inherent thermal DNA motion. This dynamics results in transient openings in the double helix and is referred to as "DNA breathing" or "DNA bubbles." The propensity to form local transient openings is important in a wide range of biological processes, such as transcription, replication, and transcription factors binding. However, the modeling and computer simulation of these phenomena, have remained a challenge due to the complex interplay of numerous factors, such as, temperature, salt content, DNA sequence, hydrogen bonding, base stacking, and others.

RESULTS

We present pyDNA-EPBD, a parallel software implementation of the Extended Peyrard-Bishop-Dauxois (EPBD) nonlinear DNA model that allows us to describe some features of DNA dynamics in detail. The pyDNA-EPBD generates genomic scale profiles of average base-pair openings, base flipping probability, DNA bubble probability, and calculations of the characteristically dynamic length indicating the number of base pairs statistically significantly affected by a single point mutation using the Markov Chain Monte Carlo algorithm.

AVAILABILITY AND IMPLEMENTATION

pyDNA-EPBD is supported across most operating systems and is freely available at https://github.com/lanl/pyDNA_EPBD. Extensive documentation can be found at https://lanl.github.io/pyDNA_EPBD/.

Backbone Conformational Equilibrium in Mismatched DNA Correlates with Enzyme Activity

T:G mismatches in mammals arise primarily from the deamination of methylated CpG sites or the incorporation of improper nucleotides. The process by which repair enzymes such as thymine DNA glycosylase (TDG) identify a canonical DNA base in the incorrect pairing context remains a mystery. However, the abundant contacts of the repair enzymes with the DNA backbone suggest a role for protein-phosphate interaction in the recognition and repair processes, where conformational properties may facilitate the proper interactions. We have previously used 31P NMR to investigate the energetics of DNA backbone BI-BII interconversion and the effect of a mismatch or lesion compared to canonical DNA and found stepwise differences in ΔG of 1-2 kcal/mol greater than equivalent steps in unmodified DNA. We have currently compared our results to substrate dependence for TDG, MBD4, M. HhaI, and CEBPβ, testing for correlations to sequence and base-pair dependence. We found strong correlations of our DNA phosphate backbone equilibrium (Keq) to different enzyme kinetics or binding parameters of these varied enzymes, suggesting that the backbone equilibrium may play an important role in mismatch recognition and/or conformational rearrangement and energetics during nucleotide flipping or other aspects of enzyme interrogation of the DNA substrate.

Examining DNA Breathing with pyDNA-EPBD

Motivation

The two strands of the DNA double helix locally and spontaneously separate and recombine in living cells due to the inherent thermal DNA motion.This dynamics results in transient openings in the double helix and is referred to as "DNA breathing" or "DNA bubbles." The propensity to form local transient openings is important in a wide range of biological processes, such as transcription, replication, and transcription factors binding. However, the modeling and computer simulation of these phenomena, have remained a challenge due to the complex interplay of numerous factors, such as, temperature, salt content, DNA sequence, hydrogen bonding, base stacking, and others.

Results

We present pyDNA-EPBD, a parallel software implementation of the Extended Peyrard-Bishop- Dauxois (EPBD) nonlinear DNA model that allows us to describe some features of DNA dynamics in detail. The pyDNA-EPBD generates genomic scale profiles of average base-pair openings, base flipping probability, DNA bubble probability, and calculations of the characteristically dynamic length indicating the number of base pairs statistically significantly affected by a single point mutation using the Markov Chain Monte Carlo (MCMC) algorithm.

Base Dynamics in the HhaI Protein Binding Site

Protein-DNA interactions play an important role in numerous biological functions within the living cell. In many of these interactions, the DNA helix is significantly distorted upon protein-DNA complex formation. The HhaI restriction-modification system is one such system, where the methylation target is flipped out of the helix when bound to the methyltransferase. However, the base flipping mechanism is not well understood. The dynamics of the binding site of the HhaI methyltransferase and endonuclease (underlined) within the DNA oligomer [d(G1A2T3A4G5C6G7C8T9A10T11C12)]2 are studied using deuterium solid-state NMR (SSNMR). SSNMR spectra obtained from DNAs deuterated on the base of nucleotides within and flanking the [5'-GCGC-3']2 sequence indicate that all of these positions are structurally flexible. Previously, conformational flexibility within the phosphodiester backbone and furanose ring within the target sequence has been observed and hypothesized to play a role in the distortion mechanism. However, whether that distortion was occurring through an active or passive mechanism remained unclear. These NMR data demonstrate that although the [5'-GCGC-3']2 sequence is dynamic, the target cytosine is not passively flipping out of the double-helix on the millisecond-picosecond time scale. Additionally, although previous studies have shown that both the furanose ring and phosphodiester backbone experience a change in dynamics upon methylation, which may play a role in recognition and cleavage by the endonuclease, our observations here indicate that methylation has no effect on the dynamics of the base itself.

Energy Landscapes for Base-Flipping in a Model DNA Duplex

We explore the process of base-flipping for four central bases, adenine, guanine, cytosine, and thymine, in a deoxyribonucleic acid (DNA) duplex using the energy landscape perspective. NMR imino-proton exchange and fluorescence correlation spectroscopy studies have been used in previous experiments to obtain lifetimes for bases in paired and extrahelical states. However, the difference of almost 4 orders of magnitude in the base-flipping rates obtained by the two methods implies that they are exploring different pathways and possibly different open states. Our results support the previous suggestion that minor groove opening may be favored by distortions in the DNA backbone and reveal links between sequence effects and the direction of opening, i.e., whether the base flips toward the major or the minor groove side. In particular, base flipping along the minor groove pathway was found to align toward the 5′ side of the backbone. We find that bases align toward the 3′ side of the backbone when flipping along the major groove pathway. However, in some cases for cytosine and thymine, the base flipping along the major groove pathway also aligns toward the 5′ side. The sequence effect may be caused by the polar interactions between the flipping-base and its neighboring bases on either of the strands. For guanine flipping toward the minor groove side, we find that the equilibrium constant for opening is large compared to flipping via the major groove. We find that the estimated rates of base opening, and hence the lifetimes of the closed state, obtained for thymine flipping through small and large angles along the major groove differ by 6 orders of magnitude, whereas for thymine flipping through small angles along the minor groove and large angles along the major groove, the rates differ by 3 orders of magnitude.

CRISPR-Cas molecular beacons as tool for studies of assembly of CRISPR-Cas effector complexes and their interactions with DNA

CRISPR-Cas systems protect prokaryotic cells from invading phages and plasmids by recognizing and cleaving foreign nucleic acid sequences specified by CRISPR RNA spacer sequences. Several CRISPR-Cas systems have been widely used as tool for genetic engineering. In DNA-targeting CRISPR-Cas nucleoprotein effector complexes, the CRISPR RNA forms a hybrid with the complementary strand of foreign DNA, displacing the noncomplementary strand to form an R-loop. The DNA interrogation and R-loop formation involve several distinct steps the molecular details of which are not fully understood. This chapter describes a recently developed fluorometric Cas beacon assay that may be used for measuring of specific affinity of various CRISPR-Cas complexes for unlabeled target DNA and model DNA probes. The Cas beacon approach also can provide a sensitive method for monitoring the kinetics of assembly of CRISPR-Cas complexes.

Evaluating the role of coherent delocalized phonon-like modes in DNA cyclization

The innate flexibility of a DNA sequence is quantified by the Jacobson-Stockmayer's J-factor, which measures the propensity for DNA loop formation. Recent studies of ultra-short DNA sequences revealed a discrepancy of up to six orders of magnitude between experimentally measured and theoretically predicted J-factors. These large differences suggest that, in addition to the elastic moduli of the double helix, other factors contribute to loop formation. Here, we develop a new theoretical model that explores how coherent delocalized phonon-like modes in DNA provide single-stranded "flexible hinges" to assist in loop formation. We combine the Czapla-Swigon-Olson structural model of DNA with our extended Peyrard-Bishop-Dauxois model and, without changing any of the parameters of the two models, apply this new computational framework to 86 experimentally characterized DNA sequences. Our results demonstrate that the new computational framework can predict J-factors within an order of magnitude of experimental measurements for most ultra-short DNA sequences, while continuing to accurately describe the J-factors of longer sequences. Further, we demonstrate that our computational framework can be used to describe the cyclization of DNA sequences that contain a base pair mismatch. Overall, our results support the conclusion that coherent delocalized phonon-like modes play an important role in DNA cyclization.

A higher degree of expression of DNA methyl transferase 1 in cervical cancer is associated with poor survival outcome

Background

Even though novel therapies based on aberrant DNA methylation could be of particular importance for the treatment of cervical cancer (CC) because the oncoproteins E6/E7 of high-risk human papillomaviruses, the causative agents for developing CC, have the capacity to bind and upregulate DNA methyltransferases (DNMTs), to our knowledge, no previous studies have evaluated the expression of this enzyme in CC in relation to survival outcomes. The purpose of the study was to evaluate the expression of DNMT1 in CC and its association with survival outcomes.

Methods

The study population consisted of 76 women treated for primary CC and followed up by the University of Alabama at Birmingham (UAB) cancer registry. The expression of DNMT1 was examined using immunohistochemistry, and the degree of expression of DNMT1 was expressed as a percentage of cells positive for DNMT1 and its intensity. Cox proportional hazards model was used to assess the relationship between the degree of expression of DNMT1 and overall survival after adjusting for relevant covariates.

Results

The expression of DNMT1 was significantly higher in CC cells compared to that in the normal cervical epithelium. A higher percentage of cells positive for DNMT1 and a higher intensity score for DNMT1 were significantly associated with poor survival outcome (hazard ratio [HR] =4.3, P=0.03 and HR =4.9, P=0.02, respectively).

Conclusion

Our findings suggested that the degree of expression of DNMT1 could be considered as a target in the epigenetic treatment of CC. Replication of our results in other study populations with CC could create the opportunity of using DNMT inhibitors to treat CC.

The impact of 6-thioguanine incorporation into DNA on the function of DNA methyltransferase Dnmt3a

The incorporation of chemotherapeutic agent 6-thioguanine (^SG) into DNA is a prerequisite for its cytotoxic action. This modification of DNA impedes the activity of enzymes involved in DNA repair and replication. Here, using hemimethylated DNA substrates we demonstrated that DNA methylation by Dnmt3a-CD is reduced if DNA is damaged by the incorporation of ^SG into one or two CpG sites separated by nine base pairs. An increase in the number of ^SG substitutions did not enhance the effect. Dnmt3a-CD binding to either of ^SG-containing DNA substrates was not distorted. Our results suggest that ^SG incorporation into DNA may influence epigenetic regulation via DNA methylation.

Somatic mutations in cancer: Stochastic versus predictable

The origins of human cancers remain unclear except for a limited number of potent environmental mutagens, such as tobacco and UV light, and in rare cases, familial germ line mutations that affect tumor suppressor genes or oncogenes. A significant component of cancer etiology has been deemed stochastic and correlated with the number of stem cells in a tissue, the number of times the stem cells divide and a low incidence of random DNA polymerase errors that occur during each cell division. While somatic mutations occur during each round of DNA replication, mutations in cancer driver genes are not stochastic. Out of a total of 2843 codons, 1031 can be changed to stop codons by a single base substitution in the tumor suppressor APC gene, which is mutated in 76% of colorectal cancers (CRC). However, the nonsense mutations, which comprise 65% of all the APC driver mutations in CRC, are not random: 43% occur at Arg CGA codons, although they represent <3% of the codons. In TP53, CGA codons comprise <3% of the total 393 codons but they account for 72% and 39% of the mutations in CRC and ovarian cancer OVC, respectively. This mutation pattern is consistent with the kinetically slow, but not stochastic, hydrolytic deamination of 5-methylcytosine residues at specific methylated CpG sites to afford T·G mismatches that lead to C→T transitions and stop codons at CGA. Analysis of nonsense mutations in CRC, OVC and a number of other cancers indicates the need to expand the predictable risk factors for cancer to include, in addition to random polymerase errors, the methylation status of gene body CGA codons in tumor suppressor genes.

Relationships between p14ARF Gene Methylation and Clinicopathological Features of Colorectal Cancer: A Meta-Analysis

We conducted a meta-analysis to explore the relationships between p14^ARF gene methylation and clinicopathological features of colorectal cancer (CRC). Databases, including Pubmed, Embase and Cochrane Library, were searched and, finally, a total of 18 eligible researches encompassing 1988 CRC patients were selected. Combined odds ratios (ORs) with 95% confidence intervals (95% CIs) were evaluated under a fixed effects model for absence of heterogeneity. Significant associations were observed between p14^ARF gene methylation and tumor location (OR = 2.35, 95% CI: 1.55–3.55, P = 0.001), microsatellite instability (MSI) status (OR = 3.28, 95% CI: 2.12–5.07, P<0.0001). However, there were no significant associations between p14^ARF gene methylation and tumor stage, tumor differentiation. We concluded that p14^ARF gene methylation may be significantly associated with tumor location, and MSI status of CRC.

Structures of Escherichia coli DNA adenine methyltransferase (Dam) in complex with a non-GATC sequence: potential implications for methylation-independent transcriptional repression

DNA adenine methyltransferase (Dam) is widespread and conserved among the γ-proteobacteria. Methylation of the Ade in GATC sequences regulates diverse bacterial cell functions, including gene expression, mismatch repair and chromosome replication. Dam also controls virulence in many pathogenic Gram-negative bacteria. An unexplained and perplexing observation about Escherichia coli Dam (EcoDam) is that there is no obvious relationship between the genes that are transcriptionally responsive to Dam and the promoter-proximal presence of GATC sequences. Here, we demonstrate that EcoDam interacts with a 5-base pair non-cognate sequence distinct from GATC. The crystal structure of a non-cognate complex allowed us to identify a DNA binding element, GTYTA/TARAC (where Y = C/T and R = A/G). This element immediately flanks GATC sites in some Dam-regulated promoters, including the Pap operon which specifies pyelonephritis-associated pili. In addition, Dam interacts with near-cognate GATC sequences (i.e. 3/4-site ATC and GAT). Taken together, these results imply that Dam, in addition to being responsible for GATC methylation, could also function as a methylation-independent transcriptional repressor.

Folate deficiency exacerbates apoptosis by inducing hypomethylation and resultant overexpression of DR4 together with altering DNMTs in Alzheimer's disease

This study was to evaluate patterns of gene expression and promoter methylation of DR4 from peripheral circulating blood lymphocytes of AD patients and folate-deficiency medium cultured neuroblast cells, and also expression levels of DNMT1, DNMT3a, and MECP2. Blood samples of 25 pairs of AD patients and age- and sex-matched controls were collected. SH-SY5Y cells were cultured with folate-deficiency medium. Bisulfite cloning coupled with sequencing was employed to analyze methylation levels of DR4 gene promoters, and quantitative real time PCR (qRT-PCR) was used to detect gene expression levels of DR4, and also DNA methyltransferase 1 (DNMT1), DNA methyltransferase 3a (DNMT3a) and methyl CpG binding protein 2 (MECP2). Folate concentration was calculated in serum of blood samples. 3-(4,5-di-methylthiazol-2-yl)-2,5-diphenyl tetrazolium (MTT) assay was used to analyze cell viability. The results showed that, the promoter of DR4 was hypomethylated in AD patients and cells cultured in folate-deficiency medium and had site-specific changes (P < 0.05), and these sites were mostly at or nearby some key transcription factor binding sites. Accordance with the hypomethylation, increased expression level of DR4 was observed (P < 0.05). DNMT1 and DNMT3a mRNA level were elevated (P < 0.05) in AD patients and folate deficient medium cultured cells compared with controls (P < 0.05), together with lower folate concentration in AD. MTT assay showed that folate deficiency inhibited cell growth. In summary, folate deficiency can induce apoptosis by increasing DR4 expression with DNA promoter hypomethylation in AD, together with upregulating DNMTs expression, which may be associated with folate deficiency-induced DNA damage.

Binding of nucleoid-associated protein fis to DNA is regulated by DNA breathing dynamics

Physicochemical properties of DNA, such as shape, affect protein-DNA recognition. However, the properties of DNA that are most relevant for predicting the binding sites of particular transcription factors (TFs) or classes of TFs have yet to be fully understood. Here, using a model that accurately captures the melting behavior and breathing dynamics (spontaneous local openings of the double helix) of double-stranded DNA, we simulated the dynamics of known binding sites of the TF and nucleoid-associated protein Fis in Escherichia coli. Our study involves simulations of breathing dynamics, analysis of large published in vitro and genomic datasets, and targeted experimental tests of our predictions. Our simulation results and available in vitro binding data indicate a strong correlation between DNA breathing dynamics and Fis binding. Indeed, we can define an average DNA breathing profile that is characteristic of Fis binding sites. This profile is significantly enriched among the identified in vivo E. coli Fis binding sites. To test our understanding of how Fis binding is influenced by DNA breathing dynamics, we designed base-pair substitutions, mismatch, and methylation modifications of DNA regions that are known to interact (or not interact) with Fis. The goal in each case was to make the local DNA breathing dynamics either closer to or farther from the breathing profile characteristic of a strong Fis binding site. For the modified DNA segments, we found that Fis-DNA binding, as assessed by gel-shift assay, changed in accordance with our expectations. We conclude that Fis binding is associated with DNA breathing dynamics, which in turn may be regulated by various nucleotide modifications.

Cellular transcription factors (TFs) are proteins that regulate gene expression, and thereby cellular activity and fate, by binding to specific DNA segments. The physicochemical determinants of protein-DNA binding specificity are not completely understood. Here, we report that the propensity of transient opening and re-closing of the double helix, resulting from thermal fluctuations, aka “DNA breathing” or “DNA bubbles,” can be associated with binding affinity in the case of Fis, a well-studied nucleoid-associated protein in Escherichia coli. We found that a particular breathing profile is characteristic of high-affinity Fis binding sites and that DNA fragments known to bind Fis in vivo are statistically enriched for this profile. Furthermore, we used simulations of DNA breathing dynamics to guide design of gel-shift experiments aimed at testing the idea that local breathing influences Fis binding. As a result, we show that via nucleotide modifications but without modifying nucleotides that directly contact Fis, we were able to transform a low-affinity Fis binding site into a high-affinity site and vice versa. The nucleotide modifications were designed only based on DNA breathing simulations. Our study suggests that strong Fis-DNA binding depends on DNA breathing - a novel physicochemical characteristic that could be used for prediction and rational design of TF binding sites.

Probing murine methyltransfease Dnmt3a interactions with benzo[a]pyrene-modified DNA by fluorescence methods

The impact of bulky carcinogen-DNA adducts positioned at or near recognition sites (CpG) of eukaryotic DNA methyltransferases on their catalytic activities is poorly understood. In the present study, we employed site-specifically modified 30-mer oligodeoxyribonucleotides containing stereoisomeric benzo[a]pyrene diol epoxide (B[a]PDE)-derived guanine (B[a]PDE-N(2)-dG) or adenine (B[a]PDE-N(6)-dA) adducts of different conformations as substrates of the catalytic domain of murine Dnmt3a (Dnmt3a-CD). The fluorescence of these lesions was used to examine interactions between Dnmt3a-CD and DNA. In B[a]PDE-DNA•Dnmt3a-CD complexes, the intensity of fluorescence of the covalently bound B[a]PDE residues is enhanced relative to the protein-free value when the B[a]PDE is positioned in the minor groove [(+)- and (-)-trans-B[a]PDE-N(2)-dG adducts in the CpG site] and when it is intercalated on the 5'-side of the CpG site [(+)-trans-B[a]PDE-N(6)-dA adduct]. The fluorescence of B[a]PDE-modified DNA•Dnmt3a-CD complexes exhibits only small changes when the B[a]PDE is intercalated with base displacement in (+)- and (-)-cis-B[a]PDE-N(2)-dG adducts and without base displacement in the (-)-trans-B[a]PDE-N(6)-dA adduct. The initial rates of methylation were significantly reduced by the minor groove trans-B[a]PDE-N(2)-dG adducts, regardless of their position in the substrate and by the intercalated cis-B[a]PDE-N(2)-dG adducts within the CpG site. The observed changes in fluorescence and methylation rates are consistent with the flipping of the target cytosine and a catalytic loop motion within the DNA•Dnmt3a-CD complexes. In the presence of the regulatory factor Dnmt3L, an enhancement of both methylation rates and fluorescence was observed, which is consistent with a Dnmt3L-mediated displacement of the catalytic loop towards the CpG site.

DNA methyltransferases: mechanistic models derived from kinetic analysis

The sequence-specific transfer of methyl groups from donor S-adenosyl-L-methionine (AdoMet) to certain positions of DNA-adenine or -cytosine residues by DNA methyltransferases (MTases) is a major form of epigenetic modification. It is virtually ubiquitous, except for some notable exceptions. Site-specific methylation can be regarded as a means to increase DNA information capacity and is involved in a large spectrum of biological processes. The importance of these functions necessitates a deeper understanding of the enzymatic mechanism(s) of DNA methylation. DNA MTases fall into one of two general classes; viz. amino-MTases and [C5-cytosine]-MTases. Amino-MTases, common in prokaryotes and lower eukaryotes, catalyze methylation of the exocyclic amino group of adenine ([N6-adenine]-MTase) or cytosine ([N4-cytosine]-MTase). In contrast, [C5-cytosine]-MTases methylate the cyclic carbon-5 atom of cytosine. Characteristics of DNA MTases are highly variable, differing in their affinity to their substrates or reaction products, their kinetic parameters, or other characteristics (order of substrate binding, rate limiting step in the overall reaction). It is not possible to present a unifying account of the published kinetic analyses of DNA methylation because different authors have used different substrate DNAs and/or reaction conditions. Nevertheless, it would be useful to describe those kinetic data and the mechanistic models that have been derived from them. Thus, this review considers in turn studies carried out with the most consistently and extensively investigated [N6-adenine]-, [N4-cytosine]- and [C5-cytosine]-DNA MTases.

Transducing methyltransferase activity into electrical signals in a carbon nanotube-DNA device()

This study creates a device where the DNA is electronically integrated to serve as both the biological target and electrical transducer in a CNT-DNA-CNT device. We detect DNA binding and methylation by the methyltransferase M.SssI at the single molecule level. We demonstrate sequence-specific, reversible binding of M.SssI and protein-catalyzed methylation that alters the protein-binding affinity of the device. This device, which relies on the exquisite electrical sensitivity of DNA, represents a unique route for the specific, single molecule detection of enzymatic activity.

A dual flip-out mechanism for 5mC recognition by the Arabidopsis SUVH5 SRA domain and its impact on DNA methylation and H3K9 dimethylation in vivo

Cytosine DNA methylation is evolutionarily ancient, and in eukaryotes this epigenetic modification is associated with gene silencing. Proteins with SRA (SET- or RING-associated) methyl-binding domains are required for the establishment and/or maintenance of DNA methylation in both plants and mammals. The 5-methyl-cytosine (5mC)-binding specificity of several SRA domains have been characterized, and each one has a preference for DNA methylation in different sequence contexts. Here we demonstrate through mobility shift assays and calorimetric measurements that the SU(VAR)3-9 HOMOLOG 5 (SUVH5) SRA domain differs from other SRA domains in that it can bind methylated DNA in all contexts to similar extents. Crystal structures of the SUVH5 SRA domain bound to 5mC-containing DNA in either the fully or hemimethylated CG context or the methylated CHH context revealed a dual flip-out mechanism where both the 5mC and a base (5mC, C, or G, respectively) from the partner strand are simultaneously extruded from the DNA duplex and positioned within binding pockets of individual SRA domains. Our structure-based in vivo studies suggest that a functional SUVH5 SRA domain is required for both DNA methylation and accumulation of the H3K9 dimethyl modification in vivo, suggesting a role for the SRA domain in recruitment of SUVH5 to genomic loci.

Determinants of precatalytic conformational transitions in the DNA cytosine methyltransferase M.HhaI

The DNA methyltransferase M.HhaI is an excellent model for understanding how recognition of a nucleic acid substrate is translated into site-specific modification. In this study, we utilize direct, real-time monitoring of the catalytic loop position via engineered tryptophan fluorescence reporters to dissect the conformational transitions that occur in both enzyme and DNA substrate prior to methylation of the target cytosine. Using nucleobase analogues in place of the target and orphan bases, the kinetics of the base flipping and catalytic loop closure rates were determined, revealing that base flipping precedes loop closure as the rate-determining step prior to methyl transfer. To determine the mechanism by which individual specific hydrogen bond contacts at the enzyme-DNA interface mediate these conformational transitions, nucleobase analogues lacking hydrogen bonding groups were incorporated into the recognition sequence to disrupt the major groove recognition elements. The consequences of binding, loop closure, and catalysis were determined for four contacts, revealing large differences in the contribution of individual hydrogen bonds to DNA recognition and conformational transitions on the path to catalysis. Our results describe how M.HhaI utilizes direct readout contacts to accelerate extrication of the target base that offer new insights into the evolutionary history of this important class of enzymes.

Direct observation of cytosine flipping and covalent catalysis in a DNA methyltransferase

Methylation of the five position of cytosine in DNA plays important roles in epigenetic regulation in diverse organisms including humans. The transfer of methyl groups from the cofactor S-adenosyl-l-methionine is carried out by methyltransferase enzymes. Using the paradigm bacterial methyltransferase M.HhaI we demonstrate, in a chemically unperturbed system, the first direct real-time analysis of the key mechanistic events—the flipping of the target cytosine base and its covalent activation; these changes were followed by monitoring the hyperchromicity in the DNA and the loss of the cytosine chromophore in the target nucleotide, respectively. Combined with studies of M.HhaI variants containing redesigned tryptophan fluorophores, we find that the target base flipping and the closure of the mobile catalytic loop occur simultaneously, and the rate of this concerted motion inversely correlates with the stability of the target base pair. Subsequently, the covalent activation of the target cytosine is closely followed by but is not coincident with the methyl group transfer from the bound cofactor. These findings provide new insights into the temporal mechanism of this physiologically important reaction and pave the way to in-depth studies of other base-flipping systems.

Theoretical study of the catalytic mechanism of DNA-(N4-cytosine)-methyltransferase from the bacterium Proteus vulgaris

In this paper the reaction mechanism for methylation of cytosine at the exocyclic N4 position catalyzed by M.PvuII has been explored by means of hybrid quantum mechanics/molecular mechanics (QM/MM) methods. A reaction model was prepared by placing a single cytosine base in the active site of the enzyme. In this model the exocyclic amino group of the base establishes hydrogen bond interactions with the hydroxyl oxygen atom of Ser53 and the carbonyl oxygen atom of Pro54. The reaction mechanism involves a direct methyl transfer from AdoMet to the N4 atom and a proton transfer from this atom to Ser53, which in turn transfers a proton to Asp96. Different timings for the proton transfers and methylation steps have been explored at the AM1/MM and B3LYP/MM levels including localization and characterization of stationary structures. At our best estimate the reaction proceeds by means of a simultaneous but asynchronous proton transfer from Ser53 to Asp96 and from N4 of cytosine to Ser53 followed by a direct methyl transfer from AdoMet to the exocyclic N4 of cytosine.

Molecular and enzymatic profiles of mammalian DNA methyltransferases: structures and targets for drugs

DNA methylation is an epigenetic event involved in a variety array of processes that may be the foundation of genetic phenomena and diseases. DNA methyltransferase is a key enzyme for cytosine methylation in DNA, and can be divided into two functional families (Dnmt1 and Dnmt3) in mammals. All mammalian DNA methyltransferases are encoded by their own single gene, and consisted of catalytic and regulatory regions (except Dnmt2). Via interactions between functional domains in the regulatory or catalytic regions and other adaptors or cofactors, DNA methyltransferases can be localized at selective areas (specific DNA/nucleotide sequence) and linked to specific chromosome status (euchromatin/heterochromatin, various histone modification status). With assistance from UHRF1 and Dnmt3L or other factors in Dnmt1 and Dnmt3a/Dnmt3b, mammalian DNA methyltransferases can be recruited, and then specifically bind to hemimethylated and unmethylated double-stranded DNA sequence to maintain and de novo setup patterns for DNA methylation. Complicated enzymatic steps catalyzed by DNA methyltransferases include methyl group transferred from cofactor Ado-Met to C5 position of the flipped-out cytosine in targeted DNA duplex. In the light of the fact that different DNA methyltransferases are divergent in both structures and functions, and use unique reprogrammed or distorted routines in development of diseases, design of new drugs targeting specific mammalian DNA methyltransferases or their adaptors in the control of key steps in either maintenance or de novo DNA methylation processes will contribute to individually treating diseases related to DNA methyltransferases.

Suicidal function of DNA methylation in age-related genome disintegration

This article is dedicated to the 60th anniversary of 5-methylcytosine discovery in DNA. Cytosine methylation can affect genetic and epigenetic processes, works as a part of the genome-defense system and has mutagenic activity; however, the biological functions of this enzymatic modification are not well understood. This review will put forward the hypothesis that the host-defense role of DNA methylation in silencing and mutational destroying of retroviruses and other intragenomic parasites was extended during evolution to most host genes that have to be inactivated in differentiated somatic cells, where it acquired a new function in age-related self-destruction of the genome. The proposed model considers DNA methylation as the generator of 5mC>T transitions that induce 40-70% of all spontaneous somatic mutations of the multiple classes at CpG and CpNpG sites and flanking nucleotides in the p53, FIX, hprt, gpt human genes and some transgenes. The accumulation of 5mC-dependent mutations explains: global changes in the structure of the vertebrate genome throughout evolution; the loss of most 5mC from the DNA of various species over their lifespan and the Hayflick limit of normal cells; the polymorphism of methylation sites, including asymmetric mCpNpN sites; cyclical changes of methylation and demethylation in genes. The suicidal function of methylation may be a special genetic mechanism for increasing DNA damage and the programmed genome disintegration responsible for cell apoptosis and organism aging and death.

Impact of 7,8-dihydro-8-oxoguanine on methylation of the CpG site by Dnmt3a

7,8-Dihydro-8-oxoguanine (8-oxoG) is a ubiquitous oxidative DNA lesion resulting from injury to DNA via reactive oxygen species. 8-oxoG lesions may play a role in the formation of aberrant DNA methylation patterns during carcinogenesis. In this study, we assessed the effects of 8-oxoG on methylation and complex formation of nine 30-mer oligodeoxynucleotide duplexes by the catalytic domain of murine Dnmt3a DNA methyltransferase (Dnmt3a-CD). The effects of 8-oxoG on the methylation rate of hemimethylated duplexes varied from a 25-fold decrease to a 1.8-fold increase, depending on the position of the lesion relative to the Dnmt3a-CD recognition site (CpG) and target cytosine (C). The most significant effect was observed when 8-oxoG replaced guanine within the recognition site immediately downstream of the target cytosine. Fluorescence polarization experiments with fluorescein-labeled duplexes revealed that two molecules of Dnmt3a-CD bind per duplex, generating sigmoid binding curves. Duplexes exhibiting the highest apparent binding cooperativity formed the least stable 1:2 complexes with Dnmt3a-CD and were methylated at the lowest rate. Kinetic analyses disclosed the formation of very stable nonproductive enzyme-substrate complexes with hemimethylated duplexes that act as suicide substrates of Dnmt3a-CD. The presence of 8-oxoG within the CpG site downstream of the target cytosine markedly diminished productive versus nonproductive binding. We propose that 8-oxoG located adjacent to the target cytosine interferes with methylation by weakening the affinity of DNA for Dnmt3a-CD, thereby favoring a nonproductive binding mode.

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.

How PspGI, catalytic domain of EcoRII and Ecl18kI acquire specificities for different DNA targets

Restriction endonucleases Ecl18kI and PspGI/catalytic domain of EcoRII recognize CCNGG and CCWGG sequences (W stands for A or T), respectively. The enzymes are structurally similar, interact identically with the palindromic CC:GG parts of their recognition sequences and flip the nucleotides at their centers. Specificity for the central nucleotides could be influenced by the strength/stability of the base pair to be disrupted and/or by direct interactions of the enzymes with the flipped bases. Here, we address the importance of these contributions. We demonstrate that wt Ecl18kI cleaves oligoduplexes containing canonical, mismatched and abasic sites in the central position of its target sequence CCNGG with equal efficiencies. In contrast, substitutions in the binding pocket for the extrahelical base alter the Ecl18kI preference for the target site: the W61Y mutant prefers only certain mismatched substrates, and the W61A variant cuts exclusively at abasic sites, suggesting that pocket interactions play a major role in base discrimination. PspGI and catalytic domain of EcoRII probe the stability of the central base pair and the identity of the flipped bases in the pockets. This ‘double check’ mechanism explains their extraordinary specificity for an A/T pair in the flipping position.

The SRA domain of UHRF1 flips 5-methylcytosine out of the DNA helix

Maintenance methylation of hemimethylated CpG dinucleotides at DNA replication forks is the key to faithful mitotic inheritance of genomic methylation patterns. UHRF1 (ubiquitin-like, containing PHD and RING finger domains 1) is required for maintenance methylation by interacting with DNA nucleotide methyltransferase 1 (DNMT1), the maintenance methyltransferase, and with hemimethylated CpG, the substrate for DNMT1 (refs 1 and 2). Here we present the crystal structure of the SET and RING-associated (SRA) domain of mouse UHRF1 in complex with DNA containing a hemimethylated CpG site. The DNA is contacted in both the major and minor grooves by two loops that penetrate into the middle of the DNA helix. The 5-methylcytosine has flipped completely out of the DNA helix and is positioned in a binding pocket with planar stacking contacts, Watson-Crick polar hydrogen bonds and van der Waals interactions specific for 5-methylcytosine. Hence, UHRF1 contains a previously unknown DNA-binding module and is the first example of a non-enzymatic, sequence-specific DNA-binding protein domain to use the base flipping mechanism to interact with DNA.

Selective recognition of pyrimidine-pyrimidine DNA mismatches by distance-constrained macrocyclic bis-intercalators

Binding of three macrocyclic bis-intercalators, derivatives of acridine and naphthalene, and two acyclic model compounds to mismatch-containing and matched duplex oligodeoxynucleotides was analyzed by thermal denaturation experiments, electrospray ionization mass spectrometry studies (ESI-MS) and fluorescent intercalator displacement (FID) titrations. The macrocyclic bis-intercalators bind to duplexes containing mismatched thymine bases with high selectivity over the fully matched ones, whereas the acyclic model compounds are much less selective and strongly bind to the matched DNA. Moreover, the results from thermal denaturation experiments are in very good agreement with the binding affinities obtained by ESI-MS and FID measurements. The FID results also demonstrate that the macrocyclic naphthalene derivative BisNP preferentially binds to pyrimidine–pyrimidine mismatches compared to all other possible base mismatches. This ligand also efficiently competes with a DNA enzyme (M.TaqI) for binding to a duplex with a TT-mismatch, as shown by competitive fluorescence titrations. Altogether, our results demonstrate that macrocyclic distance-constrained bis-intercalators are efficient and selective mismatch-binding ligands that can interfere with mismatch-binding enzymes.

13C relaxation studies of the DNA target sequence for hhai methyltransferase reveal unique motional properties

The goal of this work was to examine if sequence-dependent conformational flexibility in DNA plays a role in base extrusion, a common conformational change induced by many DNA-modifying enzymes. We studied the dynamics of the double-stranded DNA target of the HhaI methyltransferase by recording an extensive set of (13)C NMR relaxation parameters. We observe that the cytidine furanose rings experience fast (picosecond to nanosecond) motions that are not present in other nucleotides; the methylation site experiences particularly high mobility. We also observe that the bases of guanosine and cytidine residues within the HhaI recognition sequence GCGC experience motions on a much slower (1-100 micros) time scale. We compare these observations with previous solution and solid-state NMR studies of the EcoRI nuclease target sequence, and solid-state NMR studies of a similar HhaI target construct. While an increased mobility of cytidine furanose rings compared to those of other nucleotides is observed for both sequences, the slower motions are only observed in the HhaI target DNA. We propose that this inherent flexibility lowers the energetic barriers that must occur when the DNA binds to the HhaI methyltransferase and for extrusion of the cytidine prior to its methylation.

Backbone dynamics in the DNA HhaI protein binding site

The dynamics of the phosphodiester backbone in the [5'-GCGC-3'] 2 moiety of the DNA oligomer [d(G 1A 2T 3A 4 G 5 C 6 G 7 C 8T 9A 10T 11C 12)] 2 are studied using deuterium solid-state NMR (SSNMR). SSNMR spectra obtained from DNAs nonstereospecifically deuterated on the 5' methylene group of nucleotides within the [5'-GCGC-3'] 2 moiety indicated that all of these positions are structurally flexible. Previous work has shown that methylation reduces the amplitude of motion in the phosphodiester backbone and furanose ring of the same DNA, and our observations indicate that methylation perturbs backbone dynamics through not only a loss of mobility but also a change of direction of motion. These NMR data indicate that the [5'-GCGC-3'] 2 moiety is dynamic, with the largest amplitude motions occurring nearest the methylation site. The change of orientation of this moiety in DNA upon methylation may make the molecule less amenable to binding to the HhaI endonuclease.

Chemical mapping of cytosines enzymatically flipped out of the DNA helix

Haloacetaldehydes can be employed for probing unpaired DNA structures involving cytosine and adenine residues. Using an enzyme that was structurally proven to flip its target cytosine out of the DNA helix, the HhaI DNA methyltransferase (M.HhaI), we demonstrate the suitability of the chloroacetaldehyde modification for mapping extrahelical (flipped-out) cytosine bases in protein-DNA complexes. The generality of this method was verified with two other DNA cytosine-5 methyltransferases, M.AluI and M.SssI, as well as with two restriction endonucleases, R.Ecl18kI and R.PspGI, which represent a novel class of base-flipping enzymes. Our results thus offer a simple and convenient laboratory tool for detection and mapping of flipped-out cytosines in protein-DNA complexes.

Mandatory fortification with folic acid in the United States is associated with increased expression of DNA methyltransferase-1 in the cervix

OBJECTIVE

The objective of this study was to evaluate whether mandatory fortification of grain products with folic acid in the United States is associated with changes in DNA methyltransferase-1 (Dnmt-1) expression in cells involved in cervical carcinogenesis.

METHODS

Archived specimens of cervical intraepithelial neoplasia (CIN) diagnosed before (1990-1992) and after (2000-2002) mandatory folic acid fortification were used to examine the expression of Dnmt-1 in specific lesions involved in cervical carcinogenesis by immunohistochemistry. The total number of lesions examined was 101 in the prefortification period and 96 in the postfortification period. Immunohistochemical staining for Dnmt-1, its assessment, and data entry were blinded with regard to the fortification status.

RESULTS

Age- and race-adjusted mean percentage of cells positive for Dnmt-1 or the Dnmt-1 score was significantly higher in all lesion types (i.e., normal cervical epithelium, reactive cervical epithelium, metaplastic cervical epithelium, CIN, or carcinoma in situ) detected in the postfortification period compared with the prefortification period (P < 0.05, all comparisons). The degree of Dnmt-1 was significantly higher (P < 0.0001) in CIN > or = 2 lesions compared with CIN < or = 1 lesions, regardless of the fortification group.

CONCLUSION

These results suggest that mandatory fortification with folic acid in the United States seems to have resulted in a change in the degree of expression of Dnmt-1 in cells involved in cervical carcinogenesis. Because the approach we have taken to demonstrate these differences have limitations inherent to a study of this nature and this is the first study to report a folate fortification associated change in Dnmt-1, validating these results in other study populations and/or with other techniques of assessing Dnmt-1 will increase the scientific credibility of these findings.

The stereochemistry of benzo[a]pyrene-2'-deoxyguanosine adducts affects DNA methylation by SssI and HhaI DNA methyltransferases

The biologically most significant genotoxic metabolite of the environmental pollutant benzo[a]pyrene (B[a]P), (+)-7R,8S-diol 9S,10R-epoxide, reacts chemically with guanine in DNA, resulting in the predominant formation of (+)-trans-B[a]P-N(2)-dG and, to a lesser extent, (+)-cis-B[a]P-N(2)-dG adducts. Here, we compare the effects of the adduct stereochemistry and conformation on the methylation of cytosine catalyzed by two purified prokaryotic DNA methyltransferases (MTases), SssI and HhaI, with the lesions positioned within or adjacent to their CG and GCGC recognition sites, respectively. The fluorescence properties of the pyrenyl residues of the (+)-cis-B[a]P-N(2)-dG and (+)-trans-B[a]P-N(2)-dG adducts in complexes with MTases are enhanced, but to different extents, indicating that aromatic B[a]P residues are positioned in different microenvironments in the DNA-protein complexes. We have previously shown that the (+)-trans-isomeric adduct inhibits both the binding and methylating efficiencies (k(cat)) of both MTases [Subach OM, Baskunov VB, Darii MV, Maltseva DV, Alexandrov DA, Kirsanova OV, Kolbanovskiy A, Kolbanovskiy M, Johnson F, Bonala R, et al. (2006) Biochemistry45, 6142-6159]. Here we show that the stereoisomeric (+)-cis-B[a]P-N(2)-dG lesion has only a minimal effect on the binding of these MTases and on k(cat). The minor-groove (+)-trans adduct interferes with the formation of the normal DNA minor-groove contacts with the catalytic loop of the MTases. However, the intercalated base-displaced (+)-cis adduct does not interfere with the minor-groove DNA-catalytic loop contacts, allowing near-normal binding of the MTases and undiminished k(cat) values.

The role of Arg165 towards base flipping, base stabilization and catalysis in M.HhaI

Arg165 forms part of a previously identified base flipping motif in the bacterial DNA cytosine methyltransferase, M.HhaI. Replacement of Arg165 with Ala has no detectable effect on either DNA or AdoMet affinity, yet causes the base flipping and restacking transitions to be decreased approximately 16 and 190-fold respectively, thus confirming the importance of this motif. However, these kinetic changes cannot account for the mutant's observed 10(5)-fold decreased catalytic rate. The mutant enzyme/cognate DNA cocrystal structure (2.79 A resolution) shows the target cytosine to be positioned approximately 30 degrees into the major groove, which is consistent with a major groove pathway for nucleotide flipping. The pyrimidine-sugar chi angle is rotated to approximately +171 degrees, from a range of -95 degrees to -120 degrees in B DNA, and -77 degrees in the WT M.HhaI complex. Thus, Arg165 is important for maintaining the cytosine positioned for nucleophilic attack by Cys81. The cytosine sugar pucker is in the C2'-endo-C3'-exo (South conformation), in contrast to the previously reported C3'-endo (North conformation) described for the original 2.70 A resolution cocrystal structure of the WT M.HhaI/DNA complex. We determined a high resolution structure of the WT M.HhaI/DNA complex (1.96 A) to better determine the sugar pucker. This new structure is similar to the original, lower resolution WT M.HhaI complex, but shows that the sugar pucker is O4'-endo (East conformation), intermediate between the South and North conformers. In summary, Arg165 plays significant roles in base flipping, cytosine positioning, and catalysis. Furthermore, the previously proposed M.HhaI-mediated changes in sugar pucker may not be an important contributor to the base flipping mechanism. These results provide insights into the base flipping and catalytic mechanisms for bacterial and eukaryotic DNA methyltransferases.

Transition from nonspecific to specific DNA interactions along the substrate-recognition pathway of dam methyltransferase

DNA methyltransferases methylate target bases within specific nucleotide sequences. Three structures are described for bacteriophage T4 DNA-adenine methyltransferase (T4Dam) in ternary complexes with partially and fully specific DNA and a methyl-donor analog. We also report the effects of substitutions in the related Escherichia coli DNA methyltransferase (EcoDam), altering residues corresponding to those involved in specific interaction with the canonical GATC target sequence in T4Dam. We have identified two types of protein-DNA interactions: discriminatory contacts, which stabilize the transition state and accelerate methylation of the cognate site, and antidiscriminatory contacts, which do not significantly affect methylation of the cognate site but disfavor activity at noncognate sites. These structures illustrate the transition in enzyme-DNA interaction from nonspecific to specific interaction, suggesting that there is a temporal order for formation of specific contacts.

Probing a rate-limiting step by mutational perturbation of AdoMet binding in the HhaI methyltransferase

DNA methylation plays important roles via regulation of numerous cellular mechanisms in diverse organisms, including humans. The paradigm bacterial methyltransferase (MTase) HhaI (M.HhaI) catalyzes the transfer of a methyl group from the cofactor S-adenosyl-l-methionine (AdoMet) onto the target cytosine in DNA, yielding 5-methylcytosine and S-adenosyl-l-homocysteine (AdoHcy). The turnover rate (k_cat) of M.HhaI, and the other two cytosine-5 MTases examined, is limited by a step subsequent to methyl transfer; however, no such step has so far been identified. To elucidate the role of cofactor interactions during catalysis, eight mutants of Trp41, which is located in the cofactor binding pocket, were constructed and characterized. The mutants show full proficiency in DNA binding and base-flipping, and little variation is observed in the apparent methyl transfer rate k_chem as determined by rapid-quench experiments using immobilized fluorescent-labeled DNA. However, the Trp41 replacements with short side chains substantially perturb cofactor binding (100-fold higher KDAdoMet and KMAdoMet) leading to a faster turnover of the enzyme (10-fold higher k_cat). Our analysis indicates that the rate-limiting breakdown of a long-lived ternary product complex is initiated by the dissociation of AdoHcy or the opening of the catalytic loop in the enzyme.

Specificity in protein-DNA interactions: energetic recognition by the (cytosine-C5)-methyltransferase from HhaI

Sequence-specific interactions between proteins and DNA are essential for a variety of biological functions. The (cytosine-C5)-methyltransferase from HhaI (M.HhaI) specifically modifies the second base in GCGC sequences, employing a base flipping mechanism to access the target base being chemically modified. The mechanism of sequence-specific recognition of M.HhaI is not evident based on crystallographic structures, leading to the suggestion that recognition is linked to the flipping event itself, a process that may be referred to as energetic recognition. Using computational methods, it is shown that the free energy barriers to flipping are significantly higher in non-cognate versus the cognate sequence, supporting the energetic recognition mechanism. Energetic recognition is imparted by two protein "selectivity filters" that function via a "web" of protein-DNA interactions in short-lived, high energy states present along the base flipping pathway. Other sequence-specific DNA binding proteins whose function involves significant distortion of DNA's conformation may use a similar recognition mechanism.

Stopped-flow and mutational analysis of base flipping by the Escherichia coli Dam DNA-(adenine-N6)-methyltransferase

By stopped-flow kinetics using 2-aminopurine as a probe to detect base flipping, we show here that base flipping by the Escherichia coli Dam DNA-(adenine-N6)-methyltransferase (MTase) is a biphasic process: target base flipping is very fast (k(flip)>240 s(-1)), but binding of the flipped base into the active site pocket of the enzyme is slow (k=0.1-2 s(-1)). Whereas base flipping occurs in the absence of S-adenosyl-l-methionine (AdoMet), binding of the target base in the active site pocket requires AdoMet. Our data suggest that the tyrosine residue in the DPPY motif conserved in the active site of DNA-(adenine-N6)-MTases stacks to the flipped target base. Substitution of the aspartic acid residue of the DPPY motif by alanine abolished base flipping, suggesting that this residue contacts and stabilizes the flipped base. The exchange of Ser188 located in a loop next to the active center by alanine led to a seven- to eightfold reduction of k(flip), which was also reduced with substrates having altered GATC recognition sites and in the absence of AdoMet. These findings provide evidence that the enzyme actively initiates base flipping by stabilizing the transition state of the process. Reduced rates of base flipping in substrates containing the target base in a non-canonical sequence demonstrate that DNA recognition by the MTase starts before base flipping. DNA recognition, cofactor binding and base flipping are correlated and efficient base flipping takes place only if the enzyme has bound to a cognate target site and AdoMet is available.

HhaI DNA methyltransferase uses the protruding Gln237 for active flipping of its target cytosine

Access to a nucleotide by its rotation out of the DNA helix (base flipping) is used by numerous DNA modification and repair enzymes. Despite extensive studies of the paradigm HhaI methyltransferase, initial events leading to base flipping remained elusive. Here we demonstrate that the replacement of the target C:G pair with the 2-aminopurine:T pair in the DNA or shortening of the side chain of Gln237 in the protein severely perturb base flipping, but retain specific DNA binding. Kinetic analyses and molecular modeling suggest that a steric interaction between the protruding side chain of Gln237 and the target cytosine in B-DNA reduces the energy barrier for flipping by 3 kcal/mol. Subsequent stabilization of an open state by further 4 kcal/mol is achieved through specific hydrogen bonding of the side chain to the orphan guanine. Gln237 thus plays a key role in actively opening the target C:G pair by a "push-and-bind" mechanism.

Symmetry elements in DNA structure important for recognition/methylation by DNA [amino]-methyltransferases

The phage T4Dam and EcoDam DNA-[adenine-N6] methyltransferases (MTases) methylate GATC palindromic sequences, while the BamHI DNA-[cytosine-N4] MTase methylates the GGATCC palindrome (which contains GATC) at the internal cytosine residue. We compared the ability of these enzymes to interact productively with defective duplexes in which individual elements were deleted on one chain. A sharp decrease in kcat was observed for all three enzymes if a particular element of structural symmetry was disrupted. For the BamHI MTase, integrity of the ATCC was critical, while an intact GAT sequence was necessary for the activity of T4Dam, and an intact GA was necessary for EcoDam. Theoretical alignment of the region of best contacts between the protein and DNA showed that in the case of a palindromic interaction site, a zone covering the 5'-symmetric residues is located in the major groove versus a zone of contact covering the 3'-symmetric residues in the minor groove. Our data fit a simple rule of thumb that the most important contacts are aligned around the methylation target base: if the target base is in the 5' half of the palindrome, the interaction between the enzyme and the DNA occurs mainly in the major groove; if it is in the 3' half, the interaction occurs mainly in the minor groove.

Caught in the act: visualization of an intermediate in the DNA base-flipping pathway induced by HhaI methyltransferase

Rotation of a DNA or RNA nucleotide out of the double helix and into a protein pocket ('base flipping') is a mechanistic feature common to some DNA/RNA-binding proteins. Here, we report the structure of HhaI methyltransferase in complex with DNA containing a south-constrained abasic carbocyclic sugar at the target site in the presence of the methyl donor byproduct AdoHcy. Unexpectedly, the locked south pseudosugar appears to be trapped in the middle of the flipping pathway via the DNA major groove, held in place primarily through Van der Waals contacts with a set of invariant amino acids. Molecular dynamics simulations indicate that the structural stabilization observed with the south-constrained pseudosugar will not occur with a north-constrained pseudosugar, which explains its lowered binding affinity. Moreover, comparison of structural transitions of the sugar and phosphodiester backbone observed during computational studies of base flipping in the M.HhaI-DNA-AdoHcy ternary complex indicate that the south-constrained pseudosugar induces a conformation on the phosphodiester backbone that corresponds to that of a discrete intermediate of the base-flipping pathway. As previous crystal structures of M.HhaI ternary complex with DNA displayed the flipped sugar moiety in the antipodal north conformation, we suggest that conversion of the sugar pucker from south to north beyond the middle of the pathway is an essential part of the mechanism through which flipping must proceed to reach its final destination. We also discuss the possibility of the south-constrained pseudosugar mimicking a transition state in the phosphodiester and sugar moieties that occurs during DNA base flipping in the presence of M.HhaI.

Atomistic view of base flipping in DNA

Base flipping is essential for the enzyme-catalysed methylation of DNA. In our previous studies, the flipping of bases out of duplex DNA in DNA alone and when bound to the (cytosine-C5)-methyltransferase from HhaI (M.HhaI) were investigated via potential of mean force calculations. Insights into various experimental observations were obtained. In the present paper we present an overview of previous computational studies of base flipping along with new detailed structural and energetic analysis on atomic events that contribute to the free energy surfaces. The contributions from different intrinsic and environmental effects to the base-flipping process are explored, and experimental data derived from a variety of methods are reconciled. A detailed protein-facilitated base-flipping mechanism is proposed. Ground-state destabilization is achieved via disruption of the target base Watson-Crick interactions by substitution with favourable DNA-protein interactions. In addition, specific DNA-protein interactions and favourable solvation effects further promote target base flipping along the major groove through the protein matrix, and maximal interactions occur between the DNA and the protein upon reaching the fully flipped state. Other DNA binding proteins that involve distortion of DNA's conformation may use a similar mechanism to that by which M.HhaI facilitates base flipping.

DNA mismatch-specific base flipping by a bisacridine macrocycle

Most, if not all, enzymes that chemically modify nucleobases in DNA flip their target base from the inside of the double helix into an extrahelical position. This energetically unfavorable conformation is partly stabilized by specific binding of the apparent abasic site being formed. Thus, DNA base-flipping enzymes, like DNA methyltransferases and DNA glycosylases, generally bind very strongly to DNA containing abasic sites or abasic-site analogues. The macrocyclic bisacridine BisA has previously been shown to bind abasic sites. Herein we demonstrate that it is able to specifically recognize DNA base mismatches and most likely induces base flipping. Specific binding of BisA to DNA mismatches was studied by thermal denaturation experiments by using short duplex oligodeoxynucleotides containing central TT, TC, or TG mismatches or a TA match. In the presence of the macrocycle a strong increase in the melting temperature of up to 7.1 degrees C was observed for the mismatch-containing duplexes, whereas the melting temperature of the fully matched duplex was unaffected. Furthermore, BisA binding induced an enhanced reactivity of the mispaired thymine residue in the DNA toward potassium permanganate oxidation. A comparable reactivity has previously been observed for a TT target base mismatch in the presence of DNA methyltransferase M.TaqI. This similarity to a known base-flipping enzyme suggests that insertion of BisA into the DNA helix displaces the mispaired thymine residue into an extrahelical position, where it should be more prone to chemical oxidation. Thus, DNA base flipping does not appear to be limited to DNA-modifying enzymes but it is likely to also be induced by a small synthetic molecule binding to a thermodynamically weakened site in DNA.

The coupling of tight DNA binding and base flipping: identification of a conserved structural motif in base flipping enzymes

Val(121) is positioned immediately above the extrahelical cytosine in HhaI DNA C(5)-cytosine methyltransferase, and replacement with alanine dramatically interferes with base flipping and catalysis. DNA binding and k(cat) are decreased 10(5)-fold for the Val(121) --> Ala mutant that has a normal circular dichroism spectrum and AdoMet affinity. The magnitude of this loss of function is comparable with removal of the essential catalytic Cys(81). Surprisingly, DNA binding is completely recovered (increase of 10(5)-fold) with a DNA substrate lacking the target cytosine base (abasic). Thus, interfering with the base flipping transition results in a dramatic loss of binding energy. Our data support an induced fit mechanism in which tight DNA binding is coupled to both base flipping and protein loop rearrangement. The importance of the proximal protein segment (His(127)-Thr(132)) in maintaining this critical interaction between Val(121) and the flipped cytosine was probed with single site alanine substitutions. None of these mutants are significantly altered in secondary structure, AdoMet or DNA affinity, k(methylation), k(inactivation), or k(cat). Although Val(121) plays a critical role in both extrahelical base stabilization and catalysis, its position and mobility are not influenced by individual residues in the adjacent peptide region. Structural comparisons with other DNA methyltransferases and DNA repair enzymes that stabilize extrahelical nucleotides reveal a motif that includes a positively charged or polar side chain and a hydrophobic residue positioned adjacent to the target DNA base and either the 5'- or 3'-phosphate.

Mechanisms of DNA damage, DNA hypomethylation, and tumor progression in the folate/methyl-deficient rat model of hepatocarcinogenesis

Using the folate/methyl-deficient rat model of hepatocarcinogenesis, we obtained evidence that may provide new insights into a major unresolved paradox in DNA methylation and cancer research: the mechanistic basis for genome-wide hypomethylation despite an increase in DNA methyltransferase activity and gene-specific regional hypermethylation. Previous studies revealed that the methyltransferase binds with higher affinity to DNA strand breaks, gaps, abasic sites, and uracil than it does to its cognate hemimethylated CpG sites, consistent with its ancestral function as a DNA repair enzyme. These same DNA lesions are an early occurrence in models of folate and methyl deficiency and are often present in human preneoplastic cells. We hypothesized that the high-affinity binding of the maintenance DNA methyltransferase to unrepaired lesions in DNA could sequester available enzyme away from the replication fork and promote passive replication-dependent demethylation. In support of this possibility, we found that lesion-containing DNA is less efficiently methylated than lesion-free DNA from folate/methyl-deficient rats and that an increase in DNA strand breaks precedes DNA hypomethylation. Despite an adaptive increase in DNA methyltransferase activity, hemimethylated DNA from folate/methyl-deficient rats is progressively replaced by double-stranded unmethylated DNA that is resistant to remethylation with dietary methyl repletion. In promoter regions, the inappropriate binding of the DNA methyltransferase to unrepaired lesions or mispairs may promote local histone deacetylation, methylation, and regional hypermethylation associated with tumor suppressor gene silencing. These insights in an experimental model are consistent with the possibility that DNA lesions may be a necessary prerequisite for the disruption of normal DNA methylation patterns in preneoplastic and neoplastic cells.

Methods for the design and analysis of oligodeoxynucleotide-based DNA (cytosine-5) methyltransferase inhibitors

Several second-generation inhibitors of DNA (cytosine-5) methyltransferases based on studies of modified synthetic oligodeoxynucleoides have been described. As an aid to studies of these inhibitors, we present an electronic structure-based algorithm that can be used as a method for predicting the nature of the expected inhibition by any noncytosine nucleotide target. Targeting by the major human enzyme (hDnmt1) is governed by the presence of a three-nucleotide motif. In hemimethylated DNA, this motif consists of a 5-methylcytosine targeting signal that causes the enzyme to probe the opposite strand for a normally paired guanosine or inosine residue and attempt to methylate the residue 5' to that site. As a demonstration of the method, we apply these rules to the design and characterization of a novel oligodeoxynucleotide inhibitor of hDnmt1. This inhibitor takes advantage of the three-nucleotide recognition motif characteristic of hDnmt1 and shows that the enzyme is inhibited in vitro by non-CG methylation which targets the enzyme to normally basepaired but unproductive nucleotides such as dG, dA, and dT. Kinetic analysis at constant S-adenosyl-L-methionine concentration shows that representative inhibitory oligodeoxynucleotides are best viewed as weakly productive components of systems containing two DNA substrates. This model suggests that the most effective inhibitors are those with very low apparent Vmax and very low Km values. Oligodeoxynucleotides containing mispaired and unproductive targets such as dG, dA, dT, and dU are also inhibitory as secondary substrates for the human enzyme. Biologically, fail-safe mechanisms identified by the ab initio approach appear to be active in preventing potentially mutagenic deamination of dihydrocytosine and enzymatic methylation of dU.

Theoretical ab initio study of the effects of methylation on structure and stability of G:C Watson-Crick base pair

Methylation of DNA occurs most readily at N(3), N(7), and O(6) of purine bases and N(3) and O(2) of pyrimidines. Methylated bases are continuously formed through endogenous and exogenous mechanisms. The results of a theoretical ab initio study on the methylation of G:C base pair components are reported. The geometries of the local minima were optimized without symmetry restrictions by the gradient procedure at DFT level of theory and were verified by energy second derivative calculations. The standard 6-31G(d) basis set was used. The single-point calculations have been performed at the MP2/6-31G(d,p), MP2/6-31++G(d,p), and MP2/6-311++G(2d,2p) levels of theory. The geometrical parameters, relative stability and counterpoise corrected interaction energies are reported. Also, using a variation-perturbation energy decomposition scheme we have found the vital contributions to the total interaction energy.

Protein-facilitated base flipping in DNA by cytosine-5-methyltransferase

DNA methylation, various DNA repair mechanisms, and possibly early events in the opening of DNA as required for transcription and replication are initiated by flipping of a DNA base out of the DNA double helix. The energetics and structural mechanism of base flipping in the presence of the DNA-processing enzyme, cytosine 5-methyltransferase from HhaI (M.HhaI), were obtained through molecular dynamics based upon free-energy calculations. Free-energy profiles for base flipping show that, when in the closed conformation, M.HhaI lowers the free-energy barrier to flipping by 17 kcalmol and stabilizes the fully flipped state. Flipping is shown to occur via the major groove of the DNA. Structural analysis indicates that flipping is facilitated by destabilization of the DNA double-helical structure and substitution of DNA base-pairing and base-stacking interactions with DNA-protein interactions. The fully flipped state is stabilized by DNA-protein interactions that are enhanced upon binding of coenzyme. This study represents an atomic detail description of the mechanism by which a protein facilitates specific structural distortion in DNA.