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

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.

Relaxed specificity of prokaryotic DNA methyltransferases results in DNA site-specific modification of RNA/DNA heteroduplexes

RNA/DNA hybrid duplexes regularly occur in nature, for example in transcriptional R loops. Their susceptibility to modification by DNA-specific or RNA-specific enzymes is, thus, a biologically relevant question, which, in addition, has possible biotechnological implications. In this study, we investigated the activity of four isospecific DNA methyltransferases (M.EcoVIII, M.LlaCI, M.HindIII, M.BstZ1II) toward an RNA/DNA duplex carrying one 5'-AAGCUU-3'/3'-TTCGAA-5' target sequence. The analyzed enzymes belong to the β-group of adenine N6-methyltransferases and recognize the palindromic DNA sequence 5'-AAGCTT-3'/3'-TTCGAA-5'. Under standard conditions, none of these isospecific enzymes could detectibly methylate the RNA/DNA duplex. However, the addition of agents that generally relax specificity, such as dimethyl sulfoxide (DMSO) and glycerol, resulted in substantial methylation of the RNA/DNA duplex by M.EcoVIII and M.LlaCI. Only the DNA strand of the RNA/DNA duplex was methylated. The same was not observed for M.HindIII or M.BstZ1II. This is, to our knowledge, the first report that demonstrates such activity by prokaryotic DNA methyltransferases. Possible applications of these findings in a laboratory practice are also discussed.

Metadynamics simulation study on the conformational transformation of HhaI methyltransferase: an induced-fit base-flipping hypothesis

DNA methyltransferases play crucial roles in establishing and maintenance of DNA methylation, which is an important epigenetic mark. Flipping the target cytosine out of the DNA helical stack and into the active site of protein provides DNA methyltransferases with an opportunity to access and modify the genetic information hidden in DNA. To investigate the conversion process of base flipping in the HhaI methyltransferase (M.HhaI), we performed different molecular simulation approaches on M.HhaI-DNA-S-adenosylhomocysteine ternary complex. The results demonstrate that the nonspecific binding of DNA to M.HhaI is initially induced by electrostatic interactions. Differences in chemical environment between the major and minor grooves determine the orientation of DNA. Gln237 at the target recognition loop recognizes the GCGC base pair from the major groove side by hydrogen bonds. In addition, catalytic loop motion is a key factor during this process. Our study indicates that base flipping is likely to be an “induced-fit” process. This study provides a solid foundation for future studies on the discovery and development of mechanism-based DNA methyltransferases regulators.

Elucidating molecular motion through structural and dynamic filters of energy-minimized conformer ensembles

Complex RNA structures are constructed from helical segments connected by flexible loops that move spontaneously and in response to binding of small molecule ligands and proteins. Understanding the conformational variability of RNA requires the characterization of the coupled time evolution of interconnected flexible domains. To elucidate the collective molecular motions and explore the conformational landscape of the HIV-1 TAR RNA, we describe a new methodology that utilizes energy-minimized structures generated by the program “Fragment Assembly of RNA with Full-Atom Refinement (FARFAR)”. We apply structural filters in the form of experimental residual dipolar couplings (RDCs) to select a subset of discrete energy-minimized conformers and carry out principal component analyses (PCA) to corroborate the choice of the filtered subset. We use this subset of structures to calculate solution T₁ and T_1ρ relaxation times for ¹³C spins in multiple residues in different domains of the molecule using two simulation protocols that we previously published. We match the experimental T₁ times to within 2% and the T_1ρ times to within less than 10% for helical residues. These results introduce a protocol to construct viable dynamic trajectories for RNA molecules that accord well with experimental NMR data and support the notion that the motions of the helical portions of this small RNA can be described by a relatively small number of discrete conformations exchanging over time scales longer than 1 μs.

Structure of a Naegleria Tet-like dioxygenase in complex with 5-methylcytosine DNA

