Catalytic roles for carbon-oxygen hydrogen bonding in SET domain lysine methyltransferases

Product specificity and mechanism of protein lysine methyltransferases: insights from the histone lysine methyltransferase SET8

Molecular dynamics simulations employing a molecular mechanics (MM) force field and hybrid quantum mechanics (QM) and MM (QM/MM) have been carried out to investigate the product specificity and mechanism of the histone H4 lysine 20 (H4-K20) methylation by human histone lysine methyltransferase SET8. At neutral pH, the target lysine is available to only the enzyme in the protonated state. The first step in the methylation reaction must be deprotonation of the lysine target which is followed by the (+)AdoMet methylation of the neutral lysine [Enz.Lys-CH(2)-NH(3)(+).(+)AdoMet --> H(+) + Enz.Lys-CH(2)-NH(2).(+)AdoMet -->--> Enz.Lys-CH(2)-N(Me)H(2)(+).AdoHcy]. The electrostatic interactions between two positive charges on (+)AdoMet and Lys20-NH(3)(+) decrease the pK(a) of Lys20-NH(3)(+). Upon formation of Enz.Lys-NH(3)(+).(+)AdoMet, a water channel by which the proton escapes to the outer solvent phase is formed. The formation of a water channel for the escape of a proton from Lys20-N(Me)H(2)(+) in Enz.Lys20-N(Me)H(2)(+).(+)AdoMet is not formed because the methyl substituent blocks the starting of the water channel. Thus, a second methylation does not take place. The dependence of the occurrence of methyl transfer on the formation of a water channel in SET8 is in accord with our previous reports on product specificity by histone lysine monomethyltransferase SET7/9, large subunit lysine dimethyltransferase (LSMT), and viral histone lysine trimethyltransferase (vSET). The average value of the experimental DeltaG(E)() for the six lysine methyl transfer reactions catalyzed by vSET, LSMT, and SET7/9 with p53 as a substrate is 22.1 +/- 1.0 kcal/mol, and the computed average (DeltaG(C)()) is 22.2 +/- 0.8 kcal/mol. In this study, the computed free energy barrier of the methyl transfer reaction [Lys20-NH(2) + (+)AdoMet --> Lys20-N(Me)H(2)(+) + AdoHcy] catalyzed by SET8 is 20.8 kcal/mol. This is in agreement with the value of 20.6 kcal/mol calculated from the experimental rate constant (0.43 +/- 0.02 min(-1)). Our bond-order computations establish that the H4-K20 monomethylation in SET8 is a concerted linear S(N)2 displacement reaction.

Arabidopsis Histone Lysine Methyltransferases

In eukaryotes, changes in chromatin structure regulate the access of gene regulatory sequences to the transcriptional machinery and play important roles in the repression of transposable elements, thereby protecting genome integrity. Chromatin dynamics and gene expression states are highly correlated, with DNA methylation and histone post-translational modifications playing important roles in the establishment or maintenance of chromatin states in plants. Histones can be covalently modified in a variety of ways, thereby affecting nucleosome spacing and/or higher-order nucleosome interactions directly or via the recruitment of histone-binding proteins. An extremely important group of chromatin modifying enzymes are the histone lysine methyltransferases (HKMTs). These enzymes are involved in the establishment and/or maintenance of euchromatic or heterochromatic states of active or transcriptionally repressed sequences, respectively. The vast majority of HKMTs possess a SET domain named for the three Drosophila proteins that are the founding members of the family: Suppressor of variegation, Enhancer of zeste and Trithorax. It is the SET domain that is responsible for HKMT enzymatic activity. Mutation of Arabidopsis HKMT genes can result in phenotypic abnormalities due to the improper regulation of important developmental genes. Here, we review the different classes of HKMTs present in the model plant Arabidopsis thaliana and discuss what is known about their biochemical and biological functions.

Enzymatic and structural insights for substrate specificity of a family of jumonji histone lysine demethylases

Combinatorial readout of multiple covalent histone modifications is poorly understood. We provide insights into how an activating histone mark, in combination with linked repressive marks, is differentially 'read' by two related human demethylases, PHF8 and KIAA1718 (also known as JHDM1D). Both enzymes harbor a plant homeodomain (PHD) that binds Lys4-trimethylated histone 3 (H3K4me3) and a jumonji domain that demethylates either H3K9me2 or H3K27me2. The presence of H3K4me3 on the same peptide as H3K9me2 makes the doubly methylated peptide a markedly better substrate of PHF8, whereas the presence of H3K4me3 has the opposite effect, diminishing the H3K9me2 demethylase activity of KIAA1718 without adversely affecting its H3K27me2 activity. The difference in substrate specificity between the two is explained by PHF8 adopting a bent conformation, allowing each of its domains to engage its respective target, whereas KIAA1718 adopts an extended conformation, which prevents its access to H3K9me2 by its jumonji domain when its PHD engages H3K4me3.

