1. Introduction
Post-translational modifications (PTMs) are covalent modifications made via protein–protein interactions that control the activation and deactivation of cellular pathways. To date, roughly 400 distinct PTMs have been identified [
1]. PMTs have been observed on both histone and non-histone targets and are involved in a wide variety of cellular processes, including signal transduction, the control of processes such as protein synthesis and chromatin organization, and disease phenotype presentation [
2,
3]. There are three generalized categories of enzymes which catalyze PTMs, including “writers”, which are enzymes that carry out the chemical modification of a substrate, “readers”, which are enzymes that contain a recognition domain involved in identifying the modification carried out by the “writer” enzymes, leading to the downstream changes in signaling based on the PTM, and “erasers”, which are the enzymes responsible for reversing or removing the PTM [
3,
4]. An enzyme can have multiple functions as well, serving simultaneously as any of the roles mentioned above for different targets. Commonly studied PTMs include acetylation, ubiquitination, phosphorylation, and methylation [
2,
3]. Our understanding of PTMs and their impact on protein biology has extended the known functional human proteome. PTMs act as modulators of protein activity, making them a very valuable target for therapeutic development.
The rate and occurrence of PTMs change in disease phenotypes, suggesting that specific PTMs may be used as biomarkers for disease prognosis [
3]. Cells which form cancerous tumors display different growth pathways, are more resistant to apoptosis, and display higher levels of mechanical stress in comparison to their non-cancerous counterparts. This is in part due to differences in PTM-mediated regulatory pathways [
3]. As such, specific PTMs have been characterized as being involved in tumorigenesis, apoptosis resistance, and malignancy; however, only a small fraction of identified PTMs have been studied extensively in the context of disease progression [
3]. This remains an active and important area of research to better understand how disease phenotypes manifest. Additionally, targeting disease phenotype PTMs has been suggested as a promising avenue for therapeutic development [
5].
One notable PTM is protein methylation, which entails the transfer of a methyl group (-CH3) from a methyl donor, S-adenosylmethionine (SAM), to a target. Enzymes that catalyze methyl transfer are referred to as methyltransferases (MTases). Methyl transfer can occur on the side chains of amino acid residues such as lysine, arginine [
6,
7,
8], histidine [
9], glutamate, and aspartate [
10]. Additionally, the N-terminal of proteins can undergo methylation either directly on the amino group of methionine or on the amino group of the following amino acid, assuming methionine is cleaved via methionine aminopeptidases [
11,
12]. Focusing primarily on lysine and N-terminal methylation, these enzymes are designated as lysine methyltransferases (KMTs) and N-terminal methyltransferases (NTMTs), respectively [
8,
11,
12]. Lysine methyl transfer is known to be a reversible process, with demethylation catalyzed by enzymes referred to as lysine demethylases (KDMs); however, to date N-terminal methylation is thought to be an irreversible process, and no validated N-terminal demethylases have been identified [
11].
There are two main structural classes of KMTs, including the suppressor of variegation 3–9, the enhancer of zeste, and Trithorax (SET)-domain-containing KMTs and seven-beta-strand (7βS) KMTs. SET KMTs are typically responsible for histone lysine methylation events, and 7βS KMTs are currently known to be responsible for non-histone lysine targets [
13]; however, exceptions do exist for both classes as more validated substrates are documented. The SET-domain-containing class contains roughly 50 known KMTs, while the 7βS KMT class includes over 130 enzymes, making it the largest class of MTases in the human proteome [
13].
The mechanism of methyl transfer of the SET-domain-containing class has been characterized in much greater detail than that of the 7βS class [
14,
15]. Although, in general, methyl transfer is assumed to follow an SN2 reaction mechanism, the catalytic binding regions of the two classes differ greatly, suggesting that the two may use different binding mechanisms [
14]. A computational study has modeled methyl transfer in 7βS KMT DOTL1 and suggested that the mechanism of methyl transfer differs greatly from that of the SET-domain-containing KMTs in terms of SAM-binding orientation, as well as substrate coordination [
14,
16]. There is very little structure homology within the 7βS class, so making generalized claims regarding the structural classes’ mechanisms based on the computational modelling of DOTL1 is not possible [
14]. This, however, highlights the importance of additional research regarding structural interactions and determining binding mechanisms among 7βS KMTs.
The general process of methyl transfer is the same, whether catalyzed by a KMT or an NTMT. The MTase binds a substrate and SAM, the methyl donor, and catalyzes the transfer of a methyl group, resulting in the formation of S-adenosylhomocysteine (SAH), the demethylated form of SAM, and a methylated substrate [
8,
12]. Both the ε amino group of lysine and the α amino group of an N-terminal can occupy various methylation states. These are the addition of a single methyl-group (monomethylation; me1), the addition of two methyl-groups (dimethylation; me2), or the addition of three methyl-groups (trimethylation; me3). Increasing the number of methyl groups on a residue increases both the hydrophobic character of the methylated substrate and the steric bulkiness, influencing possible interactions. Trimethylation also renders the methylated amino group chemically inert by inducing a permanent positive charge and chemical saturation [
8]. It is assumed that N-terminal methylation protects against protein degradation by inhibiting N-terminal acetylation or ubiquitination [
11,
12].
Methyltransferase-like protein 13 (METTL13) is an enzyme capable of both lysine and N-terminal methylation, rendering it a dual methyltransferase, and is the subject of the remainder of this review.
2. Structure and Function of METTL13
METTL13 is a methyltransferase with two distinct seven-beta-strand (7βS) catalytic domains responsible for catalyzing the SAM-dependent methylation of both isoforms of the eukaryotic elongation factor (eEFIA), eEFIA1 and eEFIA2, at two sites: lysine 55 and the N-terminal glycine [
17]. Contextually, eEFIA1 is the more prevalent isoform, present in all cells, whereas eEFIA2 is typically only present in post-mitotic cells such as myocytes and neurons. Both isoforms, however, are active in cancerous cells [
17,
18]. Broadly, eEFIA is a GTPase responsible for delivering aminoacyl-transfer RNA molecules to the ‘A site’ of a ribosome for protein translation [
8,
17]. When METTL13 was initially identified as a novel enzyme involved in tumorigenesis, it was first given the name faint expression in normal tissues, aberrant overexpression in tumors (FEATs) [
19]. This name highlights a very important point regarding METTL13 biology in disease phenotypes: treatments targeting METTL13 that result in decreased cell proliferation in cancerous cells typically do not result in decreased cell proliferation of healthy cells [
19,
20]. Although there is currently no targeted inhibitor for METTL13, inhibitor development for numerous MTases has begun [
4,
21], with some currently in stages of clinical trials [
22,
23,
24]. These advancements, along with METTL13 expression levels in healthy versus tumor cells, highlight METTL13 as a promising therapeutic target.
The subsequent study of the structure and function of METTL13 demonstrated that each catalytic domain is responsible for just one of the currently identified methylation events [
17]. The N-terminal domain (MT13-N), which encompasses residues 46 through 160, catalyzes the dimethylation of lysine 55, and the C-terminal domain (MT13-C), encompassing residues 499 through 673, catalyzes the trimethylation of the N-terminal [
25]. These domains are structurally unique from each other, and each can function both in vitro and in vivo as orphan enzymes [
17]. It is important to note that the residues contained within each isolated domain may vary among researchers; however, the residues given above result in the retention of the enzymatic function, providing a general size of each domain.
METTL13 is known to regulate the eEFIA function. Knockout (KO) experiments have demonstrated that removing METTL13 from cells affects the occupancy levels of different amino acid codons. METTL13 KO increases the occupancy of all lysine and histidine codons, as well as [ACG] for threonine, [CCG] for proline, and [AGG] for arginine, resulting in a slower rate of translation for these codons. Conversely, METTL13 KO decreases the occupancy of all alanine codons, as well as [TGG] and [TAC] for tryptophan, [AAC] for aspartic acid, and [TCA] and [TCC] for serine. This means that these codons are translated faster in the absence of METTL13 [
17]. This clearly demonstrates that the absence of METTL13 has a direct effect on the rate of protein translation.
As stated, both METTL13 domains are 7βS methyltransferases, and METTL13, to date, is the only known member of the 7βS family, which contains two distinct catalytic domains. This classification is made based on the presence of distinct structural features, including the Rossmann fold characteristics, as well as similarities in SAM binding regions [
8]. A schematic visualization of both METTL13 domains is shown in
Figure 1A. Although many 7βS enzymes show these key structural features, there is very little structural homology in the family. Even the two subunits of METTL13 are not closely related: MT13-N most closely resembles other 7βS KMTs, including eEFIA1-KMT2, eEFIA1-KMT4, and METTL10, whereas MT13-C most closely resembles spermidine synthase, which is not a methyltransferase [
8]. Enzymes belonging to this class have seven alternating beta strands in series, separated by alpha helices [
8]. Other relevant structural features include the presence of the Post I and Post II domains, which are involved in SAM coordination and substrate recognition [
17].
In MT13-C, the Post I domain includes all residues involved in the first two beta strands, as well as sections of their downstream loop structures [
26]. The Post II domain has been identified in MT13-C as the region downstream of beta strand 4 and is thought to be involved in substrate recognition [
17]. An excess of SAH has been shown to decrease the activity of METTL13, suggesting that SAH could function as a competitive inhibitor for METTL13, as it does for several other methyltransferases [
17]. Specific amino acids identified to be potentially involved in substrate recognition include Asp575, Asp577, and Asn647 [
26]. In addition, Gly503 and Glu524 have been identified as coordinating SAH, in a region analogous to S-adenosylmethionine coordination in spermidine synthase [
8,
17]. A schematic block diagram showing the sites is shown in
Figure 1B. To date, the structural conformation of MT13-N, as well as full-length METTL13, has only been predicted using protein folding software such as AlphaFold, so the crystallization of MT13-N and full-length METTL13 in complex with SAM, SAH, or substrates is warranted to facilitate identification of characteristic structural features.
The N-terminal methylation of eEFIA has been observed in yeasts and humans, suggesting that this modification is important in eEFIA regulation [
11]. As stated above, all identified NTMTs in the human proteome have a 7βS structure: MT13-C, N-terminal methyltransferase 1 and 2 (NTMT 1 and 2 respectively), and METTL11 [
11]. Most NTMTs have a recognition motif of X-P-[K/R], resulting in N-terminal lysine or arginine methylation. The recognition motif for MT13-C has been identified as M-[GAP]-[KRFYQH]-E-[KRQHIL], suggesting 49 possible targets for MT13-C in the human proteome [
11,
17]. Subsequent in vitro study of peptides derived from the proteins containing the identified binding motif were unable to validate any novel sites of methylation catalyzed by MT13-C, indicative of a possible high level of stringent specificity by MT13-C [
17].
As stated previously, N-terminal methylation can occur on the first amino acid in any peptide sequence synthesized intracellularly, typically methionine. Alternatively, it can be removed by methionine aminopeptidases, allowing the amino group of the following residue to be targeted by NTMTs. The latter is the case for eEFIA methylation by MT13-C, as methylation occurs on the N-terminal of glycine 2 [
11]. Orphan MT13-C can methylate the N-terminal of full length eEFIA, and it has been shown that this subunit can also catalyze the methylation of peptides, including the recognition motif, suggesting that the sequence is significant in substrate docking, as opposed to the structure [
17]. The structure of MT13-C in complex with SAH has been crystalized as PDB 5WCJ (
https://www.rcsb.org/structure/5wcj, accessed on 3 June 2024) [
17], and subsequently docking experiments were carried out with the hexapeptide GKEKTH, both supporting the suggested mechanism of SAM/SAH and substrate binding [
17].
Lysine methylation is a well-defined, reversible modification observed in both histones and non-histone targets. Instances of lysine methylation are now widespread, and it is known that this modification is involved in the regulation of numerous pathways in different disease states. The N-terminal domain of METTL13, MT13-N, catalyzes the trimethylation of eEFIA lysine 55, and it is generally less understood than MT13-C. Orphan MT13-N can methylate full-length eEFIA-K55 in vitro; however, to date, in vitro methylation of a peptide form of eEFIA-K55 has not been validated, suggesting that perhaps MT13-N requires more higher level structural information to catalyze methylation compared to MT13-C [
17,
25]. It has been shown that the occupancy levels of dimethylated lysine 55 of eEFIA, hereby referred to as eEFIA-K55me2, are higher in KRAS-driven cancers, and appear to be essential to meeting the increased translational output required for promoting tumorigenesis [
25]. eEFIA1-K55 is located on the GTPase catalytic surface of eEFIA; therefore, the methylation of K55 is thought to modulate the GTPase activity directly. It has since been demonstrated that the presence of eEFIA-K55me2 increased the GTPase catalytic efficiency of eEFIA by roughly 20% [
25]. It was also shown that when METTL13 is depleted in model cell lines, global protein synthesis is decreased, indicating the eEFIA-K55me2 is directly involved in protein synthesis modulation [
25]. In general, although MT13-N bears more structural resemblance to known 7βS KMTs, its structure, function, and activity are less understood than those of MT13-C.
As previously stated, METTL13 was initially given the name FEAT, as it is faintly expressed in healthy tissues, and aberrantly expressed in tumors [
19,
26]. METTL13 expression has been identified as slightly elevated in normal testis [
26]; however, in general, it has been shown that modulating METTL13 in healthy tissues has little to no effect on cell proliferation and survival [
17,
25], suggesting that METTL13 does not have a fundamental role in normal biology in healthy cells. Despite this, the research area surrounding METTL13 remains limited due to challenges associated with protein expression and the in vitro study of the enzyme.
3. METTL13 Regulation
The expression of METTL13 is regulated at the transcriptional, translational, and post-translational levels. At the transcriptional level, cases of METTL13 regulation are present in a number of different disease states, to be discussed in a case-by-case basis later in this review.
At the translational level, METTL13 expression is regulated by microRNAs (miRNAs). METTL13 was identified as a human homolog for a suppressor of rat neuronal apoptosis [
27], and therefore, the role of METTL13 in the context of apoptosis was an area of initial study. A common proapoptotic miRNA, miR-16, was identified as an upstream regulator of METTL13 expression [
27]. Liang et al. demonstrated that METTL13 was overexpressed at the protein level in various cancer tissues relative to their non-cancerous counterparts; however, at the mRNA level, the level of overexpression was much less substantial, suggesting that METTL13 expression may be regulated at a post-transcriptional level, a function consistent with that of many miRNAs [
27]. Based on computational binding, modelled between common miRNAs and METTL13 mRNA, it was suggested that a binding interaction between miR-16 and the 3’-UTR of METTL13 mRNA was possible both through base-pairing and thermodynamically [
27]. Liang et al. subsequently demonstrated that in the context of cancer cell lines, the overexpression of miR-16 resulted in a decrease in METTL13 expression at the protein level [
27].
At the post-translational level, METTL13 forms a regulatory complex with two other known N-terminal methyltransferases, METTL11A and METTL11B, also known as NTMT1 and NTMT2, respectively [
28]. For context, METTL11A can exist as a monomer or a dimer, and it catalyzes the trimethylation of a variety of protein N-terminals, while METTL11B exists as a monomer and catalyzes the monomethylation of different protein N-terminals, following methionine cleavage. The dimerization of METTL11A provides stability to the enzyme, and the formation of a heterotrimer between a METTL11A dimer and METTL11B provides further complex stability [
28]. Aside from a strictly structural interaction, the formation of the METTL11A/METTL11B heterotrimer decreases the substrate specificity of METTL11A, allowing for more non-canonical substrates to be targeted by the enzyme. The physical interaction and resulting regulatory axes between different MTases are common within RNA and DNA methyltransferases [
28], and Parker and Tooley sought to determine if METTL13 were involved in the regulatory pathway as well. Their work demonstrated that METTL11A, regardless of METTL11B involvement, had a non-catalytic interaction with METTL13, which resulted in mutual regulation [
28].
The complexation of METTL13 with METTL11A inhibits the trimethylation activity catalyzed by the METTL11A and results in both an increase in Km and a decrease in V
max [
28]. The same complex formation has opposing regulatory effects on METTL13 activity. METTL11A binding inhibits the trimethylation of the eEF1A N-terminal, with no significant effect on Km and a decrease in V
max, suggesting that the interaction may be that of a standard non-competitive inhibitor [
28]. In contrast, METTL11A complexation with METTL13 promotes the dimethylation of eEF1A K55, the canonical target of the MT13-N [
28]. Although the results of the complex formation affect the catalytic efficiency of both enzymes, a catalytically active enzyme is not needed to elicit the regulatory effect on the other [
28].
To further study how the two catalytic domains of METTL13 are differentially involved in the METTL11A/METTL13 complex, truncated versions of the enzyme were constructed, where residues 1 through 344 were used to represent the N-terminal catalytic domain, MT13-N, and residues 345 through 699 were used to represent the C-terminal catalytic domain, MT13-C [
28]. Co-immunoprecipitation experiments conducted by Parker and Tooley demonstrated that MT13-N, which catalyzes the dimethylation of eEFIA-K55, interacts with METTL11A directly.
Finally, to determine if a METTL11A/METTL11B/METTL13 regulatory complex exists, the regulatory effect of all interactions was studied both in the presence and absence of METTL11A. This demonstrated that there is no direct interaction between METTL13 and METTL11B. Additionally, the inhibitory effect of METTL13 binding to METTL11A outcompetes the increased substrate promiscuity effects of the METTL11A/METTL11B heterotrimer [
28]. The authors do note that the cellular localization of METTL13 and METTL11B are different. METTL13 is primarily cytoplasmic, whereas METTL11B is nuclear, suggesting that the opposing regulatory effects on METTL11A of the different methyltransferases have functionalities in different organelles; however, in the case where a METTL11A/METTL11B/METTL13 complex is formed, the effects of METTL13 on the METTL11A activity are dominant [
28].