Cytosine residues in mammalian DNA occur in five forms, cytosine (C), 5-methylcytosine (5mC), 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC). The ten-eleven translocation (Tet) dioxygenases convert 5mC to 5hmC, 5fC and 5caC in three consecutive, Fe(II)- and α-ketoglutarate-dependent oxidation reactions^1–4. The Tet family of dioxygenases is widely distributed across the tree of life⁵, including the heterolobosean amoeboflagellate Naegleria gruberi. The genome of Naegleria⁶ encodes homologs of mammalian DNA methyltransferase and Tet proteins⁷. Here we study biochemically and structurally one of the Naegleria Tet-like proteins (NgTet1), which shares significant sequence conservation (approximately 14% identity or 39% similarity) with mammalian Tet1. Like mammalian Tet proteins, NgTet1 acts on 5mC and generates 5hmC, 5fC and 5caC. The crystal structure of NgTet1 complexed with DNA containing a 5mCpG site revealed that NgTet1 uses a base-flipping mechanism to access 5mC. The DNA is contacted from the minor groove and bent towards the major groove. The flipped 5mC is positioned in the active site pocket with planar stacking contacts, Watson–Crick polar hydrogen bonds and van der Waals interactions specific for 5mC. The sequence conservation between NgTet1 and mammalian Tet1, including residues involved in structural integrity and functional significance, suggests structural conservation across phyla.

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.

Lesion processing by a repair enzyme is severely curtailed by residues needed to prevent aberrant activity on undamaged DNA

DNA base excision repair is essential for maintaining genomic integrity and for active DNA demethylation, a central element of epigenetic regulation. A key player is thymine DNA glycosylase (TDG), which excises thymine from mutagenic G·T mispairs that arise by deamination of 5-methylcytosine (mC). TDG also removes 5-formylcytosine and 5-carboxylcytosine, oxidized forms of mC produced by Tet enzymes. Recent studies show that the glycosylase activity of TDG is essential for active DNA demethylation and for embryonic development. Our understanding of how repair enzymes excise modified bases without acting on undamaged DNA remains incomplete, particularly for mismatch glycosylases such as TDG. We solved a crystal structure of TDG (catalytic domain) bound to a substrate analog and characterized active-site residues by mutagenesis, kinetics, and molecular dynamics simulations. The studies reveal how TDG binds and positions the nucleophile (water) and uncover a previously unrecognized catalytic residue (Thr197). Remarkably, mutation of two active-site residues (Ala145 and His151) causes a dramatic enhancement in G·T glycosylase activity but confers even greater increases in the aberrant removal of thymine from normal A·T base pairs. The strict conservation of these residues may reflect a mechanism used to strike a tolerable balance between the requirement for efficient repair of G·T lesions and the need to minimize aberrant action on undamaged DNA, which can be mutagenic and cytotoxic. Such a compromise in G·T activity can account in part for the relatively weak G·T activity of TDG, a trait that could potentially contribute to the hypermutability of CpG sites in cancer and genetic disease.

Understanding the origins of bacterial resistance to aminoglycosides through molecular dynamics mutational study of the ribosomal A-site

Paromomycin is an aminoglycosidic antibiotic that targets the RNA of the bacterial small ribosomal subunit. It binds in the A-site, which is one of the three tRNA binding sites, and affects translational fidelity by stabilizing two adenines (A1492 and A1493) in the flipped-out state. Experiments have shown that various mutations in the A-site result in bacterial resistance to aminoglycosides. In this study, we performed multiple molecular dynamics simulations of the mutated A-site RNA fragment in explicit solvent to analyze changes in the physicochemical features of the A-site that were introduced by substitutions of specific bases. The simulations were conducted for free RNA and in complex with paromomycin. We found that the specific mutations affect the shape and dynamics of the binding cleft as well as significantly alter its electrostatic properties. The most pronounced changes were observed in the U1406C∶U1495A mutant, where important hydrogen bonds between the RNA and paromomycin were disrupted. The present study aims to clarify the underlying physicochemical mechanisms of bacterial resistance to aminoglycosides due to target mutations.

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.

Base-pair opening and bubble transport in a DNA double helix induced by a protein molecule in a viscous medium

