Structural insights of the specificity and catalysis of a viral histone H3 lysine 27 methyltransferase

Parallel functional annotation of cancer-associated missense mutations in histone methyltransferases

Using exome sequencing for biomarker discovery and precision medicine requires connecting nucleotide-level variation with functional changes in encoded proteins. However, for functionally annotating the thousands of cancer-associated missense mutations, or variants of uncertain significance (VUS), purifying variant proteins for biochemical and functional analysis is cost-prohibitive and inefficient. We describe parallel functional annotation (PFA) of large numbers of VUS using small cultures and crude extracts in 96-well plates. Using members of a histone methyltransferase family, we demonstrate high-throughput structural and functional annotation of cancer-associated mutations. By combining functional annotation of paralogs, we discovered two phylogenetic and clustering parameters that improve the accuracy of sequence-based functional predictions to over 90%. Our results demonstrate the value of PFA for defining oncogenic/tumor suppressor functions of histone methyltransferases as well as enhancing the accuracy of sequence-based algorithms in predicting the effects of cancer-associated mutations.

Histone H3K27 methylation perturbs transcriptional robustness and underpins dispensability of highly conserved genes in fungi

Epigenetic modifications are key regulators of gene expression and underpin genome integrity. Yet, how epigenetic changes affect the evolution and transcriptional robustness of genes remains largely unknown. Here, we show how the repressive histone mark H3K27me3 underpins the trajectory of highly conserved genes in fungi. We first performed transcriptomic profiling on closely related species of the plant pathogen Fusarium graminearum species complex. We determined transcriptional responsiveness of genes across environmental conditions to determine expression robustness. To infer evolutionary conservation, we used a framework of 23 species across the Fusarium genus including three species covered with histone methylation data. Gene expression variation is negatively correlated with gene conservation confirming that highly conserved genes show higher expression robustness. In contrast, genes marked by H3K27me3 do not show such associations. Furthermore, highly conserved genes marked by H3K27me3 encode smaller proteins, exhibit weaker codon usage bias, higher levels of hydrophobicity, show lower intrinsically disordered regions, and are enriched for functions related to regulation and membrane transport. The evolutionary age of conserved genes with H3K27me3 histone marks falls typically within the origins of the Fusarium genus. We show that highly conserved genes marked by H3K27me3 are more likely to be dispensable for survival during host infection. Lastly, we show that conserved genes exposed to repressive H3K27me3 marks across distantly related Fusarium fungi are associated with transcriptional perturbation at the microevolutionary scale. In conclusion, we show how repressive histone marks are entangled in the evolutionary fate of highly conserved genes across evolutionary timescales.

Chloroviruses

Chloroviruses are large dsDNA, plaque-forming viruses that infect certain chlorella-like green algae; the algae are normally mutualistic endosymbionts of protists and metazoans and are often referred to as zoochlorellae. The viruses are ubiquitous in inland aqueous environments throughout the world and occasionally single types reach titers of thousands of plaque-forming units per ml of native water. The viruses are icosahedral in shape with a spike structure located at one of the vertices. They contain an internal membrane that is required for infectivity. The viral genomes are 290 to 370 kb in size, which encode up to 16 tRNAs and 330 to ~415 proteins, including many not previously seen in viruses. Examples include genes encoding DNA restriction and modification enzymes, hyaluronan and chitin biosynthetic enzymes, polyamine biosynthetic enzymes, ion channel and transport proteins, and enzymes involved in the glycan synthesis of the virus major capsid glycoproteins. The proteins encoded by many of these viruses are often the smallest or among the smallest proteins of their class. Consequently, some of the viral proteins are the subject of intensive biochemical and structural investigation.

The Polycomb group protein CLF emerges as a specific tri-methylase of H3K27 regulating gene expression and development in Physcomitrella patens

Packaging of eukaryotic DNA largely depends on histone modifications that affect the accessibility of DNA to transcriptional regulators, thus controlling gene expression. The Polycomb group (PcG) chromatin remodeling complex deposits a methyl group on lysine 27 of histone 3 leading to repressed gene expression. Plants encode homologs of the Enhancer of zeste (E(z)), a component of the PcG complex from Drosophila, one of which is a SET domain protein designated CURLY LEAF (CLF). Although this SET domain protein exhibits a strong correlation with the presence of the H3K27me3 mark in plants, the methyl-transferase activity and specificity of its SET domain have not been directly tested in-vivo. Using the evolutionary early-diverged land plant model species Physcomitrella patens we show that abolishment of a single copy gene PpCLF, as well as an additional member of the PcG complex, FERTILIZATION-INDEPENDENT ENDOSPERM (PpFIE), results in a specific loss of tri-methylation of H3K27. Using site-directed mutagenesis of key residues, we revealed that H3K27 tri-methylation is mediated by the SET domain of the CLF protein. Moreover, the abolishment of H3K27me3 led to enhanced expression of transcription factor genes. This in turn led to the development of fertilization-independent sporophyte-like structures, as observed in PpCLF and PpFIE null mutants. Overall, our results demonstrate the role of PpCLF as a SET protein in tri-methylation of H3K27 in-vivo and the importance of this modification in regulating the expression of transcription factor genes involved in developmental programs of P. patens.

Identification of family determining residues in Jumonji-C lysine demethylases: A sequence-based, family wide classification

