Base Flipping Research Articles

Structural analysis of the BisI family of modification dependent restriction endonucleases.

The BisI family of restriction endonucleases is unique in requiring multiple methylated or hydroxymethylated cytosine residues within a short recognition sequence (GCNGC), and in cleaving directly within this sequence, rather than at a distance. Here, we report that the number of modified cytosines that are required for cleavage can be tuned by the salt concentration. We present crystal structures of two members of the BisI family, NhoI and Eco15I_Ntd (N-terminal domain of Eco15I), in the absence of DNA and in specific complexes with tetra-methylated GCNGC target DNA. The structures show that NhoI and Eco15I_Ntd sense modified cytosine bases in the context of double-stranded DNA (dsDNA) without base flipping. In the co-crystal structures of NhoI and Eco15I_Ntd with DNA, the internal methyl groups (G5mCNGC) interact with the side chains of an (H/R)(V/I/T/M) di-amino acid motif near the C-terminus of the distal enzyme subunit and arginine residue from the proximal subunit. The external methyl groups (GCNG5mC) interact with the proximal enzyme subunit, mostly through main chain contacts. Surface plasmon resonance analysis for Eco15I_Ntd shows that the internal and external methyl binding pockets contribute about equally to sensing of cytosine methyl groups.

Structural Basis of Nucleic Acid Recognition and 6mA Demethylation by Caenorhabditis elegans NMAD-1A.

N6-methyladenine (6mA) of DNA is an emerging epigenetic mark in the genomes of Chlamydomonas, Caenorhabditis elegans, and mammals recently. Levels of 6mA undergo drastic fluctuation and thus affect fertility during meiosis and early embryogenesis. Here, we showed three complex structures of 6mA demethylase C. elegans NMAD-1A, a canonical isoform of NMAD-1 (F09F7.7). Biochemical results revealed that NMAD-1A prefers 6mA Bubble or Bulge DNAs. Structural studies of NMAD-1A revealed an unexpected "stretch-out" conformation of its Flip2 region, a conserved element that is usually bent over the catalytic center to facilitate substrate base flipping in other DNA demethylases. Moreover, the wide channel between the Flip1 and Flip2 of the NMAD-1A explained the observed preference of NMAD-1A for unpairing substrates, of which the flipped 6mA was primed for catalysis. Structural analysis and mutagenesis studies confirmed that key elements such as carboxy-terminal domain (CTD) and hypothetical zinc finger domain (ZFD) critically contributed to structural integrity, catalytic activity, and nucleosome binding. Collectively, our biochemical and structural studies suggest that NMAD-1A prefers to regulate 6mA in the unpairing regions and is thus possibly associated with dynamic chromosome regulation and meiosis regulation.

Mapping the recognition pathway of cyclobutane pyrimidine dimer in DNA by Rad4/XPC.

UV radiation-induced DNA damages have adverse effects on genome integrity and cellular function. The most prevalent UV-induced DNA lesion is the cyclobutane pyrimidine dimer (CPD), which can cause skin disorders and cancers in humans. Rad4/XPC is a damage sensing protein that recognizes and repairs CPD lesions with high fidelity. However, the molecular mechanism of how Rad4/XPC interrogates CPD lesions remains elusive. Emerging viewpoints indicate that the association of Rad4/XPC with DNA, the insertion of a lesion-sensing β-hairpin of Rad4/XPC into the lesion siteand the flipping of CPD's partner bases (5'-dA and 3'-dA) are essential for damage recognition. Characterizing these slow events is challenging due to their infrequent occurrence on molecular time scales. Herein, we have used enhanced sampling and molecular dynamics simulations to investigate the mechanism and energetics of lesion recognition by Rad4/XPC, considering multiple plausible pathways between the crystal structure of the Rad4-DNA complex and nine intermediate states. Our results shed light on the most likely sequence of events, their potential coupling and energetics. Upon association, Rad4 and DNA form an encounter complex in which CPD and its partner bases remain in the duplex and the BHD3 β-hairpin is yet to be inserted into the lesion site. Subsequently, sequential base flipping occurs, with the flipping of the 5'-dA base preceding that of the 3'-dA base, followed by the insertion of the BHD3 β-hairpin into the lesion site. The results presented here have significant implications for understanding the molecular basis of UV-related skin disorders and cancers and for paving the way for novel therapeutic strategies.

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.

