The mechanism of M.HhaI DNA C5 cytosine methyltransferase enzyme: a quantum mechanics/molecular mechanics approach

Cytosine analogues as DNA methyltransferase substrates

DNA methyltransferases are drug targets for myelodysplastic syndrome (MDS), chronic myelomonocytic leukemia (CMML), acute myelogenous leukemia (AML) and possibly β-hemoglobinopathies. We characterize the interaction of nucleoside analogues in DNA with a prokaryotic CpG-specific DNA methyltransferase (M.MpeI) as a model for mammalian DNMT1 methyltransferases. We tested DNA containing 5-hydroxymethylcytosine (5hmC), 5-hydroxycytosine (5OHC), 5-methyl-2-pyrimidinone (in the ribosylated form known as 5-methylzebularine, 5mZ), 5,6-dihydro-5-azacytosine (dhaC), 5-fluorocytosine (5FC), 5-chlorocytosine (5ClC), 5-bromocytosine (5BrC) and 5-iodocytosine (5IC). Covalent complex formation was by far most efficient for 5FC. Non-covalent complexes were most abundant for dhaC and 5mZ. Surprisingly, we observed methylation of 5IC and 5BrC, and to a lesser extent 5ClC and 5FC, in the presence, but not the absence of small molecule thiol nucleophiles. For 5IC and 5BrC, we demonstrated by mass spectrometry that the reactions were due to methyltransferase driven dehalogenation, followed by methylation. Crystal structures of M.MpeI-DNA complexes capture the 'in' conformation of the active site loop for analogues with small or rotatable (5mZ) 5-substituents and its 'out' form for bulky 5-substituents. Since very similar 'in' and 'out' loop conformations were also observed for DNMT1, it is likely that our conclusions generalize to other DNA methyltransferases.

Insights into the Inhibitory Mechanisms of the Covalent Drugs for DNMT3A

The perturbations of DNA methyltransferase 3 alpha (DNMT3A) may cause uncontrolled gene expression, resulting in cancers and tumors. The DNMT inhibitors Azacytidine (AZA) and Zebularine (ZEB) inhibit the DNMT family with no specificities, and consequently would bring side effects during the treatment. Therefore, it is vital to understand the inhibitory mechanisms in DNMT3A to inform the new inhibitor design for DNMTs. Herein, we carried out molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) simulations to investigate the inhibitory mechanisms of the AZA and ZEB. The results were compared to the methyl transfer of cytosine. We showed how the AZA might stop the methyl transfer process, whereas the ZEB might be stuck in a methyl-transferred intermediate (IM3). The IM3 state then fails the elimination due to the unique protein dynamics that result in missing the catalytic water chain. Our results brought atomic-level insights into the mechanisms of the two drugs in DNMT3A, which could benefit the new generation of drug design for the DNMTs.

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.

Computational Insights in DNA Methylation: Catalytic and Mechanistic Elucidations for Forming 3-Methyl Cytosine

Methylation at C5 position of cytosine (5 mC) is the most abundantly occurring methylation process at CpG island, which has been well known as an epigenetic modification linked to many human diseases. Recently, another methylation approach has been discovered to show that DNA methyltransferases (DNMTs) promote the addition of the methyl group at position 3 to yield 3 mC. The existence of 3 mC can cause severe damages to the DNA strand, such as blocking its replication, repair, and transcription, affecting its stability, and initiating a double-strand DNA break. To gain a deeper insight into the formation of 3 mC, we have performed density functional theory (DFT) modeling studies at different levels of theory to clearly map out the mechanistic details for this new methylation approach. Our computed results are in harmony with pertinent experimental observations and shed light on a crucial off-target activity of DNMTs.

Discovery of an Unnatural DNA Modification Derived from a Natural Secondary Metabolite