Histone post-translational modifications play a critical role in the regulation of gene expression. Methylation of lysines at N-terminal tails of histones has been shown to be involved in such regulation. While this modification was long considered to be irreversible, two different classes of enzymes capable of carrying out the demethylation of histone lysines were recently identified: the oxidases, such as LSD1, and the oxygenases (JmjC-containing). Here, a family-wide analysis of the second of these classes is proposed, with over 300 proteins studied at the sequence level. We show that a correlated evolution analysis yields some position/residue pairs which are critical at comparing JmjC sequences and enables the classification of JmjC domains into five families. A few positions appear more frequently among conditions, such as positions 23 (directly C-terminal to the second iron ligand), 24, 252 and 253 (directly N-terminal to a conserved Asn). Implications of family conditions are studied in detail on PHF2, revealing the meaningfulness of the sequence-derived conditions at the structural level. These results should help obtain insights on the diversity of JmjC-containing proteins solely by considering some of the amino acids present in their JmjC domain.

SET domain-mediated lysine methylation in lower organisms regulates growth and transcription in hosts

Su(var)3-9, Enhancer-of-zeste, Trithorax (SET) domain-mediated lysine methylation, one of the major epigenetic marks, has been found to regulate chromatin-mediated gene transcription. Published studies have established further that methylation is not restricted to nuclear proteins but is involved in many cellular processes, including growth, differentiation, immune regulation, and cancer progression. The biological complexity of lysine methylation emerges from its capacity to cause gene activation or gene repression owing to the specific position of methylated-lysine moieties on the chromatin. Accumulating evidence suggests that despite the absence of chromatin, viruses and prokaryotes also express SET proteins, although their functional roles remain relatively less investigated. One possibility could be that SET proteins in lower organisms have more than one biological function, for example, in regulating growth or in manipulating host transcription machinery in order to establish infection. Thus, elucidating the role of an SET protein in host-pathogen interactions requires a thorough understanding of their functions. This review discusses the biological role of lysine methylation in prokaryotes and lower eukaryotes, as well as the underlying structural complexity and functional diversity of SET proteins.

Preparation, Biochemical Analysis, and Structure Determination of SET Domain Histone Methyltransferases

In eukaryotes, several lysine residues on histone proteins are methylated. This posttranslational modification is linked to a myriad of nuclear-based transactions such as epigenetic inheritance of heterochromatin, regulation of gene expression, DNA damage repair, and DNA replication. The majority of the enzymes responsible for writing these marks onto chromatin belong to the SET domain family of histone lysine methyltransferases. Although they often share important structural features, including a conserved catalytic domain, SET domain enzymes use different mechanisms to achieve substrate recognition, mono-, di-, or trimethylate lysine residues and some require other proteins to achieve maximal methyltransferase activity. In this chapter, we summarize our efforts to purify, crystallize, and enzymatically characterize SET domain enzymes with a specific focus on the histone H3K27 monomethyltransferase ATXR5.

[Urinary protein detection by iTRAQ® associated with renal transplant complications and its modification with therapy]

BACKGROUND

After renal transplant, surgical, infection complications, as well as graft rejection may occur; early detection through non-invasive markers is the key to change therapy and avoid biopsy.

OBJECTIVE

The aime of the study is to determine urine protein profiles in patients undergoing renal transplant with complications and detect its variation when therapy is modified.

MATERIAL AND METHODS

Urine samples were collected from patients prior the transplant and various postoperative stages. Urinary protein profiles were obtained by peptide labeling using isobaric isotopes for relative quantification (iTRAQ(®)).

RESULTS

A total of 22 patients were included, of whom 12 developed post-transplant complication: 2 with graft rejection (one male and one female) and 10 (6 males and 4 females) in the group of post-transplant infections. Using iTRAQ(®) 15/345 and 28/113 proteins were identified and fulfilled the acceptance criteria, in graft rejection and post-transplant infections group, respectively.

CONCLUSIONS

Albumin was the only protein found in both groups, the remaining proteins were different. The 5 proteins with higher scores in graft rejection were: alpha-1-microglobulin, 5'-nucleotidase cytosolic III, retinol-binding protein 4, membrane protein palmitoylated 4, and serine carboxypeptidase, while post-transplant infections were: mitochondrial acetyl-coenzyme A synthetase, putative adenosyl homocysteinase 2, zinc finger protein GLIS1, putative protein FAM157B, and zinc finger protein 615. It remains to elucidate the involvement of each of these in patients with renal transplantation.

Master manipulators: an update on Legionella pneumophila Icm/Dot translocated substrates and their host targets

Macrophages are the front line of immune defense against invading microbes. Microbes, however, have evolved numerous and diverse mechanisms to thwart these host immune defenses and thrive intracellularly. Legionella pneumophila, a Gram-negative pathogen of amoebal and mammalian phagocytes, is one such microbe. In humans, it causes a potentially fatal pneumonia referred to as Legionnaires' disease. Armed with the Icm/Dot type IV secretion system, which is required for virulence, and approximately 300 translocated proteins, Legionella is able to enter host cells, direct the biogenesis of its own vacuolar compartment, and establish a replicative niche, where it grows to high levels before lysing the host cell. Efforts to understand the pathogenesis of this bacterium have focused on characterizing the molecular activities of its many effectors. In this article, we highlight recent strides that have been made in understanding how Legionella effectors mediate host-pathogen interactions.

Biochemical reconstitution and phylogenetic comparison of human SET1 family core complexes involved in histone methylation

