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.

     

Promiscuity in alkaline phosphatase superfamily. Unraveling evolution through molecular simulations

We here present a theoretical study of the alkaline hydrolysis of a phosphodiester (methyl p-nitrophenyl phosphate or MpNPP) in the active site of Escherichia coli alkaline phosphatase (AP), a monoesterase that also presents promiscuous activity as a diesterase. The analysis of our simulations, carried out by means of molecular dynamics (MD) simulations with hybrid quantum mechanics/molecular mechanics (QM/MM) potentials, shows that the reaction takes place through a D(N)A(N) or dissociative mechanism, the same mechanism employed by AP in the hydrolysis of monoesters. The promiscuous activity observed in this superfamily can be then explained on the basis of a conserved reaction mechanism. According to our simulations the specialization in the hydrolysis of phosphomonoesters or phosphodiesters, developed in different members of the superfamily, is a consequence of the interactions established between the protein and the oxygen atoms of the phosphate group and, in particular, with the oxygen atom that bears the additional alkyl group when the substrate is a diester. A water molecule, belonging to the coordination shell of the Mg(2+) ion, and residue Lys328 seem to play decisive roles stabilizing a phosphomonoester substrate, but the latter contributes to increase the energy barrier for the hydrolysis of phosphodiesters. Then, mutations affecting the nature or positioning of Lys328 lead to an increased diesterase activity in AP. Finally, the capacity of this enzymatic family to catalyze the reaction of phosphoesters having different leaving groups, or substrate promiscuity, is explained by the ability of the enzyme to stabilize different charge distributions in the leaving group using different interactions involving either one of the zinc centers or residues placed on the outer side of the catalytic site.     

Oxidation reaction by xanthine oxidase: theoretical study of reaction mechanism

The oxidation process by molybdenum-containing enzyme, xanthine oxidase, is theoretically studied with a model complex representing the reaction center and a typical benchmark substrate, formamide. Comparisons were systematically made among reaction mechanisms proposed previously. In the concerted and stepwise mechanisms that were theoretically discussed previously, the oxidation reaction takes place with a moderate activation barrier. However, the product is less stable than the reactant complex, which indicates that these mechanisms are unlikely. Moreover, the product of the concerted mechanism is not consistent with the isotope experimental result. In addition to those mechanisms, another mechanism initiated by the deprotonation of the active site was newly investigated here. In the transition state of this reaction, the carbon atom of formamide interacts with the oxo ligand of the Mo center and the hydrogen atom is moving from the carbon atom to the thioxo ligand. This reaction takes place with a moderate activation barrier and considerably large exothermicity. Furthermore, the product by this mechanism is consistent with the isotope experimental result. Also, our computations clearly show that the deprotonation of the active site occurs with considerable exothermicity in the presence of glutamic acid and substrate. The intermediate of the stepwise mechanism could not be optimized in the case of the deprotonated active site. From all these results, it should be concluded that the one-step mechanism with the deprotonated active site is the most plausible.     

Ion chemistry of 1H-1,2,3-triazole

A combination of experimental methods, photoelectron-imaging spectroscopy, flowing afterglow-photoelectron spectroscopy and the flowing afterglow-selected ion flow tube technique, and electronic structure calculations at the B3LYP/6-311++G(d,p) level of density functional theory (DFT) have been employed to study the mechanism of the reaction of the hydroxide ion (HO-) with 1H-1,2,3-triazole. Four different product ion species have been identified experimentally, and the DFT calculations suggest that deprotonation by HO- at all sites of the triazole takes place to yield these products. Deprotonation of 1H-1,2,3-triazole at the N1-H site gives the major product ion, the 1,2,3-triazolide ion. The 335 nm photoelectron-imaging spectrum of the ion has been measured. The electron affinity (EA) of the 1,2,3-triazolyl radical has been determined to be 3.447 +/- 0.004 eV. This EA and the gas-phase acidity of 2H-1,2,3-triazole are combined in a negative ion thermochemical cycle to determine the N-H bond dissociation energy of 2H-1,2,3-triazole to be 112.2 +/- 0.6 kcal mol-1. The 363.8 nm photoelectron spectroscopic measurements have identified the other three product ions. Deprotonation of 1H-1,2,3-triazole at the C5 position initiates fragmentation of the ring structure to yield a minor product, the ketenimine anion. Another minor product, the iminodiazomethyl anion, is generated by deprotonation of 1H-1,2,3-triazole at the C4 position, followed by N1-N2 bond fission. Formation of the other minor product, the 2H-1,2,3-triazol-4-ide ion, can be rationalized by initial deprotonation of 1H-1,2,3-triazole at the N1-H site and subsequent proton exchanges within the ion-molecule complex. The EA of the 2H-1,2,3-triazol-4-yl radical is 1.865 +/- 0.004 eV.     

