Active site dynamics of the HhaI methyltransferase: insights from computer simulation

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.

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).

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.

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.

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.

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

The mechanism of DNA cytosine-5-methylation catalyzed by the bacterial M.HhaI enzyme has been considered as a stepwise nucleophilic addition of Cys-81-S- to cytosine C6 followed by C5 nucleophilic replacement of the methyl of S-adenosyl-L-methionine to produce 5-methyl-6-Cys-81-S-5,6-dihydrocytosine. In this study, we show that the reaction is concerted from a series of energy calculations by using the quantum mechanical/molecular mechanical hybrid method. Deprotonation of 5-methyl-6-Cys-81-S-5,6-dihydrocytosine and expulsion of Cys-81-S- provides the product DNA 5-methylcytosine. A required base catalyst for this deprotonation is not available as a member of the active site structure. A water channel between the active site and bulk water allows entrance of solvent to the active site. Hydroxide at 10(-7) mole fraction (pH = 7) is shown to be sufficient for the required catalysis. We also show that Glu-119-CO2H can divert the reaction by protonating cytosine N3 when Cys-81-S- attacks cytosine, to form the 6-Cys-81-S-3-hydrocytosine. The reactants and 6-Cys-81-S-3-hydrocytosine product are in rapid equilibrium, and this explains the observed hydrogen exchange of cytosine with solvent.

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.

Stochastic boundary molecular dynamics simulation of l-ribose in the active site of Actinoplanes missouriensis xylose isomerase and its Val135Asn mutant with improved reaction rate

We used molecular dynamics simulations to study how a non-natural substrate, L-ribose, interacts with the active site of Actinoplanes missouriensis xylose isomerase. The simulations showed that L-ribose does not stay liganded in the active site in the same way as D-xylose, in which the oxygens O2 and O4 are liganded to the metal M1. The oxygen O4 of L-ribose moved away from the metal M1 to an upside down position. Furthermore, the distances of the carbons C1 and C2 of L-ribose to the catalytic metal M2 were higher than in the case of D-xylose. These findings explain the extremely low reaction rate of xylose isomerase with L-ribose. The mutation V135N close to the C5-OH of the substrate increased the reaction efficiency 2- to 4-fold with L-ribose. V135N did not affect the reaction with D-xylose and L-arabinose, whereas the reaction with D-glucose was impaired, probably due to a hydrogen bond between Asn-135 and the substrate. When L-ribose was the substrate, Asn-135 formed a hydrogen bond to Glu-181. As a consequence, O4 of L-ribose stayed liganded to the metal M1 in the V135N mutant in molecular dynamics simulations. This explains the decreased K(m) of the V135N mutant with L-ribose.

Residues Distal from the Active Site that Alter Enzyme Function in M.HhaI DNA Cytosine Methyltransferase

Ten M.HhaI residues were replaced with alanine to probe the importance of distal protein elements to substrate/cofactor binding, methyl transfer, and product release. The substitutions, ranging from 6-20 A from the active site were evaluated by thermodynamic analysis, pre-steady and steady-state kinetics, to obtain Kd(AdoMet), Kd(DNA), kcat/Km(DNA), kcat, and kmethyltransfer values. For the wild-type M.HhaI, product release steps dominate catalytic turnover while the 4-fold faster internal microscopic constant kmethyltransfer presents an upper limit. The methyl transfer reaction has DeltaH and DeltaS values of 10.3 kcal/mol and -29.4 cal/(mol K), respectively, consistent with a compressed transition state similar to that observed in the gas phase. Although the ten mutants remained largely unperturbed in methyl transfer, long-range effects influencing substrate/cofactor binding and product release were observed. Positive enhancements were seen in Asp73Ala, which showed a 25-fold improvement in AdoMet affinity and in Val282Ala, which showed a 4-fold improvement in catalytic turnover. Based on an analysis of the positional probability within the C5-cytosine DNA methyltransferase family we propose that certain conserved distal residues may be important in mediating long-range effects.

Statistical coevolution analysis and molecular dynamics: identification of amino acid pairs essential for catalysis

Molecular dynamics (MD) simulations of HhaI DNA methyltransferase and statistical coupling analysis (SCA) data on the DNA cytosine methyltransferase family were combined to identify residues that are coupled by coevolution and motion. The highest ranking correlated pairs from the data matrix product (SCA.MD) are colocalized and form stabilizing interactions; the anticorrelated pairs are separated on average by 30 A and form a clear focal point centered near the active site. We suggest that these distal anti-correlated pairs are involved in mediating active-site compressions that may be important for catalysis. Mutants that disrupt the implicated interactions support the validity of our combined SCA.MD approach.

The coupling of tight DNA binding and base flipping: identification of a conserved structural motif in base flipping enzymes