Mixed lineage leukemia protein-1 (MLL1) is a member of the SET1 family of histone H3 lysine 4 (H3K4) methyltransferases that are required for metazoan development. MLL1 is the best characterized human SET1 family member, which includes MLL1-4 and SETd1A/B. MLL1 assembles with WDR5, RBBP5, ASH2L, DPY-30 (WRAD) to form the MLL1 core complex, which is required for H3K4 dimethylation and transcriptional activation. Because all SET1 family proteins interact with WRAD in vivo, it is hypothesized they are regulated by similar mechanisms. However, recent evidence suggests differences among family members that may reflect unique regulatory inputs in the cell. Missing is an understanding of the intrinsic enzymatic activities of different SET1 family complexes under standard conditions. In this investigation, we reconstituted each human SET1 family core complex and compared subunit assembly and enzymatic activities. We found that in the absence of WRAD, all but one SET domain catalyzes at least weak H3K4 monomethylation. In the presence of WRAD, all SET1 family members showed stimulated monomethyltransferase activity but differed in their di- and trimethylation activities. We found that these differences are correlated with evolutionary lineage, suggesting these enzyme complexes have evolved to accomplish unique tasks within metazoan genomes. To understand the structural basis for these differences, we employed a "phylogenetic scanning mutagenesis" assay and identified a cluster of amino acid substitutions that confer a WRAD-dependent gain-of-function dimethylation activity on complexes assembled with the MLL3 or Drosophila trithorax proteins. These results form the basis for understanding how WRAD differentially regulates SET1 family complexes in vivo.

Designing synthetic drugs against Plasmodium falciparum: a computational study of histone-lysine N-methyltransferase (PfHKMT)

Histone lysine methyltransferase (HKMT) are histone-modifying enzymes that catalyze the transfer of methyl groups to lysine and arginine residues of histone protein. HKMTs have been involved in transcriptional regulation of various proteins in organisms. Malaria parasite also has HKMT, which plays a major role in parasite development and pathogenesis and also in regulation of various biological process and pathways. Our aim is to study fundamental biology of key molecules involved in the survival of Plasmodium falciparum and use these to develop efficient synthetic peptides and chemical compounds. As a first step in this direction, we computationally predicted the three-dimensional structure of HKMT of P. falciparum (PfHKMT) by using iterative threading assembly refinement. The PfHKMT three-dimensional model was validated using PROCHECK and docked with known HKMT inhibitor Bix01294 using Autodock. Our initial results are encouraging and indicate that structural analysis of PfHKMT could be important in developing novel synthetic molecules against malaria.

De novo peptide design and experimental validation of histone methyltransferase inhibitors

Histones are small proteins critical to the efficient packaging of DNA in the nucleus. DNA–protein complexes, known as nucleosomes, are formed when the DNA winds itself around the surface of the histones. The methylation of histone residues by enhancer of zeste homolog 2 (EZH2) maintains gene repression over successive cell generations. Overexpression of EZH2 can silence important tumor suppressor genes leading to increased invasiveness of many types of cancers. This makes the inhibition of EZH2 an important target in the development of cancer therapeutics. We employed a three-stage computational de novo peptide design method to design inhibitory peptides of EZH2. The method consists of a sequence selection stage and two validation stages for fold specificity and approximate binding affinity. The sequence selection stage consists of an integer linear optimization model that was solved to produce a rank-ordered list of amino acid sequences with increased stability in the bound peptide-EZH2 structure. These sequences were validated through the calculation of the fold specificity and approximate binding affinity of the designed peptides. Here we report the discovery of novel EZH2 inhibitory peptides using the de novo peptide design method. The computationally discovered peptides were experimentally validated in vitro using dose titrations and mechanism of action enzymatic assays. The peptide with the highest in vitro response, SQ037, was validated in nucleo using quantitative mass spectrometry-based proteomics. This peptide had an IC₅₀ of 13.5 M, demonstrated greater potency as an inhibitor when compared to the native and K27A mutant control peptides, and demonstrated competitive inhibition versus the peptide substrate. Additionally, this peptide demonstrated high specificity to the EZH2 target in comparison to other histone methyltransferases. The validated peptides are the first computationally designed peptides that directly inhibit EZH2. These inhibitors should prove useful for further chromatin biology investigations.

EZH2: biology, disease, and structure-based drug discovery

EZH2 is the catalytic subunit of the polycomb repressive complex 2 (PRC2), which is a highly conserved histone methyltransferase that methylates lysine 27 of histone 3. Overexpression of EZH2 has been found in a wide range of cancers, including those of the prostate and breast. In this review, we address the current understanding of the oncogenic role of EZH2, including its PRC2-dependent transcriptional repression and PRC2-independent gene activation. We also discuss the connections between EZH2 and other silencing enzymes, such as DNA methyltransferase and histone deacetylase. We comprehensively address the architecture of the PRC2 complex and the crucial roles of each subunit. Finally, we summarize new progress in developing EZH2 inhibitors, which could be a new epigenetic therapy for cancers.

Enzyme-dependent lysine deprotonation in EZH2 catalysis

