The SET8 H4K20 protein lysine methyltransferase has a long recognition sequence covering seven amino acid residues

The SET8 histone lysine methyltransferase, which monomethylates the histone 4 lysine 20 residue plays important roles in cell cycle control and genomic stability. By employing peptide arrays we have shown that it has a long recognition sequence motif covering seven amino acid residues, viz. R¹⁷–H¹⁸–(R¹⁹KY)–K²⁰–(V²¹ILFY)–(L²²FY)–R²³. Celluspots peptide array methylation studies confirmed specific monomethylation of H4K20 and revealed that the symmetric and asymmetric methylation on R¹⁷ of the H4 tail inhibits methylation on H4K20. Similarly, dimethylation of the R located at the −3 position also reduced methylation of p53 K382 which had been shown previously to be methylated by SET8. Based on the derived specificity profile, we identified 4 potential non-histone substrate proteins. After relaxing the specificity profile, we identified several more candidate substrates and showed efficient methylation of 20 novel non-histone peptides by SET8. However, apart from H4 and p53 none of the identified novel peptide targets was methylated at the protein level. Since H4 and p53 both contain the target lysine in an unstructured part of the protein, we conclude that the long recognition sequence of SET8 makes it difficult to methylate a lysine in a folded region of a protein, because amino acid side chains essential for recognition will be buried.     

The Legionella pneumophila Methyltransferase RomA Methylates Also Non-histone Proteins during Infection

RomA is a SET-domain containing protein lysine methyltransferase encoded by the Gram-negative bacterium Legionella pneumophila. It is exported into human host cells during infection and has been previously shown to methylate histone H3 at lysine 14 [Rolando et al. (2013), Cell Host Microbe, 13, 395–405]. Here, we investigated the substrate specificity of RomA on peptide arrays showing that it mainly recognizes a G-K-X-(PA) sequence embedded in a basic amino acid sequence context. Based on the specificity profile, we searched for possible additional RomA substrates in the human proteome and identified 34 novel peptide substrates. For nine of these, the corresponding full-length protein or protein domains could be cloned and purified. Using radioactive and antibody-based methylation assays, we showed that seven of them are methylated by RomA, four of them strongly, one moderately, and two weakly. Mutagenesis confirmed for the seven methylated proteins that methylation occurs at target lysine residues fitting to the specificity profile. Methylation of one novel substrate (AROS) was investigated in HEK293 cells overexpressing RomA and during infection with L. pneumophila. Methylation could be detected in both conditions, confirming that RomA methylates non-histone proteins in human cells. Our data show that the bacterial methyltransferase RomA methylates also human non-histone proteins suggesting a multifaceted role in the infection process.     

Structural basis for the methylation site specificity of SET7/9

Human SET7/9 is a protein lysine methyltransferase (PKMT) that methylates histone H3, the tumor suppressor p53 and the TBP-associated factor TAF10. To elucidate the determinants of its substrate specificity, we have solved the enzyme's structure bound to a TAF10 peptide and examined its ability to methylate histone H3, TAF10 and p53 substrates bearing either mutations or covalent modifications within their respective methylation sites. Collectively, our data reveal that SET7/9 recognizes a conserved K/R-S/T/A motif preceding the lysine substrate and has a propensity to bind aspartates and asparagines on the C-terminal side of the lysine target. We then used a sequence-based approach with this motif to identify novel substrates for this PKMT. Among the putative targets is TAF7, which is methylated at Lys5 by the enzyme in vitro. These results demonstrate the predictive value of the consensus motif in identifying novel substrates for SET7/9.     

On your histone mark,SET, methylate!

Lysine methylation of histones and non-histone proteins has emerged in recent years as a posttranslational modification with wide-ranging cellular implications beyond epigenetic regulation. The molecular interactions between lysine methyltransferases and their substrates appear to be regulated by posttranslational modifications surrounding the lysine methyl acceptor. Two very interesting examples of this cross-talk between methyl-lysine sites are found in the SET (Su(var)3–9, Enhancer-of-zeste, Trithorax) domain-containing lysine methyltransferases SET7 and SETDB1, whereby the histone H3 trimethylated on lysine 4 (H3K4^me3) modification prevents methylation by SETDB1 on H3 lysine 9 (H3K9) and the histone H3 trimethylated on lysine 9 (H3K9^me3) modification prevents methylation by SET7 on H3K4. A similar cross-talk between posttranslational modifications regulates the functions of non-histone proteins such as the tumor suppressor p53 and the DNA methyltransferase DNMT1. Herein, in cis effects of acetylation, phosphorylation, as well as arginine and lysine methylation on lysine methylation events will be discussed.     

