Determinants of sequence-specific DNA methylation: target recognition and catalysis are coupled in M.HhaI

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.

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.

Unveiling the Methyl Transfer Mechanisms in the Epigenetic Machinery DNMT3A-3L: A Comprehensive Study Integrating Assembly Dynamics with Catalytic Reactions

In epigenetic mechanisms, DNA methyltransferase 3 alpha (DNMT3A) acts as an initiator for DNA methylation and prevents the downstream genes from expressing. Perturbations of DNMT3A functions may cause uncontrolled gene expression, resulting in pathogenic consequences such as cancers. It is, therefore, vitally important to understand the catalytic process of DNMT3A in its biological macromolecule assembly, viz., heterotetramer: (DNMT3A-3 L)^dimer. In this study, we utilized molecular dynamics (MD) simulations, Markov State Models (MSM), and quantum mechanics/molecular mechanics simulations (QM/MM) to investigate the de novo methyl transfer process. We identified the dynamics of the key residues relevant to the insertion of the target cytosine (dC) into the catalytic domain of DNMT3A, and the detailed potential energy surface of the seven-step reaction referring to methyl transfer. Our calculated potential energy barrier (22.51 kcal/mol) approximates the former experimental data (23.12 kcal/mol). The conformational change of the 5-methyl-cytosine (5mC) intermediate was found necessary in forming a four-water chain for the elimination step, which is unique to the other DNMTs. The biological assembly facilitates the creation of such a water chain, and the elimination occurs in an asynchronized mechanism in the two catalytic pockets. We anticipate the findings can enable a better understanding of the general mechanisms of the de novo methyl transfer for fulfilling the key enzymatic functions in epigenetics. And the unique elimination of DNMT3A might ignite novel methods for designing anti-cancer and tumor inhibitors of DNMTs.

Role of Conformational Fluctuations of Protein toward Methylation in DNA by Cytosine-5-methyltransferase

Methylation of cytosine is the common epigenetic modification in genomes ranging from bacteria to mammals, and aberrant methylation leads to human diseases including cancer. Recognition of a cognate DNA sequence by DNA methyltransferases and flipping of a target base into the enzyme active site pocket are the key steps in DNA methylation. Using molecular dynamics simulations and enhanced sampling techniques here we elucidate the role of conformational fluctuations of protein and active or passive involvement of protein elements that mediate base flipping and formation of the closed catalytic complex. The free energy profiles for the flipping of target cytosine into the enzyme active site support the major groove base eversion pathway; and the results show that the closed state of enzyme increases the free energy barrier, whereas the open state reduces it. We found that the interactions of the key loop residues of protein with cognate DNA altered the protein motions, and modulation of protein fluctuations relates to the closed catalytic complex formation. Methylation of cytosine in the active site of the closed complex destabilizes the interactions of catalytic loop residues with cognate DNA and reduces the stability of the closed state. Our study provides microscopic insights on the base flipping mechanism coupled with enzyme's loop motions and provides evidence for the role of conformational fluctuations of protein in the enzyme-catalyzed DNA processing mechanism.

The highly specific, cell cycle-regulated methyltransferase from Caulobacter crescentus relies on a novel DNA recognition mechanism

Two DNA methyltransferases, Dam and β-class cell cycle-regulated DNA methyltransferase (CcrM), are key mediators of bacterial epigenetics. CcrM from the bacterium Caulobacter crescentus (CcrM C. crescentus, methylates adenine at 5'-GANTC-3') displays 10⁵-10⁷-fold sequence discrimination against noncognate sequences. However, the underlying recognition mechanism is unclear. Here, CcrM C. crescentus activity was either improved or mildly attenuated with substrates having one to three mismatched bp within or adjacent to the recognition site, but only if the strand undergoing methylation is left unchanged. By comparison, single-mismatched substrates resulted in up to 10⁶-fold losses of activity with α (Dam) and γ-class (M.HhaI) DNA methyltransferases. We found that CcrM C. crescentus has a greatly expanded DNA-interaction surface, covering six nucleotides on the 5' side and eight nucleotides on the 3' side of its recognition site. Such a large interface may contribute to the enzyme's high sequence fidelity. CcrM C. crescentus displayed the same sequence discrimination with single-stranded substrates, and a surprisingly large (>10⁷-fold) discrimination against ssRNA was largely due to the presence of two or more riboses within the cognate (DNA) site but not outside the site. Results from C-terminal truncations and point mutants supported our hypothesis that the recently identified C-terminal, 80-residue segment is essential for dsDNA recognition but is not required for single-stranded substrates. CcrM orthologs from Agrobacterium tumefaciens and Brucella abortus share some of these newly discovered features of the C. crescentus enzyme, suggesting that the recognition mechanism is conserved. In summary, CcrM C. crescentus uses a previously unknown DNA recognition mechanism.