Acyl migration from N to C in aziridine-2-carboxylic esters

Acyl group migration from N to C in aziridine-2-carboxylates takes place in deprotonation reactions and, as a result, aziridine-2,2-dicarboxylates are formed. Mechanistic studies proved that the observed migration is an intramolecular reaction.     

Pyridine ring opening at room temperature at a rhenium tricarbonyl bipyridine complex

Pyridine ring opening occurs in the reaction of [Re(CO)3(MeIm)(bipy)]OTf with KN(SiMe3)2 followed by double methylation with methyl trifluoromethanesulfonate. Analogues of the neutral product of the initial deprotonation and of the product of the first methylation were isolated by using mesitylimidazole (MesIm) in place of methylimidazole (MeIm) and/or 1,10-phenanthroline (phen) instead of 2,2'-bipyridine (bipy).     

Ionization mechanism of oligonucleotides in matrix-assisted laser desorption/ionization time-of-flight mass spectrometry

The ionization of nucleosides in matrix-assisted laser desorption/ionization time-of-flight mass spectrometry was systematically investigated using adenine (A), thymine (T), guanine (G) and cytosine (C) with several common matrices. Experimental results of the protonation and deprotonation of the bases of A, T, G and C in the matrices 2,5-dihydroxybenzoic acid (2,5-DHB), alpha-cyano-4-hydroxycinnamic acid (alpha-CHCA) and 3-hydroxypicolinic acid (3-HPA) provide an insight into the ionization mechanism of oligonucleotides in MALDI. It was found that the low ion signal from DNA in poly-G in MALDI as reported in earlier work could be attributed to the fact that the base of G is difficult to ionize. Our results suggest that the ionization of DNA in MALDI is dominated by the protonation and deprotonation of bases and it is basically independent of the backbone of DNA. Both the protonation and deprotonation are strongly structure dependent. The protonation is dominated by pre-protonation before laser ablation, while the deprotonation is controlled by the thermal reaction.     

Mechanism of epoxide hydrolysis in microsolvated nucleotide bases adenine, guanine and cytosine: a DFT study

Six water molecules have been used for microsolvation to outline a hydrogen bonded network around complexes of ethylene epoxide with nucleotide bases adenine (EAw), guanine (EGw) and cytosine (ECw). These models have been developed with the MPWB1K-PCM/6-311++G(3df,2p)//MPWB1K/6-31+G(d,p) level of DFT method and calculated S(N)2 type ring opening of the epoxide due to amino group of the nucleotide bases, viz. the N6 position of adenine, N2 position of guanine and N4 position of cytosine. Activation energy (E(act)) for the ring opening was found to be 28.06, 28.64, and 28.37 kcal mol(-1) respectively for EAw, EGw and ECw. If water molecules were not used, the reactions occurred at considerably high value of E(act), viz. 53.51 kcal mol(-1) for EA, 55.76 kcal mol(-1) for EG and 56.93 kcal mol(-1) for EC. The ring opening led to accumulation of negative charge on the developing alkoxide moiety and the water molecules around the charge localized regions showed strong hydrogen bond interactions to provide stability to the intermediate systems EAw-1, EGw-1 and ECw-1. This led to an easy migration of a proton from an activated water molecule to the alkoxide moiety to generate a hydroxide. Almost simultaneously, a proton transfer chain reaction occurred through the hydrogen bonded network of water molecules and resulted in the rupture of one of the N-H bonds of the quaternized amino group. The highest value of E(act) for the proton transfer step of the reaction was 2.17 kcal mol(-1) for EAw, 2.93 kcal mol(-1) for EGw and 0.02 kcal mol(-1) for ECw. Further, the overall reaction was exothermic by 17.99, 22.49 and 13.18 kcal mol(-1) for EAw, EGw and ECw, respectively, suggesting that the reaction is irreversible. Based on geometric features of the epoxide-nucleotide base complexes and the energetics, the highest reactivity is assigned for adenine followed by cytosine and guanine. Epoxide-mediated damage of DNA is reported in the literature and the present results suggest that hydrated DNA bases become highly S(N)2 active on epoxide systems and the occurrence of such reactions can inflict permanent damage to the DNA.     

Electron transfer in amino acid.nucleic acid base complexes: EPR, ENDOR, and DFT study of X-irradiated N-formylglycine.cytosine complex crystals

