Molecular recognition of H3/H4 histone tails by the tudor domains of JMJD2A: a comparative molecular dynamics simulations study

The emerging roles of lysine-specific demethylase 4A in cancer: Implications in tumorigenesis and therapeutic opportunities

Lysine-specific demethylase 4 A (KDM4A, also named JMJD2A, KIA0677, or JHDM3A) is a demethylase that can remove methyl groups from histones H3K9me2/3, H3K36me2/3, and H1.4K26me2/me3. Accumulating evidence suggests that KDM4A is not only involved in body homeostasis (such as cell proliferation, migration and differentiation, and tissue development) but also associated with multiple human diseases, especially cancers. Recently, an increasing number of studies have shown that pharmacological inhibition of KDM4A significantly attenuates tumor progression in vitro and in vivo in a range of solid tumors and acute myeloid leukemia. Although there are several reviews on the roles of the KDM4 subfamily in cancer development and therapy, all of them only briefly introduce the roles of KDM4A in cancer without systematically summarizing the specific mechanisms of KDM4A in various physiological and pathological processes, especially in tumorigenesis, which greatly limits advances in the understanding of the roles of KDM4A in a variety of cancers, discovering targeted selective KDM4A inhibitors, and exploring the adaptive profiles of KDM4A antagonists. Herein, we present the structure and functions of KDM4A, simply outline the functions of KDM4A in homeostasis and non-cancer diseases, summarize the role of KDM4A and its distinct target genes in the development of a variety of cancers, systematically classify KDM4A inhibitors, summarize the difficulties encountered in the research of KDM4A and the discovery of related drugs, and provide the corresponding solutions, which would contribute to understanding the recent research trends on KDM4A and advancing the progression of KDM4A as a drug target in cancer therapy.

The role of lysine-specific demethylase 6A (KDM6A) in tumorigenesis and its therapeutic potentials in cancer therapy

Histone demethylation is a key post-translational modification of chromatin, and its dysregulation affects a wide array of nuclear activities including the maintenance of genome integrity, transcriptional regulation, and epigenetic inheritance. Lysine specific demethylase 6A (KDM6A, also known as UTX) is an Fe²⁺- and α-ketoglutarate- dependent oxidase which belongs to KDM6 Jumonji histone demethylase subfamily, and it can remove mono-, di- and tri-methyl groups from methylated lysine 27 of histone H3 (H3K27me1/2/3). Mounting studies indicate that KDM6A is responsible for driving multiple human diseases, particularly cancers and pharmacological inhibition of KDM6A is an effective strategy to treat varieties of KDM6A-amplified cancers in cellulo and in vivo. Although there are several reviews on the roles of KDM6 subfamily in cancer development and therapy, all of them only simply introduce the roles of KDM6A in cancer without systematically summarizing the specific mechanisms of KDM6A in tumorigenesis, which greatly limits the advances on the understanding of roles KDM6A in varieties of cancers, discovering targeting selective KDM6A inhibitors, and exploring the adaptive profiles of KDM6A antagonists. Herein, we present the structure and functions of KDM6A, simply outline the functions of KDM6A in homeostasis and non-cancer diseases, summarize the role of KDM6A and its distinct target genes/ligand proteins in development of varieties of cancers, systematically classify KDM6A inhibitors, sum up the difficulties encountered in the research of KDM6A and the discovery of related drugs, and provide the corresponding solutions, which will contribute to understanding the roles of KDM6A in carcinogenesis and advancing the progression of KDM6A as a drug target in cancer therapy.

Mechanistic aspects of reversible methylation modifications of arginine and lysine of nuclear histones and their roles in human colon cancer