Despite widespread interest for understanding how modified bases have evolved their contemporary functions, limited experimental evidence exists for measuring how close an organism is to accidentally creating a new, modified base within the framework of its existing genome. Here, we describe the biochemical and structural basis for how a single-point mutation in E. coli's naturally occurring cytosine methyltransferase can surprisingly endow a neomorphic ability to create the unnatural DNA base, 5-carboxymethylcytosine (5cxmC), in vivo. Mass spectrometry, bacterial genetics, and structure-guided biochemistry reveal this base to be exclusively derived from the natural but sparse secondary metabolite carboxy-S-adenosyl-L-methionine (CxSAM). Our discovery of a new, unnatural DNA modification reveals insights into the substrate selectivity of DNA methyltransferase enzymes, offers a promising new biotechnological tool for the characterization of the mammalian epigenome, and provides an unexpected model for how neomorphic bases could arise in nature from repurposed host metabolites.

Structural insights into CpG-specific DNA methylation by human DNA methyltransferase 3B

DNA methyltransferases are primary enzymes for cytosine methylation at CpG sites of epigenetic gene regulation in mammals. De novo methyltransferases DNMT3A and DNMT3B create DNA methylation patterns during development, but how they differentially implement genomic DNA methylation patterns is poorly understood. Here, we report crystal structures of the catalytic domain of human DNMT3B-3L complex, noncovalently bound with and without DNA of different sequences. Human DNMT3B uses two flexible loops to enclose DNA and employs its catalytic loop to flip out the cytosine base. As opposed to DNMT3A, DNMT3B specifically recognizes DNA with CpGpG sites via residues Asn779 and Lys777 in its more stable and well-ordered target recognition domain loop to facilitate processive methylation of tandemly repeated CpG sites. We also identify a proton wire water channel for the final deprotonation step, revealing the complete working mechanism for cytosine methylation by DNMT3B and providing the structural basis for DNMT3B mutation-induced hypomethylation in immunodeficiency, centromere instability and facial anomalies syndrome.

RNA (C5-cytosine) Methyltransferases

The review summarizes the data on pro- and eukaryotic RNA (C5-cytosine) methyltransferases. The structure, intracellular location, RNA targets, and catalytic mechanisms of these enzymes, as well as the functional role of methylated cytosine residues in RNA are presented. The functions of RNA (C5-cytosine) methyltransferases unassociated with their methylation activity are discussed. Special attention is given to the similarities and differences in the structures and mechanisms of action of RNA and DNA methyltransferases. The data on the association of mutations in the RNA (C5-cytosine) methyltransferases genes and human diseases are presented.

Understanding the R882H mutation effects of DNA methyltransferase DNMT3A: a combination of molecular dynamics simulations and QM/MM calculations

The AML-related high-frequent R882H mutation of DNA (cytosine-5)-methyltransferase 3A (DNMT3A), a key enzyme forde novoepigenetic methylation in human beings, was characterized by a disturbing conformation ofS-adenosylmethionine (SAM).

Structural insights into the specificity and catalytic mechanism of mycobacterial nucleotide pool sanitizing enzyme MutT2

Mis-incorporation of modified nucleotides, such as 5-methyl-dCTP or 8-oxo-dGTP, in DNA can be detrimental to genomic integrity. MutT proteins are sanitization enzymes which function by hydrolyzing such nucleotides and regulating the pool of free nucleotides in the cytoplasm. Mycobacterial genomes have a set of four MutT homologs, namely, MutT1, MutT2, MutT3 and MutT4. Mycobacterial MutT2 hydrolyzes 5 m-dCTP and 8-oxo-dGTP to their respective monophosphate products. Additionally, it can hydrolyze canonical nucleotides dCTP and CTP, with a suggested role in sustaining their optimal levels in the nucleotide pool. The structures of M. smegmatis MutT2 and its complexes with cytosine derivatives have been determined at resolutions ranging from 1.10 Å to 1.73 Å. The apo enzyme and its complexes with products (dCMP, CMP and 5 m-dCMP) crystallize in space group P2₁2₁2, while those involving substrates (dCTP, CTP and 5 m-dCTP) crystallize in space group P2_1. The molecule takes an α/β/α sandwich fold arrangement, as observed in other MutT homologs. The nucleoside moiety of the ligands is similarly located in all the complexes, while the location of the remaining tail exhibits variability. This is the first report of a MutT2-type protein in complex with ligands. A critical interaction involving Asp116 confers the specificity of the enzyme towards cytosine moieties. A conserved set of enzyme-ligand interactions along with concerted movements of important water molecules provide insights into the mechanism of action.