Register-Shifted Structures in Base-Flipped Uracil-Damaged DNA.

We report the occurrence of register-shifted structures in simulations of uracil-containing dsDNA. These occur when the 3' base vicinal to uracil is thymine in U:A base-paired DNA. Upon base flipping of uracil, this 3' thymine hydrogen bonds with the adenine across the uracil instead of its complementary base. The register-shifted structure is persistent and sterically blocks re-entry of uracil into the helix stack. Register shifting might be important for DNA repair since the longer exposure of the lesion in register-shifted structures could facilitate enzymatic recognition and repair.

Evaluating orthogonality between ion-pair reversed phase, anion exchange, and hydrophilic interaction liquid chromatography for the separation of synthetic oligonucleotides

The orthogonality of separation between ion-pair reversed phase (IP-RP), anion exchange (AEX), and hydrophilic interaction liquid chromatography (HILIC) was evaluated for oligonucleotides. A polythymidine standard ladder was first used to evaluate the three methods and showed zero orthogonality, where retention and selectivity were based on oligonucleotide charge/size under all three conditions. Next, a model 23-mer synthetic oligonucleotide containing 4 phosphorothioate bonds with 2′ fluoro and 2′-O-methyl ribose modifications typical of small interfering RNA was used for evaluating orthogonality. The resolution and orthogonality were evaluated between the three modes of chromatography in terms of selectivity differences for nine common impurities, including truncations (n-1, n-2), addition (n + 1), oxidation, and de-fluorination. We first evaluated different ion-pairing reagents that provided the best separation of the key impurities while suppressing diastereomer separation due to phosphorothioate linkages. Although different ion-pairing reagents affected resolution, very little orthogonality was observed. We then compared the retention times between IP-RP, HILIC, and AEX for each impurity of the model oligonucleotide and observed various selectivity changes. The results suggest that coupling HILIC with either AEX or IP-RP provide the highest degree of orthogonality due to the differences in retention for hydrophilic nucleobases and modifications under HILIC conditions. IP-RP provided the highest overall resolution for the impurity mixture, whereas more co-elution was observed with HILIC and AEX. The unique selectivity patterns offered by HILIC provides an interesting alternative to IP-RP or AEX, in addition to the potential for coupling with multidimensional separations. Future work should explore orthogonality for oligonucleotides with subtle sequence differences such as nucleobase modifications and base flip isomers, longer strands such as guide RNA and messenger RNA, and other biotherapeutic modalities such as peptides, antibodies, and antibody-drug-conjugates.

Systematic assessment of the flexibility of uracil damaged DNA

Uracil is a common DNA lesion which is recognized and removed by uracil DNA-glycosylase (UDG) as a part of the base excision repair pathway. Excision proceeds by base flipping, and UDG efficiency is thought to depend on the ease of deformability of the bases neighboring the lesion. We used molecular dynamics simulations to assess the flexibility of a large library of dsDNA strands, containing all tetranucleotide motifs with U:A, U:G, T:A or C:G base pairs. Our study demonstrates that uracil damaged DNA largely follows trends in flexibility of undamaged DNA. Measured bending persistence lengths, groove widths, step parameters and base flipping propensities demonstrate that uracil increases the flexibility of DNA, and that U:G base paired strands are more flexible than U:A strands. Certain sequence contexts are more deformable than others, with a key role for the 3′ base next to uracil. Flexibilities are large when this base is an A or G, and repressed for a C or T. A 5′ T adjacent to the uracil strongly promotes flexibility, but other 5′ bases are less influential. DNA bending is correlated to step deformations and base flipping, and bending aids flipping. Our study implies that the link between substrate flexibility and UDG efficiency is widely valid, helps explain why UDG prefers to bind U:G base paired strands, and suggests that the DNA bending angle of the UDG-substrate complex is optimal for base flipping. Communicated by Ramaswamy H. Sarma

Inhibitors of UHRF1 base flipping activity showing cytotoxicity against cancer cells