Developmental proceedings and maintenance of cellular homeostasis are regulated by the precise orchestration of a series of epigenetic events that eventually control gene expression. DNA methylation and post-translational modifications (PTMs) of histones are well-characterized epigenetic events responsible for fine-tuning gene expression. PTMs of histones bear molecular logic of gene expression at chromosomal territory and have become a fascinating field of epigenetics. Nowadays, reversible methylation on histone arginine and lysine is gaining increasing attention as a significant PTM related to reorganizing local nucleosomal structure, chromatin dynamics, and transcriptional regulation. It is now well-accepted and reported that histone marks play crucial roles in colon cancer initiation and progression by encouraging abnormal epigenomic reprogramming. It is becoming increasingly clear that multiple PTM marks at the N-terminal tails of the core histones cross-talk with one another to intricately regulate DNA-templated biological processes such as replication, transcription, recombination, and damage repair in several malignancies, including colon cancer. These functional cross-talks provide an additional layer of message, which spatiotemporally fine-tunes the overall gene expression regulation. Nowadays, it is evident that several PTMs instigate colon cancer development. How colon cancer-specific PTM patterns or codes are generated and how they affect downstream molecular events are uncovered to some extent. Future studies would address more about epigenetic communication, and the relationship between histone modification marks to define cellular functions in depth. This chapter will comprehensively highlight the importance of histone arginine and lysine-based methylation modifications and their functional cross-talk with other histone marks from the perspective of colon cancer development.

Engineering human JMJD2A tudor domains for an improved understanding of histone peptide recognition

JMJD2A is a histone lysine demethylase which recognizes and demethylates H3K9me3 and H3K36me3 residues and is overexpressed in various cancers. It utilizes a tandem tudor domain to facilitate its own recruitment to histone sites, recognizing various di- and tri-methyl lysine residues with moderate affinity. In this study, we successfully engineered the tudor domain of JMJD2A to specifically bind to H4K20me3 with a 20-fold increase of affinity and improved selectivity. To reveal the molecular basis, we performed molecular dynamics and free energy decomposition analysis on the human JMJD2A tandem tudor domains bound to H4K20me2, H4K20me3, and H3K23me3 peptides to uncover the residues and conformational changes important for the enhanced binding affinity and selectivity toward H4K20me2/3. These analyses revealed new insights into understanding chromatin reader domains recognizing histone modifications and improving binding affinity and selectivity of these domains. Furthermore, we showed that the tight binding of JMJD2A to H4K20me2/3 is not sufficient to improve the efficiency of CRISPR-CAS9 mediated homology directed repair (HDR), suggesting a complicated relationship between JMJD2A and the DNA damage response beyond binding affinity toward the H4K20me2/3 mark.

KDM4 Involvement in Breast Cancer and Possible Therapeutic Approaches

Breast cancer (BC) is the second leading cause of cancer death in women, although recent scientific and technological achievements have led to significant improvements in progression-free disease and overall survival of patients. Genetic mutations and epigenetic modifications play a critical role in deregulating gene expression, leading to uncontrolled cell proliferation and cancer progression. Aberrant histone modifications are one of the most frequent epigenetic mechanisms occurring in cancer. In particular, methylation and demethylation of specific lysine residues alter gene accessibility via histone lysine methyltransferases (KMTs) and histone lysine demethylases (KDMs). The KDM family includes more than 30 members, grouped into six subfamilies and two classes based on their sequency homology and catalytic mechanisms, respectively. Specifically, the KDM4 gene family comprises six members, KDM4A-F, which are associated with oncogene activation, tumor suppressor silencing, alteration of hormone receptor downstream signaling, and chromosomal instability. Blocking the activity of KDM4 enzymes renders them "druggable" targets with therapeutic effects. Several KDM4 inhibitors have already been identified as anticancer drugs in vitro in BC cells. However, no KDM4 inhibitors have as yet entered clinical trials due to a number of issues, including structural similarities between KDM4 members and conservation of the active domain, which makes the discovery of selective inhibitors challenging. Here, we summarize our current knowledge of the molecular functions of KDM4 members in BC, describe currently available KDM4 inhibitors, and discuss their potential use in BC therapy.

The COMPASS-like complex modulates fungal development and pathogenesis by regulating H3K4me3-mediated targeted gene expression in Magnaporthe oryzae