Isospecific adenine DNA methyltransferases show distinct preferences towards DNA substrates

Here, we report results on systematic analysis of DNA substrate preferences of three N6-adenine β-class DNA methyltransferases that are part of the type II restriction-modification systems. The studied enzymes were: M.EcoVIII, M.HindIII and M.LlaCI, which although found in phylogenetically distant bacteria (γ-proteobacteria and low-GC Gram-positive bacteria), recognize the same palindromic specific sequence 5′-AAGCTT-3′ and catalyze formation of N6-methyladenine at the first A-residue. As expected overall the enzymes share the most analyzed features, but they show also some distinct differences in substrate recognition. Therefore DNA methylation reactions were carried out not only under standard, but also under relaxed conditions using DMSO or glycerol. We found that all of these enzymes preferred DNA containing a hemimethylated target site, but differ in modification of ssDNA, especially more pronounced for M.EcoVIII under relaxed conditions. In these conditions they also have shown varied preferences toward secondary sites, which differ by one nucleotide from specific sequence. They preferred sequences with substitutions at the 1^st (A¹ → G/C) and at the 2^nd position (A² → C), while sites with substitutions at the 3^rd position (G³ → A/C) were modified less efficiently. Kinetic parameters of the methylation reaction carried out by M.EcoVIII were determined. Methylation efficiency (k_cat/K_m) of secondary sites was 4.5–10 times lower when compared to the unmethylated specific sequences, whilst efficiency observed for the hemimethylated substrate was almost 4.5 times greater. We also observed a distinct effect of analyzed enzymes on unspecific interaction with DNA phosphate backbone. We concluded that for all three enzymes the most critical is the phosphodiester bond between G³-C⁴ nucleotides at the center of the target site.

Caulobacter crescentus Cell Cycle-Regulated DNA Methyltransferase Uses a Novel Mechanism for Substrate Recognition

Caulobacter crescentus relies on DNA methylation by the cell cycle-regulated methyltransferase (CcrM) in addition to key transcription factors to control the cell cycle and direct cellular differentiation. CcrM is shown here to efficiently methylate its cognate recognition site 5'-GANTC-3' in single-stranded and hemimethylated double-stranded DNA. We report the K_m, k_cat, k_methylation, and K_d for single-stranded and hemimethylated substrates, revealing discrimination of 10⁷-fold for noncognate sequences. The enzyme also shows a similar discrimination against single-stranded RNA. Two independent assays clearly show that CcrM is highly processive with single-stranded and hemimethylated DNA. Collectively, the data provide evidence that CcrM and other DNA-modifying enzymes may use a new mechanism to recognize DNA in a key epigenetic process.

Fidelity index determination of DNA methyltransferases

DNA methylation is the most frequent form of epigenetic modification in the cell, which involves gene regulation in eukaryotes and protection against restriction enzymes in prokaryotes. Even though many methyltransferases exclusively modify their cognate sites, there have been reports of those that exhibit promiscuity. Previous experimental approaches used to characterize these methyltransferases do not provide the exact concentration at which off-target methylation occurs. Here, we present the first reported fidelity index (FI) for a number of DNA methyltransferases. We define the FI as the ratio of the highest amount of methyltransferase that exhibits no star activity (off-target effects) to the lowest amount that exhibits complete modification of the cognate site. Of the methyltransferases assayed, M.MspI and M.AluI exhibited the highest fidelity of ≥250 and ≥500, respectively, and do not show star activity even at very high concentrations. In contrast, M.HaeIII, M.EcoKDam and M.BamHI have the lowest fidelity of 4, 4 and 2, respectively, and exhibit star activity at concentrations close to complete methylation of the cognate site. The fidelity indexes provide vital information on the usage of methyltransferases and are especially important in applications where site specific methylation is required.

Enzyme-promoted base flipping controls DNA methylation fidelity

A quantitative understanding of how conformational transitions contribute to enzyme catalysis and specificity remains a fundamental challenge. A suite of biophysical approaches was used to reveal several transient states of the enzyme-substrate complexes of the model DNA cytosine methyltransferase M.HhaI. Multidimensional, transverse relaxation-optimized nuclear magnetic resonance (NMR) experiments show that M.HhaI has the same conformation with noncognate and cognate DNA sequences. The high-affinity cognatelike mode requires the formation of a subset of protein-DNA interactions that drive the flipping of the target base from the helix to the active site. Noncognate substrates lacking these interactions undergo slow base flipping, and fluorescence tracking of the catalytic loop corroborates the NMR evidence of a loose, nonspecific binding mode prior to base flipping and subsequent closure of the catalytic loop. This slow flipping transition defines the rate-limiting step for the methylation of noncognate sequences. Additionally, we present spectroscopic evidence of an intermediate along the base flipping pathway that has been predicted but never previously observed. These findings provide important details of how conformational rearrangements are used to balance specificity with catalytic efficiency.

