Literature DB >> 30344952

Prominent role of histone lysine demethylases in cancer epigenetics and therapy.

Avilala Janardhan1,2, Chandrasekhar Kathera2, Amrutha Darsi2, Wajid Ali2, Lingfeng He2, Yanhua Yang1, Libo Luo1, Zhigang Guo2.   

Abstract

Protein methylation has an important role in the regulation of chromatin, gene expression and regulation. The protein methyl transferases are genetically altered in various human cancers. The enzymes that remove histone methylation have led to increased awareness of protein interactions as potential drug targets. Specifically, Lysine Specific Demethylases (LSD) removes methylated histone H3 lysine 4 (H3K4) and H3 lysine 9 (H3K9) through formaldehyde-generating oxidation. It has been reported that LSD1 and its downstream targets are involved in tumor-cell growth and metastasis. Functional studies of LSD1 indicate that it regulates activation and inhibition of gene transcription in the nucleus. Here we made a discussion about the summary of histone lysine demethylase and their functions in various human cancers.

Entities:  

Keywords:  epigenetics; histone modifications; lysine specific demethylases; protein methylation; therapeutics

Year:  2018        PMID: 30344952      PMCID: PMC6188137          DOI: 10.18632/oncotarget.24319

Source DB:  PubMed          Journal:  Oncotarget        ISSN: 1949-2553


INTRODUCTION

Epigenetics refers to heritable changes regulating gene expression that are not a result of changes in the primary DNA sequence. In cancer, aberrant epigenetic silencing of tumour suppressor genes is a common occurrence that is associated with abnormal DNA methylation patterns and changes in covalent histone modifications. Histone proteins are subjected to post-translational modifications, like acetylation, methylation, phosphorylation, and ubiquitination, and these histone modifications act as the molecular switches that alter the state of compaction of chromatin to allow gene activation or repression [1-3]. Some histone modifications like acetylation and phosphorylation are highly dynamic, whereas the methylation have been regarded as ‘‘permanent’’ chromatin marks. Methylation of lysine (K) residues at the ε-amino group are one of the most common post-translational modifications of histone N-terminal tails [4, 5]. Lysine residues can be mono-, di- and tri-methylated. Histone lysine methylation plays critical roles in virtually all chromatin-templated biological processes, such as transcription and DNA repair [5, 6]. Cumulative evidence indicates that histone lysine methylation exerts diverse biological functions in a site- and methylation state-dependent manner [5, 6]. Like other histone modifications, histone lysine methylation is reversible and dynamically regulated through the balancing activities of histone lysine methyltransferases and demethylases [7-9]. To date, two families of proteins have been found to possess lysine demethylase activity. Members of one family are amine oxidases that can catalyze demethylation of mono- and di-methylated lysine residues using flavin adenine dinucleotide (FAD) as a cofactor [10]. Members of the other family are the JmjC family proteins that catalyze oxidative demethylation reactions with iron and α-ketoglutarate as cofactors [11-13]. In this review we made a generalized description about the histone lysine demethylases- which includes the two known lysine demethylase families, their role in human tumorigenesis and the potential therapeutic targets of lysine demethylases.

Historical development of protein methylation

The research on ‘protein methylation’ were first started in the 1960s. In 1962, first Huang and Bonner's [14] investigations stated that the histones might function as gene regulators. Later in 1964 Allfrey et al. [15]. observed that the amino acids of histone might be acetylated or methylated. In 1967, Paik et al. [16]. and in 1968 Hempel et al. [17]. reported the presence of N-dimethyl-lysine in histone hydrolysates in addition to the ε-N-monomethyl-lysine. In 1968 it was the reason for the identification and purification of protein methylase I [18]. The discovery of the arginine methyltransferase was by Baldwin and Carnegie [19], and Brostoff and Eylar [20], independently and they finally reported that Arg107 is methylated at the guanidino group. In 1967 Liss and Edelstein [21] identified another enzyme called protein methylase II, that is capable of esterifying the dicarboxylic amino acid residues of proteins. It was the second enzyme discovered that transfers a methyl group to a protein side chain. In 1970 protein lysine methyltransferase, which methylates the ε-amino group of lysine residues in a protein from calf thymus was characterized and Identified as protein methylase III [22, 23]. Now, these enzymes are known as a family of lysine specific histone methyltransferases [3]. In the early 1980s, it was known that the specific enzymes methyltransferases were responsible for methylation of lysine, arginine, histidine and dicarboxylic amino acids. In the 1990s, a vast research in protein methylation occurred and now it is clear that protein methylation carries many important biological functions, including gene regulation and signal transduction. The protein methylase I have two subtypes: (i) histone or heterogenous nuclear ribonucleoprotein (hnRNP)-specific and (ii) myelin basic protein (MBP)-specific. In the mid-1990s the knowledge of protein- arginine methylation and its applications were vastly spread. It is now recognized that Protein methyl transferase (PRMT) exists as a family with several subtypes. The PRMT family has been shown to include at least nine methyltransferases, designated as PRMT 1–9 based on differences in primary sequence and substrate specificity [24-28].

Types of protein-methylation

The protein-methylation was by specific enzymes called protein methyl transferases. These protein methyl transferases can be classified into two major groups: The first group modifies carboxyl groups to form methyl esters. These are reversible reactions and can regulate the activities of the methylatable protein. The second group of protein-methylation reactions is irreversible methyl transfer to sulfur and nitrogen atoms. The function of these reactions is much less clear but appears to be the generation of a variety of new amino acids for specialized cellular roles.

Glutamate methylation in bacteria

Glutamate methylation was observed in bacteria but not in eukaryotic cells, particularly the chemotactic bacteria glutamate methylation was common but Robert Sprung first time reported glutamate methylation in bacteria [29]. The enzyme recognizes the sequence (S/A)-(S/A)-X-X-(F/Q)-(F/Q)-X-A-A (here X represents a variety of amino acids) and then methylate's a glutamate residue in the sixth position [30, 31]. The enzyme responsible for the demethylation reaction was the CheB methylesterase, can also recognize glutamine residues and catalyzes their deamidation as well as the demethylation of glutamate γ-methyl esters. Both the glutamate methyl ester and glutamine residues are uncharged at neutral pH, while the glutamate residue has a net negative charge. Thus, the methylation/demethylation system can convert receptors with glutamine/γ-glutamate methyl ester residues with one type of signaling function to receptors with glutamate residues and the second type of signaling function [32]. Initially, these bacterial enzymes were the examples of a family of glutamate methyltransferases/methylesterases that could reversibly modify other bacterial and animal cell proteins and regulate their activity. However, glutamate methylation has not found to be utilized in eukaryotic cells, or even in non-chemotactic systems in bacteria.