Histone-3-lysine-4 (H3K4) methylation is catalysed by the multiprotein complex known as the Set1/COMPASS or MLL/COMPASS-like complex, an element that is highly evolutionarily conserved from yeast to humans. However, the components and mechanisms by which the COMPASS-like complex targets the H3K4 methylation of plant-pathogenic genes in fungi remain elusive. Here we present a comprehensive analysis combining biochemical, molecular, and genome-wide approaches to characterize the roles of the COMPASS-like family in the rice blast fungus Magnaporthe oryzae, a model plant pathogen. We purified and identified six conserved subunits of COMPASS from M. oryzae: MoBre2 (Cps60/ASH2L), MoSpp1 (Cps40/Cfp1), MoSwd2 (Cps35), MoSdc1 (Cps25/DPY30), MoSet1 (MLL/ALL), and MoRbBP5 (Cps50), using an affinity tag on MoBre2. We determined the sequence repeat in dual-specificity kinase splA and ryanodine receptors domain of MoBre2 can interact directly with the DPY30 domain of MoSdc1 in vitro. Furthermore, we found that deletion of the genes encoding COMPASS subunits of MoBre2, MoSPP1, and MoSwd2 caused similar defects regarding invasive hyphal development and pathogenicity. Genome-wide profiling of H3K4me3 revealed that it has remarkable co-occupancy at the transcription start site regions of target genes. Significantly, these target genes are often involved in spore germination and pathogenesis. Decreased gene expression caused by the deletion of MoBre2, MoSwd2, or MoSpp1 was highly correlated with a decrease in H3K4me3. These results suggest that MoBre2, MoSpp1, and MoSwd2 function as a whole COMPASS complex, contributing to fungal development and pathogenesis by regulating H3K4me3-targeted genes in M. oryzae.

Unresolved Complexity in the Gene Regulatory Network Underlying EMT

Epithelial to mesenchymal transition (EMT) is the process whereby a polarized epithelial cell ceases to maintain cell-cell contacts, loses expression of characteristic epithelial cell markers, and acquires mesenchymal cell markers and properties such as motility, contractile ability, and invasiveness. A complex process that occurs during development and many disease states, EMT involves a plethora of transcription factors (TFs) and signaling pathways. Whilst great advances have been made in both our understanding of the progressive cell-fate changes during EMT and the gene regulatory networks that drive this process, there are still gaps in our knowledge. Epigenetic modifications are dynamic, chromatin modifying enzymes are vast and varied, transcription factors are pleiotropic, and signaling pathways are multifaceted and rarely act alone. Therefore, it is of great importance that we decipher and understand each intricate step of the process and how these players at different levels crosstalk with each other to successfully orchestrate EMT. A delicate balance and fine-tuned cooperation of gene regulatory mechanisms is required for EMT to occur successfully, and until we resolve the unknowns in this network, we cannot hope to develop effective therapies against diseases that involve aberrant EMT such as cancer. In this review, we focus on data that challenge these unknown entities underlying EMT, starting with EMT stimuli followed by intracellular signaling through to epigenetic mechanisms and chromatin remodeling.

Understanding the Histone DNA Repair Code: H4K20me2 Makes Its Mark

Chromatin is a highly compact structure that must be rapidly rearranged in order for DNA repair proteins to access sites of damage and facilitate timely and efficient repair. Chromatin plasticity is achieved through multiple processes, including the posttranslational modification of histone tails. In recent years, the impact of histone posttranslational modification on the DNA damage response has become increasingly well recognized, and chromatin plasticity has been firmly linked to efficient DNA repair. One particularly important histone posttranslational modification process is methylation. Here, we focus on the regulation and function of H4K20 methylation (H4K20me) in the DNA damage response and describe the writers, erasers, and readers of this important chromatin mark as well as the combinatorial histone posttranslational modifications that modulate H4K20me recognition. Finally, we discuss the central role of H4K20me in determining if DNA double-strand breaks (DSB) are repaired by the error-prone, nonhomologous DNA end joining pathway or the error-free, homologous recombination pathway. This review article discusses the regulation and function of H4K20me2 in DNA DSB repair and outlines the components and modifications that modulate this important chromatin mark and its fundamental impact on DSB repair pathway choice. Mol Cancer Res; 16(9); 1335-45. ©2018 AACR.