Protein methyltransferases as a target class for drug discovery

The protein methyltransferases (PMTs) - which methylate protein lysine and arginine residues and have crucial roles in gene transcription - are emerging as an important group of enzymes that play key parts in normal physiology and human diseases. The collection of human PMTs is a large and diverse group of enzymes that have a common mechanism of catalysis. Here, we review the biological, biochemical and structural data that together present PMTs as a novel, chemically tractable target class for drug discovery.

Structural basis for G9a-like protein lysine methyltransferase inhibition by BIX-01294

Histone lysine methylation is an important epigenetic mark that regulates gene expression and chromatin organization. G9a and G9a-like protein (GLP) are euchromatin-associated methyltransferases that repress transcription by methylating histone H3 Lys9. BIX-01294 was originally identified as a G9a inhibitor during a chemical library screen of small molecules and has previously been used in the generation of induced pluripotent stem cells. Here we present the crystal structure of the catalytic SET domain of GLP in complex with BIX-01294 and S-adenosyl-L-homocysteine. The inhibitor is bound in the substrate peptide groove at the location where the histone H3 residues N-terminal to the target lysine lie in the previously solved structure of the complex with histone peptide. The inhibitor resembles the bound conformation of histone H3 Lys4 to Arg8, and is positioned in place by residues specific for G9a and GLP through specific interactions.

Rubisco in complex with Rubisco large subunit methyltransferase

SET domain protein lysine methyltransferases (PKMT) are a structurally unique class of enzymes that catalyze the specific methylation of lysine residues in a number of different substrates. Especially histone-specific SET domain PKMTs have received widespread attention because of their roles in the regulation of epigenetic gene expression and the development of some cancers. Rubisco large subunit methyltransferase (RLSMT) is a chloroplast-localized SET domain PKMT responsible for the formation of trimethyl-lysine-14 in the large subunit of Rubisco, an essential photosynthetic enzyme. Here, we have used cryoelectron microscopy to produce an 11-A density map of the Rubisco-RLSMT complex. The atomic model of the complex, obtained by fitting crystal structures of Rubisco and RLSMT into the density map, shows that the extensive contact regions between the 2 proteins are mainly mediated by hydrophobic residues and leucine-rich repeats. It further provides insights into potential conformational changes that may occur during substrate binding and catalysis. This study presents the first structural analysis of a SET domain PKMT in complex with its intact polypeptide substrate.

Identification of stable S-adenosylmethionine (SAM) analogues derivatised with bioorthogonal tags: effect of ligands on the affinity of the E. coli methionine repressor, MetJ, for its operator DNA

The efficient synthesis of a range of stable SAM mimetics, and their ability to promote the binding of the E. coli methionine repressor (MetJ) to its operator DNA, is described. Active analogues functionalised with bioorthogonal tags were identified.

The efficient synthesis of a range of stable SAM mimetics, and their ability to promote the binding of the E. coli methionine repressor (MetJ) to its operator DNA, is described.

Structural origins for the product specificity of SET domain protein methyltransferases

SET domain protein lysine methyltransferases (PKMTs) regulate transcription and other cellular functions through site-specific methylation of histones and other substrates. PKMTs catalyze the formation of monomethylated, dimethylated, or trimethylated products, establishing an additional hierarchy with respect to methyllysine recognition in signaling. Biochemical studies of PKMTs have identified a conserved position within their active sites, the Phe/Tyr switch, that governs their respective product specificities. To elucidate the mechanism underlying this switch, we have characterized a Phe/Tyr switch mutant of the histone H4 Lys-20 (H4K20) methyltransferase SET8, which alters its specificity from a monomethyltransferase to a dimethyltransferase. The crystal structures of the SET8 Y334F mutant bound to histone H4 peptides bearing unmodified, monomethyl, and dimethyl Lys-20 reveal that the phenylalanine substitution attenuates hydrogen bonding to a structurally conserved water molecule adjacent to the Phe/Tyr switch, facilitating its dissociation. The additional space generated by the solvent's dissociation enables the monomethyllysyl side chain to adopt a conformation that is catalytically competent for dimethylation and furnishes sufficient volume to accommodate the dimethyl epsilon-ammonium product. Collectively, these results indicate that the Phe/Tyr switch regulates product specificity through altering the affinity of an active-site water molecule whose dissociation is required for lysine multiple methylation.