Cysteine and leucine methylation

Cysteine and leucine methylation was observed in eukaryotic cells but not in prokaryotic cells. It was observed in carboxyl-terminal of isoprenyl cysteine and leucine residues. This type of methylation may control the assembly and disassembly of nuclear lamins [33, 34]. Evidence has also been presented for a role of this type of methylation in signal transduction reactions involving G proteins. Inhibitor studies have suggested more specific roles of the methylation reaction in platelet stimulation leading to their aggregation [35] and in the response of macrophages [36], neutrophils [37] and HL-60 granulocytes [38] to chemotactic stimulation.

Histidine methylation

Methylation Of histidine residues has been described in actin, myosin, histones and rhodopsin [39]. The role of actin methylation has been tested by creating actin variants that contain an arginine or tyrosine residue in place of histidine-73 [40]. Several properties of these variants, including their ability to become incorporated into the cytoskeleton of a non-muscle cell, were found to be identical to those of normally methylated actin. However, actin interacts with a very large number of proteins and the conserved methylhistidine residue may be crucial for one or more of these interactions.

Asparagine methylation

The side-chain amide nitrogen of the asparagine-72 residue of the g-subunit of the C- and R-phycocyanins and the B-, C- and R-phycoerythrins (photosynthetic accessory proteins present in cyanobacteria and red algae) is methylated. On the C-phycocyanin, this methylated residue is located near one of the chromatophoric groups, and evidence has been presented that methylation serves to enhance the efficiency of energy transfer in the light-harvesting complexes [41, 42].

Arginine methylation

Arginine methylation has been identified down to the earliest branches of eukaryotes [43, 44]. It can be either activatory or repressive for transcription. Methylation of arginine is regulated by arginine methyltransferases (PRMTs), but there are no enzymes yet identified that can reverse arginine methylation [45]. Arginine methylation was initially thought to be a rather permanent modification [46], but there are several analyses supporting a reversible nature of methyl arginine [47]. The estrogen receptor α (ERα) is methylated in an RGG motif at R260 by PRMT1 [48]. Methylation of arginine is catalyzed by a group of nine protein arginine methyltransferases [49]. The PRMT family of enzymes is well conserved within multicellular organisms ranging from cnidarians to humans [50]. Methylated arginine residues have been found in both histone and non-histone proteins as myelin basic protein, myosin, heat shock proteins, nuclear proteins, ribosomal proteins, and tooth matrix proteins. These methylated arginine residues involved in functions such as RNA processing, DNA repair, and transcription [51]. In histones, the best-characterized methylated-arginine residues include R2, R8, R17, and R26 of histone H3 and R3 in histone H4 and H2A [52]. Due to the methylation state of individual arginine residues, the adjacent chromatin region is either transcriptional active or repressed. Asymmetric di-methylation of H3R2, for instance, has been linked to inactive chromatin regions. In contrast, mono-methylation or symmetric di-methylation of H3R2 is associated with an active chromatin state [53, 54].

Lysine methylation

The discovery of LSD1, another protein family with more than 30 histone demethylases structurally different from LSD1 was described, all of which share a motif designated the Jumonji C (JmjC) domain and revealing substrate specificity. Identification of these enzymes opened a new era in understanding how chromatin dynamic is regulated and further understanding of the regulation of these enzymes will provide significant insight into fundamental mechanisms of many biological processes and human diseases. Six lysine (K) residues on histone H3 and H4 (H3K4, H3K9, H3K27, H3K36, H3K79 and H4K20) are subjected to mono, di and tri methylation. Importantly each methylation state represents a specific epigenetic mark with a precise biological meaning and well defined chromatin localization [2]. H3K4, H3K36, and H3K79 are implicated in activation of transcription, whereas H3K9, H3K27 and H4K20 are connected to transcriptional repression.

Histone lysine demethylation

Histones may be methylated on either lysine (K) or arginine (R) residues. Methylation of histone does not affect the chromatin structure but they create a binding sites for regulatory proteins [55]. Arginine side chains may be mono or di-methylated as symmetrically or asymmetrical, whereas the lysine side chains can be mono, di, or tri-methylated. As there are number of distinct lysine and arginine residues present in the N-terminal tails of histones the methylation possibilities were more at this sites. Lysine methylation was limited to certain proteins that recognize the modifications. For instance, hetero chromatin protein HP1 binds to H3k9 where as chromodomain helicase DNA binding protein CHD1 recognizes H3K4 [2]. On the other side p53-binding protein 1 (p53BP1) recognizes methylated H3K79 [56], and the WDR5 a vertebrate transcriptional activator recognize both di and tri methylated H3K79 [57]. All methylated residues play a key role in transcription and DNA damage [57]. Methylation of lysine is by histone methyltransferases (HMTs) and as there are two classes of lysine methyl transferases have been identified, the proteins with amine oxidases LSD1 and LSD2. Both LSD1 and LSD2 use FAD-cofactor and releases demethylated lysine, FADH2 and formaldehyde. The JmjC-domain, another class of lysine demethylases uses iron and α-ketoglutarate as cofactors leads to release succinate, carbon dioxide, and formaldehyde along with demethylated lysine. The histone lysine methylation regulation landscape has been changed considerably in the recent years but the biological functions of lysine demethylases are uncovered. Lysine methylation was associated with different types of cancers. The mutations or overexpression of several lysine demethylases associated with cancer has been raising the possibility of therapeutic potential. Past evidence suggests that HDMs participate in nuclear receptor [58]. For instance the H3K4-Me associated with prostate cancer, H3K4-Me2 was associated with lung, kidney, prostate and breast cancers, H3k9-Me2 with gastro adenocarcinoma and H3K27-Me3 with breast and ovarian cancer [55]. Transcriptional activation by nuclear receptor is a multistep process and involves primary and secondary coactivators. The binding of hormone to the receptor induces a structural reorganization which allows the nuclear receptor to recognize specific DNA elements at promoter regions [59]. These changes will result in recruitment of coactivators, leads to intrinsic enzymatic activities, such as demethylation. At the same time these proteins may change the local chromatin structure and subsequent recruitment of RNA polymerase II transcriptional machinery [60]. The histone demethylases were the key players in the regulation of androgen receptor (AR) -mediated transcription. Both LSD1 and KDM4C interact with AR and induce demethylation of H3K9 at the prostate-specific antigen (PSA) promoter [61, 62]. Similarly, transcriptional activation of thyroid receptor leads to demethylation of H3K9 by a JmjC-domain containing family member [63]. Although studies associated with arginine methylation to estrogen receptor (ER) function, very little has been done in elucidating the role of lysine methylation in regulating ER-dependent transcription. Several H3K9 histone methyl transferases have been involved in repressing ER targets pS2 and GREB1 in the absence of estrogen [64]. List of targets that corresponds to the methylation/demethylation sites and their role in transcription was tabulated in Table 1.
Table 1