Methyltransferase-Directed Labeling of Biomolecules and its Applications

Methyltransferases (MTases) form a large family of enzymes that methylate a diverse set of targets, ranging from the three major biopolymers to small molecules. Most of these MTases use the cofactor S-adenosyl-l-Methionine (AdoMet) as a methyl source. In recent years, there have been significant efforts toward the development of AdoMet analogues with the aim of transferring moieties other than simple methyl groups. Two major classes of AdoMet analogues currently exist: doubly-activated molecules and aziridine based molecules, each of which employs a different approach to achieve transalkylation rather than transmethylation. In this review, we discuss the various strategies for labelling and functionalizing biomolecules using AdoMet-dependent MTases and AdoMet analogues. We cover the synthetic routes to AdoMet analogues, their stability in biological environments and their application in transalkylation reactions. Finally, some perspectives are presented for the potential use of AdoMet analogues in biology research, (epi)genetics and nanotechnology.

Structure and dynamics of H. pylori 98-10 C5-cytosine specific DNA methyltransferase in complex with S-adenosyl-l-methionine and DNA

Helicobacter pylori is a Gram-negative bacterium that inhabits the human gastrointestinal tract, and some strains of this bacterium cause gastric ulcers and cancer. DNA methyltransferases (MTases) are promising drug targets for the treatment of cancer and other diseases that are also caused by epigenetic alternations of the genome. The C5-cytosine specific DNA methyltransferase from H. pylori (M. Hpy C5mC) catalyzes the transfer of the methyl group from the cofactor S-adenosyl-l-methionine (AdoMet) to the flipped cytosine of the substrate DNA. Herein we report the sequence analyses, 3-D structure modeling and molecular dynamics simulations of M. Hpy C5mC, when complexed with AdoMet as well as DNA. We analyzed the protein-DNA interactions prominently established by the flipped cytosine and the interactions between the protein and cofactor in the active site. We propose that the contacts made by cytosine O2 with Arg155 and Arg157, and the water-mediated interactions with cytosine N3 may be essential for the activity of methyl transfer as well as the deprotonation at the C5 position in our C5mC model. Specific recognition of DNA was mediated mainly by residues from Ser221-Arg229 and Ser243-Gln246 of the target recognition domain (TRD) and some residues of the loop Ser75-Lys83 from the large domain. These findings are further supported by alanine scanning mutagenesis studies. The results reported here explain the sequence, structure and binding features necessary for the recognition between the cofactor and the substrate by the key epigenetic enzyme, M. Hpy C5mC.

Conserved motif VIII of murine DNA methyltransferase Dnmt3a is essential for methylation activity

BACKGROUND

Dnmt3a is a DNA methyltransferase that establishes de novo DNA methylation in mammals. The structure of the Dnmt3a C-terminal domain is similar to the bacterial M. HhaI enzyme, a well-studied prokaryotic DNA methyltransferase. No X-ray structure is available for the complex of Dnmt3a with DNA and the mechanistic details of DNA recognition and catalysis by mammalian Dnmts are not completely understood.

RESULTS

Mutant variants of the catalytic domain of the murine Dnmt3a carrying substitutions of highly conserved N167, R200, and R202 have been generated by site directed mutagenesis and purified. Their methylation activity, DNA binding affinity, ability to flip the target cytosine out of the DNA double helix and covalent complex formation with DNA have been examined. Substitutions of N167 lead to reduced catalytic activity and reduced base flipping. Catalytic activity, base flipping, and covalent conjugate formation were almost completely abolished for the mutant enzymes with substitutions of R200 or R202.

CONCLUSIONS