Single crystals of the 1:1 complex of the nucleic acid base cytosine and the dipeptide N-formylglycine (C.NFG) have been irradiated at 10 and 273 K to doses of about 70 kGy and studied at temperatures between 10 and 293 K using 24 GHz (K-band) and 9.5 GHz (X-band) electron paramagnetic resonance (EPR), electron nuclear double resonance (ENDOR), and ENDOR-induced EPR (EIE) spectroscopy. In this complex, the cytosine base is hydrogen bonded at positions N3 and N4 to the carboxylic group of the dipeptide, and the N3 position of cytosine has become protonated by the carboxylic group. At 10 K, two major radicals were characterized and identified. One of these (R1) is ascribed to the decarboxylated N-formylglycine one-electron oxidized species. The other (R2) is the N3-protonated cytosine one-electron reduced species. A third minority species (R3) appears to be a different conformation or protonation state of the one-electron reduced cytosine radical. Upon warming, the R2 and R3 radicals decay at about 100 K, and at 295 K, the only cytosine-centered radicals present are the C5 and C6 H-addition radicals (R5, R6). The R1 radical decays at about 150 K, and a glycine backbone radical (R4) grows in slowly. Thus, in the complex, a complete separation of initial oxidation and reduction events occurs, with oxidation localized at the dipeptide moiety, whereas reduction occurs at the nucleic acid base moiety. DFT calculations indicate that this separation is driven by large differences in electron affinities and ionization potentials between the two constituents of the complex. Once the initial oxidation and reduction products are trapped, no further electron transfer between the two constituents of the complex takes place.     

Assessing the Brønsted Basicity of Diaminoboryl Anions: Reactivity toward Methylated Benzenes and Dihydrogen

Treatment of toluene or p‐xylene with diaminoboryllithium results in consecutive reactions, involving boryl‐anion‐mediated deprotonation at the benzylic position followed by nucleophilic substitution at the boron center, producing benzylborane species and LiH. Diaminoboryllithium also cleaves H₂ heterolytically affording diaminohydroborane and LiH, while the reaction of lithium diaminoboryl(bromo)cuprate with H₂ takes place accompanied by reduction of Cu^I to give diaminohydroborane, LiH, and Cu⁰.     

Glutathione peroxidase-like antioxidant activity of diaryl diselenides: a mechanistic study.

The synthesis, structure, and thiol peroxidase-like antioxidant activities of several diaryl diselenides having intramolecularly coordinating amino groups are described. The diselenides derived from enantiomerically pure R-(+)- and S-(-)-N,N-dimethyl(1-ferrocenylethyl)amine show excellent peroxidase activity. To investigate the mechanistic role of various organoselenium intermediates, a detailed in situ characterization of the intermediates has been carried out by (77)Se NMR spectroscopy. While most of the diselenides exert their peroxidase activity via selenol, selenenic acid, and selenenyl sulfide intermediates, the differences in the relative activities of the diselenides are due to the varying degree of intramolecular Se.N interaction. The diselenides having strong Se.N interactions are found to be inactive due to the ability of their selenenyl sulfide derivatives to enhance the reverse GPx cycle (RSeSR + H(2)O(2) = RSeOH). In these cases, the nucleophilic attack of thiol takes place preferentially at selenium rather than sulfur and this reduces the formation of selenol by terminating the forward reaction. On the other hand, the diselenides having weak Se.N interactions are found to be more active due to the fast reaction of the selenenyl sulfide derivatives with thiol to produce diphenyl disulfide and the expected selenol (RSeSR + PhSH = PhSSPh + RSeH). The unsubstituted diaryl diselenides are found to be less active due to the slow reactions of these diselenides with thiol and hydrogen peroxide and also due to the instability of the intermediates. The catalytic cycles of 18 and 19 strongly resemble the mechanism by which the natural enzyme, glutathione peroxidase, catalyzes the reduction of hydroperoxides.     

Mechanistic study on DNA mutation of the cytosine methylation reaction at C5 position

DNA methylation is a crucial epigenetic mark connected to conventionally changing the DNA bases, typically by adding methyl groups into DNA bases. Methylation of cytosine at the C5 position (5-methylcytosine) occurs mostly in the context of cytosine-phosphate-guanine dinucleotides, the methylation of which has important impacts on gene regulation and expression. However, the mechanistic details of this reaction are still debatable concerning the concertedness of the key reaction steps and the roles played by the base that abstracts the proton in the β-elimination and water molecules at the active site. To gain a deeper insight into the formation of 5-mehtylcytosine, an extensive density functional theory (DFT) study was performed with the B3LYP functional in conjunction with different basis sets. Our study has clearly established the mechanistic details of this methylation approach, based on which the roles of conserved active site residues, such as glutamic acid and waters, are well understood. Our results show that the reaction of 5-methylcytosine follows a concerted mechanism in which water molecules are critically involved. Moreover, arginine and alanine give more significant catalytic effects than glutamic acid on the 5-methylcytosine process. Considering the effect of Alanine, Arginine, and one water bridging molecule, the activation energy is 31 kJ mol^?1 calculated at B3LYP/6-31G(d) level of theory.     