Ubiquitin-like containing PHD and RING finger domain 1 (UHRF1) is a nuclear multi-domain protein overexpressed in numerous human cancer types. We previously disclosed the anthraquinone derivative UM63 that inhibits UHRF1-SRA domain base-flipping activity, although having DNA intercalating properties. Herein, based on the UM63 structure, new UHRF1-SRA inhibitors were identified through a multidisciplinary approach, combining molecular modelling, biophysical assays, molecular and cell biology experiments. We identified AMSA2 and MPB7, that inhibit UHRF1-SRA mediated base flipping at low micromolar concentrations, but do not intercalate into DNA, which is a key advantage over UM63. These molecules prevent UHRF1/DNMT1 interaction at replication forks and decrease the overall DNA methylation in cells. Moreover, both compounds specifically induce cell death in numerous cancer cell lines, displaying marginal effect on non-cancer cells, as they preferentially affect cells with high level of UHRF1. Overall, these two compounds are promising leads for the development of anti-cancer drugs targeting UHRF1.

Uracil-DNA glycosylase efficiency is modulated by substrate rigidity

Uracil DNA-glycosylase (UNG) is a DNA repair enzyme that removes the highly mutagenic uracil lesion from DNA using a base flipping mechanism. Although this enzyme has evolved to remove uracil from diverse sequence contexts, UNG excision efficiency depends on DNA sequence. To provide the molecular basis for rationalizing UNG substrate preferences, we used time-resolved fluorescence spectroscopy, NMR imino proton exchange measurements, and molecular dynamics simulations to measure UNG specificity constants (kcat/KM) and DNA flexibilities for DNA substrates containing central AUT, TUA, AUA, and TUT motifs. Our study shows that UNG efficiency is dictated by the intrinsic deformability around the lesion, establishes a direct relationship between substrate flexibility modes and UNG efficiency, and shows that bases immediately adjacent to the uracil are allosterically coupled and have the greatest impact on substrate flexibility and UNG activity. The finding that substrate flexibility controls UNG efficiency is likely significant for other repair enzymes and has major implications for the understanding of mutation hotspot genesis, molecular evolution, and base editing.

Uracil-DNA glycosylase efficiency is modulated by substrate rigidity.

DNA sequence effects have been identified for many glycosylases and alkyl-transferases that repair base mismatches, uracil bases, alkylated bases, etc. Here, we focused on understanding how DNA sequence impacts the repair of uracil, a highly mutagenic and common lesion in DNA. Uracil DNA-glycosylase (UNG) is a base excision repair enzyme that removes uracil from DNA by a base flipping mechanism. UNG excision efficiency depends on DNA sequence. We show that UNG efficiency is dictated by the intrinsic local deformability of the substrate sequence around the uracil. UNG specificity constants (kcat/KM) and DNA flexibilities were measured for an engineered set of DNA substrates. Time-resolved fluorescence spectroscopy, NMR imino proton exchange measurements, and molecular dynamics simulations of the bare DNA indicated significant differences in substrate flexibilities. We observed a strong correlation between UNG efficiency and substrate flexibility, with higher kcat/KM values measured for more flexible strands. DNA bending and base flipping were observed in simulations, with more frequent uracil flipping observed for the more bendable sequences. Experiments show that bases immediately adjacent to the uracil are allosterically coupled and have the greatest impact on substrate flexibility and resultant UNG activity. A better understanding of the fundamental principles that dictate UNG activity has broad implications in diverse fields, from cancer to evolution.

Structural and functional determinants of the archaeal 8-oxoguanine-DNA glycosylase AGOG for DNA damage recognition and processing.