We conclude that R202 plays a similar role in catalysis in Dnmt3a-CD as R232 in M.SssI and R165 in M.HhaI, which could be positioning of the cytosine for nucleophilic attack by a conserved Cys. R200 of Dnmt3a-CD is important in both catalysis and cytosine flipping. Both conserved R200 and R202 are involved in creating and stabilizing of the transient covalent intermediate of the methylation reaction. N167 might contribute to the positioning of the residues from the motif VI, but does not play a direct role in catalysis.

Human DNMT1 transition state structure

Human DNA methyltransferase 1 (DNMT1) maintains the epigenetic state of DNA by replicating CpG methylation signatures from parent to daughter strands, producing heritable methylation patterns through cell divisions. The proposed catalytic mechanism of DNMT1 involves nucleophilic attack of Cys(1226) to cytosine (Cyt) C6, methyl transfer from S-adenosyl-l-methionine (SAM) to Cyt C5, and proton abstraction from C5 to form methylated CpG in DNA. Here, we report the subangstrom geometric and electrostatic structure of the major transition state (TS) of the reaction catalyzed by human DNMT1. Experimental kinetic isotope effects were used to guide quantum mechanical calculations to solve the TS structure. Methyl transfer occurs after Cys(1226) attack to Cyt C6, and the methyl transfer step is chemically rate-limiting for DNMT1. Electrostatic potential maps were compared for the TS and ground states, providing the electronic basis for interactions between the protein and reactants at the TS. Understanding the TS of DNMT1 demonstrates the possibility of using similar analysis to gain subangstrom geometric insight into the complex reactions of epigenetic modifications.

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.

DNA cytosine methylation: structural and thermodynamic characterization of the epigenetic marking mechanism

DNA cytosine methyltransferases regulate the expression of the genome through the precise epigenetic marking of certain cytosines with a methyl group, and aberrant methylation is a hallmark of human diseases including cancer. Targeting these enzymes for drug design is currently a high priority. We have utilized ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations to investigate extensively the reaction mechanism of the representative DNA methyltransferase HhaI (M.HhaI) from prokaryotes, whose overall mechanism is shared with the mammalian enzymes. We obtain for the first time full free energy profiles for the complete reaction, together with reaction dynamics in atomistic detail. Our results show an energetically preferred mechanism in which nucleophilic attack of cytosine C5 on the S-adenosyl-L-methionine (AdoMet) methyl group is concerted with formation of the Michael adduct between a conserved Cys in the active site with cytosine C6. Spontaneous and reversible proton transfer between a conserved Glu in the active site and cytosine N3 at the transition state was observed in our simulations, revealing the chemical participation of this Glu residue in the catalytic mechanism. Subsequently, the β-elimination of the C5 proton utilizes as base an OH(-) derived from a conserved crystal water that is part of a proton wire water channel, and this syn β-elimination reaction is the rate-limiting step. Design of novel cytosine methylation inhibitors would be advanced by our structural and thermodynamic characterization of the reaction mechanism.

Structural origins of DNA target selection and nucleobase extrusion by a DNA cytosine methyltransferase

BACKGROUND

How DNA 5-cytosine methyltransferases (DCMTases) select their substrate nucleobase for extrusion from DNA duplex is poorly understood.

RESULTS

The crystal structure of a pre-extrusion M.HaeIII DCMTase-substrate DNA complex is reported here.

CONCLUSION

M.HaeIII selects its substrate cytosine for extrusion by selectively interfering with its stacking and hydrogen bonding interactions within the DNA duplex.

SIGNIFICANCE

This is the first structural elucidation of the target cytosine selection by a DCMTase. Epigenetic methylation of cytosine residues in DNA is an essential element of genome maintenance and function in organisms ranging from bacteria to humans. DNA 5-cytosine methyltransferase enzymes (DCMTases) catalyze cytosine methylation via reaction intermediates in which the DNA is drastically remodeled, with the target cytosine residue extruded from the DNA helix and plunged into the active site pocket of the enzyme. We have determined a crystal structure of M.HaeIII DCMTase in complex with its DNA substrate at a previously unobserved state, prior to extrusion of the target cytosine and frameshifting of the DNA recognition sequence. The structure reveals that M.HaeIII selects the target cytosine and destabilizes its base-pairing through a precise, focused, and coordinated assault on the duplex DNA, which isolates the target cytosine from its nearest neighbors and thereby facilitates its extrusion from DNA.