The protein-DNA interaction dynamics is studied by modeling the DNA bases as classical spins in a coupled spin system, which are bosonized and coupled to thermal phonons and longitudinal motion of the protein molecule in the nonviscous limit. The nonlinear dynamics of this protein-DNA complex molecular system is governed by the completely integrable nonlinear Schrödinger (NLS) equation which admits N -soliton solutions. The soliton excitations of the DNA bases in the two strands make localized base-pair opening and travel along the DNA chain in the form of a bubble. This may characterize the bubble generated during the transcription process, when an RNA polymerase binds to a promoter site in the DNA double helical chain. When the protein-DNA molecular system interacts with the surrounding viscous solvating water medium, the dynamics is governed by a perturbed NLS equation. This equation is solved using a multiple scale perturbation analysis, by treating the viscous effect as a weak perturbation, and the results show that the viscosity of the solvent medium damps out the soliton as time progresses.

In the absence of writhe, DNA relieves torsional stress with localized, sequence-dependent structural failure to preserve B-form

To understand how underwinding and overwinding the DNA helix affects its structure, we simulated 19 independent DNA systems with fixed degrees of twist using molecular dynamics in a system that does not allow writhe. Underwinding DNA induced spontaneous, sequence-dependent base flipping and local denaturation, while overwinding DNA induced the formation of Pauling-like DNA (P-DNA). The winding resulted in a bimodal state simultaneously including local structural failure and B-form DNA for both underwinding and extreme overwinding. Our simulations suggest that base flipping and local denaturation may provide a landscape influencing protein recognition of DNA sequence to affect, for examples, replication, transcription and recombination. Additionally, our findings help explain results from single-molecule experiments and demonstrate that elastic rod models are strictly valid on average only for unstressed or overwound DNA up to P-DNA formation. Finally, our data support a model in which base flipping can result from torsional stress.

Characterizing structural transitions using localized free energy landscape analysis

BACKGROUND

Structural changes in molecules are frequently observed during biological processes like replication, transcription and translation. These structural changes can usually be traced to specific distortions in the backbones of the macromolecules involved. Quantitative energetic characterization of such distortions can greatly advance the atomic-level understanding of the dynamic character of these biological processes.

METHODOLOGY/PRINCIPAL FINDINGS

Molecular dynamics simulations combined with a variation of the Weighted Histogram Analysis Method for potential of mean force determination are applied to characterize localized structural changes for the test case of cytosine (underlined) base flipping in a GTCAGCGCATGG DNA duplex. Free energy landscapes for backbone torsion and sugar pucker degrees of freedom in the DNA are used to understand their behavior in response to the base flipping perturbation. By simplifying the base flipping structural change into a two-state model, a free energy difference of upto 14 kcal/mol can be attributed to the flipped state relative to the stacked Watson-Crick base paired state. This two-state classification allows precise evaluation of the effect of base flipping on local backbone degrees of freedom.

CONCLUSIONS/SIGNIFICANCE

The calculated free energy landscapes of individual backbone and sugar degrees of freedom expectedly show the greatest change in the vicinity of the flipping base itself, but specific delocalized effects can be discerned upto four nucleotide positions away in both 5' and 3' directions. Free energy landscape analysis thus provides a quantitative method to pinpoint the determinants of structural change on the atomic scale and also delineate the extent of propagation of the perturbation along the molecule. In addition to nucleic acids, this methodology is anticipated to be useful for studying conformational changes in all macromolecules, including carbohydrates, lipids, and proteins.

Molecular dynamics of a kappaB DNA element: base flipping via cross-strand intercalative stacking in a microsecond-scale simulation