Targeting protein lysine methylation and demethylation in cancers

During the last decade, we saw an explosion of studies investigating the role of lysine methylation/demethylation of histones and non-histone proteins, such as p53, NF-kappaB, and E2F1. These 'Ying-Yang' post-translational modifications are important to fine-tuning the activity of these proteins. Lysine methylation and demethylation are catalyzed by protein lysine methyltransferases (PKMTs) and protein lysine demethylases (PKDMs). PKMTs, PKDMs, and their substrates have been shown to play important roles in cancers. Although the underlying mechanisms of tumorigenesis are still largely unknown, growing evidence is starting to link aberrant regulation of methylation to tumorigenesis. This review focuses on summarizing the recent progress in understanding of the function of protein lysine methylation, and in the discovery of small molecule inhibitors for PKMTs and PKDMs. We also discuss the potential and the caveats of targeting protein lysine methylation for the treatment of cancer.     

Epigenetic regulation: methylation of histone and non-histone proteins

Histone methylation is believed to play important roles in epigenetic memory in various biological processes. However, questions like whether the methylation marks themselves are faithfully transmitted into daughter cells and through what mechanisms are currently under active investigation. Previously, methylation was considered to be irreversible, but the recent discovery of histone lysine demethylases revealed a dynamic nature of histone methylation regulation on four of the main sites of methylation on histone H3 and H4 tails (H3K4, H3K9, H3K27 and H3K36). Even so, it is still unclear whether demethylases specific for the remaining two sites, H3K79 and H4K20, exist. Furthermore, besides histone proteins, the lysine methylation and demethylation also occur on non-histone proteins, which are probably subjected to similar regulation as histones. This review discusses recent progresses in protein lysine methylation regulation focusing on the above topics, while referring readers to a number of recent reviews for the biochemistry and biology of these enzymes     

A novel non-SET domain multi-subunit methyltransferase required for sequential nucleosomal histone H3 methylation by the mixed lineage leukemia protein-1 (MLL1) core complex

Gene expression within the context of eukaryotic chromatin is regulated by enzymes that catalyze histone lysine methylation. Histone lysine methyltransferases that have been identified to date possess the evolutionarily conserved SET or Dot1-like domains. We previously reported the identification of a new multi-subunit histone H3 lysine 4 methyltransferase lacking homology to the SET or Dot1 family of histone lysine methyltransferases. This enzymatic activity requires a complex that includes WRAD (WDR5, RbBP5, Ash2L, and DPY-30), a complex that is part of the MLL1 (mixed lineage leukemia protein-1) core complex but that also exists independently of MLL1 in the cell. Here, we report that the minimal complex required for WRAD enzymatic activity includes WDR5, RbBP5, and Ash2L and that DPY-30, although not required for enzymatic activity, increases the histone substrate specificity of the WRAD complex. We also show that WRAD requires zinc for catalytic activity, displays Michaelis-Menten kinetics, and is inhibited by S-adenosyl-homocysteine. In addition, we demonstrate that WRAD preferentially methylates lysine 4 of histone H3 within the context of the H3/H4 tetramer but does not methylate nucleosomal histone H3 on its own. In contrast, we find that MLL1 and WRAD are required for nucleosomal histone H3 methylation, and we provide evidence suggesting that each plays distinct structural and catalytic roles in the recognition and methylation of a nucleosome substrate. Our results indicate that WRAD is a new H3K4 methyltransferase with functions that include regulating the substrate and product specificities of the MLL1 core complex.     

Trimethylation of histone H3 lysine 4 impairs methylation of histone H3 lysine 9: Regulation of lysine methyltransferases by physical interaction with their substrates