Genome stability of Lyme disease spirochetes: comparative genomics of Borrelia burgdorferi plasmids

Lyme disease is the most common tick-borne human illness in North America. In order to understand the molecular pathogenesis, natural diversity, population structure and epizootic spread of the North American Lyme agent, Borrelia burgdorferi sensu stricto, a much better understanding of the natural diversity of its genome will be required. Towards this end we present a comparative analysis of the nucleotide sequences of the numerous plasmids of B. burgdorferi isolates B31, N40, JD1 and 297. These strains were chosen because they include the three most commonly studied laboratory strains, and because they represent different major genetic lineages and so are informative regarding the genetic diversity and evolution of this organism. A unique feature of Borrelia genomes is that they carry a large number of linear and circular plasmids, and this work shows that strains N40, JD1, 297 and B31 carry related but non-identical sets of 16, 20, 19 and 21 plasmids, respectively, that comprise 33–40% of their genomes. We deduce that there are at least 28 plasmid compatibility types among the four strains. The B. burgdorferi ∼900 Kbp linear chromosomes are evolutionarily exceptionally stable, except for a short ≤20 Kbp plasmid-like section at the right end. A few of the plasmids, including the linear lp54 and circular cp26, are also very stable. We show here that the other plasmids, especially the linear ones, are considerably more variable. Nearly all of the linear plasmids have undergone one or more substantial inter-plasmid rearrangements since their last common ancestor. In spite of these rearrangements and differences in plasmid contents, the overall gene complement of the different isolates has remained relatively constant.

Characterization of DNA-binding activity in the N-terminal domain of the DNA methyltransferase Dnmt3a

The Dnmt3a gene, which encodes de novo-type DNA methyltransferase, encodes two isoforms, full-length Dnmt3a and Dnmt3a2, which lacks the N-terminal 219 amino acid residues. We found that Dnmt3a showed higher DNA-binding and DNA-methylation activities than Dnmt3a2. The N-terminal sequence from residues 1 to 211 was able to bind to DNA, but could not distinguish methylated and unmethylated CpG. Its binding to DNA was inhibited by a major groove binder. Four basic amino acid residues, Lys51, Lys53, Arg177 and Arg179, in the N-terminal region were crucial for the DNA-binding activity. The ectopically expressed N-terminal sequence (residues 1-211) was localized in nuclei, whereas that harbouring mutations at the four basic amino acid residues was also detected in the cytoplasm. The DNA-methylation activity of Dnmt3a with the mutations was suppressed under physiological salt conditions, which is similar that of Dnmt3a2. In addition, ectopically expressed Dnmt3a with mutations, as well as Dnmt3a2, could not be retained efficiently in nuclei on salt extraction. We conclude that the DNA-binding activity of the N-terminal domain contributes to the DNA-methyltransferase activity via anchoring of the whole molecule to DNA under physiological salt conditions.

Introduction--Epiphanies in epigenetics

The combinatorial pattern of DNA and histone modifications and their associated histone variants constitute an epigenetic code that shapes gene expression patterns by increasing or decreasing the transcriptional potential of genomic domains. The epigenetic coding status, at any given chromosomal location, is subject to modulation by noncoding RNAs and remodeling complexes. DNA methylation is associated with histone modifications, particularly the absence of histone H3 lysine 4 methylation (H3K4me0) and the presence of histone H3 lysine 9 methylation (H3K9m). We briefly discuss four protein domains (ADD, CXXC, MBD, and SRA), and the functional implications of their architecture in linking histone methylation to that of DNA in mammalian cells. We also consider the domain structure of the DNA methyltransferase DNMT1, its accessory protein UHRF1, and their associated proteins. Finally, we discuss a mechanism by which methylation of DNA and of histones may be coordinately maintained during mitotic cell division, allowing for the transmission of parental methylation patterns to newly replicated chromatin.