List of targets that corresponds to the methylation/demethylation sites and their role in transcription

Methylation site to be modified by KMT/KDMTargetsRole in transcription
H3k4MLL, ALL1,ALR1, ASH1, Set 1, Set 9/7, NSD1Activation
H3K9SUV39H, G9a, EU-HMTase/GLP, ESET/SETBD1, CLL*, Spclr4Activation/Repression
H3K27EZH2Repression
H3K36HYPB, Symd2, NSD1, Set2Elongation
H3K79Dot1Activation
H3K20PrSet7/8, SUV420H, SpSet9Silencing

Structure of histone lysine demethylases

Histone lysine demethylases are most frequent targets for epigenetic drugs and are gaining more interest in research. A lysine in histone proteins can be mono-, di-, and trimethylated. Every modification on the each/same amino acid can have different biological effects [65]. The recent studies on histone lysine demethylases reveal two types of enzyme mechanisms [66]. The iron-dependent enzymes, which can demethylate side chains of lysine in all three methylation states. With reference to this, the oxidative chemistry that involves in the function of flavin-dependent histone demethylases makes it impossible for these enzymes to act on trimethylated lysine and restricts their activity to mono- and dimethylated substrates [67]. The two main flavoenzyme demethylases are LSD1 and LSD2 also known as KDM1A and KDM1B [10, 68].

LSD1

The first FAD-dependent histone lysine demethylase family identified was LSD1 [10] and its homolog LSD2. The LSD1 structure consists of three domains, the amine oxidase domain, FAD binding domain and SWIRM domain. The FAD dependent amine oxidase domain consists of an active site to bind substrates. The FAD binding domain consists of a tower domain which interacts with CoREST and the SWIRM domain function in LSD1 was unknown but the SWIRM domain in other histone protein modifications binds to DNA as it is having solvent exposed positive charged patch. The tower domain has an insert of 102 amino acids which seperates N and C terimus of amine oxidase domain. The tower domian is linked with the corepressor protein CoREST is essential for demethylation and removes methyl groups from mono- and dimethyl Lys4 of histone H3, a gene activation mark [69]. These enzymes are often deregulated in human diseases, it is essential to understand their physiological and pathological/biological functions by elucidating their exact structure, regulation and substrate specificity. The enzyme is an attractive target for epigenetic drugs as it is overexpressed in solid tumors [70], its involvement in various infection, differentiation processes [71, 72].

LSD2

LSD2 is also like LSD1 which has specificity in removing methyl groups from mono and dimethyl Lys4 of H3. But, the biology of LSD2 is different from that of LSD1 as LSD2 does not bind CoREST [73, 74]. The catalytic domains of these LSD1 and LSD2 proteins shows 45% identity and is structurally homologous with the amine oxidases that act on biogenic amines [10, 67]. Among all these proteins, a vast research for more than 50 years was on human monoamine oxidases (MAOs) A and B which led to the development of a multitude of inhibitors including antidepressive and antiparkinson drugs [75]. The LSD1 has been vastly studied [10, 76] when compared to LSD2, the enzyme activity and biological function of LSD2 have become more prominent. Several studies [68, 77] reported that LSD2 particularly demethylates mono- and di-methylated H3K4, but van Essen et al. [78] reported that LSD2 also demethylates dimethylated H3K9. The LSD2 requires a long stretch of H3 tail sequence for demethylation of H3K4 [68]. Biochemical studies of LSD2 indicate that it is associated with transcriptional elongation factors and phosphorylated RNA polymerase II [77, 79] (Figure 1). LSD2 is enriched in the coding regions of actively transcribed genes and has a role in maintaining a repressive chromatin environment within gene bodies through its H3K4 demethylase activity [77]. The H3K4 demethylation by LSD2 indicates that it is linked with the establishment of maternal imprinting in mouse oocytes [73]. As in LSD1 the LSD2 also contains a N-terminal SWIRM domain and a C-terminal amine oxidase domain. But the LSD2 does not have tower domain which is a particular insert present in the amine oxidase domain of LSD1 that is required for LSD1 binding of CoREST and LSD1's enzymatic activity [69, 80, 81]. The LSD2 has a N-terminal region that is absent in LSD1, which forms a CW-type zinc-finger domain [74]. The CW domain is present in several chromatin-remodeling proteins and has capable to bind methylated histones [82, 83]. Deletion mutant experiments explained that the N-terminal zinc finger domain and the SWIRM domain are essential for LSD2 demethylase activity [74].
Figure 1

Regulation of transcriptional elongation: LSD2 associated with transcriptional elongation factors and phosphorylated RNA polymerase II

Crystal structure of LSD1