Chromatin is broadly compartmentalized in two defined states: euchromatin and heterochromatin. Generally, euchromatin is trimethylated on histone H3 lysine 4 (H3K4^me3) while heterochromatin contains the H3K9^me3 mark. The H3K9^me3 modification is added by lysine methyltransferases (KMTs) such as SETDB1. Herein, we show that SETDB1 interacts with its substrate H3, but only in the absence of the euchromatic mark H3K4^me3. In addition, we show that SETDB1 fails to methylate substrates containing the H3K4^me3 mark. Likewise, the functionally related H3K9 KMTs G9A, GLP and SUV39H1 also fail to bind and to methylate H3K4^me3 substrates. Accordingly, we provide in vivo evidence that H3K9^me2-enriched histones are devoid of H3K4^me2/3 and that histones depleted of H3K4^me2/3 have elevated H3K9^me2/3. The correlation between the loss of interaction of these KMTs with H3K4^me3 and concomitant methylation impairment leads to the postulate that at least these four KMTs require stable interaction with their respective substrates for optimal activity. Thus, novel substrates could be discovered via the identification of KMT interacting proteins. Indeed, we find that SETDB1 binds to and methylates a novel substrate, the inhibitor of growth protein ING2, while SUV39H1 binds to and methylates the heterochromatin protein HP1α. Thus, our observations suggest a mechanism of post-translational regulation of lysine methylation and propose a potential mechanism for the segregation of the biologically opposing marks, H3K4^me3 and H3K9^me3. Furthermore, the correlation between H3-KMTs interaction and substrate methylation highlights that the identification of novel KMT substrates may be facilitated by the identification of interaction partners.Key words: histone methylation, lysine methyltransferase, H3K4^me3, H3K9^me3, SETDB1, G9A, ING2     

Trimethylation of histone H3 lysine 4 impairs methylation of histone H3 lysine 9

Chromatin is broadly compartmentalized in two defined states: euchromatin and heterochromatin. Generally, euchromatin is trimethylated on histone H3 lysine 4 (H3K4^me3) while heterochromatin contains the H3K9^me3 marks. The H3K9^me3 modification is added by lysine methyltransferases (KMTs) such as SETDB1. Herein, we show that SETDB1 interacts with its substrate H3, but only in the absence of the euchromatic mark H3K4^me3. In addition, we show that SETDB1 fails to methylate substrates containing the H3K4^me3 mark. Likewise, the functionally related H3K9 KMTs G9A, GLP, and SUV39H1 also fail to bind and to methylate H3K4^me3 substrates. Accordingly, we provide in vivo evidence that H3K9^me2-enriched histones are devoid of H3K4^me2/3 and that histones depleted of H3K4^me2/3 have elevated H3K9^me2/3. The correlation between the loss of interaction of these KMTs with H3K4^me3 and concomitant methylation impairment leads to the postulate that, at least these four KMTs, require stable interaction with their respective substrates for optimal activity. Thus, novel substrates could be discovered via the identification of KMT interacting proteins. Indeed, we find that SETDB1 binds to and methylates a novel substrate, the inhibitor of growth protein ING2, while SUV39H1 binds to and methylates the heterochromatin protein HP1α. Thus, our observations suggest a mechanism of post-translational regulation of lysine methylation and propose a potential mechanism for the segregation of the biologically opposing marks, H3K4^me3 and H3K9^me3. Furthermore, the correlation between H3-KMTs interaction and substrate methylation highlights that the identification of novel KMT substrates may be facilitated by the identification of interaction partners.     

A Proteomic Strategy Identifies Lysine Methylation of Splicing Factor snRNP70 by the SETMAR Enzyme

The lysine methyltransferase (KMT) SETMAR is implicated in the response to and repair of DNA damage, but its molecular function is not clear. SETMAR has been associated with dimethylation of histone H3 lysine 36 (H3K36) at sites of DNA damage. However, SETMAR does not methylate H3K36 in vitro. This and the observation that SETMAR is not active on nucleosomes suggest that H3K36 methylation is not a physiologically relevant activity. To identify potential non-histone substrates, we utilized a strategy on the basis of quantitative proteomic analysis of methylated lysine. Our approach identified lysine 130 of the mRNA splicing factor snRNP70 as a SETMAR substrate in vitro, and we show that the enzyme primarily generates monomethylation at this position. Furthermore, we show that SETMAR methylates snRNP70 Lys-130 in cells. Because snRNP70 is a key early regulator of 5′ splice site selection, our results suggest a model in which methylation of snRNP70 by SETMAR regulates constitutive and/or alternative splicing. In addition, the proteomic strategy described here is broadly applicable and is a promising route for large-scale mapping of KMT substrates.     

KDAC8 substrate specificity quantified by a biologically relevant,label‐free deacetylation assay