APP-dependent up-regulation of Ptch1 underlies proliferation impairment of neural precursors in Down syndrome

Mental retardation in Down syndrome (DS) appears to be related to severe neurogenesis impairment during critical phases of brain development. Recent lines of evidence in the cerebellum of a mouse model for DS (the Ts65Dn mouse) have shown a defective responsiveness to Sonic Hedgehog (Shh), a potent mitogen that controls cell division during brain development, suggesting involvement of the Shh pathway in the neurogenesis defects of DS. Based on these premises, we sought to identify the molecular mechanisms underlying derangement of the Shh pathway in neural precursor cells (NPCs) from Ts65Dn mice. By using an in vitro model of NPCs obtained from the subventricular zone and hippocampus, we found that trisomic NPCs had an increased expression of the Shh receptor Patched1 (Ptch1), a membrane protein that suppresses the action of a second receptor, Smoothened (Smo), thereby maintaining the pathway in a repressed state. Partial silencing of Ptch1 expression in trisomic NPCs restored cell proliferation, indicating that proliferation impairment was due to Ptch1 overexpression. The overexpression of Ptch1 in trisomic NPCs resulted from increased levels of AICD [a transcription-promoting fragment of amyloid precursor protein (APP)] and increased AICD binding to the Ptch1 promoter. Our data provide novel evidence that Ptch1 overexpression underlies derangement of the Shh pathway in trisomic NPCs with consequent proliferation impairment. The demonstration that Ptch1 overexpression in trisomic NPCs is due to an APP fragment provides a link between this trisomic gene and the defective neuronal production that characterizes the DS brain.

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.

Mechanism of DNA methylation: the double role of DNA as a substrate and as a cofactor

Methylation of cytosine residues in the DNA is one of the most important epigenetic marks central to the control of differential expression of genes. We perform quantum mechanical calculations to investigate the catalytic mechanism of the bacterial HhaI DNA methyltransferase. We find that the enzyme nucleophile, Cys81, can attack C6 of cytosine only after it is deprotonated by the DNA phosphate group, a reaction facilitated by a bridging water molecule. This finding, which indicates that the DNA acts as both the substrate and the cofactor, can explain the total loss of activity observed in an analogous enzyme, thymidylate synthase, when the phosphate group of the substrate was removed. Furthermore, our results displaying the inability of the phosphate group to deprotonate the side chain of serine is in agreement with the total, or the large extent of, inactivity observed for the C81S mutant. In contrast to results from previous calculations, we find that the active site conserved residues, Glu119, Arg163, and Arg165, are crucial for catalysis. In addition, the enzyme-DNA adduct formation and the methyl transfer from the cofactor S-adenosyl-L-methionine are not concerted but proceed via stepwise mechanism. In many of the different steps of this methylation reaction, the transfer of a proton is found to be necessary. To render these processes possible, we find that several water molecules, found in the crystal structure, play an important role, acting as a bridge between the donating and accepting proton groups.

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

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

Homology modeling and molecular dynamics simulations of HgiDII methyltransferase in complex with DNA and S-adenosyl-methionine: catalytic mechanism and interactions with DNA

M.HgiDII is a methyltransferase (MTase) from Herpetosiphon giganteus that recognizes the sequence GTCGAC. This enzyme belongs to a group of MTases that share a high degree of amino acid similarity, albeit none of them has been thoroughly characterized. To study the catalytic mechanism of M.HgiDII and its interactions with DNA, we performed molecular dynamics simulations with a homology model of M.HgiDII complexed with DNA and S-adenosyl-methionine. Our results indicate that M.HgiDII may not rely only on Glu119 to activate the cytosine ring, which is an early step in the catalysis of cytosine methylation; apparently, Arg160 and Arg162 may also participate in the activation by interacting with cytosine O2. Another residue from the catalytic site, Val118, also played a relevant role in the catalysis of M.HgiDII. Val118 interacted with the target cytosine and kept water molecules from accessing the region of the catalytic pocket where Cys79 interacts with cytosine, thus preventing water-mediated disruption of interactions in the catalytic site. Specific recognition of DNA was mediated mainly by amino acids of the target recognition domain, although some amino acids (loop 80-88) of the catalytic domain may also contribute to DNA recognition. These interactions involved direct contacts between M.HgiDII and DNA, as well as indirect contacts through water bridges. Additionally, analysis of sequence alignments with closely related MTases helped us to identify a motif in the TRD of M.HgiDII that may be relevant to specific DNA recognition.