Small molecule KDM4s inhibitors as anti-cancer agents

Histone demethylation is a vital process in epigenetic regulation of gene expression. A number of histone demethylases are present to control the methylated states of histone. Among these enzymes, KDM4s are one subfamily of JmjC KDMs and play important roles in both normal and cancer cells. The discovery of KDM4s inhibitors is a potential therapeutic strategy against different diseases including cancer. Here, we summarize the development of KDM4s inhibitors and some related pharmaceutical information to provide an update of recent progress in KDM4s inhibitors.

Computational Analysis and Predictive Cheminformatics Modeling of Small Molecule Inhibitors of Epigenetic Modifiers

BACKGROUND

The dynamic and differential regulation and expression of genes is majorly governed by the complex interactions of a subset of biomolecules in the cell operating at multiple levels starting from genome organisation to protein post-translational regulation. The regulatory layer contributed by the epigenetic layer has been one of the favourite areas of interest recently. This layer of regulation as we know today largely comprises of DNA modifications, histone modifications and noncoding RNA regulation and the interplay between each of these major components. Epigenetic regulation has been recently shown to be central to development of a number of disease processes. The availability of datasets of high-throughput screens for molecules for biological properties offer a new opportunity to develop computational methodologies which would enable in-silico screening of large molecular libraries.

METHODS

In the present study, we have used data from high throughput screens for the inhibitors of epigenetic modifiers. Computational predictive models were constructed based on the molecular descriptors. Machine learning algorithms for supervised training, Naive Bayes and Random Forest, were used to generate predictive models for the small molecule inhibitors of histone methyl-transferases and demethylases. Random forest, with the accuracy of 80%, was identified as the most accurate classifier. Further we complemented the study with substructure search approach filtering out the probable pharmacophores from the active molecules leading to drug molecules.

RESULTS

We show that effective use of appropriate computational algorithms could be used to learn molecular and structural correlates of biological activities of small molecules. The computational models developed could be potentially used to screen and identify potential new biological activities of molecules from large molecular libraries and prioritise them for in-depth biological assays. To the best of our knowledge, this is the first and most comprehensive computational analysis towards understanding activities of small molecules inhibitors of epigenetic modifiers.

Molecular features of interaction between VEGFA and anti-angiogenic drugs used in retinal diseases: a computational approach

Anti-angiogenic agents are biological drugs used for treatment of retinal neovascular degenerative diseases. In this study, we aimed at in silico analysis of interaction of vascular endothelial growth factor A (VEGFA), the main mediator of angiogenesis, with binding domains of anti-angiogenic agents used for treatment of retinal diseases, such as ranibizumab, bevacizumab and aflibercept. The analysis of anti-VEGF/VEGFA complexes was carried out by means of protein-protein docking and molecular dynamics (MD) coupled to molecular mechanics-Poisson Boltzmann Surface Area (MM-PBSA) calculation. Molecular dynamics simulation was further analyzed by protein contact networks. Rough energetic evaluation with protein-protein docking scores revealed that aflibercept/VEGFA complex was characterized by electrostatic stabilization, whereas ranibizumab and bevacizumab complexes were stabilized by Van der Waals (VdW) energy term; these results were confirmed by MM-PBSA. Comparison of MM-PBSA predicted energy terms with experimental binding parameters reported in literature indicated that the high association rate (Kon) of aflibercept to VEGFA was consistent with high stabilizing electrostatic energy. On the other hand, the relatively low experimental dissociation rate (Koff) of ranibizumab may be attributed to lower conformational fluctuations of the ranibizumab/VEGFA complex, higher number of contacts and hydrogen bonds in comparison to bevacizumab and aflibercept. Thus, the anti-angiogenic agents have been found to be considerably different both in terms of molecular interactions and stabilizing energy. Characterization of such features can improve the design of novel biological drugs potentially useful in clinical practice.