8-Oxoguanine (GO) is a major purine oxidation product in DNA. Because of its highly mutagenic properties, GO absolutely must be eliminated from DNA. To do this, aerobic and anaerobic organisms from the three kingdoms of life have evolved repair mechanisms to prevent its deleterious effect on genetic integrity. The major way to remove GO is the base excision repair pathway, usually initiated by a GO-DNA glycosylase. First identified in bacteria (Fpg) and eukaryotes (OGG1), GO-DNA glycosylases were more recently identified in archaea (OGG2 and AGOG). AGOG is the less documented enzyme and its mode of damage recognition and removing remains to be clarified at the molecular and atomic levels. This study presents a complete structural characterisation of apo AGOGs from Pyrococcus abyssi (Pab) and Thermococcus gammatolerans (Tga) and the first structure of Pab-AGOG bound to lesion-containing single- or double-stranded DNA. By combining X-ray structure analysis, site directed mutagenesis and biochemistry experiments, we identified key amino acid residues of AGOGs responsible for the specific recognition of the lesion and the base opposite the lesion and for catalysis. Moreover, a unique binding mode of GO, involving double base flipping, never observed for any other DNA glycosylases, is revealed. In addition to unravelling the properties of AGOGs, our study, through comparative biochemical and structural analysis, offers new insights into the evolutionary plasticity of DNA glycosylases across all three kingdoms of life.

Extension of the CHARMM Classical Drude Polarizable Force Field to N- and O-Linked Glycopeptides and Glycoproteins.

Molecular dynamic simulations are an effective tool to study complex molecular systems and are contingent upon the availability of an accurate and reliable molecular mechanics force field. The Drude polarizable force field, which allows for the explicit treatment of electronic polarization in a computationally efficient fashion, has been shown to reproduce experimental properties that were difficult or impossible to reproduce with the CHARMM additive force field, including peptide folding cooperativity, RNA hairpin structures, and DNA base flipping. Glycoproteins are essential components of glycoconjugate vaccines, antibodies, and many pharmaceutically important molecules, and an accurate polarizable force field that includes compatibility between the protein and carbohydrate aspect of the force field is essential to study these types of systems. In this work, we present an extension of the Drude polarizable force field to glycoproteins, including both N- and O-linked species. Parameter optimization focused on the dihedral terms using a reweighting protocol targeting NMR solution J-coupling data for model glycopeptides. Validation of the model include eight model glycopeptides and four glycoproteins with multiple N- and O-linked glycosylations. The new glycoprotein carbohydrate force field can be used in conjunction with the remainder of Drude polarizable force field through a variety of MD simulation programs including GROMACS, OPENMM, NAMD, and CHARMM and may be accessed through the Drude Prepper module in the CHARMM-GUI.

Markov state models elucidate the stability of DNA influenced by the chiral 5S-Tg base.

The static and dynamic structures of DNA duplexes affected by 5S-Tg (Tg, Thymine glycol) epimers were studied using MD simulations and Markov State Models (MSMs) analysis. The results show that the 5S,6S-Tg base caused little perturbation to the helix, and the base-flipping barrier was determined to be 4.4 kcal mol−1 through the use of enhanced sampling meta-eABF calculations, comparable to 5.4 kcal mol−1 of the corresponding thymine flipping. Two conformations with the different hydrogen bond structures between 5S,6R-Tg and A19 were identified in several independent MD trajectories. The 5S,6R-Tg:O6HO6•••N1:A19 hydrogen bond is present in the high-energy conformation displaying a clear helical distortion, and near barrier-free Tg base flipping. The low-energy conformation always maintains Watson–Crick base pairing between 5S,6R-Tg and A19, and 5S-Tg base flipping is accompanied by a small barrier of ca. 2.0 KBT (T = 298 K). The same conformations are observed in the MSMs analysis. Moreover, the transition path and metastable structures of the damaged base flipping are for the first time verified through MSMs analysis. The data clearly show that the epimers have completely different influence on the stability of the DNA duplex, thus implying different enzymatic mechanisms for DNA repair.

A comprehensive thermodynamic model for RNA binding by the Saccharomyces cerevisiae Pumilio protein PUF4

Genomic methods have been valuable for identifying RNA-binding proteins (RBPs) and the genes, pathways, and processes they regulate. Nevertheless, standard motif descriptions cannot be used to predict all RNA targets or test quantitative models for cellular interactions and regulation. We present a complete thermodynamic model for RNA binding to the S. cerevisiae Pumilio protein PUF4 derived from direct binding data for 6180 RNAs measured using the RNA on a massively parallel array (RNA-MaP) platform. The PUF4 model is highly similar to that of the related RBPs, human PUM2 and PUM1, with one marked exception: a single favorable site of base flipping for PUF4, such that PUF4 preferentially binds to a non-contiguous series of residues. These results are foundational for developing and testing cellular models of RNA-RBP interactions and function, for engineering RBPs, for understanding the biophysical nature of RBP binding and the evolutionary landscape of RNAs and RBPs.