Analysis of the human proteome has identified thousands of unique protein sequences that contain acetylated lysine residues in vivo. These modifications regulate a variety of biological processes and are reversed by the lysine deacetylase (KDAC) family of enzymes. Despite the known prevalence and importance of acetylation, the details of KDAC substrate recognition are not well understood. While several methods have been developed to monitor protein deacetylation, none are particularly suited for identifying enzyme‐substrate pairs of label‐free substrates across the entire family of lysine deacetylases. Here, we present a fluorescamine‐based assay which is more biologically relevant than existing methods and amenable to probing substrate specificity. Using this assay, we evaluated the activity of KDAC8 and other lysine deacetylases, including a sirtuin, for several peptides derived from known acetylated proteins. KDAC8 showed clear preferences for some peptides over others, indicating that the residues immediately surrounding the acetylated lysine play an important role in substrate specificity. Steady‐state kinetics suggest that the sequence surrounding the acetylated lysine affects binding affinity and catalytic rate independently. Our results provide direct evidence that potential KDAC8 substrates previously identified through cell based experiments can be directly deacetylated by KDAC8. Conversely, the data from this assay did not correlate well with predictions from previous screens for KDAC8 substrates using less biologically relevant substrates and assay conditions. Combining results from our assay with mass spectrometry‐based experiments and cell‐based experiments will allow the identification of specific KDAC‐substrate pairs and lead to a better understanding of the biological consequences of these interactions.     

Catalytic properties and kinetic mechanism of human recombinant Lys-9 histone H3 methyltransferase SUV39H1: participation of the chromodomain in enzymatic catalysis

Histone H3 lysine 9 (H3K9) methylation is a major component of gene regulation and chromatin organization. SUV39H1 methylates H3K9 at the pericentric heterochromatin region and participates in the maintenance of genome stability. In this study, a recombinant purified SUV39H1 is used for substrate specificity and steady-state kinetic analysis with peptides representing the un- or dimethylated lysine 9 histone H3 tail or full-length human recombinant H3 (rH3). Recombinant SUV39H1 methylated its substrate via a nonprocessive mechanism. Binding of either peptide or AdoMet first to the enzyme made a catalytically competent binary complex. Product inhibition studies with SUV39H1 showed that S-adenosyl-l-homocysteine is a competitive inhibitor of S-adenosyl-l-methionine and a mixed inhibitor of substrate peptide. Similarly, the methylated peptide was a competitive inhibitor of the unmethylated peptide and a mixed inhibitor of AdoMet, suggesting a random mechanism in a bi-bi reaction for recombinant SUV39H1 in which either substrate can bind to the enzyme first and either product can release first. The turnover numbers (k(cat)) for the H3 tail peptide and rH3 were comparable (12 and 8 h(-)(1), respectively) compared to the value of 1.5 h(-)(1) for an identical dimethylated lysine 9 H3 tail peptide. The Michaelis constant for the methylated peptide (K(m)(pep)) was 13-fold lower compared to that of the unmethylated peptide. The Michaelis constants for AdoMet (K(m)(AdoMet)) were 12 and 6 microM for the unmethylated peptide substrate and rH3, respectively. A reduction in the level of methylation was observed at high concentrations of rH3, implying substrate inhibition. Deletion of the chromodomain or point mutation of the conserved amino acids, W64A or W67A, of SUV39H1 impaired enzyme activity despite the presence of an intact catalytic SET domain. Thus, SUV39H1 utilizes both the chromodomain and the SET domain for catalysis.     

Synthetic lethal screens identify gene silencing processes in yeast and implicate the acetylated amino terminus of Sir3 in recognition of the nucleosome core

Dot1 methylates histone H3 lysine 79 (H3K79) on the nucleosome core and is involved in Sir protein-mediated silencing. Previous studies suggested that H3K79 methylation within euchromatin prevents nonspecific binding of the Sir proteins, which in turn facilitates binding of the Sir proteins in unmethylated silent chromatin. However, the mechanism by which the Sir protein binding is influenced by this modification is unclear. We performed genome-wide synthetic genetic array (SGA) analysis and identified interactions of DOT1 with SIR1 and POL32. The synthetic growth defects found by SGA analysis were attributed to the loss of mating type identity caused by a synthetic silencing defect. By using epistasis analysis, DOT1, SIR1, and POL32 could be placed in different pathways of silencing. Dot1 shared its silencing phenotypes with the NatA N-terminal acetyltransferase complex and the conserved N-terminal bromo adjacent homology (BAH) domain of Sir3 (a substrate of NatA). We classified all of these as affecting a common silencing process, and we show that mutations in this process lead to nonspecific binding of Sir3 to chromatin. Our results suggest that the BAH domain of Sir3 binds to histone H3K79 and that acetylation of the BAH domain is required for the binding specificity of Sir3 for nucleosomes unmethylated at H3K79.     

