Structural basis of arginine asymmetrical dimethylation by PRMT6
PRMT6 is a type I protein arginine methyltransferase, generating the asymmetric dimethylarginine mark on proteins such as histone H3R2. Asymmetric dimethylation of histone H3R2 by PRMT6 acts as a repressive mark that antagonizes trimethylation of H3 lysine 4 by the MLL histone H3K4 methyltransferase. PRMT6 is overexpressed in several cancer types, including prostate, bladder and lung cancers; therefore, it is of great interest to develop potent and selective inhibitors for PRMT6. Here we report the synthesis of a potent bi-substrate inhibitor GMS (6\xe2\x80\xb2-methyleneamine sinefungin, an analogue of sinefungin), and the crystal structures of human PRMT6 in complex respectively with SAH and the bi-substrate inhibitor GMS that shed light on the significantly improved inhibition effect of GMS on methylation activity of PRMT6 compared to SAH and a SAM competitive methyltransferase inhibitor sinefungin (SNF). In addition, we also crystallized PRMT6 in complex with SAH and a short arginine containing peptide. Based on the structural information here and available in the PDB database, we propose a mechanism that can rationalize the distinctive arginine methylation product specificity of different types of arginine methyltransferases and pinpoint the structural determinant of such a specificity.
Introduction<!>Cloning, expression and purification<!>Inhibition assay<!>Methylarginine quantitation<!>Crystallization<!>Data Collection and Structure Determination<!>Crystal structure of PRMT6 in complex with SAH<!>Synthesis of a bi-substrate compound as a potent PRMT6 inhibitor<!>Crystal structure of PRMT6 in complex with the bi-substrate inhibitor GMS<!>Structural basis of asymmetrical dimethylation ability of PRMT6
<p>Arginine methylation is an abundant covalent post-translational modification. In two separate accounts, it was reported that about 2% of arginine residues in the total protein extracts from rat liver nuclei [1], and over 10% of proteins encoded in the T. brucei genome are arginine methylated [2]. Arginine methylation exists in three forms, i.e., mono-methylation (Rme1), asymmetric di-methylation (Rme2a) and symmetric di-methylation (Rme2s), which are carried out by three different types of protein arginine methyltransferases (PRMTs). PRMT1, 2, 3, 4, 6, and 8 (type I PRMT) generate asymmetric di-methylarginine modifications, PRMT5 and 9 (type II PRMT) generate symmetric di-methylarginine modifications, and PRMT7 (type III) is an arginine mono-methyltransferase [3, 4]. Arginine methylation is involved in a variety of cellular processes, such as transcriptional regulation, RNA processing, signal transduction, DNA repair and genomic stability [5].</p><p>Arginine methylation executes its functions normally through regulating protein-protein interactions either positively or negatively. For instance, the PRMT5-containing methylosome modifies some specific arginine sites to symmetric di-methylation in several spliceosomal Sm proteins, which are recognized by the Tudor domain of the SMN protein [6–8]. The SMN complex binds both arginine-methylated Sm proteins and snRNAs, bringing them together and facilitating Sm core assembly [9]. The same PRMT5-containing methyltransferase complex also methylates the Piwi proteins (PIWIL1-4), which contain multiple arginine-glycine (RG) and arginine-alanine (RA) repeats at their N-termini. Methylation of the Piwi proteins is required for their interaction with the TDRD group of germline-enriched Tudor domain proteins, such as TDRD1, TDRD2, TDRD4-9 and TDRD12, which leads to subsequent localization of these proteins to the meiotic nuage [10–13]. CARM1 (PRMT4) is a transcriptional coactivator that asymmetrically dimethylates histone H3R17. The histone H3R17me2a mark is recognized by the Tudor domain of TDRD3, consistent with the observation that TDRD3 is recruited to an estrogen-responsive element in a CARM1-dependent manner and promotes transcription by binding methyl-arginine marks on histone tails [8, 14]. PRMT6 is a histone H3R2 methyltransferase. Asymmetric dimethylation of histone H3R2 by PRMT6 acts as a repressive mark that antagonizes trimethylation of H3 lysine 4 by the MLL histone H3K4 methyltransferase [15–18]. Asymmetric dimethylation of histone H3R2 diminishes its binding to WDR5 and impedes the recruitment of WDR5 to euchromatic regions [19, 20]. WDR5 is a common component of the SET1/MLL family of histone H3K4 methyltransferases, and has been shown to bind different arginine containing peptides, including histone H3R2 [20–24]. Strikingly, when histone H3R2 is symmetrically dimethylated (H3R2me2s) by PRMT5, its binding to WDR5 is enhanced and the target genes are poised for transcriptional activation [19].</p><p>PRMT6 has been shown to modify histone H3R2 and a few other substrates, and it was reported recently that PRMT6 is overexpressed in several cancer types, and knockdown of PRMT6 significantly suppresses growth of bladder and lung cancer cells [25]. It was further shown that the tumor suppressor genes p21 and p27, two members of the CIP/KIP family of cyclin-dependent kinase (CDK) inhibitors, are direct targets of PRMT6 by methylating histone H3R2 in their promoters. Knockdown of PRMT6 leads to up-regulation of p21 and p27 and cellular senescence [26–28]. Hence, PRMT6 promotes cell growth and prevents senescence, making it an attractive therapeutic target for various types of cancer. In addition, PRMT6 is also implicated in the regulation of gene expression of TSP-1, a potent natural inhibitor of angiogenesis [29, 30], and the RUNX1 target genes, a group of genes implicated in differentiation of hematopoietic stem/progenitor cells [31], and genes involved in maintaining embryonic stem cell identity [32]. PRMT6 also acts as a restriction factor for viral replication in human immunodeficiency virus pathogenesis by methylating TAT and other HIV proteins [33].</p><p>So far, crystal structures for several arginine methyltransferases have been determined, but it is still unclear how different types of arginine methyltransferases achieve the distinctive arginine methylation product specificity. In addition, considering the importance of PRMT6 in gene regulation and its dysregulation implicated in many different cancers, it is of great interest to develop potent inhibitors for PRMT6. In this study, we synthesized a potent bi-substrate inhibitor GMS for PRMT6, obtained its crystal structure in complex with PRMT6, and demonstrated that GMS exhibits better inhibition than sinefungin (SNF), a natural SAM (S-adenosyl-L-methionine) analog and potent inhibitor of methyltransferases. Furthermore, we crystallized PRMT6 in complex with SAH and a short peptide of a sequence GR(me1)G. By comparing these structures to those available crystal structures of other arginine methyltransferases, we are able to explain why different types of arginine methyltransferases generate the distinctive arginine methylation products specificity.</p><!><p>A DNA fragment encoding full-length human PRMT6 (residues 1–375) was cloned into a baculovirus expression vector pFBOH-MHL (http://www.thesgc.org/sites/default/files/toronto_vectors/pFBOH-MHL.pdf). The protein was expressed in Sf9 cells (Invitrogen) with addition of 18 amino acid residues including a hexa-His tag followed by a TEV cleavage site (MHHHHHHSSGRENLYFQG) at the N-terminus. The harvested cells were resuspended in lysis buffer containing 20 mM Tris-HCl, pH 8.0, 500 mM NaCl, 5 mM imidazole, 2 mM β-mercaptoethanol, 5% glycerol, 0.6% NP-40, protease inhibitor cocktail (Roche), 3000 U of benzonase (Novagen). Cells were lysed by brief sonication. The clarified lysate was loaded onto a 2-mL TALON column (Clonetech). The column was washed with 50 column volumes of 20 mM Tris-HCl buffer, pH 8.0, containing 500 mM NaCl, 5% glycerol and 5 mM imidazole. The bound protein was eluted with elution buffer containing 20 mM Tris-HCl, pH 8.0, 500 mM NaCl, 5% glycerol and 250 mM imidazole. The eluted protein was loaded onto a Superdex200 column (GE Healthcare), equilibrated with 20 mM Tris-HCl buffer, pH 8.0 and 150 mM NaCl. Pooled fractions containing PRMT6 were subjected to TEV treatment to remove the His-tag. The protein was further purified to homogeneity by ion-exchange chromatography.</p><!><p>A scintillation proximity assay (SPA) method was used to determine the activity of PRMT6 as described previously [34, 35]. In brief, the reaction mixture (20 μL volume) was composed of 20 mM bis-tris-propane (pH 7.5), 0.01% Tween-20, 0.5% DMSO, 10 mM DTT, 50 nM PRMT6, 0.6 μM biotinylated H4(1–24) peptide (Tufts University Peptide Synthesis Core Facility, Boston, Ma). To start the reaction SAM was added at final concentration of 2.3 μM which include 57% 3H-SAM (PerkinElmer Life Sciences, cat#NET155V001MC; specific activity range 12–18 Ci/mmol) and 43% SAM (AK Scientific, Union City, CA). Due to highly acidic nature of the 3H-SAM solution, non-tritiated (cold) SAM was used to supplement the reactions when necessary. For IC50 determination, different concentrations of compounds were added to the reaction mixture. The reactions were started by adding the substrate. Reaction mixtures were incubated at 23 °C for 20 minutes and were quenched by addition of 20 μL of 7.5 M guanidinium hydrochloride (GuHCl) followed by the addition of 180 μL of 20 mM Tris-HCl, pH 8.0. The quenched reactions were then transferred to the wells of a streptavidin and scintillant-coated microplate (FlashPlate® PLUS; PerkinElmer Life Sciences). The amount of the methylated peptide was quantified by tracing the radioactivity (CPM; counts per minute) as measured by a TopCount NXT™ Microplate Scintillation Counter (PerkinElmer Life Sciences). The IC50 values were calculated using SigmaPlot® software (SYSTAT Software Inc., CA, USA). To determine the optimal pH, the reactions (with no compound) were carried out in bis-tris-propane (20 mM) at pH values between 6–9.5 (Supplementary Fig. 4). The potency of the compounds were determined against PRMT1, PRMT3, CARM1 (PRMT4), PRMT5, PRMT6, and PRMT8 as well as few lysine methyltransferases (DOT1L, G9a, SETD7, and PRDM9) using a radioactivity based method, as described previously [34, 36, 37].</p><!><p>The reactions were run overnight at 23 °C in a reaction mixture containing 20 mM bis-tris-propane (pH 7.5), 0.01% Tween-20, 0.5% DMSO, 10 mM DTT, 1 μM PRMT6, 8 μM peptide, and 5 μM of 3H-SAM (Perkin Elmer Life Sciences, specific activity range 12–18 Ci/mmol). The samples were spin-filtered at 14,000 g at 4 °C using 10-kDa molecular weight cut-off filters (Nanosep, OD010C34) for 15 min to separate the enzyme from the peptide substrate. Sample eluates were transferred into 300 μL glass tubes and dried using a Thermo Savant SC110A speed vacuum. The dried reactions were hydrolyzed with 200 μL 6 N HCl at 110 °C for 24 h in vacuo. The dried samples were reconstituted in 100 μL of 0.5% acetic acid and 0.01% trifluoroacetic acid (TFA) mixture. The samples were quantified by LC-MS/MS with previously described buffer conditions and multiple reaction monitoring protocols [38].</p><!><p>Purified PRMT6 (5 mg/mL) was incubated with SAH (S-adenosyl-L-homocysteine, Sigma), GSM and SAH/peptide at 1:5 molar ratio of protein:ligand and crystallized using the hanging drop vapor diffusion method at 20 °C by mixing 1 μL of the protein solution with 1 μL of the reservoir solution. The complex crystals were obtained in solution containing 15% PEG 3350, 0.1 M succinate acid, pH 7.0 for PRMT6-SAH, 25% PEG 3350, 0.2 M KSCN for PRMT6-GSM and 10% PEG 8000, 0.1 M Tris-HCl, pH 8.5 for PRMT6-SAH-GR(me)G peptide, respectively.</p><!><p>X-ray diffraction data for PRMT6-SAH, PRMT6-GMS and PRMT6-SAH-GR(me)G peptide complexes were collected at 100K at beam line 19ID-D of Advanced Photon Source (APS), Argonne National Laboratory. Data sets were processed using the HKL-3000 suite [39]. All the structures were solved by molecular replacement using MOLREP [40] with PDB entry 1G6Q as the search template. REFMAC [41] was used for structure refinement. The graphics program COOT [42] was used for model building and visualization. Molprobity [43] was used for structure validation.</p><!><p>The crystal structure of full-length human PRMT6 in complex with SAH has been determined at a resolution of 1.97 Å (Fig. 1). The detailed crystal diffraction data and refinement statistics are summarized in Table 1. In this high-resolution crystal structure, the N-terminal 47 residues of PRMT6 are not visible. Although all the PRMT proteins share a conserved catalytic domain (Fig. 2), they have a variable N-terminal fragment that has been proposed to regulate methyl transfer activity and substrate specificity [44]. Consistently, the N-terminal fragment of PRMT6 has been reported to be necessary and sufficient for its association with other binding partners and plays a role in substrate specificity [45].</p><p>Similar to other solved type I PRMT structures, including PRMT1 [46, 47], PRMT3 [48, 49], CARM1/PRMT4 [50, 51], mouse PRMT6 [52] and Trypanosoma brucei PRMT6 (TbPRMT6) [53], the catalytic domain of human PRMT6 consists of two domains, i.e., the N-terminal Rossmann fold, and the C-terminal β-barrel domain with a dimerization arm embedded within it (Fig. 1). Like other type I PRMT proteins, PRMT6 also exists as a dimer, with the dimerization arm (helixes α4–6, colored in green) from one monomer packing against helixes αY/Z and α1/2 of the Rossmann domain from the other monomer to form a ring-like dimer architecture (Fig. 1B). The SAH molecule is bound in an extended conformation in a pocket formed by the Rossmann fold domain (Fig. 1B and 3A). The adenine ring, ribose moiety and homocysteine carboxylate of SAH are bound by PRMT6 through a series of hydrogen bonds and salt bridges, which are highly conserved in the PRMT family (Fig. 2 and 3A).</p><p>Structural comparisons of PRMT6 to the other type I enzymes reveals two interesting structural features of PRMT6 (Fig. 3C). First, the conserved Y(F/Y)xxY motif in the N-terminus of the Rossmann fold domain points away from the SAH binding pocket. This motif for CARM1 is disordered in the absence of SAH, but becomes an ordered helix and acts as a lid to cover the adenosine part of SAH in the CARM1-SAH binary complex (Fig. 3B and 3C). As a result, the SAH molecule is almost completely buried with the methyl group of the methyl donor cofactor SAM being only accessible to the substrate arginine via a narrow opening in the active site (Fig. 3D, PDB: 3B3F) [51]. The same phenomenon was also observed in the SAH-free and SAH bound PRMT3 structures [45, 49]. On the other hand, this motif in PRMT1 was found to be disordered or point away from SAH in either the presence or absence of SAH [46, 47], which leaves the SAH molecule widely exposed to the solvent, similar to what we have observed in the PRMT6-SAH structure (Fig. 3E). It has been shown previously that Y(F/Y)xxY plays an important role in SAH binding and catalysis in PRMT1, because deletion of this motif diminishes SAH binding and abolishes enzymatic activity [47]. Hence, we believe that this aromatic residue rich motif in the highly conserved catalytic domain of the type I PRMT proteins would adopt different conformations to open and close the SAH binding site to allow exit of SAH during catalysis. Second, although all the type I PRMT proteins have a dimerization arm, structural alignment of all the available PRMT structures reveal that the dimerization arm in PRMT6 has a different conformation, which in turn leads to a flat dimer ring structure with a wide central cavity, whereas the two monomers in the PRMT1, PRMT3 and CARM1 dimers form a concave surface (Supplementary Fig. 1). The substrates access the active site of methyl transfer through the central cavity and the size and shape of the substrate binding groove determines the substrate specificity of different PRMTs [51, 54]; therefore, the unique dimer ring structure of PRMT6 may be a structural determinant for its substrate preference.</p><!><p>The crystal structures of PRMT6 and other PRMTs in complex with SAH or sinefungin further inspired us to develop sinefungin analogues containing the features of substrates as bi-substrate inhibitors of PRMTs, such as 6′-methyleneamine sinefungin (GMS, 1), which was synthesized from a known oxazolidinone imide precursor 2 [55] (Fig. 4A). The reduction of the chiral auxiliary imide of 2 led to the primary alcohol 3 with the desired steric center maintained at the C6 position of D-ribose ring. The approach was then implemented to introduce the amine functionality via 4 by Dess-Martin oxidation of the primary alcohol 3 and then reductive amination with benzyl amine, followed by carbobenzyloxy (Cbz) protection. The terminal alkene of the intermediate 4 was readily converted into the protected amino acid 5 as described previously [56]. After acidic cleavage of dihydropyrazine and N-trifuoroacetylation (TFA), the protecting groups on glycoside were exchanged with acetyl to facilitate the following adenosylation. The β-ribosyl adenine in 6 was installed through Vorbrüggen glycosylation of bis-silyl-N-benzoyladenine [56]. Here the amino acid in 6 was masked with N-trifuoroacetyl group rather than the previously-reported Cbz or acetyl group [57] for its orthogonality to hydrogenolysis and ease of removal. The hydrogenolysis condition set the primary amine on 6′ position free for further guanidylation. Trifluoroethanol was used as solvent to avoid N-alkylation in alcoholic solvent [58]. At last, global deprotection of base-labile blocking groups in the present of lithium hydroxide, followed by acidic removal of t-butoxycarbonyl completed the synthesis of the target compound GMS (Refer to the Supporting Methods for the detailed synthesis method).</p><p>We then applied the scintillation proximity assay (SPA) method to determine the inhibition activity of this compound against PRMT6, which showed an inhibition IC50 of 90 nM. That is significantly more potent than the SAM reaction product SAH and the SAH analog sinefungin (Fig. 4B). We also determined the specificity of GMS over other PRMTs and some selected lysine methyltransferases (Supplementary table 1). These data indicated that GMS preferentially inhabited some type 1 PRMTs, such as PRMT6, PRMT8 and CRAM1 (Supplementary table 1).</p><!><p>In order to understand the enhanced inhibition effect of GMS against PRMT6, we have determined the crystal structure of PRMT6 in complex with the bi-substrate inhibitor GMS at a resolution of 1.88 Å (Fig. 5A). The PRMT6-GMS binary structure is almost identical to that of the PRMT6-SAH complex, and the Cα of these two structures can be superimposed with a RMSD of just 0.25 Å. When these two complex structures are overlaid together, the SAH molecule in the PRMT6-SAH structure can perfectly superimpose with the corresponding moiety of the GMS molecule in the PRMT6-GMS complex structure. All the interactions between PRMT6 and SAH in the PRMT6-SAH binary structure are conserved in the GMS complex structure. However, GMS makes two extra hydrogen bonds through its guanidino moiety with PRMT6 (Fig. 5A). Namely, the terminal guanidino nitrogen atom in the guanidino moiety forms hydrogen bonds with the side chain carboxylate oxygen of E155 and the main chain carbonyl oxygen of M157, respectively (Fig. 5A). The extra interactions between GMS and PRMT6 presumably account for the significantly improved inhibition effect of GMS on the methylation activity of PRMT6, compared to SNF and SAH. When we superimposed the PRMT6-GMS complex structure with the PRMT1 ternary complex structure (PDB code: 1OR8 [47]), we found that the SAH molecule in the PRMT1 complex structure can also perfectly align with the SAH counterpart of the GMS molecule (Fig. 5B). Interestingly, the substrate arginine in the PRMT1 ternary complex is located in a similar position to the guanidino moiety of GMS in the PRMT6 structure (Fig. 5B). Structural comparison of the PRMT6-GMS complex with those of the recently reported other PRMT6 inhibitor complexes revealed that those inhibitors, such as MS023 [34], EPZ0204111 [59] and fragment 7 [35] only occupy the substrate arginine binding pocket, different from the GMS molecule, which occupies both substrate and cofactor binding pockets. Hence, the GMS inhibitor acts as a bi-substrate inhibitor by capturing the interacting features of both the cofactor and the substrate.</p><!><p>There are three types of protein arginine methyltransferases, which produce asymmetrical di-methylated, symmetrical di-methylated, and mono-methylated arginine modifications, respectively. But, it is still not very clear how these arginine methyltransferases generate different arginine modification marks. In order to gain insights into the different product specificities of arginine methyltransferases, it is essential to obtain the complex structures of these different types of arginine methyltransferases with their substrates. So far, the substrate complex structures have been reported for the arginine mono-methyltransferase Trypanosoma brucei PRMT7 (TbPRMT7) [60, 61], and the human arginine symmetrical di-methyltransferase PRMT5 [54]. Arginine asymmetrical di-methyltransferase PRMT1 has been crystallized with its substrate peptide, and the electron density map revealed just an arginine residue in the substrate binding groove [47]. In this study, we crystallized PRMT6 with its cofactor SAH and a short peptide of GR(me1)G, and obtained a high-resolution crystal structure. From the electron density map, some extra electron density was found near the SAM binding pocket (Supplementary Fig. 2A), and when we superimposed it with the PRMT1-SAH-GR(me1)G ternary complex structure, we found that the side chain of the substrate mono-methyl arginine residue could be fitted into the extra electron density perfectly (Supplementary Fig. 2A). In addition, as we mentioned previously that the guanidino moiety of GMS also resides in the same position as the side chain of the substrate arginine in the PRMT1-SAH-arginine complex structure. Furthermore, during preparation of the manuscript, a CARM1-SNF-peptide complex structure (PDB: 5DWQ) was published, and the target arginine interaction is consistent with what we observed in this study (Supplementary Fig. 3) [62]. Therefore, it is very likely that the extra electron density could be attributed to the mono-methyl arginine residue of the substrate peptide. Nevertheless, we could not exclude that the extra electron density might arise from some other small molecule, such as the Tris (Tris(hydroxymethyl)aminomethane) buffer molecule (Supplementary Fig. 2B), which we used in our purification buffer.</p><p>Interestingly, the putative PRMT6 ternary complex could provide us with insights into how different arginine methylation products are generated. In the PRMT6-SAH-GR(me1)G ternary complex structure (Fig. 5C), the mono-methyl arginine residue forms a salt bridge with the highly conserved E155, and this interaction is also found in the TbPRMT7-SAH-H4R3 peptide complex structure (Fig. 5E, PDB: 4M38) [60] and the PRMT5-SAH-H4R3 peptide complex structure (Fig. 5F, PDB: 4GQB) [54]. The other glutamic acid residue E164 in the double-E loop in both of our PRMT6 structures presented here points away from the substrate arginine, but in a recently reported PRMT6 structure in complex with SAH and an inhibitor EPZ0204111, E164 points towards the inhibitor although it does not form a hydrogen bond with the inhibitor (Fig. 5D, PDB: 4Y30) [59]. When we superimposed E164 from the PRMT6-SAH-EPZ0204111 structure into our PRMT6-SAH-GR(me1)G complex, it could also form a hydrogen bond with the substrate arginine (Fig. 5D). The glutamic acid residue corresponding to E164 of PRMT6 could also form one or two hydrogen bond-mediated salt bridges with the target arginine residue in the TbPRMT7 and PRMT5 complex structures, respectively (Fig. 5E and 5F). Therefore, the two glutamic acid residues in the double-E loop are engaged in the substrate target arginine recognition in all types of arginine methyltransferases.</p><p>The two glutamic acid residues in the double-E loop have been found to be critical for methyl transfer activity in different PRMTs [47, 48, 51, 63]. In both CARM1 and PRMT3 structures, the side chains of these two highly conserved glutamic acid residues point towards the pocket, and were proposed to deprotonate the target arginine residue for nucleophilic attack on the methyl donor SAM. In our PRMT6 complex structures, the side chain of E155 points to the arginine binding pocket, but the side chain of the other glutamic acid E164 points away from the pocket (Fig. 5B). The same conformation is observed in the rat PRMT1 structure [47]. The rat PRMT1 was crystallized at a low pH (~4.7), and the protein is inactive under that condition [47]. However, the corresponding glutamic acid in the yeast PRMT1 homolog also adopts the same conformation even when it was crystallized at a high pH (7.5) [46]. In addition, all the PRMT6 structures reported in this study were crystallized at around pH 7.5, which has the optimal enzymatic activity (Supplementary Fig. 4). On the other hand, in the PRMT6 structure in complex with SAH and an inhibitor EPZ0204111, E164 points to the substrate arginine-binding pocket (Fig. 5D, PDB: 4Y30). We then examined all the available structures of PRMT5 and PRMT7 in the PDB database, and found that in all of these available structures, both glutamic acid residues in the double-E loop for these type II and III enzymes always point towards the arginine-binding pocket regardless of the substrate presence in the crystal structures [54, 60, 64, 65]. Taken together, the second glutamic acid in the double-E loop (E164 in the case of PRMT6) is very dynamic in the type I PRMT proteins, and the dynamic nature may be important for catalysis [51].</p><p>The substrate mono-methyl arginine in the PRMT6-SAH-GR(me1)G complex structure also forms a hydrogen bond with the side chain of H317 in the PRMT6 structure (Fig. 5C). This histidine residue is conserved in all the other type I arginine methyltransferases, which corresponds to a glutamine residue in both human and Trypanosoma brucei PRMT7 (Q329 in TbPRMT7) (Fig. 2). In the TbPRMT7-SAH-H4R3 ternary structure, Q329 also forms a hydrogen bond with the terminal guanidino nitrogen (Fig. 5E). However, this hydrogen bond is absent in the type II arginine methyltransferase PRMT5, in which the corresponding residue S578 is far away from the substrate arginine (Fig. 5F). Based on the available structural information and previous catalytic studies revealing that arginine di-methylation is processive [66–68], possible methylation mechanisms could be proposed for these different types of arginine methyltransferases. All the arginine methyltransferases utilize the two glutamic acid residues in the double-E loop to deprotonate the Nη1 atom of the target arginine residue, which leads to a methyl group transfer to the Nη1 atom (Fig. 2, 5C, 5E and 5F). Once the first methyl group is transferred, for PRMT7, due to its limited space around the Nη1 atom, an additional methyl group cannot be added to the Nη1 atom to produce asymmetrical di-methylation (Fig. 5E) [60]. The importance of this limited space in determining the product specificity was confirmed by a recent study, in which a point mutation E181D in TbPTMT7, which creates an enlarging space around the Nη1 atom, converts TbPTMT7 from a type III PRMT to a type I PRMT [61]. On the other hand, the hydrogen bond interactions between Q329 or E172, and Nη2 of the substrate arginine would deter methyl-Nη1 and Nη2 from swapping positions because methyl-Nη1 would cause steric clashes with Q329 and E172 when it swaps positions with Nη2 (Fig. 5E). That explains why the type III arginine methyltransferases can only carry out mono-methylation. For the type II arginine methyltransferases, based on the PRMT5-substrate complex structure, S578 does not contact the substrate arginine. This setting would create enough room to accommodate an extra methyl group on the Nη2 side (Fig. 5F). Therefore, once a methyl group is attached to the Nη1 atom, because the second glutamic acid in the double-E loop forms double hydrogen bonds with the substrate arginine, which limits this nitrogen atom to take the second methyl group, the methyl-Nη1 would then rotate to swap positions with Nη2, and Nη2 would be deprotonated again to accept another methyl group to form symmetrical dimethylation. This rotation is necessary because Nη2 is far from the reactive methyl group of the methyl donor SAM to accept the methyl transfer directly (Fig. 5E). For the type I arginine methyltransferases, because the second glutamic acid is very dynamic (Fig. 5C and 5D), it should take a different conformation once it aids the first glutamic acid to deprotonate the target arginine; therefore, two methyl groups could be transferred to the Nη1 atom to form asymmetrical di-methylation.</p><p>Because the type I specific histidine residue (H317 in PRMT6) is a critical structural determinant in the asymmetrical di-methylation ability of the type I arginine methyltransferases, mutating it to a type II specific serine residue could potentially alter its product specificity. We, therefore, made a H317S PRMT6 mutant. Our enzymatic assay results revealed that H317S mutant only displayed a slightly reduced activity, but it could not produce any detectable symmetrical di-methylation (Fig. 6A). Structural comparison of the PRMT6, PRMT5 and TbPRMT7 substrate complexes revealed that the type I specific histidine residue (H317 in PRMT6), the type II specific serine residue (S578 in PRMT5) and the type III specific glutamine (Q329 in TbPRMT7) are located in a structural motif, which is highly conserved in their own subfamily, but diverse among these different types of arginine methyltransferases (Fig. 2 and Fig. 6B). The serine residue (S578 in PRMT5) is significantly farther to the substrate arginine residue than its corresponding residues in type I and III arginine methyltransferases, which would create big enough room to accommodate a methyl group (Fig. 6B). Therefore, the structural element harboring the type-specific residue (H317 in PRMT6) as a whole is the structural determinant in conferring the arginine methylation product specificity of different types of arginine methyltransferases.</p>