Protein lysine methyltransferases (PKMTs) are key players in epigenetic regulation and have been associated with a variety of diseases, including cancers. The catalytic subunit of Polycomb Repressive Complex 2, EZH2 (EC 2.1.1.43), is a PKMT and a member of a family of SET domain lysine methyltransferases that catalyze the transfer of a methyl group from S-adenosyl-l-methionine to lysine 27 of histone 3 (H3K27). Wild-type (WT) EZH2 primarily catalyzes the mono- and dimethylation of H3K27; however, a clinically relevant active site mutation (Y641F) has been shown to alter the reaction specificity, dominantly catalyzing trimethylation of H3K27, and has been linked to tumor genesis and maintenance. Herein, we explore the chemical mechanism of methyl transfer by EZH2 and its Y641F mutant with pH-rate profiles and solvent kinetic isotope effects (sKIEs) using a short peptide derived from histone H3 [H3(21-44)]. A key component of the chemical reaction is the essential deprotonation of the ε-NH3(+) group of lysine to accommodate subsequent methylation. This deprotonation has been suggested by independent studies (1) to occur prior to binding to the enzyme (by bulk solvent) or (2) to be facilitated within the active site following binding, either (a) by the enzyme itself or (b) by a water molecule with access to the binding pocket. Our pH-rate and sKIE data best support a model in which lysine deprotonation is enzyme-dependent and at least partially rate-limiting. Furthermore, our experimental data are in agreement with prior computational models involving enzyme-dependent solvent deprotonation through a channel providing bulk solvent access to the active site. The mechanism of deprotonation and the rate-limiting catalytic steps appear to be unchanged between the WT and Y641F mutant enzymes, despite their activities being highly dependent on different substrate methylation states, suggesting determinants of substrate and product specificity in EZH2 are independent of catalytic events limiting the steady-state rate.

Convergent evolution of chromatin modification by structurally distinct enzymes: comparative enzymology of histone H3 Lys²⁷ methylation by human polycomb repressive complex 2 and vSET

H3K27 (histone H3 Lys27) methylation is an important epigenetic modification that regulates gene transcription. In humans, EZH (enhancer of zeste homologue) 1 and EZH2 are the only enzymes capable of catalysing methylation of H3K27. There is great interest in understanding structure-function relationships for EZH2, as genetic alterations in this enzyme are thought to play a causal role in a number of human cancers. EZH2 is challenging to study because it is only active in the context of the multi-subunit PRC2 (polycomb repressive complex 2). vSET is a viral lysine methyltransferase that represents the smallest protein unit capable of catalysing H3K27 methylation. The crystal structure of this minimal catalytic protein has been solved and researchers have suggested that vSET might prove useful as an EZH2 surrogate for the development of active site-directed inhibitors. To test this proposition, we conducted comparative enzymatic analysis of human EZH2 and vSET and report that, although both enzymes share similar preferences for methylation of H3K27, they diverge in terms of their permissiveness for catalysing methylation of alternative histone lysine sites, their relative preferences for utilization of multimeric macromolecular substrates, their active site primary sequences and, most importantly, their sensitivity to inhibition by drug-like small molecules. The cumulative data led us to suggest that EZH2 and vSET have very distinct active site structures, despite the commonality of the reaction catalysed by the two enzymes. Hence, the EZH2 and vSET pair of enzymes represent an example of convergent evolution in which distinct structural solutions have developed to solve a common catalytic need.

The promise and failures of epigenetic therapies for cancer treatment

Genetic mutations and gross structural defects in the DNA sequence permanently alter genetic loci in ways that significantly disrupt gene function. In sharp contrast, genes modified by aberrant epigenetic modifications remain structurally intact and are subject to partial or complete reversal of modifications that restore the original (i.e. non-diseased) state. Such reversibility makes epigenetic modifications ideal targets for therapeutic intervention. The epigenome of cancer cells is extensively modified by specific hypermethylation of the promoters of tumor suppressor genes relative to the extensive hypomethylation of repetitive sequences, overall loss of acetylation, and loss of repressive marks at microsatellite/repeat regions. In this review, we discuss emerging therapies targeting specific epigenetic modifications or epigenetic modifying enzymes either alone or in combination with other treatment regimens. The limitations posed by cancer treatments elicit unintended epigenetic modifications that result in exacerbation of tumor progression are also discussed. Lastly, a brief discussion of the specificity restrictions posed by epigenetic therapies and ways to address such limitations is presented.

Inner workings and regulatory inputs that control Polycomb repressive complex 2

Polycomb repressive complex 2 (PRC2) is a conserved multisubunit enzyme that methylates histone H3 on lysine-27. This chromatin modification is a hallmark of target genes transcriptionally silenced by the Polycomb system. At its core, PRC2 activity depends upon the SET domain active site of its catalytic subunit, EZH2, as well as critical stimulatory inputs from noncatalytic subunits, especially EED and SU(Z)12. We review recent progress on this core PRC2 machinery, including key features of the active site, control mechanisms that operate via EZH2 phosphorylation, and subunit elements and architectures that influence PRC2 function. Among these, we highlight work identifying an EED regulatory site that enables PRC2 to bind pre-existing methylated H3-K27 and stimulate enzyme output. These advances illuminate basic inner workings of PRC2 and also provide insights that could aid design of PRC2 inhibitors. The chromatin landscape that PRC2 encounters in vivo is decorated with many histone modifications that accompany active transcription, such as H3-K4 methylation. It has long been assumed that these "active" modifications oppose PRC2 at some level but, until recently, mechanisms of this antagonistic cross-talk have been elusive. We discuss new findings that illuminate how H3-K4 and H3-K36 methylation, H3-K27 acetylation, and H3-S28 phosphorylation each exert a negative impact on PRC2 function. The emerging picture presents PRC2 as a cooperative multipart machine, intricately outfitted to sense and respond to the local chromatin environment and other cues. This PRC2 design ensures flexibility and fine tuning of its fundamental gene silencing roles in diverse biological contexts.