Identification of allosteric hotspots regulating the ribosomal RNA binding by antibiotic resistance-conferring Erm methyltransferases

Antibiotic resistance via epigenetic methylation of ribosomal RNA is one of the most prevalent strategies adopted by multidrug resistant pathogens. The erythromycin-resistance methyltransferase (Erm) methylates rRNA at the conserved A2058 position and imparts resistance to macrolides such as erythromycin. However, the precise mechanism adopted by Erm methyltransferases for locating the target base within a complicated rRNA scaffold remains unclear. Here, we show that a conserved RNA architecture, including specific bulge sites, present more than 15 Å from the reaction center, is key to methylation at the pathogenic site. Using a set of RNA sequences site-specifically labeled by fluorescent nucleotide surrogates, we show that base flipping is a prerequisite for effective methylation and that distal bases assist in the recognition and flipping at the reaction center. The Erm–RNA complex model revealed that intrinsically flipped-out bases in the RNA serve as a putative anchor point for the Erm. Molecular dynamic simulation studies demonstrated the RNA undergoes a substantial change in conformation to facilitate an effective protein–rRNA handshake. This study highlights the importance of unique architectural features exploited by RNA to impart fidelity to RNA methyltransferases via enabling allosteric crosstalk. Moreover, the distal trigger sites identified here serve as attractive hotspots for the development of combination drug therapy aimed at reversing resistance.

DNA Deformation Exerted by Regulatory DNA-Binding Motifs in Human Alkyladenine DNA Glycosylase Promotes Base Flipping.

Human alkyladenine DNA glycosylase (AAG) is a key enzyme that corrects a broad range of alkylated and deaminated nucleobases to maintain genomic integrity. When encountering the lesions, AAG adopts a base-flipping strategy to extrude the target base from the DNA duplex to its active site, thereby cleaving the glycosidic bond. Despite its functional importance, the detailed mechanism of such base extrusion and how AAG distinguishes the lesions from an excess of normal bases both remain elusive. Here, through the Markov state model constructed on extensive all-atom molecular dynamics simulations, we find that the alkylated nucleobase (N3-methyladenine, 3MeA) everts through the DNA major groove. Two key AAG motifs, the intercalation and E131-N146 motifs, play active roles in bending/pressing the DNA backbone and widening the DNA minor groove during 3MeA eversion. In particular, the intercalated residue Y162 is involved in buckling the target site at the early stage of 3MeA eversion. Our traveling-salesman based automated path searching algorithm further revealed that a non-target normal adenine tends to be trapped in an exo site near the active site, which however barely exists for a target base 3MeA. Collectively, these results suggest that the Markov state model combined with traveling-salesman based automated path searching acts as a promising approach for studying complex conformational changes of biomolecules and dissecting the elaborate mechanism of target recognition by this unique enzyme.

Thienoguanosine, a unique non-perturbing reporter for investigating rotational dynamics of DNA duplexes and their complexes with proteins

Time-resolved fluorescence anisotropy (TRFA) provides key information on the dynamics of biomolecules and their interaction with ligands. However, since natural nucleosides are almost non-fluorescent, its application to DNA duplexes (dsDNA) requires fluorescent labels, which can alter dsDNA stability, hinder protein binding, and complicate interpretation of TRFA experiments due to their local motion. As shown here, thienoguanosine (thG), a fluorescent analogue of guanosine, overcomes all these limitations. We recorded the TRFA decays of thG-labelled dsDNA of different lengths. thG behaved as a rigid, non-perturbing reporter, since no fast correlation time was recorded for any tested dsDNA. Due to its perfect stacking, only two correlation times, instead of the typical three, describe thG-labelled dsDNA rotational dynamics. Thanks to these features, we provided a complete description of the tumbling of the different dsDNA and their complexes with the Set and Ring Associated (SRA) domain of UHRF1, a key epigenetic regulator, obtaining values in excellent agreement with theoretical predictions. Moreover, thG was also found sensitive to SRA-induced base flipping of neighboring nucleobases. In the DNA label toolbox, thG thus stands out as a unique reporter for investigating the rotational dynamics of dsDNA and protein/dsDNA complexes.