High-throughput RNA sequencing reveals structural differences of orthologous brain-expressed genes between western lowland gorillas and humans

The human brain and human cognitive abilities are strikingly different from those of other great apes despite relatively modest genome sequence divergence. However, little is presently known about the interspecies divergence in gene structure and transcription that might contribute to these phenotypic differences. To date, most comparative studies of gene structure in the brain have examined humans, chimpanzees, and macaque monkeys. To add to this body of knowledge, we analyze here the brain transcriptome of the western lowland gorilla (Gorilla gorilla gorilla), an African great ape species that is phylogenetically closely related to humans, but with a brain that is approximately one-third the size. Manual transcriptome curation from a sample of the planum temporale region of the neocortex revealed 12 protein-coding genes and one noncoding-RNA gene with exons in the gorilla unmatched by public transcriptome data from the orthologous human loci. These interspecies gene structure differences accounted for a total of 134 amino acids in proteins found in the gorilla that were absent from protein products of the orthologous human genes. Proteins varying in structure between human and gorilla were involved in immunity and energy metabolism, suggesting their relevance to phenotypic differences. This gorilla neocortical transcriptome comprises an empirical, not homology- or prediction-driven, resource for orthologous gene comparisons between human and gorilla. These findings provide a unique repository of the sequences and structures of thousands of genes transcribed in the gorilla brain, pointing to candidate genes that may contribute to the traits distinguishing humans from other closely related great apes.

Molecular dynamics for computational proteomics of methylated histone H3

BACKGROUND

Post-translational modifications of histones, and in particular of their disordered N-terminal tails, play a major role in epigenetic regulation. The identification of proteins and proteic domains that specifically bind modified histones is therefore of paramount importance to understand the molecular mechanisms of epigenetics.

METHODS

We performed an energetic analysis using the MM/PBSA method in order to study known complexes between methylated histone H3 and effector domains of the PHD family. We then developed a simple molecular dynamics based predictive model based on our analysis.

RESULTS

We present a thorough validation of our procedure, followed by the computational predictions of new PHD domains specific for binding histone H3 methylated on lysine 4 (K4).

CONCLUSIONS

PHD domains recognize methylated K4 on histone H3 in the context of a linear interaction motif (LIM) formed by the first four amino acids of histone H3 as opposed to recognition of a single methylated site. PHD domains with different sequences find chemically equivalent solutions for stabilizing the histone LIM and these can be identified from energetic analysis. This analysis, in turn, allows for the identification of new PHD domains that bind methylated H3K4 using information that cannot be retrieved from sequence comparison alone.

GENERAL SIGNIFICANCE

Molecular dynamics simulations can be used to devise computational proteomics protocols that are both easy to implement and interpret, and that yield reliable predictions that compare favorably to and complement experimental proteomics methods. This article is part of a Special Issue entitled Recent developments of molecular dynamics.

The oncogenic potential of Jumonji D2 (JMJD2/KDM4) histone demethylase overexpression

The Jumonji D2 proteins (JMJD2/KDM4) function to demethylate di- and trimethylated (me2/3) histone 3 lysine 9 (H3K9me2/3) and H3K36me3. Knockout mouse models for Kdm4b and Kdm4d have not resulted in gross abnormalities, while mouse models for Kdm4a and Kdm4c have not been reported. However, the KDM4 subfamily of demethylases are overexpressed in several tumor types. Overexpression of KDM4 proteins alters transcription and chromatin remodeling, driving cellular proliferation, anchorage-independent growth, invasion, and migration. Increased proliferation occurs through KDM4-mediated modification of cell cycle timing, as well as through increased numbers of replication forks. Recent evidence also suggests that KDM4C overexpression contributes to the maintenance of a pluripotent state. Together these data suggest that overexpression of KDM4 proteins induces numerous oncogenic effects.