An investigation of the catalytic mechanism of S-adenosylmethionine synthetase by QM/MM calculations

Catalysis by S-adenosylmethionine synthetase has been investigated by quantum mechanical/molecular mechanical calculations, exploiting structures of the active crystalline enzyme. The transition state energy of +19.1 kcal/mol computed for a nucleophilic attack of the methionyl sulfur on carbon-5' of the nucleotide was indistinguishable from the experimental (solution) value when the QM residues were an uncharged histidine that hydrogen bonds to the leaving oxygen-5' and an aspartate that chelates a Mg2+ ion, and was similar (+18.8 kcal/mol) when the QM region also included the active site arginine and lysines. The computed energy difference between reactant and product was also consistent with their equimolar abundance in co-crystals. The calculated geometrical changes support catalysis of a S(N)2 reaction through hydrogen bonding of the liberated oxygen-5' to the histidine, charge neutralization by the two Mg2+ ions, and stabilization of the product sulfonium cation through a close, non-bonded, contact between the sulfur and the ribose oxygen-4'.

QM/MM methods for biomolecular systems

Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.

Crystal structure of the Escherichia coli 23S rRNA:m5C methyltransferase RlmI (YccW) reveals evolutionary links between RNA modification enzymes

Methylation is the most common RNA modification in the three domains of life. Transfer of the methyl group from S-adenosyl-l-methionine (AdoMet) to specific atoms of RNA nucleotides is catalyzed by methyltransferase (MTase) enzymes. The rRNA MTase RlmI (rRNA large subunit methyltransferase gene I; previously known as YccW) specifically modifies Escherichia coli 23S rRNA at nucleotide C1962 to form 5-methylcytosine. Here, we report the crystal structure of RlmI refined at 2 A to a final R-factor of 0.194 (R(free)=0.242). The RlmI molecule comprises three domains: the N-terminal PUA domain; the central domain, which resembles a domain previously found in RNA:5-methyluridine MTases; and the C-terminal catalytic domain, which contains the AdoMet-binding site. The central and C-terminal domains are linked by a beta-hairpin structure that has previously been observed in several MTases acting on nucleic acids or proteins. Based on bioinformatics analyses, we propose a model for the RlmI-AdoMet-RNA complex. Comparative structural analyses of RlmI and its homologs provide insight into the potential function of several structures that have been solved by structural genomics groups and furthermore indicate that the evolutionary paths of RNA and DNA 5-methyluridine and 5-methylcytosine MTases have been closely intertwined.

Mammalian DNA methyltransferases: a structural perspective

The methylation of mammalian DNA, primarily at CpG dinucleotides, has long been recognized to play a major role in controlling gene expression, among other functions. Given their importance, it is surprising how many basic questions remain to be answered about the proteins responsible for this methylation and for coordination with the parallel chromatin-marking system that operates at the level of histone modification. This article reviews recent studies on, and discusses the resulting biochemical and structural insights into, the DNA nucleotide methyltransferase (Dnmt) proteins 1, 3a, 3a2, 3b, and 3L.

AdoMet-dependent Methyl-transfer: Glu119 Is Essential for DNA C5-Cytosine Methyltransferase M.HhaI