A bacteria colony-based screen for optimal linker combinations in genetically encoded biosensors

Background

Fluorescent protein (FP)-based biosensors based on the principle of intramolecular Förster resonance energy transfer (FRET) enable the visualization of a variety of biochemical events in living cells. The construction of these biosensors requires the genetic insertion of a judiciously chosen molecular recognition element between two distinct hues of FP. When the molecular recognition element interacts with the analyte of interest and undergoes a conformational change, the ratiometric emission of the construct is altered due to a change in the FRET efficiency. The sensitivity of such biosensors is proportional to the change in ratiometric emission, and so there is a pressing need for methods to maximize the ratiometric change of existing biosensor constructs in order to increase the breadth of their utility.

Results

To accelerate the development and optimization of improved FRET-based biosensors, we have developed a method for function-based high-throughput screening of biosensor variants in colonies of Escherichia coli. We have demonstrated this technology by undertaking the optimization of a biosensor for detection of methylation of lysine 27 of histone H3 (H3K27). This effort involved the construction and screening of 3 distinct libraries: a domain library that included several engineered binding domains isolated by phage-display; a lower-resolution linker library; and a higher-resolution linker library.

Conclusion

Application of this library screening methodology led to the identification of an optimized H3K27-trimethylation biosensor that exhibited an emission ratio change (66%) that was 2.3 × improved relative to that of the initially constructed biosensor (29%).

Structural insights into the regulation and the recognition of histone marks by the SET domain of NSD1

The development of epigenetic therapies fuels cancer hope. DNA-methylation inhibitors, histone-deacetylase and histone-methyltransferase (HMTase) inhibitors are being developed as the utilization of epigenetic targets is emerging as an effective and valuable approach to chemotherapy as well as chemoprevention of cancer. The nuclear receptor binding SET domain (NSD) protein is a family of three HMTases, NSD1, NSD2/MMSET/WHSC1, and NSD3/WHSC1L1 that are critical in maintaining the chromatin integrity. A growing number of studies have reported alterations or amplifications of NSD1, NSD2, or NSD3 in numerous carcinogenic events. Reducing NSDs activity through specific lysine-HMTase inhibitors appears promising to help suppressing cancer growth. However, little is known about the NSD pathways and our understanding of the histone lysine-HMTase mechanism is partial. To shed some light on both the recognition and the regulation of epigenetic marks by the SET domain of the NSD family, we investigate the structural mechanisms of the docking of the histone-H4 tail on the SET domain of NSD1. Our finding exposes a key regulatory and recognition mechanism driven by the flexibility of a loop at the interface of the SET and postSET region. Finally, we prospect the special value of this regulatory region for developing specific and selective NSD inhibitors for the epigenetic therapy of cancers.

Cancers and the NSD family of histone lysine methyltransferases

Both genetic and epigenetic alterations are responsible for the stepwise initiation and progression of cancers. Only epigenetic aberrations can be reversible, allowing the malignant cell population to revert to a more benign phenotype. The epigenetic therapy of cancers is emerging as an effective and valuable approach to both the chemotherapy and the chemoprevention of cancer. The utilization of epigenetic targets that include histone methyltransferase (HMTase), Histone deacetylatase, and DNA methyltransferase, are emerging as key therapeutic targets. The nuclear receptor binding SET domain (NSD) protein is a family of three HMTases, NSD1, NSD2/MMSET/WHSC1, and NSD3/WHSC1L1, and plays a critical part in chromatin integrity as evidenced by a growing number of conditions linked to the alterations and/or amplification of NSD1, NSD2, and/or NSD3. NSD1, NSD2 and NSD3 are associated with multiple cancers. The amplification of either NSD1 or NSD2 triggers the cellular transformation and thus is key in the early carcinogenesis events. In most cases, reducing the levels of NSD proteins would suppress cancer growth. NSD1 and NSD2 were isolated as genes linked to developmental diseases, such as Sotos syndrome and Wolf-Hirschhorn syndrome, respectively, implying versatile aspects of the NSD proteins. The NSD pathways, however, are not well understood. It is noteworthy that the NSD family is phylogenetically distinct compared to other known lysine-HMTases, Here, we review the current knowledge on NSD1/NSD2/NSD3 in tumorigenesis and prospect their special value for developing novel anticancer drugs.

Structural and biochemical studies of human lysine methyltransferase Smyd3 reveal the important functional roles of its post-SET and TPR domains and the regulation of its activity by DNA binding

The SET- and MYND-domain containing (Smyd) proteins constitute a special subfamily of the SET-containing lysine methyltransferases. Here we present the structure of full-length human Smyd3 in complex with S-adenosyl-l-homocysteine at 2.8 Å resolution. Smyd3 affords the first example that other region(s) besides the SET domain and its flanking regions participate in the formation of the active site. Structural analysis shows that the previously uncharacterized C-terminal domain of Smyd3 contains a tetratrico-peptide repeat (TPR) domain which together with the SET and post-SET domains forms a deep, narrow substrate binding pocket. Our data demonstrate the important roles of both TPR and post-SET domains in the histone lysine methyltransferase (HKMT) activity of Smyd3, and show that the hydroxyl group of Tyr239 is critical for the enzymatic activity. The characteristic MYND domain is located nearby to the substrate binding pocket and exhibits a largely positively charged surface. Further biochemical assays show that DNA binding of Smyd3 can stimulate its HKMT activity and the process may be mediated via the MYND domain through direct DNA binding.