The sequence-dependent structural variability and conformational dynamics of DNA play pivotal roles in many biological milieus, such as in the site-specific binding of transcription factors to target regulatory elements. To better understand DNA structure, function, and dynamics in general, and protein···DNA recognition in the ‘κB’ family of genetic regulatory elements in particular, we performed molecular dynamics simulations of a 20-bp DNA encompassing a cognate κB site recognized by the proto-oncogenic ‘c-Rel’ subfamily of NF-κB transcription factors. Simulations of the κB DNA in explicit water were extended to microsecond duration, providing a broad, atomically detailed glimpse into the structural and dynamical behavior of double helical DNA over many timescales. Of particular note, novel (and structurally plausible) conformations of DNA developed only at the long times sampled in this simulation—including a peculiar state arising at ≈0.7 μs and characterized by cross-strand intercalative stacking of nucleotides within a longitudinally sheared base pair, followed (at ≈1 μs) by spontaneous base flipping of a neighboring thymine within the A-rich duplex. Results and predictions from the microsecond-scale simulation include implications for a dynamical NF-κB recognition motif, and are amenable to testing and further exploration via specific experimental approaches that are suggested herein.

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.

Solid-state nuclear magnetic resonance spectroscopy studies of furanose ring dynamics in the DNA HhaI binding site

The dynamics of the furanose rings in the GCGC 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 by using deuterium solid-state NMR (SSNMR). SSNMR spectra obtained from DNAs selectively deuterated on the furanose rings of nucleotides within the 5'-GCGC-3' moiety indicated that all of these positions are structurally flexible. The furanose ring within the deoxycytidine that is the methylation target displays the largest-amplitude structural changes according to the observed deuterium NMR line shapes, whereas the furanose rings of nucleotides more remote from the methylation site have less-mobile furanose rings (i.e., with puckering amplitudes < 0.3 A). Previous work has shown that methylation reduces the amplitude of motion in the phosphodiester backbone of the same DNA, and our observations indicate that methylation perturbs backbone dynamics through the furanose ring. These NMR data indicate that the 5'-GCGC-3' is dynamic, with the largest-amplitude motions occurring nearest the methylation site. The inherent flexibility of this moiety in DNA makes the molecule more amenable to the large-amplitude structural rearrangements that must occur when the DNA binds to the HhaI methyltransferase.

A molecular dynamics study of slow base flipping in DNA using conformational flooding

Individual DNA bases are known to be able to flip out of the helical stack, providing enzymes with access to the genetic information otherwise hidden inside the helix. Consequently, base flipping is a necessary first step to many more complex biological processes such as DNA transcription or replication. Much remains unknown about this elementary step, despite a wealth of experimental and theoretical studies. From the theoretical point of view, the involved timescale of milliseconds or longer requires the use of enhanced sampling techniques. In contrast to previous theoretical studies employing umbrella sampling along a predefined flipping coordinate, this study attempts to induce flipping without prior knowledge of the pathway, using information from a molecular dynamics simulation of a B-DNA fragment and the conformational flooding method. The relevance to base flipping of the principal components of the simulation is assayed, and a combination of modes optimally related to the flipping of the base through either helical groove is derived for each of the two bases of the central guanine-cytosine basepair. By applying an artificial flooding potential along these collective coordinates, the flipping mechanism is accelerated to within the scope of molecular dynamics simulations. The associated free energy surface is found to feature local minima corresponding to partially flipped states, particularly relevant to flipping in isolated DNA; further transitions from these minima to the fully flipped conformation are accelerated by additional flooding potentials. The associated free energy profiles feature similar barrier heights for both bases and pathways; the flipped state beyond is a broad and rugged attraction basin, only a few kcal/mol higher in energy than the closed conformation. This result diverges from previous works but echoes some aspects of recent experimental findings, justifying the need for novel approaches to this difficult problem: this contribution represents a first step in this direction. Important structural factors involved in flipping, both local (sugar-phosphate backbone dihedral angles) and global (helical axis bend), are also identified.

Biophysical studies of DNA modified with conformationally constrained nucleotides: comparison of 2'-exo (north) and 3'-exo (south) 'locked' templates