N4‐Methylation of Cytosine Drastically Favors the Formation of (6‐4) Photoproducts in a TCG Context

Methylation of cytosine is a common biological process both in prokaryotic and eukaryotic cells. In addition to 5‐methylcytosine (5mC), some bacterial species contain in their genome N⁴‐methylcytosine (N4mC). Methylation at C5 has been shown to enhance the formation of pyrimidine dimeric photoproducts but nothing is known of the effect of N4 methylation on UV‐induced DNA damage. In the present work, we compared the yield and the nature of bipyrimidine photoproducts induced in a series of trinucleotides exhibiting a TXG sequence where X is either T, C, 5mC or N4mC. HPLC associated to tandem mass spectrometry was used to quantify cyclobutane pyrimidine dimers (CPD), (6‐4) photoproducts (64PP) and their Dewar valence isomer. Methylation at position N4 was found to drastically increase the reactivity of C upon exposure to both UVC and UVB and to favor the formation of 64PP. In contrast methylation at C5 increased the yield of CPD at the expense of 64PP. In addition, enhancement of photoreactivity by C5 methylation was much higher in the UVB than in the UVC range. These results show the drastic effect of the methylation site on the photochemistry of cytosine.     

Development and validation of a gas chromatography/mass spectrometry method for the assessment of genomic DNA methylation

A method for the determination of DNA global methylation, taken as the ratio (%) of 5‐methylcytosine (5mCyt) versus the sum of cytosine (Cyt) and 5mCyt, via gas chromatography/mass spectrometry (GC/MS), was developed and validated. DNA (2.5 µg) was hydrolyzed with aqueous formic acid 88%, spiked with cytosine‐2,4‐¹³C₂,¹⁵N₃ and 5‐methyl‐²H₃‐cytosine‐6‐²H₁ as internal standards, and derivatized with N‐methyl‐N‐(tert‐butyldimethylsilyl)trifluoroacetamide and 1% tert‐butyldimethylchlorosilane, in the presence of acetonitrile and pyridine. GC/MS, operating in single ion monitoring mode, separated and specifically detected all nucleobases as tert‐butyldimethylsilyl derivatives, without interferences, with the exception of guanosine. The method was linear throughout the range of clinical interest and had good sensitivity, with a limit of quantification of 3.2 pmol for Cyt and 0.056 pmol for 5mCyt, the latter corresponding to a methylation level of 0.41%. Intra‐ and inter‐day precision and accuracy were below 4.0% for both analytes and methylation. The matrix absolute effect, process efficiency and coefficient of variation ranged from 96.5 to 101.2%. The matrix relative effect was below 1%. The method was applied to the analysis of different human DNAs, including: nonmethylated DNA from PCR (methylation 0.00%), hypermethylated DNA prepared using M.SssI CpG methyltransferase (methylation 18.05%), DNA from peripheral blood leukocytes of healthy subjects (N = 6, median methylation 5.45%), DNA from bone marrow of leukemia patients (N = 5, 3.58%) and DNA from myeloma cell lines (N = 4, 2.74%). Copyright © 2009 John Wiley & Sons, Ltd.     

Enzymatic effects on reactant and transition states. The case of chalcone isomerase

Chalcone isomerase catalyzes the transformation of chalcone to naringerin as a part of flavonoid biosynthetic pathways. The global reaction takes place through a conformational change of the substrate followed by chemical reaction, being thus an excellent example to analyze current theories about enzyme catalysis. We here present a detailed theoretical study of the enzymatic action on the conformational pre-equilibria and on the chemical steps for two different substrates of this enzyme. Free-energy profiles are obtained in terms of potentials of mean force using hybrid quantum mechanics/molecular mechanics potentials. The role of the enzyme becomes clear when compared to the counterpart equilibria and reactions in aqueous solution. The enzyme does not only favor the chemical reaction lowering the corresponding activation free energy but also displaces the conformational equilibria of the substrates toward the reactive form. These results, which can be rationalized in terms of the electrostatic interactions established in the active site between the substrate and the environment, agree with a more general picture of enzyme catalysis. According to this, an active site designed to accommodate the transition state of the reaction would also have consequences on the reactant state, stabilizing those forms which are geometrically and/or electronically closer to the transition structure.     