Crystal structure of LSD1 represents a number of transcriptional corepressor complexes like CoREST, CtBP, and HDAC complexes which plays a major role in silencing neuronal-specific genes in nonneuronal cells [84-89]. The C-terminal FAD dependent amine oxidase with a tower insert occupies major part in the LSD1 structure [90-92]. The N-terminal of LSD1 contains SWIRM domain, which was often found in chromatin remodeling and modifying complexes with unknown function [93]. The SWIRM and FAD dependent amine oxidase domains pack together through extensive forces and form a globular structure. The SWIRM domain contains mostly the α-helices (six-helical bundle structures characterized as a long helix and in the center it was surrounded by five other helices) and the amine oxidase domain exhibits topology with several Flavin dependent oxidases. In addition to this the SWIRM domain has two stranded β sheets at the C-terminal end which helps to anchor the interactions between the SWIRM and AOL domains. The amine oxidase domain has two sub domains, a FAD binding sub-domain and a substrate binding subdomain. The FAD binding sub-domain has a mixed structures with both α-helices and β-sheets. The substrate binding domain has six stranded mixed β sheet flanked by six α helices. The tower insert of amine oxidase domain consists of a pair of long helices that adopts a typical antiparallel coiled confirmation [94]. Demethylation of mono or di-methylated lysine on histone is catalyzed by LSD1. This enzyme has an amine oxidase containing flavin adenine dinucleotide (FAD) as an electron acceptor to first oxidize the lysine N-methyl amine to lysine N-methylimine. FADH2 is re-oxidized to FAD by molecular oxygen producing hydrogen peroxide. The N-methylimine is non-enzymatically hydrolyzed to a carbinolamine which spontaneously dissociates to a demethylated lysine and formaldehyde. If the substrate is a dimethylated lysine the enzyme performs sequential removal of both methyl groups by the same mechanism. LSD1 demethylates histone H3 as part of a multimeric protein complex in which it directly interacts with the co-repressor of the repressor element 1 silencing transcription factor, CoREST. Interaction of these two proteins has been studied in detail and the crystal structure of the complex containing an analog of the histone H3 peptide was determined and reported [81, 95].

LSD1 mechanism and function

As LSD1 uses FAD as cofactor, it is a member of the FAD-dependent amine oxidase class lysine methylases. Demethylation by LSD1 occurs by oxidizing the carbon-nitrogen bond between the methyl group of lysine residues, and produce demethylated lysine, FADH2 and formaldehyde. LSD1 can only demethylate mono- and dimethylated lysines, but not trimethylated lysines, as they need a lone pair of electrons on the amine group to from imine. Initially LSD1 was reported as gene repressor as it demethylate H3K4 and repress the transcription of neuron [85], but the LSD1 was also able to activate the gene as it demethylate H3K9 which leads to upregulate the expression of gene [61]. Garcia Barcet etal reported that 80% of the promoters occupied by LSD1 were bound to RNA polymerase II, suggesting that LSD1 was liked more in active genes rather than the inactive genes [64]. LSD1 has been involved in the maintenance of diseases such as neuroblastoma [70] and high-risk prostate cancer [96]. This enzyme plays important regulatory roles in cell differentiation [97] and it is also required for transcriptional repression of the hTERT gene in both normal and cancerous cells [98]. Other than histone substrates, LSD1 was shown to demethylate K370 of the tumor suppressor p53, repressing its function by interacting with 53BP1 (Figure 2), a DNA damage checkpoint [99]. The SNAG domain of Snail1 seems a H3-like structure and functions as a molecular hook for recruitment of LSD1 to repress gene expression in metastasis [100]. Further, LSD1 is highly expressed in breast cancer cells and it is an essential mediator of the inter chromosomal interactions necessary for estrogen-dependent transcription [101, 102]. Downstream target genes regulated by lysine-specifc demethylase 1 (LSD1) in colon cancer cells was investigated by Jiang Chen [103]. They report that total 3633 and 4642 differentially expressed genes were signifcantly upregulated and downregulated respectively in LSD1-silenced SW620 cells.
Figure 2

Regulation of p53: LSD1 demethylate lysine 370 of the tumor suppressor p53, repressing its function by interacting with 53BP1

Biochemical characteristics of lysine methylation

The protein lysine methyl transferases and protein lysine demethyal transferases have a key role in epigenetic regulation. The protein lysine methyl transferase methylates histone H3 at lysine 27 (polycomb repressive complex) is over expressed in different types of cancer [104]. Polycomb-repressive complexes has important role in controlling the expression of downstream genes and it also promotes cell cycle progression in physiological and pathological conditions [105]. Similar reports on protein lysine methyl transferases and protein lysine demethyl transferases promote malignant transformation via histone methylation dependent transcriptional regulation [106-110].

Lysine methyltransferases

Mono, di and trimethylation of the lysine ε-amino group are catalysed by lysine methyltransferases in an S-adenosyl methionine dependent manner [111]. Histone methyl transferases (HMT) consists of two main classes, the SET domain family and DOT1 family. Most lysine methyltransferases contain the SET domain which has N-methyl transferase activity, but some non-SET-domain proteins like DOT1L methyltransferase like 10 (METTL10) and METTL21A are also known to have lysine N-methyltransferase activity [112, 113]. DOT1 methyl transferases lacks SET domain and methylate lysine 79 with in the core domain of H3. DOT1 methyl transferases can only methylate nucleosomal substrates but not free histones. The SET-domain methyl transferases catalyze a sequential bi-bi kinetic mechanism in which both substrate association and product release occur in a random manner [111]. The methylation increases the hydrophobic and basic nature of the lysine residue, which allows other proteins to recognize methylated lysine.

Lysine demethylases

The methylation of lysine residues is by lysine demethylases and this is an irreversible reaction. Lysine demethylases are a group of proteins that are categorized into two functional enzymatic families. The first family includes LSD1 and LSD. LSD1 is a flavin-dependent monoamine oxidase and the LSD2 is the only homologue of LSD1 in the human genome [77]. These amine oxidases can only demethylate mono and dimethyl lysine residues, but not trimethyl lysine residues because they require a lone pair of electrons that are not present on trimethyl lysine residues. The second family of lysine demethylases is JmjC domain containing proteins which are metalloenzymes that can use an oxygenase mechanism to demethylate mono, dimethyl and trimethyl lysine residues.

Roles of KMTs and KDMs in human tumorigenesis