The methyl-CpG-binding protein MeCP2 and neurological disease

The methyl-CpG-binding protein MeCP2 was discovered over 15 years ago as part of a search for proteins that selectively bind methylated DNA. It is a nuclear protein that is largely chromatin-bound and has a strong preference for binding to methylated DNA sequences in vivo. Evidence from model systems shows that MeCP2 can recruit the Sin3a co-repressor complex to promoters leading to transcriptional repression, therefore suggesting that MeCP2 can interpret the DNA methylation signal to bring about gene silencing. Mutations in the human MECP2 gene cause the autism spectrum disorder Rett Syndrome. MeCP2 is most highly expressed in neurons, and mice lacking this protein show symptoms that strikingly parallel those of Rett patients. Surprisingly, these symptoms are efficiently reversed by delayed activation of a ‘stopped’ Mecp2 gene, raising hopes that human Rett syndrome may also be reversible. Future studies of MeCP2 promise to shed light upon brain function, neurological disease and the biology of DNA methylation.

Chemical mechanisms of histone lysine and arginine modifications

Histone lysine and arginine residues are subject to a wide array of post-translational modifications including methylation, citrullination, acetylation, ubiquitination, and sumoylation. The combinatorial action of these modifications regulates critical DNA processes including replication, repair, and transcription. In addition, enzymes that modify histone lysine and arginine residues have been correlated with a variety of human diseases including arthritis, cancer, heart disease, diabetes, and neurodegenerative disorders. Thus, it is important to fully understand the detailed kinetic and chemical mechanisms of these enzymes. Here, we review recent progress towards determining the mechanisms of histone lysine and arginine modifying enzymes. In particular, the mechanisms of S-adenosyl-methionine (AdoMet) dependent methyltransferases, FAD-dependent demethylases, iron dependent demethylases, acetyl-CoA dependent acetyltransferases, zinc dependent deacetylases, NAD(+) dependent deacetylases, and protein arginine deiminases are covered. Particular attention is paid to the conserved active-site residues necessary for catalysis and the individual chemical steps along the catalytic pathway. When appropriate, areas requiring further work are discussed.

MeCP2 binding to DNA depends upon hydration at methyl-CpG

MeCP2 is an essential transcriptional repressor that mediates gene silencing through binding to methylated DNA. Binding specificity has been thought to depend on hydrophobic interactions between cytosine methyl groups and a hydrophobic patch within the methyl-CpG-binding domain (MBD). X-ray analysis of a methylated DNA-MBD cocrystal reveals, however, that the methyl groups make contact with a predominantly hydrophilic surface that includes tightly bound water molecules. This suggests that MeCP2 recognizes hydration of the major groove of methylated DNA rather than cytosine methylation per se. The MeCP2-DNA complex also identifies a unique structural role for T158, the residue most commonly mutated in Rett syndrome.

Enzymatic mechanism and product specificity of SET-domain protein lysine methyltransferases

Molecular dynamics and hybrid quantum mechanics/molecular mechanics have been used to investigate the mechanisms of (+)AdoMet methylation of protein-Lys-NH(2) catalyzed by the lysine methyltransferase enzymes: histone lysine monomethyltransferase SET7/9, Rubisco large-subunit dimethyltransferase, viral histone lysine trimethyltransferase, and the Tyr245Phe mutation of SET7/9. At neutrality in aqueous solution, primary amines are protonated. The enzyme reacts with Lys-NH(3)(+) and (+)AdoMet species to provide an Enz.Lys-NH(3)(+).(+)AdoMet complex. The close positioning of two positive charges lowers the pK(a) of the Lys-NH(3)(+) entity, a water channel appears, and the proton escapes to the aqueous solvent; then the reaction Enz.Lys-NH(2).(+)AdoMet --> Enz.Lys-N(Me)H(2)(+).AdoHcy occurs. Repeat of the sequence provides dimethylated lysine, and another repeat yields a trimethylated lysine. The sequence is halted at monomethylation when the conformation of the Enz.Lys-N(Me)H(2)(+).(+)AdoMet has the methyl positioned to block formation of a water channel. The sequence of reactions stops at dimethylation if the conformation of Enz.Lys-N(Me)(2)H(+).(+)AdoMet has a methyl in position, which forbids the formation of the water channel.