Structure and function of histone H3 lysine 9 methyltransferases and demethylases

Histone lysine methylation is a dynamic chromatin modification that plays key regulatory roles in gene expression and other genomic functions. Methylation of Lys9 in histone H3 (H3K9) is a prominent modification that has been implicated in diverse processes, including transcriptional silencing, heterochromatin formation, and DNA methylation. In this review, we summarize recent advances in understanding the structure and substrate specificity of the H3K9-specific methyltransferases G9A and GLP and explore current efforts to develop inhibitors of these enzymes. In addition, we discuss the structure and specificity of the recently discovered PHF8 family of histone demethylases that target H3K9 as well as other methylation sites in histones H3 and H4. Finally, we conclude by comparing the H3K9 binding modes displayed by these enzymes and examine the relevance of these studies to their biological functions and to structure-based inhibitor design.

An Ash2L/RbBP5 heterodimer stimulates the MLL1 methyltransferase activity through coordinated substrate interactions with the MLL1 SET domain

Histone H3 lysine 4 (K4) methylation is a prevalent mark associated with transcription activation and is mainly catalyzed by the MLL/SET1 family histone methyltransferases. A common feature of the mammalian MLL/SET1 complexes is the presence of three core components (RbBP5, Ash2L and WDR5) and a catalytic subunit containing a SET domain. Unlike most other histone lysine methyltransferases, all four proteins are required for efficient H3 K4 methylation. Despite extensive efforts, mechanisms for how three core components regulate MLL/SET1 methyltransferase activity remain elusive. Here we show that a heterodimer of Ash2L and RbBP5 has intrinsic histone methyltransferase activity. This activity requires the highly conserved SPRY domain of Ash2L and a short peptide of RbBP5. We demonstrate that both Ash2L and the MLL1 SET domain are capable of binding to S-adenosyl-L- [methyl-³H] methionine in the MLL1 core complex. Mutations in the MLL1 SET domain that fail to support overall H3 K4 methylation also compromise SAM binding by Ash2L. Taken together, our results show that the Ash2L/RbBP5 heterodimer plays a critical role in the overall catalysis of MLL1 mediated H3 K4 methylation. The results we describe here provide mechanistic insights for unique regulation of the MLL1 methyltransferase activity. It suggests that both Ash2L/RbBP5 and the MLL1 SET domain make direct contacts with the substrates and contribute to the formation of a joint catalytic center. Given the shared core configuration among all MLL/SET1 family HMTs, it will be interesting to test whether the mechanism we describe here can be generalized to other MLL/SET1 family members in the future.

Dimerization of a viral SET protein endows its function

Histone modifications are regarded as the most indispensible phenomena in epigenetics. Of these modifications, lysine methylation is of the greatest complexity and importance as site- and state-specific lysine methylation exerts a plethora of effects on chromatin structure and gene transcription. Notably, paramecium bursaria chlorella viruses encode a conserved SET domain methyltransferase, termed vSET, that functions to suppress host transcription by methylating histone H3 at lysine 27 (H3K27), a mark for eukaryotic gene silencing. Unlike mammalian lysine methyltransferases (KMTs), vSET functions only as a dimer, but the underlying mechanism has remained elusive. In this study, we demonstrate that dimeric vSET operates with negative cooperativity between the two active sites and engages in H3K27 methylation one site at a time. New atomic structures of vSET in the free form and a ternary complex with S-adenosyl homocysteine and a histone H3 peptide and biochemical analyses reveal the molecular origin for the negative cooperativity and explain the substrate specificity of H3K27 methyltransferases. Our study suggests a "walking" mechanism, by which vSET acts all by itself to globally methylate host H3K27, which is accomplished by the mammalian EZH2 KMT only in the context of the Polycomb repressive complex.

SET7/9 catalytic mutants reveal the role of active site water molecules in lysine multiple methylation

SET domain lysine methyltransferases (KMTs) methylate specific lysine residues in histone and non-histone substrates. These enzymes also display product specificity by catalyzing distinct degrees of methylation of the lysine ε-amino group. To elucidate the molecular mechanism underlying this specificity, we have characterized the Y245A and Y305F mutants of the human KMT SET7/9 (also known as KMT7) that alter its product specificity from a monomethyltransferase to a di- and a trimethyltransferase, respectively. Crystal structures of these mutants in complex with peptides bearing unmodified, mono-, di-, and trimethylated lysines illustrate the roles of active site water molecules in aligning the lysine ε-amino group for methyl transfer with S-adenosylmethionine. Displacement or dissociation of these solvent molecules enlarges the diameter of the active site, accommodating the increasing size of the methylated ε-amino group during successive methyl transfer reactions. Together, these results furnish new insights into the roles of active site water molecules in modulating lysine multiple methylation by SET domain KMTs and provide the first molecular snapshots of the mono-, di-, and trimethyl transfer reactions catalyzed by these enzymes.

Structure-function relationships of wheat flavone O-methyltransferase: Homology modeling and site-directed mutagenesis

BACKGROUND