Lysine methylases and demethylases are the factors for development of many diseases, especially cancer. LSD1 is over expressed in various types of cancer like bladder and colorectal cancer, prostate cancer and oestrogen-receptor-negative breast cancer [114, 115]. Over expression of LSD1 was found to prevent the differentiation process and hold the malignant phenotype in neuroblastoma. The p53 activity was decreased when LSD1 was directly interacted with p53, which is associated with tumorigenesis [99]. In addition, lysine demethylases was not found to be altered in many cancers though it may functions as ubiquitin ligase. In the other side recurrent mutations in genes that encodes histone demethylases have been discovered in the gene encoding ubiquitously transcribed TPR protein on the X chromosome, where inactivating mutations are present in various types of tumors [116, 117], indicating the role of lysine methylases in human tumorigenesis. KDM6A catalyses the demethylation of H3K27me2 and H3K27me3, and a mutation in KDM6A is therefore to lead to an increase in H3K27 methylation, which is functionally equivalent to the enhanced EZH2 activity found in different human tumours [118, 119]. Recent studies on LSD1 activity have therapeutic potential in cancer, which includes the finding of LSD1 that is necessary for the maintenance of acute myeloid leukemia containing mixed lineage leukemia translocations [120] and inhibition of LSD1 can reactivate the all-trans-retinoic acid differentiation pathway in acute myeloid leukemia [121]. In addition to providing genetic evidence for a role of LSD1 in a mouse model of mixed lineage leukemia AF9 induced leukemia, Harris et al. [120]. showed that the MAO inhibitor tranylcypromine significantly inhibited the colony-forming ability of AML cells. The JMJD2 subfamily member's demethylases H3K9me3 and H3K9me2 as well as H3K36me3 and H3K36me2 which are overexpressed in cancer. The locus containing the JMJD2C (GASC1) gene is genomically amplified in squamous cell carcinoma, medulloblastoma and breast cancer [122, 123]. The exhaustion of JMJD2B slows the growth of various cancer cell lines [124-127]. The development of cancer in different mouse models was linked with the loss of H3K9 trimethylation. Here, loss of the two specific H3K9me3 histone methyltransferases SUV39H1 and SUV39H2 in double-knockout mice results in genomic instability and the growth of cancer [128]. The decrease in H3K9me3 levels can also increase carcinogenesis in a T cell lymphoma mouse model [129]. The over expression of H3K9me3 and H3K9me2 by JMJD2 demethylases would show same effects, and targets the enzyme activity of the JMJD2 subgroup of Jumonji proteins could therefore have therapeutic potential. as discussed above, proof of concept has been achieved for JMJD2B, and different studies have also shown that JMJD2C is required for growth of various cancer cell lines like prostate carcinoma [62], breast carcinoma [122], squamous cell carcinoma [130] and diffuse large B cell lymphoma [131]. The H3K4me3 and H3K4me2 demethylase by JARID1B binds to H3K4me3-enriched promoters and it is required for normal development [132, 133]. This protein is overexpressed in various cancers like prostate cancer, breast cancer and bladder cancer [134, 135]. It was reported that JARID1B is necessary for the growth of a breast cancer cell line [136] and for the continuous growth of melanoma [137]. In addition to JMJC domain a multifunctional protein KDM2B also catalyzes the demethylation of H3K36me2 and H3K36me1 [138]. The catalytic demethylase activity of KDM2B is required for the proliferation of acute myeloid leukemia and pancreatic ductal adenocarcinoma, and it also found to be overexpressed [139, 140]. The previous research explains that KDM2B regulates the expression of polycomb target genes [138, 141] suggesting that KDM2B might contribute to tumorigenesis through the regulation of these genes.

Functions of protein lysine demethylases

Protein lysine methylation appears to play a key role in gene repression, gene activation and in gene differentiation process. Several posttranslational modifications of histones on specific amino acids like altering chromatin structure or by recruiting enzyme regulatory complexes or the combined effects determine the outcome of transcription, either activating or suppressing gene expression [45, 90]. The general functions of protein lysine demethylases in are explained below.

Gene repression

The LSD1-CoREST-HDAC core is involved in various biological process like in like hematopoiesis where it interacts with growth factor independent1 transcription repressor and repress the target genes. It also involved in silencing mature B-cell genes [142] by interacting directly with the transcriptional repressor B lymphocyte-induced maturation protein-1 (Blimp-1). This LSD1-CoREST-HDAC also interacts with constitutive transrepressor TLX and forms a complex to repress PTEN gene and inhibiting cell proliferation [143]. LSD1 can also interact directly with p53 to show p53-mediated transcriptional repression and repress the alpha-fetoprotein. The well-known p53 target gene p21 can transcribe actively without LSD1.

Gene activation

The activation of androgen and estrogen receptors requires LSD1 dependent histone H3K9 demethylation [61, 64]. AR and LSD1 colocalize on promoters and stimulate H3K9 demethylation without altering the H3K4 methylation status and promote ligand dependent transcription of AR target genes resulting in enhanced tumor cell growth (Figure 3). LSD1 knock-down decrease the activation of AR-responsive promoters. Genome-wide analysis of LSD1 promoter occupancy following estrogen treatment of MCF7 cells has showed better results regarding the activatory role of LSD1 [64]. The activation markers like H3K4-me2 and acetyl-H3K9 suggesting that LSD1 is involved in gene activation rather than repression. The two roles of LSD1 in gene repression and activation is also explained by the fine regulation of growth hormone expression during pituitary development [144].
Figure 3

Coactivator function: Androgen and LSD1 colocalize on promoters and stimulate H3K9 demethylation without altering the H3K4 methylation status and promote ligand dependent transcription of AR target genes resulting in enhanced tumor cell growth

Development and differentation

The functional role of LSD1 has been actively investigated and it appears to be key in development and differentiation. Knockout of LSD1results in mouse embryonic lethality at or before embryonic day 5.5 [97, 145]. The zygotic LSD1 expression initially appeared at the morular stage and became wide in postimplantation embryos. ES cells derived from LSD1 knockdown mouse showed severe growth impairment, it may be due to increased cell death, impaired cell cycle progression, and defects in differentiation. Conditional knockout of LSD1 showed defects in pituitary gland development and the Notch signalling pathway [97]. The interference RNA inhibition of LSD1 in several mammalian haematopoietic lineages resulted in impairment of differentiation in vitro [146].

Functions of lysine methylation

Effect on other protein modifications

According to the previous reports the protein lysine methylation affects other post translational modifications by directly through competitive inhibition or indirectly. The lysine methylase EZH2 methylates histone H2B at lysine 120, which is a site of ubiquitylation, inhibits ubiquitylation and suppresses the transcription of genes involved in cancer [147]. Besides to this, lysine methylation also affects phosphorylation and acetylation of neighbouring distinct amino acids. So, methylated lysine residues change the affinity of enzymes such as kinases or phosphatases for their substrates, which may alter post translational modifications at other sites on the other substrates [148, 149].

Interacting protein partners