How do SET-domain protein lysine methyltransferases achieve the methylation state specificity? Revisited by Ab initio QM/MM molecular dynamics simulations

A distinct protein lysine methyltransferase (PKMT) only transfers a certain number of methyl group(s) to its target lysine residue in spite of the fact that a lysine residue can be either mono-, di-, or tri-methylated. In order to elucidate how such a remarkable product specificity is achieved, we have carried out ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations on two SET-domain PKMTs: SET7/9 and Rubisco large subunit methyltransferase (LSMT). The results indicate that the methylation state specificity is mainly controlled by the methyl-transfer reaction step, and confirm that SET7/9 is a mono-methyltransferase while LSMT has both mono-and di-methylation activities. It is found that the binding of the methylated lysine substrate in the active site of SET7/ 9 opens up the cofactor AdoMet binding channel so that solvent water molecules get access to the active site. This disrupts the catalytic machinery of SET7/9 for the di-methylation reaction, which leads to a higher activation barrier, whereas for the LSMT, its active site is more spacious than that of SET7/9, so that the methylated lysine substrate can be accommodated without interfering with its catalytic power. These detailed insights take account of protein dynamics and are consistent with available experimental results as well as recent theoretical findings regarding the catalytic power of SET7/9.

Mechanism of product specificity of AdoMet methylation catalyzed by lysine methyltransferases: transcriptional factor p53 methylation by histone lysine methyltransferase SET7/9

The catalysis by SET7/9 histone lysine methyltransferase of AdoMet N-methylation of the transcriptional factor p53-Lys4-NH 2 has been investigated with particular attention paid to the means of product specificity. After formation of the SET7/9.p53-Lys4-NH 3 (+).AdoMet complex, the following events occur: (i) the appearance of a water channel, (ii) a depronation of p53-Lys4-NH 3 (+) via this water channel into the aqueous solvent, and (iii) AdoMet methylation of p53-Lys4-NH 2 to form p53-Lys4-N(Me)H 2 (+). The formation of a water channel does not occur on formation of the SET7/9.p53-Lys4-NH 3 (+), SET7/9.p53-Lys4-N(Me)H 2 (+).AdoHcy, or SET7/9.p53-Lys4-N(Me)H 2 (+).AdoMet complex. Without a water channel, the substrate p53-Lys4-N(Me)H is not available because the proton dissociation p53-Lys4-N(Me)H 2 (+) --> p53-Lys4-N(Me)H + H (+) does not occur. The lack of formation of a water channel is due to the positioning of the methyl substituent of the SET7/9.p53-Lys4-N(Me)H 2 (+).AdoMet complex. By quantum mechanics/molecular mechanics, the computed free energy barrier of the methyl transfer reaction [p53-Lys4-NH 2 + AdoMet --> p53-Lys4-N(Me)H 2 (+) + AdoHcy] in the SET7/9 complex is Delta G (++) = 20.1 +/- 2.9 kcal/mol. This Delta G (++) is in agreement with the value of 20.9 kcal/mol calculated from the experimental rate constant (1.2 +/- 0.1 min (-1)). Our bond-order computations establish that the methyl transfer reaction in protein lysine methyltransferases occurs via a linear S N2 associative reaction with bond making of approximately 50%.

A quantum mechanics/molecular mechanics study of the catalytic mechanism and product specificity of viral histone lysine methyltransferase

There are three reaction steps in the S-adenosylmethionine (AdoMet) methylation of lysine-NH2 catalyzed by a methyltransferase. They are (i) combination of enzyme.Lys-NH3+ with AdoMet, (ii) substrate ionization to provide enzyme.AdoMet.Lys-NH2, and (iii) methyl transfer providing enzyme.AdoHcy.Lys-N(Me)H2+ and the dissociation of AdoHcy. In this study of the viral histone methyltransferase (vSET), we find that substrate ionization of vSET.Lys27-NH3+, vSET.Lys27-N(Me)H2+, and vSET.Lys27-N(Me)2H+ takes place upon combination with AdoMet. The presence of a water channel allows dissociation of a proton to the solvent. There is no water channel in the absence of AdoMet. That the formation of a water channel is combined with AdoMet binding was first discovered in our investigation of Rubisco large subunit methyltransferase. Via a quantum mechanics/molecular mechanics (QM/MM) approach, the calculated free energy barrier (DeltaG++) of the first methyl transfer reaction catalyzed by vSET [Lys27-NH2 + AdoMet --> Lys27-N(Me)H2+ + AdoHcy] equals 22.5 +/- 4.3 kcal/mol, which is in excellent agreement with the free energy barrier (21.7 kcal/mol) calculated from the experimental rate constant (0.047 min-1). The calculated DeltaG++ of the second methyl transfer reaction [AdoMet + Lys27-N(Me)H --> AdoHcy + Lys27-N(Me)2H+] at the QM/MM level is 22.6 +/- 3.6 kcal/mol, which is in agreement with the value of 22.4 kcal/mol determined from the experimental rate constant (0.015 min-1). The third methylation [Lys27-N(Me)2 + AdoMet --> Lys27-N(Me)3+ + AdoHcy] is associated with a DeltaG++ of 23.1 +/- 4.0 kcal/mol, which is in agreement with the value of 23.0 kcal/mol determined from the experimental rate constant (0.005 min-1). Our computations establish that the first, second, and third methyl transfer steps catalyzed by vSET are linear SN2 reactions with the bond making being approximately 50% associative.