Computational characterization of substrate and product specificities,and functionality of S‐adenosylmethionine binding pocket in histone lysine methyltransferases from Arabidopsis,rice and maize

Histone lysine methylation by histone lysine methyltransferases (HKMTs) has been implicated in regulation of gene expression. While significant progress has been made to understand the roles and mechanisms of animal HKMT functions, only a few plant HKMTs are functionally characterized. To unravel histone substrate specificity, degree of methylation and catalytic activity, we analyzed Arabidopsis Trithorax‐like protein (ATX), Su (var)3‐9 h omologs protein (SUVH), Su(var)3‐9 related protein (SUVR), ATXR5, ATXR6, and E(Z) HKMTs of Arabidopsis, maize and rice through sequence and structure comparison. We show that ATXs may exhibit methyltransferase specificity toward histone 3 lysine 4 (H3K4) and might catalyse the trimethylation. Our analyses also indicate that most SUVH proteins of Arabidopsis may bind histone H3 lysine 9 (H3K9). We also predict that SUVH7, SUVH8, SUVR1, SUVR3, ZmSET20 and ZmSET22 catalyse monomethylation or dimethylation of H3K9. Except for SDG728, which may trimethylate H3K9, all SUVH paralogs in rice may catalyse monomethylation or dimethylation. ZmSET11, ZmSET31, SDG713, SDG715, and SDG726 proteins are predicted to be catalytically inactive because of an incomplete S‐adenosylmethionine (SAM) binding pocket and a post‐SET domain. E(Z) homologs can trimethylate H3K27 substrate, which is similar to the Enhancer of Zeste homolog 2 of humans. Our comparative sequence analyses reveal that ATXR5 and ATXR6 lack motifs/domains required for protein‐protein interaction and polycomb repressive complex 2 complex formation. We propose that subtle variations of key residues at substrate or SAM binding pocket, around the catalytic pocket, or presence of pre‐SET and post‐SET domains in HKMTs of the aforementioned plant species lead to variations in class‐specific HKMT functions and further determine their substrate specificity, the degree of methylation and catalytic activity.     

Structural Basis of a Histone H3 Lysine 4 Demethylase Required for Stem Elongation in Rice

Histone lysine methylation is an important epigenetic modification in regulating chromatin structure and gene expression. Histone H3 lysine 4 methylation (H3K4me), which can be in a mono-, di-, or trimethylated state, has been shown to play an important role in gene expression involved in plant developmental control and stress adaptation. However, the resetting mechanism of this epigenetic modification is not yet fully understood. In this work, we identified a JmjC domain-containing protein, JMJ703, as a histone lysine demethylase that specifically reverses all three forms of H3K4me in rice. Loss-of-function mutation of the gene affected stem elongation and plant growth, which may be related to increased expression of cytokinin oxidase genes in the mutant. Analysis of crystal structure of the catalytic core domain (c-JMJ703) of the protein revealed a general structural similarity with mammalian and yeast JMJD2 proteins that are H3K9 and H3K36 demethylases. However, several specific features were observed in the structure of c-JMJ703. Key residues that interact with cofactors Fe(II) and N-oxalylglycine and the methylated H3K4 substrate peptide were identified and were shown to be essential for the demethylase activity in vivo. Several key residues are specifically conserved in known H3K4 demethylases, suggesting that they may be involved in the specificity for H3K4 demethylation.     

植物组蛋白赖氨酸甲基化建立过程及其遗传性研究进展