Lysine methylation can regulate interactions with other proteins. Biochemically, methylation of lysine residues doesn't change or neutralize its charge and has a small change in size, but it changes the hydration energies and hydrogen bonding potential of side chains in lysine residues. There are some methyl lysine binding proteins which have a specific motif like chromodomain for recognizing the methylation states of lysine, called the ‘aromatic cage’ which contains a collection of aromatic proteins, often accompanied by one or more neighbouring anionic residues [150]. The combination of favorable cation-π, electrostatic, and van der Waals interactions, as well as size matching, gives these proteins a high degree of specificity for the methylation state [151]. By contrast, the PP2A phosphatase complex (serine/threonine protein phosphatase 2A), a key negative regulator of the MAP Kinase pathway, binds to MAP3K2 and this interaction is blocked by methylation, implying that methyl lysine can also block protein-protein interactions [149]. All the protein interactions can influence the functions lysine demethylases and its involvement in DNA methylation. LSD1 interacting partners and their functions are represented in the Table 2. A picture was also generated by STRING software (Figure 4).
Table 2

List of interacting proteins with lysine specific demethylases with their molecular weights and related pathways

KDM1A interacting protein nameMolecular weight (Kda)Function/related pathway
MTA181Activated PKN1 stimulates transcription of AR (androgen receptor) regulated genes KLK2 and KLK3 and Validated targets of C-MYC transcriptional activation.
RBBP747.8Activated PKN1 stimulates transcription of AR (androgen receptor) regulated genes KLK2 and KLK3 and Cellular Senescence
RBBP447.6
CTBP149Notch signaling pathway (KEGG) and Signaling by Wnt.
HDAC155Activated PKN1 stimulates transcription of AR (androgen receptor) regulated genes KLK2 and KLK3 and PEDF Induced Signaling
HDAC255
XRCC582DNA Double-Strand Break Repair and DNA damage_NHEJ mechanisms of DSBs repair
XRCC669.8
PRKDC461
RHOA22Bisphosphonate Pathway, Pharmacodynamics and Development VEGF signaling via VEGFR2 - generic cascades.
SIN3B133PEDF Induced Signaling and Development NOTCH1-mediated pathway for NF-KB activity modulation
KLK328.7Adhesion and Transcription Androgen Receptor nuclear signaling.
HIST2H2BE14.2DNA Double-Strand Break Repair and Activated PKN1 stimulates transcription of AR (androgen receptor) regulated genes KLK2 and KLK3.
HIST1H2BK13.8
HIST1H2BJ13.9
HMG20B33.9Activated PKN1 stimulates transcription of AR (androgen receptor) regulated genes KLK2 and KLK3 and Factors involved in megakaryocyte development and platelet production.
NCOA2159CCR5 Pathway in Macrophages and Circadian rythm related genes.
PHF21A74.8Factors involved in megakaryocyte development and platelet production and Response to elevated platelet cytosolic Ca2+.
AR110CCR5 Pathway in Macrophages and Transcription_Role of VDR in regulation of genes involved in osteoporosis.
RCOR174.8fMLP Pathway and Respiratory electron transport, ATP synthesis by chemiosmotic coupling, and heat production by uncoupling proteins.
PKN1103.9Signaling mediated by p38-gamma and p38-delta and Translation Translation regulation by Alpha-1 adrenergic receptors.
PRKDC469DNA Double-Strand Break Repair and DNA damage_NHEJ mechanisms of DSBs repair.
RAD50146
MRE11A81DNA Double-Strand Break Repair and Resolution of D-loop Structures through Synthesis-Dependent Strand Annealing (SDSA).
KDM1B interacting proteinsMolecular weight (Kda)Function/Related Pathway
ASXL1165.4Deubiquitination and Metabolism of proteins
ASXL298.3
MBD5160
MBD6101
FOXk278
HCFC1208Mitochondrial Gene Expression
BAP191DNA Double-Strand Break Repair and Deubiquitination.
HIST2H2AA314Activated PKN1 stimulates transcription of AR (androgen receptor) regulated genes KLK2 and KLK3 and Meiosis.
HIST2H2AC14
HIST2H3A16.9
HIST3H2BB13.9DNA Double-Strand Break Repair and Activated PKN1 stimulates transcription of AR (androgen receptor) regulated genes KLK2 and KLK3.
HIST1H2BA14.16
HIST1H2BN13.9
HIST2H2BF13.9
HIST2H2BE13.9
HIST1H2BD13.9
HIST1H2BB13.8
HIST1H2BO13.9
HIST1H2BK13.9
HIST1H2BJ13.9
Figure 4

Interacting proteins map for LSD1: A group of proteins that interact with LSD1, picture was generated by STRING software

Protein stability

The lysine methylation block the action of ubiquitin ligase, preventing proteins from degradation via the ubiquitin/proteasome system and hence the lysine methylation may increase the stability of proteins by preventing polyubiquitylation. In Saccharomyces cerevisiae, lysine methylated proteins show a significant longer half-life than proteins which no methylation was found. Furthermore 43% of methylated lysine sites were predicted to be amenable to ubiquitylation, suggesting that methylated lysine residues might block the action of ubiquitin ligase [152]. However, methylation-dependent ubiquitination is carried out by damage-specific DNA binding protein 1 (DDB1)/cullin4 (CUL4) E3 ubiquitin ligase complex and a DDB1-CUL4-associated factor 1 (DCAF1) adaptor, which recognizes monomethylated lysine residues and promotes polyubiquitylation of the other lysine residues on substrates such as retinoic acid related orphan nuclear receptor α (RORα) [153], stating that lysine methylation may also destabilize target proteins through regulation of distant polyubiquitylation.

Subcellular localization

A nuclear localization signal comprises of some short sequences of positively charged lysine or arginine residues exposed on the protein surface [154]. As both lysine and arginine residues are critical for the nuclear localization of proteins, one could mediate that methylation of lysine or arginine residues may influence subcellular localization. Infact, some lysine methylated proteins like heat shock protein 70 and p53 are commonly localized in the nucleus, whereas the unmodified versions of these same proteins are localized in both the cytoplasm and the nucleus [155, 156].

Promoter binding

The methylation of lysine regulates the binding affinity of transcription factors for promoters, that changes the transcription levels of target genes [157]. For instance, lysine methylation of p53 by the SETD7 and the nuclear factor κB by the methyl transferase nuclear receptor-binding SET domain-containing protein 1 distinctly increase their ability to bind the promoter and activation of downstream genes [156, 158]. In a structural modelling analysis of RELA(a subunit of NF-κB) with DNA complex, two methylated lysine's on RELA were interact with DNA through hydrophobic contacts [159]. This report indicates that increased hydrophobicity of lysine residues by methylation can increase the promoter binding ability of transcription factors.