Specificity and mechanism of JMJD2A, a trimethyllysine-specific histone demethylase

JMJD2A is a JmjC histone demethylase (HDM) that catalyzes the demethylation of di- and trimethylated Lys9 and Lys36 in histone H3 (H3K9me2/3 and H3K36me2/3). Here we present the crystal structures of the JMJD2A catalytic domain in complex with H3K9me3, H3K36me2 and H3K36me3 peptides. The structures reveal that histone substrates are recognized through a network of backbone hydrogen bonds and hydrophobic interactions that deposit the trimethyllysine into the active site. The trimethylated epsilon-ammonium cation is coordinated within a methylammonium-binding pocket through carbon-oxygen (CH...O) hydrogen bonds that position one of the zeta-methyl groups adjacent to the Fe(II) center for hydroxylation and demethylation. Mutations of the residues comprising this pocket abrogate demethylation by JMJD2A, with the exception of an S288A substitution, which augments activity, particularly toward H3K9me2. We propose that this residue modulates the methylation-state specificities of JMJD2 enzymes and other trimethyllysine-specific JmjC HDMs.

Mechanism of histone methylation catalyzed by protein lysine methyltransferase SET7/9 and origin of product specificity

Methylation of certain lysine residues in the N-terminal tails of core histone proteins in nucleosome is of fundamental importance in the regulation of chromatin structure and gene expression. Such histone modification is catalyzed by protein lysine methyltransferases (PKMTs). PKMTs contain a conserved SET domain in almost all of the cases and may transfer one to three methyl groups from S-adenosyl-L-methionine (AdoMet) to the epsilon-amino group of the target lysine residue. Here, quantum mechanical/molecular mechanical molecular dynamics and free-energy simulations are performed on human PKMT SET7/9 and its mutants to understand two outstanding questions for the reaction catalyzed by PKMTs: the mechanism for deprotonation of positively charged methyl lysine (lysine) and origin of product specificity. The results of the simulations suggest that Tyr-335 (an absolute conserved residue in PKMTs) may play the role as the general base for the deprotonation after dissociation of AdoHcy (S-adenosyl-L-homocysteine) and before binding of AdoMet. It is shown that conformational changes could bring Y335 to the target methyl lysine (lysine) for proton abstraction. This mechanism provides an explanation why methyl transfers could be catalyzed by PKMTs processively. The free-energy profiles for methyl transfers are reported and analyzed for wild type and certain mutants (Y305F and Y335F) and the active-site interactions that are of importance for the enzyme's function are discussed. The results of the simulations provide important insights into the catalytic process and lead to a better understanding of experimental observations concerning the origin of product specificity for PKMTs.

Catalytic mechanism and product specificity of rubisco large subunit methyltransferase: QM/MM and MD investigations