Val(121) is positioned immediately above the extrahelical cytosine in HhaI DNA C(5)-cytosine methyltransferase, and replacement with alanine dramatically interferes with base flipping and catalysis. DNA binding and k(cat) are decreased 10(5)-fold for the Val(121) --> Ala mutant that has a normal circular dichroism spectrum and AdoMet affinity. The magnitude of this loss of function is comparable with removal of the essential catalytic Cys(81). Surprisingly, DNA binding is completely recovered (increase of 10(5)-fold) with a DNA substrate lacking the target cytosine base (abasic). Thus, interfering with the base flipping transition results in a dramatic loss of binding energy. Our data support an induced fit mechanism in which tight DNA binding is coupled to both base flipping and protein loop rearrangement. The importance of the proximal protein segment (His(127)-Thr(132)) in maintaining this critical interaction between Val(121) and the flipped cytosine was probed with single site alanine substitutions. None of these mutants are significantly altered in secondary structure, AdoMet or DNA affinity, k(methylation), k(inactivation), or k(cat). Although Val(121) plays a critical role in both extrahelical base stabilization and catalysis, its position and mobility are not influenced by individual residues in the adjacent peptide region. Structural comparisons with other DNA methyltransferases and DNA repair enzymes that stabilize extrahelical nucleotides reveal a motif that includes a positively charged or polar side chain and a hydrophobic residue positioned adjacent to the target DNA base and either the 5'- or 3'-phosphate.

Sequence-structure-function studies of tRNA:m5C methyltransferase Trm4p and its relationship to DNA:m5C and RNA:m5U methyltransferases

Three types of methyltransferases (MTases) generate 5-methylpyrimidine in nucleic acids, forming m5U in RNA, m5C in RNA and m5C in DNA. The DNA:m5C MTases have been extensively studied by crystallographic, biophysical, biochemical and computational methods. On the other hand, the sequence-structure-function relationships of RNA:m5C MTases remain obscure, as do the potential evolutionary relationships between the three types of 5-methylpyrimidine-generating enzymes. Sequence analyses and homology modeling of the yeast tRNA:m5C MTase Trm4p (also called Ncl1p) provided a structural and evolutionary platform for identification of catalytic residues and modeling of the architecture of the RNA:m5C MTase active site. The analysis led to the identification of two invariant residues that are important for Trm4p activity in addition to the conserved Cys residues in motif IV and motif VI that were previously found to be critical. The newly identified residues include a Lys residue in motif I and an Asp in motif IV. A conserved Gln found in motif X was found to be dispensable for MTase activity. Locations of essential residues in the model of Trm4p are in very good agreement with the X-ray structure of an RNA:m5C MTase homolog PH1374. Theoretical and experimental analyses revealed that RNA:m5C MTases share a number of features with either RNA:m5U MTases or DNA:m5C MTases, which suggested a tentative phylogenetic model of relationships between these three classes of 5-methylpyrimidine MTases. We infer that RNA:m5C MTases evolved from RNA:m5U MTases by acquiring an additional Cys residue in motif IV, which was adapted to function as the nucleophilic catalyst only later in DNA:m5C MTases, accompanied by loss of the original Cys from motif VI, transfer of a conserved carboxylate from motif IV to motif VI and sequence permutation.

Preorganization and protein dynamics in enzyme catalysis

Recently, an alternative has been offered to the concept of transition state (TS) stabilization as an explanation for rate enhancements in enzyme-catalyzed reactions. Instead, most of the rate increase has been ascribed to preorganization of the enzyme active site to bind substrates in a geometry close to that of the TS, which then transit the activation barrier impelled by motions along the reaction coordinate. The question as to how an enzyme achieves such preorganization and concomitant TS stabilization as well as potential coupled motions along the reaction coordinate leads directly to the role of protein dynamic motion. Dihydrofolate reductase (DHFR) is a paradigm in which the role of dynamics in catalysis continues to be unraveled by a wealth of kinetic, structural, and computational studies. DHFR has flexible loop regions adjacent to the active site whose motions modulate passage through the kinetically preferred pathway. The participation of residues distant from the DHFR active site in enhancing the rate of hydride transfer, however, is unanticipated and may signify the importance of long range protein motions. The general significance of protein dynamics in understanding other biological processes is briefly discussed.

A theoretical examination of the factors controlling the catalytic efficiency of the DNA-(adenine-N6)-methyltransferase from Thermus aquaticus

Ab initio and density functional calculations have been carried out to more fully understand the factors controlling the catalytic activity of the Thermus aquaticus DNA methyltransferase (MTaqI) in the N-methylation at the N(6) of an adenine nucleobase. The noncatalyzed reaction was modeled as a methyl transfer from trimethylsulfonium to the N(6) of adenine. Activation barriers of 32.0 kcal/mol and 24.0 kcal/mol were predicted for the noncatalyzed reaction in the gas phase by MP2/6-31+G(d,p)//HF/6-31+G(d,p) and B3LYP/6-31+G(d,p) calculations, respectively. Calculations performed to evaluate the effect of substrate positioning in the active site of MTaqI on the reaction determine the barrier to be 23.4 kcal/mol and 17.3 kcal/mol for the MP2/6-31+G(d,p)//HF/6-31+G(d,p) and B3LYP/6-31+G(d,p) gas phase calculations, respectively. The effect of hydrogen bonding between the N(6) of adenine and the terminal oxygen of Asn-105 on the activation barrier was also studied. A formamide molecule was modeled into the system to mimic the function of active site residue Asn-105. The activation barrier for this reaction was found to be 21.8 kcal/mol and 15.8 kcal/mol as determined from the MP2/6-31+G(d,p)//HF/6-31+G(d,p) and B3LYP/6-31+G(d,p) calculations, respectively. This result predicts a contribution of less than 2 kcal/mol to the lowering of the activation barrier from amide hydrogen bonding between formamide and N(6) of adenine. Comparison of the reaction coordinates suggest that it is not the hydrogen bonding of the Asn-105 that lends to the catalytic prowess of the enzyme since the organization of the substrates in the active site of the enzyme has a far greater effect on reducing the activation barrier. The results also suggest a stepwise mechanism for the removal of the hydrogen from the N(6) of adenine as opposed to a concerted reaction in which a proton is abstracted simultaneously with the transfer of the methyl group. The hydrogen on the N(6) of the intermediate methyl adenine product is far more acidic than in the reactant complex and may be subsequently abstracted by basic groups in the active site that are too weak to abstract the proton before the full sp(3) hybridization of the attacking nitrogen.