Lysine methylation as a therapeutic target

It is clear that the existing evidence showed that the inhibition of lysine methylase is key role for the development of therapeutic agents. The clinical use of lysine methylase inhibitors gives a better explanation of antitumor effects particularly in animal models of cancer. The fungal metabolite Chaetocin, the first described lysine methylase inhibitor was specific for the methyltransferase SU(VAR)3-9 both in vitro and in vivo and therefore it was used to study heterochromatin-mediated gene repression [160]. Later, in 2011 Daigle et al. [161]. reported the clear evidence of PKMT inhibition by DOT1L inhibitor could lead to selective cancer cell killing and resultant in-vivo efficacy in an animal model of cancer. After that, Jo et al. [162]. showed that DOT1L has a critical role in hematopoiesis using a postnatal conditional knockout technique. Another inhibitor BIX-01294, which is an inhibitor of the PKMT EHMT2 [163], was shown to effectively suppress the growth of cancer cell lines [108]. Later two anticancer drugs were reported that target EZH2, both inhibit EZH2 carrying a cancer-specific mutation. The first one is GSK126 that targets Y641 and A677 mutated EZH2 for the treatment of DLBCL [164]. GSK126 is a potent, highly selective, S-adenosyl-methionine-competitive, small-molecule inhibitor of EZH2 methyltransferase activity, which effectively inhibits the proliferation of EZH2 mutant DLBCL cell lines and markedly inhibits the growth of EZH2 mutant DLBCL xenografts in mice [164]. The second anticancer drug is EPZ005687, it also targets Y641 and A677 mutated EZH2 found in nonHodgkin lymphoma where it reduces H3K27 methylation in lymphoma cells [165]. Recently, Knutson, S. K. et al. [166]. reported that the selective inhibition of EZH2 by EPZ-6438 leads to potent antitumor activity in EZH2 mutant non-Hodgkin lymphoma. EPZ5676, a potent and selective aminonucleoside inhibitor of the PKMT DOT1L [167, 168], has also been studied for the treatment of patients with acute leukaemia in which the MLLgene has undergone rearrangement or tandem duplication.

LSD1 inhibitors

As the importance of LSD function and their role in different diseases was widely spreaded the recent research was focused on LSD inhibitiors. The inhibitors targeting LSD1 have been actively studied for biological and biochemical charaterizations which help for drug development. Initially the catalytic domains of the LSD show homology with MAOA and MAOB. Some MAO inhibitors including tranylcypromine and phenelzine were found to inhibit nucleosomal demethylation of H3K4 [169]. The phenelzine was more potent against LSD1, whereas the tranylcypromine acts as an inhibitor in cellular environments [169, 170]. The tranylcypromine inhibitor forms a covalent adduct with FAD factor in amine oxidase domain [171]. The wide spread knowledge in H3K4 methylation were observed in embryonal carcinoma cells after treatment with tranylcypromine [169]. Thus the derivatives of tranylcypromine with enhanced potency and target selectivity for LSD were identified and shown to induce differentiation of promyelocytic leukemia cells and slow the growth of prostate cancer cell lines [170, 172]. Pargyline a clinicially useful molecule inhibits MAO B, has also been reported to inhibit LSD1 at H3K9 in the presence of the androgen receptor [61]. The propargylamine or aziridine derived from pargyline was proposed to inhibit LSD1 activity [173]. The aziridine peptide exhibited as a reversible competitive inhibitor, but the propargylamine peptide functions as time-dependent inactivation of the LSD1 by modification of the FAD. Other polyamine derivatives includes bisguanidine and biguanide have been shown to be effective noncompetitive LSD1 inhibitors and treatment results in re-expression of silenced tumor suppressor genes in colon cancer [174, 175]. Recently these polyamine derivatives have been shown to modulate gene expression in breast cancer cells [176].

JmjC demethylases inhibitors

In addition to the LSD family proteins, JmjC domain containing proteins also have lysine demethylase activity. JMJC demethylases share some similarities with HDACs, and their substrates are also present a lysine side chain with its N-terminus in the proximity of a divalent metal ion [177]. HDAC inhibitors are poor inhibitors of the JMJD2 family of demethylases but few have sufficient potency and target selectivity for LSD were selected and considered for drug development [178]. Kruidenier et al. [179]. reported a selective JmjC H3K27 demethylase inhibitor which reduces lipopolysaccharide-induced proinflammatory cytokine production by human primary macrophages. The inhibition of JmjC demethylase proteins may have broad therapeutic application [179], several JmjC containing demethylase have been reported as candidates for anticancer therapy. Inhibitors targeting the Jumonji type demethylase activity showed antitumor effects even though no drugs have yet been evaluated in a clinical trial [180]. Some potential lead compounds specifically GSK-J1 is an inhibitor of the JMJD3 subfamily [179].

CONCLUSIONS

The present review was focused on the histone lysine demethylases and their role in carcinogenesis. In the recent years, research has been focused in the involvement of histone lysine demethylases in epigenetic regulation of cell differentiation and in tumor growth, specifically for the LSD1 and Jmjc domain containing families. The structure and functions for these families were finely studied to display the enzyme catalytic mechanisms and to develop the drug targets for the tumor growth. But the role/functions of interacting proteins in tumor development of these families mechanisms are unclear. So, the future research has to focus on new manner to find the exact mechanism of LSD when it interacts with other proteins. It's an alternative way to develop new therapeutics. It is also suggested that the finding of new histone-regulating enzymes related to histone demethylases will be well beneficial to the field.
  178 in total

Review 1.  Biochemistry and regulation of the protein arginine methyltransferases (PRMTs).

Authors:  Yalemi Morales; Tamar Cáceres; Kyle May; Joan M Hevel
Journal:  Arch Biochem Biophys       Date:  2015-12-02       Impact factor: 4.013

2.  PLU-1 is an H3K4 demethylase involved in transcriptional repression and breast cancer cell proliferation.

Authors:  Kenichi Yamane; Keisuke Tateishi; Robert J Klose; Jia Fang; Laura A Fabrizio; Hediye Erdjument-Bromage; Joyce Taylor-Papadimitriou; Paul Tempst; Yi Zhang
Journal:  Mol Cell       Date:  2007-03-15       Impact factor: 17.970

3.  Transcription regulation through promoter-proximal pausing of RNA polymerase II.

Authors:  Leighton J Core; John T Lis
Journal:  Science       Date:  2008-03-28       Impact factor: 47.728