The role of Glu119 in S-adenosyl-L-methionine-dependent DNA methyltransferase M.HhaI-catalyzed DNA methylation was studied. Glu119 belongs to the highly conserved Glu/Asn/Val motif found in all DNA C5-cytosine methyltransferases, and its importance for M.HhaI function remains untested. We show that formation of the covalent intermediate between Cys81 and the target cytosine requires Glu119, since conversion to Ala, Asp or Gln lowers the rate of methyl transfer 10(2)-10(6) fold. Further, unlike the wild-type M.HhaI, these mutants are not trapped by the substrate in which the target cytosine is replaced with the mechanism-based inhibitor 5-fluorocytosine. The DNA binding affinity for the Glu119Asp mutant is decreased 10(3)-fold. Thus, the ability of the enzyme to stabilize the extrahelical cytosine is coupled directly to tight DNA binding. The structures of the ternary protein/DNA/AdoHcy complexes for both the Glu119Ala and Glu119Gln mutants (2.70 A and 2.75 A, respectively) show that the flipped base is positioned nearly identically with that observed in the wild-type M.HhaI complex. A single water molecule in the Glu119Ala structure between Ala119 and the extrahelical cytosine N3 is lacking in the Glu119Gln and wild-type M.HhaI structures, and most likely accounts for this mutant's partial activity. Glu119 has essential roles in activating the target cytosine for nucleophilic attack and contributes to tight DNA binding.

Observing an Induced-fit Mechanism during Sequence-specific DNA Methylation

The characterization of conformational changes that drive induced-fit mechanisms and their quantitative importance to enzyme specificity are essential for a full understanding of enzyme function. Here, we report on M.HhaI, a sequence-specific DNA cytosine C(5) methyltransferase that reorganizes a flexible loop (residues 80-100) upon binding cognate DNA as part of an induced-fit mechanism. To directly observe this approximately 26A conformational rearrangement and provide a basis for understanding its importance to specificity, we replaced loop residues Lys-91 and Glu-94 with tryptophans. The double mutants W41F/K91W and W41F/E94W are relatively unperturbed in kinetic and thermodynamic properties. W41F/E94W shows DNA sequence-dependent changes in fluorescence: significant changes in equilibrium and transient state fluorescence that occur when the enzyme binds cognate DNA are absent with nonspecific DNA. These real-time, solution-based results provide direct evidence that binding to cognate DNA induces loop reorganization into the closed conformer, resulting in the correct assembly of the active site. We propose that M.HhaI scans nonspecific DNA in the loop-open conformer and rearranges to the closed form once the cognate site is recognized. The fluorescence data exclude mechanisms in which loop motion precedes base flipping, and we show loop rearrangements are directly coupled to base flipping, because the sequential removal of single hydrogen bonds within the target guanosine:cytosine base pair results in corresponding changes in loop motion.

Reaction mechanism of guanidinoacetate methyltransferase, concerted or step-wise

We describe a quantum mechanics/molecular mechanics investigation of the guanidinoacetate methyltransferase catalyzed reaction, which shows that proton transfer from guanidinoacetate (GAA) to Asp-134 and methyl transfer from S-adenosyl-L-methionine (AdoMet) to GAA are concerted. By self-consistent-charge density functional tight binding/molecular mechanics, the bond lengths in the concerted mechanism's transition state are 1.26 A for both the OD1 (Asp-134)-H(E) (GAA) and H(E) (GAA)-N(E) (GAA) bonds, and 2.47 and 2.03 A for the S8 (AdoMet)-C9 (AdoMet) and C9 (AdoMet)-N(E) (GAA) bonds, respectively. The potential-energy barrier (DeltaE++) determined by single-point B3LYP/6-31+G*//MM is 18.9 kcal/mol. The contributions of the entropy (-TDeltaS++) and zero-point energy corrections Delta(ZPE)++ by normal mode analysis are 2.3 kcal/mol and -1.7 kcal/mol, respectively. Thus, the activation enthalpy of this concerted mechanism is predicted to be DeltaH++ = DeltaE++ plus Delta(ZPE)++ = 17.2 kcal/mol. The calculated free-energy barrier for the concerted mechanism is DeltaG++ = 19.5 kcal/mol, which is in excellent agreement with the value of 19.0 kcal/mol calculated from the experimental rate constant (3.8 +/- 0.2.min(-1)).