The biophysical properties of oligodeoxyribonucleotides (ODNs) selectively modified with conformationally ‘locked’ bicyclo[3.1.0]hexane pseudosugars (Maier,M.A., Choi,Y., Gaus,H., Barchi,J.J. Jr, Marquez,V.E., Manoharan,M. (2004) Synthesis and characterization of oligonucleotides containing conformationally constrained bicyclo[3.1.0]hexane pseudosugar analogs Nucleic Acids Res., 32, 3642–3650) have been studied by various techniques. Six separate synthetic ODNs based on the Dickerson Drew dodecamer sequence (CGCGAAT*T*CGCG) were examined where each one (or both) of the thymidines (T*) were substituted with a bicyclic pseudosugar locked in either a North (2′-exo) or South (3′-exo) ring pucker. Circular dichroism spectroscopy, differential scanning calorimetry and ¹H NMR spectroscopy were used to examine the duplex stability and conformational properties of the ODNs. Replacement of one or both thymidines with North-locked sugars (RNA-like) into the dodecamer did not greatly affect duplex formation or melt temperatures but distinct differences in thermodynamic parameters were observed. In contrast, incorporation of South-locked sugar derivatives that were predicted to stabilize this standard B-DNA, had the unexpected effect of causing a conformational equilibrium between different duplex forms at specific strand and salt concentrations. Our data and those of others suggest that although DNA can tolerate modifications with RNA-like (North) nucleotides, a more complicated spectrum of changes emerges with modifications restricted to South (DNA-like) puckers.

Engineered Extrahelical Base Destabilization Enhances Sequence Discrimination of DNA Methyltransferase M.HhaI

Improved sequence specificity of the DNA cytosine methyltransferase HhaI was achieved by disrupting interactions at a hydrophobic interface between the active site of the enzyme and a highly conserved flexible loop. Transient fluorescence experiments show that mutations disrupting this interface destabilize the positioning of the extrahelical, "flipped" cytosine base within the active site. The ternary crystal structure of the F124A M.HhaI bound to cognate DNA and the cofactor analogue S-adenosyl-l-homocysteine shows an increase in cavity volume between the flexible loop and the core of the enzyme. This cavity disrupts the interface between the loop and the active site, thereby destabilizing the extrahelical target base. The favored partitioning of the base-flipped enzyme-DNA complex back to the base-stacked intermediate results in the mutant enzyme discriminating better than the wild-type enzyme against non-cognate sites. Building upon the concepts of kinetic proofreading and our understanding of M.HhaI, we describe how a 16-fold specificity enhancement achieved with a double mutation at the loop/active site interface is acquired through destabilization of intermediates prior to methyltransfer rather than disruption of direct interactions between the enzyme and the substrate for M.HhaI.

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.

Conserved patterns in backbone torsional changes allow for single base flipping from duplex DNA with minimal distortion of the double helix

Base flipping is a structural mechanism common to many DNA processing and repair enzymes. Changes in the local backbone torsions that occur during base flipping and the effect of environment on their behavior are of particular interest in understanding different base flipping mechanisms. In the present study, structures sampled during umbrella sampling molecular dynamics (MD) simulations of base flipping in aqueous and protein-bound environments, carried out with two different MD simulation strategies, are analyzed to find the most significant phosphodiester backbone distortions in the vicinity of the flipping base. Torsional sampling on the 5' side of the flipping base during flipping through the major groove shows similarities to the torsional sampling on the 3' side during flipping through the minor groove and vice versa. In differing environments, this behavior varies only marginally. These compensating torsional changes in the DNA backbone on 5' and 3' sides of the flipping base limit overall distortion of the DNA double helix during single base flipping. Rotameric intermediate states observed during base flipping are identified and postulated to be metastable states implicated in both large-scale structural changes and functional effects of chemical modifications in DNA.

Mutational Analysis of the Catalytic Domain of the Murine Dnmt3a DNA-(cytosine C5)-methyltransferase