Exploring PHD fingers and H3K4me0 interactions with molecular dynamics simulations and binding free energy calculations: AIRE-PHD1, a comparative study

PHD fingers represent one of the largest families of epigenetic readers capable of decoding post-translationally modified or unmodified histone H3 tails. Because of their direct involvement in human pathologies they are increasingly considered as a potential therapeutic target. Several PHD/histone-peptide structures have been determined, however relatively little information is available on their dynamics. Studies aiming to characterize the dynamic and energetic determinants driving histone peptide recognition by epigenetic readers would strongly benefit from computational studies. Herein we focus on the dynamic and energetic characterization of the PHD finger subclass specialized in the recognition of histone H3 peptides unmodified in position K4 (H3K4me0). As a case study we focused on the first PHD finger of autoimmune regulator protein (AIRE-PHD1) in complex with H3K4me0. PCA analysis of the covariance matrix of free AIRE-PHD1 highlights the presence of a "flapping" movement, which is blocked in an open conformation upon binding to H3K4me0. Moreover, binding free energy calculations obtained through Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) methodology are in good qualitative agreement with experiments and allow dissection of the energetic terms associated with native and alanine mutants of AIRE-PHD1/H3K4me0 complexes. MM/PBSA calculations have also been applied to the energetic analysis of other PHD fingers recognizing H3K4me0. In this case we observe excellent correlation between computed and experimental binding free energies. Overall calculations show that H3K4me0 recognition by PHD fingers relies on compensation of the electrostatic and polar solvation energy terms and is stabilized by non-polar interactions.

Epigenetic control and cancer: the potential of histone demethylases as therapeutic targets

The development of cancer involves an immense number of factors at the molecular level. These factors are associated principally with alterations in the epigenetic mechanisms that regulate gene expression profiles. Studying the effects of chromatin structure alterations, which are caused by the addition/removal of functional groups to specific histone residues, are of great interest as a promising way to identify markers for cancer diagnosis, classify the disease and determine its prognosis, and these markers could be potential targets for the treatment of this disease in its different forms. This manuscript presents the current point of view regarding members of the recently described family of proteins that exhibit histone demethylase activity; histone demethylases are genetic regulators that play a fundamental role in both the activation and repression of genes and whose expression has been observed to increase in many types of cancer. Some fundamental aspects of their association with the development of cancer and their relevance as potential targets for the development of new therapeutic strategies at the epigenetic level are discussed in the following manuscript.

Molecular mechanisms and potential functions of histone demethylases

Histone modifications are thought to regulate chromatin structure, transcription and other nuclear processes. Histone methylation was originally believed to be an irreversible modification that could only be removed by histone eviction or by dilution during DNA replication. However, the isolation of two families of enzymes that can demethylate histones has changed this notion. The biochemical activities of these histone demethylases towards specific Lys residues on histones, and in some cases non-histone substrates, have highlighted their importance in developmental control, cell-fate decisions and disease. Their ability to be regulated through protein-targeting complexes and post-translational modifications is also beginning to shed light on how they provide dynamic control during transcription.

Comparative molecular dynamics simulation study of crystal environment effect on protein structure

Crystal structures of proteins are under the influence from the crystal environment. In this study, we used molecular dynamics (MD) simulations to explore the possibility of eliminating the effect of the crystal packing and recovering the structure in solution. Ten representative proteins were chosen from the Protein Structural Change Database as the target systems, and 50 ns MD stimulations starting from two crystal structures having different domain arrangements were performed for each. The MD trajectories of the relaxation processes upon the release from the crystal environment revealed that the behaviors of the proteins were classified into three groups: "single domain linker", "harmonic motion", and "large barrier". We discuss the structural features common to the proteins in each group.

The Histone Database: an integrated resource for histones and histone fold-containing proteins