4.  Phycobiliprotein methylation. Effect of the gamma-N-methylasparagine residue on energy transfer in phycocyanin and the phycobilisome.

Authors:  R V Swanson; A N Glazer
Journal:  J Mol Biol       Date:  1990-08-05       Impact factor: 5.469

5.  Protein methylation.

Authors:  W K Paik; S Kim
Journal:  Science       Date:  1971-10-08       Impact factor: 47.728

6.  The histone demethylase JMJD2B is regulated by estrogen receptor alpha and hypoxia, and is a key mediator of estrogen induced growth.

Authors:  Jun Yang; Adrian M Jubb; Luke Pike; Francesca M Buffa; Helen Turley; Dilair Baban; Russell Leek; Kevin C Gatter; Jiannis Ragoussis; Adrian L Harris
Journal:  Cancer Res       Date:  2010-08-03       Impact factor: 12.701

7.  Androgen receptor coactivators lysine-specific histone demethylase 1 and four and a half LIM domain protein 2 predict risk of prostate cancer recurrence.

Authors:  Philip Kahl; Lucia Gullotti; Lukas Carl Heukamp; Susanne Wolf; Nicolaus Friedrichs; Roland Vorreuther; Gerold Solleder; Patrick J Bastian; Jörg Ellinger; Eric Metzger; Roland Schüle; Reinhard Buettner
Journal:  Cancer Res       Date:  2006-12-01       Impact factor: 12.701

8.  Comprehensive genomic analysis of desmoplastic medulloblastomas: identification of novel amplified genes and separate evaluation of the different histological components.

Authors:  A Ehrbrecht; U Müller; M Wolter; A Hoischen; A Koch; B Radlwimmer; B Actor; A Mincheva; T Pietsch; P Lichter; G Reifenberger; R G Weber
Journal:  J Pathol       Date:  2006-03       Impact factor: 7.996

9.  The functional role of the CARM1-SNF5 complex and its associated HMT activity in transcriptional activation by thyroid hormone receptor.

Authors:  Hyo-Kyoung Choi; Kyung-Chul Choi; So-Young Oh; Hee-Bum Kang; Yoo-Hyun Lee; Seungjoo Haam; Yong-Ho Ahn; Kyung-Sup Kim; Kunhong Kim; Ho-Geun Yoon
Journal:  Exp Mol Med       Date:  2007-08-31       Impact factor: 8.718

10.  A selective jumonji H3K27 demethylase inhibitor modulates the proinflammatory macrophage response.

Authors:  Laurens Kruidenier; Chun-wa Chung; Zhongjun Cheng; John Liddle; KaHing Che; Gerard Joberty; Marcus Bantscheff; Chas Bountra; Angela Bridges; Hawa Diallo; Dirk Eberhard; Sue Hutchinson; Emma Jones; Roy Katso; Melanie Leveridge; Palwinder K Mander; Julie Mosley; Cesar Ramirez-Molina; Paul Rowland; Christopher J Schofield; Robert J Sheppard; Julia E Smith; Catherine Swales; Robert Tanner; Pamela Thomas; Anthony Tumber; Gerard Drewes; Udo Oppermann; Dinshaw J Patel; Kevin Lee; David M Wilson
Journal:  Nature       Date:  2012-08-16       Impact factor: 49.962

View more
  8 in total

1.  Dual inhibitors of LSD1 and spermine oxidase.

Authors:  Steven Holshouser; Matthew Dunworth; Tracy Murray-Stewart; Yuri K Peterson; Pieter Burger; Joy Kirkpatrick; Huan-Huan Chen; Robert A Casero; Patrick M Woster
Journal:  Medchemcomm       Date:  2019-02-08       Impact factor: 3.597

2.  Deficiency of G9a Inhibits Cell Proliferation and Activates Autophagy via Transcriptionally Regulating c-Myc Expression in Glioblastoma.

Authors:  Xiao Xue Ke; Rui Zhang; Xi Zhong; Lei Zhang; Hongjuan Cui
Journal:  Front Cell Dev Biol       Date:  2020-11-27

Review 3.  The emerging role of KDM5A in human cancer.

Authors:  Guan-Jun Yang; Ming-Hui Zhu; Xin-Jiang Lu; Yan-Jun Liu; Jian-Fei Lu; Chung-Hang Leung; Dik-Lung Ma; Jiong Chen
Journal:  J Hematol Oncol       Date:  2021-02-17       Impact factor: 17.388

Review 4.  The Role of H3K4 Trimethylation in CpG Islands Hypermethylation in Cancer.

Authors:  Giuseppe Zardo
Journal:  Biomolecules       Date:  2021-01-22

5.  Systematic Review and Meta-Analysis of Lysine-Specific Demethylase 1 Expression as a Prognostic Biomarker of Cancer Survival and Disease Progression.

Authors:  Clement Agboyibor; Jianshu Dong; Clement Y Effah; Emmanuel K Drokow; Waqar Pervaiz; Dié Li; Lei Kang; Xinli Ma; Jian Li; Zhenzhen Liu; Hong-Min Liu
Journal:  Cancer Control       Date:  2021 Jan-Dec       Impact factor: 3.302

6.  Sanguinarine, identified as a natural alkaloid LSD1 inhibitor, suppresses lung cancer cell growth and migration.

Authors:  Ting-Ting Qin; Zhong-Hua Li; Li-Xin Li; Kun Du; Ji-Ge Yang; Zhen-Qiang Zhang; Xiang-Xiang Wu; Jin-Lian Ma
Journal:  Iran J Basic Med Sci       Date:  2022-06       Impact factor: 2.532

Review 7.  Advances of Epigenetic Biomarkers and Epigenome Editing for Early Diagnosis in Breast Cancer.

Authors:  Pourya Sarvari; Pouya Sarvari; Ivonne Ramírez-Díaz; Frouzandeh Mahjoubi; Karla Rubio
Journal:  Int J Mol Sci       Date:  2022-08-23       Impact factor: 6.208

Review 8.  Salvia miltiorrhiza in cancer: Potential role in regulating MicroRNAs and epigenetic enzymes.

Authors:  Meng Lu; Xintian Lan; Xi Wu; Xiaoxue Fang; Yegang Zhang; Haoming Luo; Wenyi Gao; Donglu Wu
Journal:  Front Pharmacol       Date:  2022-09-12       Impact factor: 5.988

  8 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.