On the basis of amino acid sequence alignments and structural data of related enzymes, we have performed a mutational analysis of 14 amino acid residues in the catalytic domain of the murine Dnmt3a DNA-(cytosine C5)-methyltransferase. The target residues are located within the ten conserved amino acid sequence motifs characteristic for cytosine-C5 methyltransferases and in the putative DNA recognition domain of the enzyme (TRD). Mutant proteins were purified and tested for their catalytic properties and their abilities to bind DNA and AdoMet. We prepared a structural model of Dnmt3a to interpret our results. We demonstrate that Phe50 (motif I) and Glu74 (motif II) are important for AdoMet binding and catalysis. D96A (motif III) showed reduced AdoMet binding but increased activity under conditions of saturation with S-adenosyl-L-methionine (AdoMet), indicating that the contact of Asp96 to AdoMet is not required for catalysis. R130A (following motif IV), R241A and R246A (in the TRD), R292A, and R297A (both located in front of motif X) showed reduced DNA binding. R130A displayed a strong reduction in catalytic activity and a complete change in flanking sequence preferences, indicating that Arg130 has an important role in the DNA interaction of Dnmt3a. R292A also displayed reduced activity and changes in the flanking sequence preferences, indicating a potential role in DNA contacts farther away from the CG target site. N167A (motif VI) and R202A (motif VIII) have normal AdoMet and DNA binding but reduced catalytic activity. While Asn167 might contribute to the positioning of residues from motif VI, according to structural data Arg202 has a role in catalysis of cytosine-C5 methyltransferases. The R295A variant was catalytically inactive most likely because of destabilization of the hinge sub-domain of the protein.

Structure and substrate recognition of the Escherichia coli DNA adenine methyltransferase