Enzymatic deamination of the epigenetic base N-6-methyladenine

Two enzymes of unknown function from the amidohydrolase superfamily were discovered to catalyze the deamination of N-6-methyladenine to hypoxanthine and methyl amine. The methylation of adenine in bacterial DNA is a common modification for the protection of host DNA against restriction endonucleases. The enzyme from Bacillus halodurans, Bh0637, catalyzes the deamination of N-6-methyladenine with a k(cat) of 185 s(-1) and a k(cat)/K(m) of 2.5 × 10(6) M(-1) s(-1). Bh0637 catalyzes the deamination of N-6-methyladenine 2 orders of magnitude faster than adenine. A comparative model of Bh0637 was computed using the three-dimensional structure of Atu4426 (PDB code: 3NQB) as a structural template and computational docking was used to rationalize the preferential utilization of N-6-methyladenine over adenine. This is the first identification of an N-6-methyladenine deaminase (6-MAD).     

Methyl transfer in glycine N-methyltransferase. A theoretical study

Density functional theory calculations using the hybrid functional B3LYP have been performed to study the methyl transfer step in glycine N-methyltransferase (GNMT). This enzyme catalyzes the S-adenosyl-L-methionine (SAM)-dependent methylation of glycine to form sarcosine. The starting point for the calculations is the recent X-ray crystal structure of GNMT complexed with SAM and acetate. Several quantum chemical models with different sizes, employing up to 98 atoms, were used. The calculations demonstrate that the suggested mechanism, where the methyl group is transferred in a single S(N)2 step, is thermodynamically plausible. By adding or eliminating various groups at the active site, it was furthermore demonstrated that hydrogen bonds to the amino group of the glycine substrate lower the reaction barrier, while hydrogen bonds to the carboxylate group raise the barrier.     

Activation barriers for DNA alkylation by carcinogenic methane diazonium ions

Methylation reactions of the DNA bases with the methane diazonium ion, which is the reactive intermediate formed from several carcinogenic methylating agents, were examined. The SN2 transition states of the methylation reactions at N7, N3, and O6 of guanine; N7, N3, and N1 of adenine; N3 and O2 of cytosine; and O2 and O4 of thymine were calculated using the B3LYP density functional method. Solvation effects were examined using the conductor-like polarizable continuum method and the combined discrete/SCRF method. The transition states for reactions at guanine N3, adenine N7, and adenine N1 are influenced by steric interactions between the methane diazonium ion and exocyclic amino groups. Both in the gas phase and in aqueous solution, the methylation reactions at N atoms have transition states that are looser, and generally occur earlier along the reaction pathways than reactions at O atoms. The forming bonds in the transition states in water are 0.03 to 0.13 A shorter than those observed in the gas phase, and the activation energies are 13 to 35 kcal/mol higher. The combined discrete/SCRF solvation energy calculations using base-water complexes with three water molecules yield base solvation energies that are larger than those obtained from the CPCM continuum method, especially for cytosine. Reactivities calculated using barriers obtained with the discrete/SCRF method are consistent with the experimentally observed high reactivity at N7 of guanine.     

3-deoxy-3,3-difluoro-D-arabinofuranose: first stereoselective synthesis and application in preparation of gem-difluorinated sugar nucleosides

The design and synthesis of gem-difluorinated sugar nucleosides were described. The key intermediate, 3-deoxy-3,3-difluoro-d-arabinofuranose 9, was first stereoselectively prepared from the chiral gem-difluorohomoallyl alcohol 12. The kinetic formation of single anti-14 in the benzylation of 12 could be accomplished by controlling the amount of sodium hydride used. The dihydroxylation of 14 (a mixture of anti and syn isomers) followed by deprotection and oxidation stereoselectively afforded furanose 9 with the arabino configuration at the C2 position. N(1)-(3-Deoxy-3,3-difluoro-beta-D-arabinofuranosyl)cytosine 6 was prepared from 9 by the glycosylation reaction. 4'-Thiofuranose 25 was easily synthesized from 9. The oxidation of 25 followed by the condensation with silylated N(4)-benzoylcytosine (Pummerer reaction) failed to give our desired protected nucleoside l-3'-deoxy-3',3'-difluoro- 4'-thiocytidine 27', but the regioisomer 27 was obtained. The regiochemistry of the Pummerer reaction was determined by the kinetic acidity of the alpha-proton of 4'-thiofuranose 25.