表观遗传学主要包括DNA甲基化、组蛋白修饰和非编码RNA,组蛋白甲基化作为组蛋白修饰中的一种重要修饰,在植物体的发育和环境适应中发挥着重要作用。组蛋白甲基化主要发生在赖氨酸残基上,同时根据不同的赖氨酸位点和每个赖氨酸位点甲基化程度的不同,形成了不同的赖氨酸甲基化修饰。根据对基因的不同功能,通常将组蛋白赖氨酸甲基化修饰分为2大类:(1)能够促进基因表达的,如H3K4me3和H3K36me3;(2)能够抑制基因表达的,如H3K9me2和H3K27me3。不同的组蛋白赖氨酸甲基化去甲基化过程需要相应的阅读(reader)、书写(writer)和擦除(eraser)3种蛋白。同时,组蛋白赖氨酸甲基化的遗传性质目前还不是很清楚。综述了植物中组蛋白赖氨酸甲基化建立与去除过程,以及对组蛋白赖氨酸甲基化可遗传性的探讨。     

组蛋白甲基转移酶G9a在表观遗传调控中的作用

组蛋白赖氨酸甲基化在表观遗传调控中起着关键作用。组蛋白甲基转移酶G9a(又称作常染色质组蛋白赖氨酸N-甲基转移酶2(euchromatic histone-lysine N-methyltransferase 2,EHMT2))含经典的SET结构域,是常染色质主要的甲基转移酶之一,可以甲基化组蛋白H3K9、H3K27和H1bK26等。此外,G9a也可以直接甲基化一些非组蛋白,并与DNA甲基化密切相关。G9a功能紊乱可以导致胚胎发育异常、免疫系统及神经系统发育障碍、甚至癌症的发生发展。     

Molecular and structural insight into lysine selection on substrate and ubiquitin lysine 48 by the ubiquitin-conjugating enzyme Cdc34

The attachment of ubiquitin (Ub) to lysines on substrates or itself by ubiquitin-conjugating (E2) and ubiquitin ligase (E3) enzymes results in protein ubiquitination. Lysine selection is important for generating diverse substrate-Ub structures and targeting proteins to different fates; however, the mechanisms of lysine selection are not clearly understood. The positioning of lysine(s) toward the E2/E3 active site and residues proximal to lysines are critical in their selection. We investigated determinants of lysine specificity of the ubiquitin-conjugating enzyme Cdc34, toward substrate and Ub lysines. Evaluation of the relative importance of different residues positioned −2, −1, +1 and +2 toward ubiquitination of its substrate, Sic1, on lysine 50 showed that charged residues in the −1 and −2 positions negatively impact on ubiquitination. Modeling suggests that charged residues at these positions alter the native salt-bridge interactions in Ub and Cdc34, resulting in misplacement of Sic1 lysine 50 in the Cdc34 catalytic cleft. During polyubiquitination, Cdc34 showed a strong preference for Ub lysine 48 (K48), with lower activity towards lysine 11 (K11) and lysine 63 (K63). Mutating the −2, −1, +1 and +2 sites surrounding K11 and K63 to mimic those surrounding K48 did not improve their ubiquitination, indicating that further determinants are important for Ub K48 specificity. Modeling the ternary structure of acceptor Ub with the Cdc34~Ub complex as well as in vitro ubiquitination assays unveiled the importance of K6 and Q62 of acceptor Ub for Ub K48 polyubiquitination. These findings provide molecular and structural insight into substrate lysine and Ub K48 specificity by Cdc34.     

QM/MM MD and free energy simulations of G9a-like protein (GLP) and its mutants: understanding the factors that determine the product specificity

Certain lysine residues on histone tails could be methylated by protein lysine methyltransferases (PKMTs) using S-adenosyl-L-methionine (AdoMet) as the methyl donor. Since the methylation states of the target lysines play a fundamental role in the regulation of chromatin structure and gene expression, it is important to study the property of PKMTs that allows a specific number of methyl groups (one, two or three) to be added (termed as product specificity). It has been shown that the product specificity of PKMTs may be controlled in part by the existence of specific residues at the active site. One of the best examples is a Phe/Tyr switch found in many PKMTs. Here quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) and free energy simulations are performed on wild type G9a-like protein (GLP) and its F1209Y and Y1124F mutants for understanding the energetic origin of the product specificity and the reasons for the change of product specificity as a result of single-residue mutations at the Phe/Tyr switch as well as other positions. The free energy barriers of the methyl transfer processes calculated from our simulations are consistent with experimental data, supporting the suggestion that the relative free energy barriers may determine, at least in part, the product specificity of PKMTs. The changes of the free energy barriers as a result of the mutations are also discussed based on the structural information obtained from the simulations. The results suggest that the space and active-site interactions around the ε-amino group of the target lysine available for methyl addition appear to among the key structural factors in controlling the product specificity and activity of PKMTs.