Wheat (Triticum aestivum L.) O-methyltransferase (TaOMT2) catalyzes the sequential methylation of the flavone, tricetin, to its 3'-methyl- (selgin), 3',5'-dimethyl- (tricin) and 3',4',5'-trimethyl ether derivatives. Tricin, a potential multifunctional nutraceutical, is the major enzyme reaction product. These successive methylations raised the question as to whether they take place in one, or different active sites. We constructed a 3-D model of this protein using the crystal structure of the highly homologous Medicago sativa caffeic acid/5-hydroxyferulic acid O-methyltransferase (MsCOMT) as a template with the aim of proposing a mechanism for multiple methyl transfer reactions in wheat.

RESULTS

This model revealed unique structural features of TaOMT2 which permit the stepwise methylation of tricetin. Substrate binding is mediated by an extensive network of H-bonds and van der Waals interactions. Mutational analysis of structurally guided active site residues identified those involved in binding and catalysis. The partly buried tricetin active site, as well as proximity and orientation effects ensured sequential methylation of the substrate within the same pocket. Stepwise methylation of tricetin involves deprotonation of its hydroxyl groups by a His262-Asp263 pair followed by nucleophilic attack of SAM-methyl groups. We also demonstrate that Val309, which is conserved in a number of graminaceous flavone OMTs, defines the preference of TaOMT2 for tricetin as the substrate.

CONCLUSIONS

We propose a mechanism for the sequential methylation of tricetin, and discuss the potential application of TaOMT2 to increase the production of tricin as a nutraceutical. The single amino acid residue in TaOMT2, Val309, determines its preference for tricetin as the substrate, and may define the evolutionary differences between the two closely related proteins, COMT and flavone OMT.

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.

Mixed lineage leukemia: a structure-function perspective of the MLL1 protein

Several acute lymphoblastic and myelogenous leukemias are correlated with alterations in the human mixed lineage leukemia protein-1 (MLL1) gene. MLL1 is a member of the evolutionarily conserved SET1 family of histone H3 lysine 4 (H3K4) methyltransferases, which are required for the regulation of distinct groups of developmentally regulated genes in metazoans. Despite the important biological role of SET1 family enzymes and their involvement in human leukemias, relatively little is understood about how these enzymes work. Here we review several recent structural and biochemical studies that are beginning to shed light on the molecular mechanisms for the regulation of H3K4 methylation by the human MLL1 enzyme.

Structural biology of human H3K9 methyltransferases

SET domain methyltransferases deposit methyl marks on specific histone tail lysine residues and play a major role in epigenetic regulation of gene transcription. We solved the structures of the catalytic domains of GLP, G9a, Suv39H2 and PRDM2, four of the eight known human H3K9 methyltransferases in their apo conformation or in complex with the methyl donating cofactor, and peptide substrates. We analyzed the structural determinants for methylation state specificity, and designed a G9a mutant able to tri-methylate H3K9. We show that the I-SET domain acts as a rigid docking platform, while induced-fit of the Post-SET domain is necessary to achieve a catalytically competent conformation. We also propose a model where long-range electrostatics bring enzyme and histone substrate together, while the presence of an arginine upstream of the target lysine is critical for binding and specificity.

Enhanced version

This article can also be viewed as an enhanced version in which the text of the article is integrated with interactive 3D representations and animated transitions. Please note that a web plugin is required to access this enhanced functionality. Instructions for the installation and use of the web plugin are available in Text S1.

Viral-encoded enzymes that target host chromatin functions

Ever since their existence, there has been an everlasting arms race between viruses and their host cells. Host cells have developed numerous strategies to silence viral gene expression whereas viruses always find their ways to overcome these obstacles. Recent studies show that viruses have also evolved to take full advantage of existing cellular chromatin components to activate or repress its own genes when needed. While in most cases viruses encode certain proteins to recruit or inhibit cellular factors through physical interactions, growing examples show that viral-encoded enzymes affect host chromatin structure through post-translationally modifying histones or other cellular proteins important for chromatin function. The most well-studied example is vSET encoded by paramecium bursaria chlorella virus 1. vSET specifically methylates histone H3 at lysine 27, causing genome-wide silencing of Polycomb target genes upon infection, thus mimicking the function of Polycomb repressive complex 2 (PRC2) in eukaryotes. Other examples include BGLF4 from Epstein-Barr virus that affects both condensin and topoisomerase II activity and Us3 from Herpes Simplex virus 1 that inhibits HDAC1 function through phosphorylation.

On the mechanism of multiple lysine methylation by the human mixed lineage leukemia protein-1 (MLL1) core complex

Transcription in eukaryotic genomes depends on enzymes that regulate the degree of histone H3 lysine 4 (H3K4) methylation. The mixed lineage leukemia protein-1 (MLL1) is a member of the SET1 family of H3K4 methyltransferases and is frequently rearranged in acute leukemias. Despite sequence comparisons that predict that SET1 family enzymes should only monomethylate their substrates, mono-, di-, and trimethylation of H3K4 has been attributed to SET1 family complexes in vivo and in vitro. To better understand this paradox, we have biochemically reconstituted and characterized a five-component 200-kDa MLL1 core complex containing human MLL1, WDR5, RbBP5, Ash2L, and DPY-30. We demonstrate that the isolated MLL1 SET domain is a slow monomethyltransferase and that tyrosine 3942 of MLL1 prevents di- and trimethylation of H3K4. In contrast, a complex containing the MLL1 SET domain, WDR5, RbBP5, Ash2L, and DPY-30, displays a marked approximately 600-fold increase in enzymatic activity but only to the dimethyl form of H3K4. Single turnover kinetic experiments reveal that the reaction leading to H3K4 dimethylation involves the transient accumulation of a monomethylated species, suggesting that the MLL1 core complex uses a non-processive mechanism to catalyze multiple lysine methylation. We have also discovered that the non-SET domain components of the MLL1 core complex possess a previously unrecognized methyltransferase activity that catalyzes H3K4 dimethylation within the MLL1 core complex. Our results suggest that the mechanism of multiple lysine methylation by the MLL1 core complex involves the sequential addition of two methyl groups at two distinct active sites within the complex.

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.