Eukaryotic chromatin is composed of DNA and protein components-core histones-that act to compactly pack the DNA into nucleosomes, the fundamental building blocks of chromatin. These nucleosomes are connected to adjacent nucleosomes by linker histones. Nucleosomes are highly dynamic and, through various core histone post-translational modifications and incorporation of diverse histone variants, can serve as epigenetic marks to control processes such as gene expression and recombination. The Histone Sequence Database is a curated collection of sequences and structures of histones and non-histone proteins containing histone folds, assembled from major public databases. Here, we report a substantial increase in the number of sequences and taxonomic coverage for histone and histone fold-containing proteins available in the database. Additionally, the database now contains an expanded dataset that includes archaeal histone sequences. The database also provides comprehensive multiple sequence alignments for each of the four core histones (H2A, H2B, H3 and H4), the linker histones (H1/H5) and the archaeal histones. The database also includes current information on solved histone fold-containing structures. The Histone Sequence Database is an inclusive resource for the analysis of chromatin structure and function focused on histones and histone fold-containing proteins.

Structural cooperativity in histone H3 tail modifications

Post-translational modifications of histone H3 tails have crucial roles in regulation of cellular processes. There is cross-regulation between the modifications of K4, K9, and K14 residues. The modifications on these residues drastically promote or inhibit each other. In this work, we studied the structural changes of the histone H3 tail originating from the three most important modifications; tri-methylation of K4 and K9, and acetylation of K14. We performed extensive molecular dynamics simulations of four types of H3 tails: (i) the unmodified H3 tail having no chemical modification on the residues, (ii) the tri-methylated lysine 4 and lysine 9 H3 tail (K4me3K9me3), (iii) the tri-methylated lysine 4 and acetylated lysine 14 H3 tail (K4me3K14ace), and (iv) tri-methylated lysine 9 and acetylated lysine 14 H3 tail (K9me3K14ace). Here, we report the effects of K4, K9, and K14 modifications on the backbone torsion angles and relate these changes to the recognition and binding of histone modifying enzymes. According to the Ramachandran plot analysis; (i) the dihedral angles of K4 residue are significantly affected by the addition of three methyl groups on this residue regardless of the second modification, (ii) the dihedral angle values of K9 residue are similarly altered majorly by the tri-methylation of K4 residue, (iii) different combinations of modifications (tri-methylation of K4 and K9, and acetylation of K14) have different influences on phi and psi values of K14 residue. Finally, we discuss the consequences of these results on the binding modes and specificity of the histone modifying enzymes such as DIM-5, GCN5, and JMJD2A.

A comparative molecular dynamics study of methylation state specificity of JMJD2A

Histone modifications have great importance in epigenetic regulation. JMJD2A is a histone demethylase which is selective for di- and trimethyl forms of residues Lys9 and Lys36 of Histone 3 tail (H3K9 and H3K36). We present a molecular dynamics simulations of mono-, di- and trimethylated histone tails in complex with JMJD2A catalytic domain to gain insight into how JMJD2A discriminates between the methylation states of H3K9. The methyl groups are located at specific distances and orientations with respect to Fe(II) in methylammonium binding pocket. For the trimethyllysine the mechanism which provides the effectual orientation of methyl groups is the symmetry, whereas for the dimethyllysine case the determining factors are the interactions between methyllysine head and its environment and subsequently the restriction on angular motion. The occurrence frequency of methyl groups in a certain proximity of Fe(II) comes out as the explanation of the enzyme activity difference on di- and tri-methylated peptides. Energy analysis suggests that recognition is mostly driven by van der Waals and followed by Coulombic interactions in the enzyme-substrate interface. The number (mono, di or tri) and orientations of methyl groups and water molecules significantly affect the extent of van der Waals interaction strengths. Hydrogen bonding analysis suggests that the interaction between JMJD2A and its substrates mainly comes from main chain-side chain interactions. Binding free energy analysis points out Arg8 as an important residue forming an intra-substrate hydrogen bond with tri and dimethylated Lys9 of the H3 chain. Our study provides new insights into how JMJD2A discriminates between its substrates from both a structural and dynamical point of view.