CRISPR-Cas9 bends and twists DNA to read its sequence.

In bacterial defense and genome editing applications, the CRISPR-associated protein Cas9 searches millions of DNA base pairs to locate a 20-nucleotide, guide-RNA-complementary target sequence that abuts a protospacer-adjacent motif (PAM). Target capture requires Cas9 to unwind DNA at candidate sequences using an unknown ATP-independent mechanism. Here we show that Cas9 sharply bends and undertwists DNA upon PAM binding, thereby flipping DNA nucleotides out of the duplex and toward the guide RNA for sequence interrogation. Cryo-electron-microscopy (EM) structures of Cas9:RNA:DNA complexes trapped at different states of the interrogation pathway, together with solution conformational probing, reveal that global protein rearrangement accompanies formation of an unstacked DNA hinge. Bend-induced base flipping explains how Cas9 “reads” snippets of DNA to locate target sites within a vast excess of non-target DNA, a process crucial to both bacterial antiviral immunity and genome editing. This mechanism establishes a physical solution to the problem of complementarity-guided DNA search and shows how interrogation speed and local DNA geometry may influence genome editing efficiency.

50S subunit recognition and modification by the Mycobacterium tuberculosis ribosomal RNA methyltransferase TlyA

Changes in bacterial ribosomal RNA (rRNA) methylation status can alter the activity of diverse groups of ribosome-targeting antibiotics. These modifications are typically incorporated by a single methyltransferase that acts on one nucleotide target and rRNA methylation directly prevents drug binding, thereby conferring drug resistance. Loss of intrinsic methylation can also result in antibiotic resistance. For example, Mycobacterium tuberculosis becomes sensitized to tuberactinomycin antibiotics, such as capreomycin and viomycin, due to the action of the intrinsic methyltransferase TlyA. TlyA is unique among antibiotic resistance-associated methyltransferases as it has dual 16S and 23S rRNA substrate specificity and can incorporate cytidine-2′-O-methylations within two structurally distinct contexts. Here, we report the structure of a mycobacterial 50S subunit-TlyA complex trapped in a postcatalytic state with a S-adenosyl-L-methionine analog using single-particle cryogenic electron microscopy. Together with complementary functional analyses, this structure reveals critical roles in 23S rRNA substrate recognition for conserved residues across an interaction surface that spans both TlyA domains. These interactions position the TlyA active site over the target nucleotide C2144, which is flipped from 23S Helix 69 in a process stabilized by stacking of TlyA residue Phe157 on the adjacent A2143. Base flipping may thus be a common strategy among rRNA methyltransferase enzymes, even in cases where the target site is accessible without such structural reorganization. Finally, functional studies with 30S subunit suggest that the same TlyA interaction surface is employed to recognize this second substrate, but with distinct dependencies on essential conserved residues.

Structural insights into DNMT5-mediated ATP-dependent high-fidelity epigenome maintenance

Epigenetic evolution occurs over million-year timescales in Cryptococcus neoformans and is mediated by DNMT5, the first maintenance type cytosine methyltransferase identified in the fungal or protist kingdoms, the first dependent on adenosine triphosphate (ATP), and the most hemimethyl-DNA-specific enzyme known. To understand these novel properties, we solved cryo-EM structures of CnDNMT5 in three states. These studies reveal an elaborate allosteric cascade in which hemimethylated DNA binding first activates the SNF2 ATPase domain by a large rigid body rotation while the target cytosine partially flips out of the DNA duplex. ATP binding then triggers striking structural reconfigurations of the methyltransferase catalytic pocket to enable cofactor binding, completion of base flipping, and catalysis. Bound unmethylated DNA does not open the catalytic pocket and is instead ejected upon ATP binding, driving high fidelity. This unprecedented chaperone-like, enzyme-remodeling role of the SNF2 ATPase domain illuminates how energy is used to enable faithful epigenetic memory.