Epigenetic transcriptional repression of cellular genes by a viral SET protein

Viruses recruit host proteins to secure viral genome maintenance and replication. However, whether they modify host histones directly to interfere with chromatin-based transcription is unknown. Here we report that Paramecium bursaria chlorella virus 1 (PBCV-1) encodes a functional SET domain histone Lys methyltransferase (HKMTase) termed vSET, which is linked to rapid inhibition of host transcription after viral infection. We show that vSET is packaged in the PBCV-1 virion, and that it contains a nuclear localization signal and probably represses host transcription by methylating histone H3 at Lys 27 (H3K27), a modification known to trigger gene silencing in eukaryotes. We also show that vSET induces cell accumulation at the G2/M phase by recruiting the Polycomb repressive complex CBX8 to the methylated H3K27 site in a heterologous system, vSET-like proteins that have H3K27 methylation activity are conserved in chlorella viruses. Our findings suggest a viral mechanism to repress gene transcription by direct modification of chromatin by PBCV-1 vSET.

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.

Dominant alleles identify SET domain residues required for histone methyltransferase of Polycomb repressive complex 2

Polycomb gene silencing requires histone methyltransferase activity of Polycomb repressive complex 2 (PRC2), which methylates lysine 27 of histone H3. Information on how PRC2 works is limited by lack of structural data on the catalytic subunit, Enhancer of zeste (E(Z)), and the paucity of E(z) mutant alleles that alter its SET domain. Here we analyze missense alleles of Drosophila E(z), selected for molecular study because of their dominant genetic effects. Four missense alleles identify key E(Z) SET domain residues, and a fifth is located in the adjacent CXC domain. Analysis of mutant PRC2 complexes in vitro, and H3-K27 methylation in vivo, shows that each SET domain mutation disrupts PRC2 histone methyltransferase. Based on known SET domain structures, the mutations likely affect either the lysine-substrate binding pocket, the binding site for the adenosylmethionine methyl donor, or a critical tyrosine predicted to interact with the substrate lysine epsilon-amino group. In contrast, the CXC mutant retains catalytic activity, Lys-27 specificity, and trimethylation capacity. Deletion analysis also reveals a functional requirement for a conserved E(Z) domain N-terminal to CXC and SET. These results identify critical SET domain residues needed for PRC2 enzyme function, and they also emphasize functional inputs from outside the SET domain.

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

Histone lysine methyltransferase SET7/9: formation of a water channel precedes each methyl transfer

Molecular dynamics (MD) simulations and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations have been carried out in an investigation of histone lysine methyltransferase (SET7/9). Proton dissociation (SET7/9.Lys4-NH3+.AdoMet --> SET7/9.Lys4-NH2.AdoMet + H+) must be prior to the methylation by S-adenosylmethionine (AdoMet). We find that a water channel is formed to allow escape of the proton to solvent. The water channel appears in the presence of AdoMet, but is not present in the species SET7/9.Lys4-NH3+ or SET7/9.Lys4-N(Me)H2+.AdoHcy. A water channel is not formed in the ground state of SET7/9.Lys4-N(Me)H2+.AdoMet, and the second methyl transfer does not occur. The structure of SET7/9.Lys4-N(Me)H2+.AdoMet includes a greater distance (6.1 +/- 0.3 A) between Cgamma(AdoMet) and N(MeLys4) than is present in SET7/9.Lys4-NH3+.AdoMet (5.7 +/- 0.2 A). The electrostatic interactions between the positive charges on AdoMet and SET7/9.Lys4-NH3+ decrease the pKa of the latter from 10.9 +/- 0.4 to 8.2 +/- 0.6, and this is not seen in the SET7/9.Lys4-N(Me)H2+.AdoMet species. The formation, or not, of a water channel, the distance between Sdelta(AdoMet) and N(Lys4), and the angle Sdelta(AdoMet)-Cgamma(AdoMet)-N(Lys4) determine whether methyl transfer can occur. By QM/MM, the calculated free energy barrier of the methyl transfer reaction in the SET7/9 [Lys4-NH2 + AdoMet --> Lys4-N(Me)H2+ + AdoHcy] complex is DeltaG++ = 19.0 +/- 1.6 kcal/mol. This DeltaG++ is in agreement with the value of 20.9 kcal/mol calculated from the experimental rate constant (0.24 min(-1)).

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.

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.

Histone-modifying enzymes: encrypting an enigmatic epigenetic code

Histone-modifying enzymes catalyze a diverse array of post-translational modifications of core and linker histones within chromatin. These modifications govern a multitude of genomic functions, particularly gene expression, and are believed to constitute an epigenetic code. Histone-modifying enzymes inscribe this code by catalyzing site-selective modifications, which are subsequently interpreted by effector proteins that recognize specific covalent marks. The substrate specificity of these enzymes is of fundamental biological importance because it underpins this epigenetic code. Recently, the structural basis of this specificity has been examined with regards to recently determined structures of GCN5 acetyltransferases and SET domain methyltransferases in complex with their cognate histone substrates.