Distal structural elements coordinate a conserved base flipping network

One of the most dramatic illustrations of enzymatic promotion of a high-energy intermediate is observed in DNA modification and repair enzymes where an individual base is rotated (flipped) 180° around the deoxyribose-phosphate backbone and into the active site. While the end states have been extensively characterized, experimental techniques have yet to yield a full description of the base flipping process and the role played by the enzyme. The C5 cytosine methyltransferase M.HhaI coordinates an ensemble of reciprocal DNA and enzyme rearrangements to efficiently flip the target cytosine from the DNA helix. We sought to understand the role of individual amino acids during base flipping. Our results demonstrate that M.HhaI initiates base flipping before closure of the catalytic loop and utilizes the conserved serine 85 in the catalytic loop to accelerate flipping and maintain distortion of the DNA backbone. Serine 87, which forms specific contacts within the DNA helix after base flipping, is not involved in the flipping process or in maintaining the catalytically competent complex. At the base of the catalytic loop, glycine 98 acts as a hinge to allow conformational dynamism of the loop and mutation to alanine inhibits stabilization of the closed loop. Our results illustrate how an enzyme utilizes numerous, distal residues in concert to transform substrate recognition into catalysis.

Evolutionary transitions to new DNA methyltransferases through target site expansion and shrinkage

DNA-binding and modifying proteins show high specificity but also exhibit a certain level of promiscuity. Such latent promiscuous activities comprise the starting points for new protein functions, but this hypothesis presents a paradox: a new activity can only evolve if it already exists. How then, do novel activities evolve? DNA methyltransferases, for example, are highly divergent in their target sites, but how transitions toward novel sites occur remains unknown. We performed laboratory evolution of the DNA methyltransferase M.HaeIII. We found that new target sites emerged primarily through expansion of the original site, GGCC, and the subsequent shrinkage of evolved expanded sites. Variants evolved for sites that are promiscuously methylated by M.HaeIII [GG(^A/_T)CC and GGCGCC] carried mutations in ‘gate-keeper’ residues. They could thereby methylate novel target sites such as GCGC and GGATCC that were neither selected for nor present in M.HaeIII. These ‘generalist’ intermediates were further evolved to obtain variants with novel target specificities. Our results demonstrate the ease by which new DNA-binding and modifying specificities evolve and the mechanism by which they occur at both the protein and DNA levels.

Molecular drivers of base flipping during sequence-specific DNA methylation

One step at a time: Substrates containing nucleotide analogues lacking sequence-specific contacts to the C5 methyltransferase M.HhaI were used to probe the role of individual interactions in effecting conformational transitions during base flipping. A segregation of duties, that is, specific recognition and chemomechanical force for base flipping and active site assembly, within the enzyme is confirmed.

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.

A directed evolution design of a GCG-specific DNA hemimethylase

DNA cytosine-5 methyltransferases (C5-MTases) are valuable models to study sequence-specific modification of DNA and are becoming increasingly important tools for biotechnology. Here we describe a structure-guided rational protein design combined with random mutagenesis and selection to change the specificity of the HhaI C5-MTase from GCGC to GCG. The specificity change was brought about by a five-residue deletion and introduction of two arginine residues within and nearby one of the target recognizing loops. DNA protection assays, bisulfite sequencing and enzyme kinetics showed that the best selected variant is comparable to wild-type M.HhaI in terms of sequence fidelity and methylation efficiency, and supersedes the parent enzyme in transalkylation of DNA using synthetic cofactor analogs. The designed C5-MTase can be used to produce hemimethylated CpG sites in DNA, which are valuable substrates for studies of mammalian maintenance MTases.

Coupling sequence-specific recognition to DNA modification

Enzymes that modify DNA are faced with significant challenges in specificity for both substrate binding and catalysis. We describe how single hydrogen bonds between M.HhaI, a DNA cytosine methyltransferase, and its DNA substrate regulate the positioning of a peptide loop which is approximately 28 A away. Stopped-flow fluorescence measurements of a tryptophan inserted into the loop provide real-time observations of conformational rearrangements. These long-range interactions that correlate with substrate binding and critically, enzyme turnover, will have broad application to enzyme specificity and drug design for this medically relevant class of 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.