Molecular dynamics (MD) simulations and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations have been carried out in an investigation of Rubisco large subunit methyltransferase (LSMT). It was found that the appearance of a water channel is required for the stepwise methylation by S-adenosylmethionine (AdoMet). The water channel appears in the presence of AdoMet (LSMT.Lys-NH3+.AdoMet), but is not present immediately after methyl transfer (LSMT.Lys-N(Me)H2+.AdoHcy). The water channel allows proton dissociation from both LSMT.AdoMet.Lys-NH3+ and LSMT.AdoMet.Lys-N(Me)H2+. The water channel does not appear for proton dissociation from LSMT.AdoMet.Lys-N(Me)2H+, and a third methyl transfer does not occur. By QM/MM, the calculated free energy barrier of the first methyl transfer reaction catalyzed by LSMT (Lys-NH2 + AdoMet --> Lys-N(Me)H2+ + AdoHcy) is DeltaG++ = 22.8 +/- 3.3 kcal/mol. This DeltaG++ is in remarkable agreement with the value 23.0 kcal/mol calculated from the experimental rate constant (6.2 x 10-5 s-1). The calculated DeltaG++ of the second methyl transfer reaction (AdoMet + Lys-N(Me)H --> AdoHcy + Lys-N(Me)2H+) at the QM/MM level is 20.5 +/- 3.6 kcal/mol, which is in agreement with the value 22.0 kcal/mol calculated from the experimental rate constant (2.5 x 10-4 s-1). The third methyl transfer (Lys-N(Me)2 + AdoMet --> Lys-N(Me)3+ + AdoHcy) is associated with an allowed DeltaG++ of 25.9 +/- 3.2 kcal/mol. However, this reaction does not occur because a water channel does not form to allow the proton dissociation of Lys-N(Me)2H+. Future studies will determine whether the product specificity of lysine (mono, di, and tri) methyltransferases is determined by the formation of water channels.

Ab initio quantum mechanical/molecular mechanical molecular dynamics simulation of enzyme catalysis: the case of histone lysine methyltransferase SET7/9

To elucidate enzyme catalysis through computer simulation, a prerequisite is to reliably compute free energy barriers for both enzyme and solution reactions. By employing on-the-fly Born-Oppenheimer molecular dynamics simulations with the ab initio quantum mechanical/molecular mechanical approach and the umbrella sampling method, we have determined free energy profiles for the methyl-transfer reaction catalyzed by the histone lysine methyltransferase SET7/9 and its corresponding uncatalyzed reaction in aqueous solution, respectively. Our calculated activation free energy barrier for the enzyme catalyzed reaction is 22.5 kcal/mol, which agrees very well with the experimental value of 20.9 kcal/mol. The difference in potential of mean force between a corresponding prereaction state and the transition state for the solution reaction is computed to be 30.9 kcal/mol. Thus, our simulations indicate that the enzyme SET7/9 plays an essential catalytic role in significantly lowering the barrier for the methyl-transfer reaction step. For the reaction in solution, it is found that the hydrogen bond network near the reaction center undergoes a significant change, and there is a strong shift in electrostatic field from the prereaction state to the transition state, whereas for the enzyme reaction, such an effect is much smaller and the enzyme SET7/9 is found to provide a preorganized electrostatic environment to facilitate the methyl-transfer reaction. Meanwhile, we find that the transition state in the enzyme reaction is a little more dissociative than that in solution.

Kinetic manifestation of processivity during multiple methylations catalyzed by SET domain protein methyltransferases

Processive versus distributive methyl group transfer was assessed for pea Rubisco large subunit methyltransferase, a SET domain protein lysine methyltransferase catalyzing the formation of trimethyllysine-14 in the large subunit of Rubisco. Catalytically competent complexes between an immobilized form of des(methyl) Rubisco and Rubisco large subunit methyltransferase were used to demonstrate enzyme release that was co-incident with and dependent on formation of trimethyllysine. Catalytic rate constants determined for formation of trimethyllysine were considerably lower ( approximately 10-fold) than rate constants determined for total radiolabel incorporation from [3H-methyl]-S-adenosylmethionine. Double-reciprocal velocity plots under catalytic conditions favoring monomethyllysine indicated a random or ordered reaction mechanism, while conditions favoring trimethyllysine suggested a hybrid ping-pong mechanism. These results were compared with double-reciprocal velocity plots and product analyses obtained for HsSET7/9 (a monomethyltransferase) and SpCLR4 (a dimethyltransferase) and suggest a predictive ability of double-reciprocal velocity plots for single versus multiple methyl group transfers by SET domain protein lysine methyltransferases. A model is proposed for SET domain protein lysine methyltransferases in which initial binding of polypeptide substrate and S-adenosylmethionine is random, with polypeptide binding followed by deprotonation of the epsilon-amine of the target lysyl residue and subsequent methylation. Following methyl group transfer, S-adenosylhomocysteine and monomethylated polypeptide dissociate from monomethyltransferases, but di- and trimethyltransferases begin a successive and catalytically obligatory deprotonation of enzyme-bound methylated lysyl intermediates, which along with binding and release of S-adenosylmethionine and S-adenosylhomocysteine is manifested as a hybrid ping-pong-like reaction mechanism.