Water-assisted dual mode cofactor recognition by HhaI DNA methyltransferase

Energetically competent binary recognition of the cofactor S-adenosyl-L-methionine (AdoMet) and the product S-adenosyl-L-homocysteine (AdoHcy) by the DNA (cytosine C-5) methyltransferase (M.HhaI) is demonstrated herein. Titration calorimetry reveals a dual mode, involving a primary dominant exothermic reaction followed by a weaker endothermic one, for the recognition of AdoMet and AdoHcy by M.HhaI. Conservation of the bimodal recognition in W41I and W41Y mutants of M.HhaI excludes the cation-pi interaction between the methylsulfonium group of AdoMet and the pi face of the Trp(41) indole ring from a role in its origin. Small magnitude of temperature-independent heat capacity changes upon AdoMet or AdoHcy binding by M.HhaI preclude appreciable conformational alterations in the reacting species. Coupled osmotic-calorimetric analyses of AdoMet and AdoHcy binding by M.HhaI indicate that a net uptake of nearly eight and 10 water molecules, respectively, assists their primary recognition. A change in water activity at constant temperature and pH is sufficient to engender and conserve enthalpy-entropy compensation, consistent with a true osmotic effect. The results implicate solvent reorganization in providing the major contribution to the origin of this isoequilibrium phenomenon in AdoMet and AdoHcy recognition by M.HhaI. The observations provide unequivocal evidence for the binding of AdoMet as well as AdoHcy to M.HhaI in solution state. Isotope partitioning analysis and preincubation studies favor a random mechanism for M.HhaI-catalyzed reaction. Taken together, the results clearly resolve the issue of cofactor recognition by free M.HhaI, specifically in the absence of DNA, leading to the formation of an energetically and catalytically competent binary complex.

Simulating enzyme reactions: challenges and perspectives

Elucidating how enzymes enhance the rates of the reactions that they catalyze is a major goal of contemporary biochemistry, and it is an area in which computational and theoretical techniques can make a major contribution. This article outlines some of the processes that need to be investigated if enzyme catalysis is to be understood, reviews the current state-of-the-art in enzyme simulation work, and highlights challenges for the future.

AdoMet-dependent methylation, DNA methyltransferases and base flipping

Twenty AdoMet-dependent methyltransferases (MTases) have been characterized structurally by X-ray crystallography and NMR. These include seven DNA MTases, five RNA MTases, four protein MTases and four small molecule MTases acting on the carbon, oxygen or nitrogen atoms of their substrates. The MTases share a common core structure of a mixed seven-stranded beta-sheet (6 downward arrow 7 upward arrow 5 downward arrow 4 downward arrow 1 downward arrow 2 downward arrow 3 downward arrow) referred to as an 'AdoMet-dependent MTase fold', with the exception of a protein arginine MTase which contains a compact consensus fold lacking the antiparallel hairpin strands (6 downward arrow 7 upward arrow). The consensus fold is useful to identify hypothetical MTases during structural proteomics efforts on unannotated proteins. The same core structure works for very different classes of MTase including those that act on substrates differing in size from small molecules (catechol or glycine) to macromolecules (DNA, RNA and protein). DNA MTases use a 'base flipping' mechanism to deliver a specific base within a DNA molecule into a typically concave catalytic pocket. Base flipping involves rotation of backbone bonds in double-stranded DNA to expose an out-of-stack nucleotide, which can then be a substrate for an enzyme-catalyzed chemical reaction. The phenomenon is fully established for DNA MTases and for DNA base excision repair enzymes, and is likely to prove general for enzymes that require access to unpaired, mismatched or damaged nucleotides within base-paired regions in DNA and RNA. Several newly discovered MTase families in eukaryotes (DNA 5mC MTases and protein arginine and lysine MTases) offer new challenges in the MTase field.

Computational enzymology

Recent advances in computational methods and the availability of fast, affordable computers have made the modeling of enzymatic reactions practical. The remaining challenges include achieving the accuracy level at which thermodynamic parameters and kinetic constants for different substrates, mutant enzymes, or in the presence of allosteric effectors can be predicted quantitatively.