The structure of the Escherichia coli Dam DNA-(adenine-N6)-methyltransferase in complex with cognate DNA was determined at 1.89 A resolution in the presence of S-adenosyl-L-homocysteine. DNA recognition and the dynamics of base-flipping were studied by site-directed mutagenesis, DNA methylation kinetics and fluorescence stopped-flow experiments. Our data illustrate the mechanism of coupling of DNA recognition and base-flipping. Contacts to the non-target strand in the second (3') half of the GATC site are established by R124 to the fourth base-pair, and by L122 and P134 to the third base-pair. The aromatic ring of Y119 intercalates into the DNA between the second and third base-pairs, which is essential for base-flipping to occur. Compared to previous published structures of bacteriophage T4 Dam, three major new observations are made in E.coli Dam. (1) The first Gua is recognized by K9, removal of which abrogates the first base-pair recognition. (2) The flipped target Ade binds to the surface of EcoDam in the absence of S-adenosyl-L-methionine, which illustrates a possible intermediate in the base-flipping pathway. (3) The orphaned Thy residue displays structural flexibility by adopting an extrahelical or intrahelical position where it is in contact to N120.

Assessment of adenyl residue reactivity within model nucleic acids by surface enhanced Raman spectroscopy

We rank the reactivity of the adenyl residues (A) of model DNA and RNA molecules with electropositive subnano size [Ag]n+ sites as a function of nucleic acid primary sequences and secondary structures and in the presence of biological amounts of Cl- and Na+ or Mg2+ ions. In these conditions A is markedly more reactive than any other nucleic acid bases. A reactivity is higher in ribo (r) than in deoxyribo (d) species [pA>pdA and (pA)n>>(pdA)n]. Base pairing decreases A reactivity in corresponding duplexes but much less in r than in d. In linear single and paired dCAG or dGAC loci, base stacking inhibits A reactivity even if A is bulged or mispaired (A.A). dA tracts are highly reactive only when dilution prevents self-association and duplex structures. In d hairpins the solvent-exposed A residues are reactive in CAG and GAC triloops and even more in ATC loops. Among the eight rG1N2R3A4 loops, those bearing a single A (A4) are the least reactive. The solvent-exposed A2 is reactive, but synergistic structural transitions make the initially stacked A residues of any rGNAA loop much more reactive. Mg2+ cross-bridging single strands via phosphates may screen A reactivity. In contrast d duplexes cross-bridging enables "A flipping" much more in rA.U pairs than in dA.T. Mg2+ promotes A reactivity in unpaired strands. For hairpins Mg2+ binding stabilizes the stems, but according to A position in the loops, A reactivity may be abolished, reduced, or enhanced. It is emphasized that not only accessibility but also local flexibility, concerted docking, and cation and anion binding control A reactivity.

Low-frequency normal mode in DNA HhaI methyltransferase and motions of residues involved in the base flipping

The results of normal-mode analyses are in accord with the proposal that a low-frequency motion of the HhaI methyltransferase enzyme is responsible for base flipping in bound DNA. The vectors of the low-frequency normal mode of residues Ser-85 and Ile-86 point directly to the phosphate and ribose moieties of the DNA backbone near the target base in position to rotate the dihedral angles and flip the base out of the DNA duplex. The vector of residue Gln-237 on the major groove is in the proper orientation to assist base separation. Our results favor the major groove pathway and the protein active process in base flipping.

Structural evidence of a passive base-flipping mechanism for AGT, an unusual GT-B glycosyltransferase

The Escherichia coli T4 bacteriophage uses two glycosyltransferases to glucosylate and thus protect its DNA: the retaining alpha-glucosyltransferase (AGT) and the inverting beta-glucosyltransferase (BGT). They glucosylate 5-hydroxymethyl cytosine (5-HMC) bases of duplex DNA using UDP-glucose as the sugar donor to form an alpha-glucosidic linkage and a beta-glucosidic linkage, respectively. Five structures of AGT have been determined: a binary complex with the UDP product and four ternary complexes with UDP or UDP-glucose and oligonucleotides containing an A:G, HMU:G (hydroxymethyl uracyl) or AP:G (apurinic/apyrimidinic) mismatch at the target base-pair. AGT adopts the GT-B fold, one of the two folds known for GTs. However, while the sugar donor binding mode is classical for a GT-B enzyme, the sugar acceptor binding mode is unexpected and breaks the established consensus: AGT is the first GT-B enzyme that predominantly binds both the sugar donor and acceptor to the C-terminal domain. Its active site pocket is highly similar to four retaining GT-B glycosyltransferases (trehalose-6-phosphate synthase, glycogen synthase, glycogen and maltodextrin phosphorylases) strongly suggesting a common evolutionary origin and catalytic mechanism for these enzymes. Structure-guided mutagenesis and kinetic analysis do not permit identification of a nucleophile residue responsible for a glycosyl-enzyme intermediate for the classical double displacement mechanism. Interestingly, the DNA structures reveal partially flipped-out bases. They provide evidence for a passive role of AGT in the base-flipping mechanism and for its specific recognition of the acceptor base.

Eukaryotic cytosine methyltransferases

Large-genome eukaryotes use heritable cytosine methylation to silence promoters, especially those associated with transposons and imprinted genes. Cytosine methylation does not reinforce or replace ancestral gene regulation pathways but instead endows methylated genomes with the ability to repress specific promoters in a manner that is buffered against changes in the internal and external environment. Recent studies have shown that the targeting of de novo methylation depends on multiple inputs; these include the interaction of repeated sequences, local states of histone lysine methylation, small RNAs and components of the RNAi pathway, and divergent and catalytically inert cytosine methyltransferase homologues that have acquired regulatory roles. There are multiple families of DNA (cytosine-5) methyltransferases in eukaryotes, and each family appears to be controlled by different regulatory inputs. Sequence-specific DNA-binding proteins, which regulate most aspects of gene expression, do not appear to be involved in the establishment or maintenance of genomic methylation patterns.

Molecular dynamics simulations and free energy calculations of base flipping in dsRNA

The family of adenosine deaminases acting on RNA (ADARs) targets adenosines in RNA that is mainly double stranded. Some substrates are promiscuously deaminated whereas others, such as the mammalian glutamate receptor B (gluR-B) pre-mRNA, are more selectively deaminated. Many DNA/RNA-base modification enzymes use a base flipping mechanism to be able to reach their target base and it is believed that ADARs function in a similar way. In this study we used molecular dynamics (MD) simulations to describe two sites on the gluR-B pre-mRNA, the selectively targeted R/G site and the nontargeted 46 site, in an attempt to explain the substrate specificity. We used regular MD and also a forced base flipping method with umbrella sampling to calculate the free energy of base opening. Spontaneous opening of the mismatched adenosine was observed for the R/G site but not for the 46 site.

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.